日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

(小费马定理降幂)Sum

發布時間:2025/3/12 编程问答 17 豆豆
生活随笔 收集整理的這篇文章主要介紹了 (小费马定理降幂)Sum 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

題目:

分析與解答:

參考思路:
https://www.cnblogs.com/stepping/p/7144512.html
https://blog.csdn.net/strangedbly/article/details/50996908
根據隔板定理,把N分成一份的分法數為C(1,n-1),
把N分成兩份的分法數為C(2,n-1),
把N分成三份的分法數為C(3,n-1),…. ,
把N分成N份的分法數為C(n-1,n-1)。
設sum=S(1)+S(2)+…+S(N),根據組合數求和公式,sum=2^(n-1)。
原式可化為(2^(N-1))mod(10^9+7)
但是現在n太大了,需要考慮簡化。

對于(a^N)mod(p)
使用費馬小定理的前提條件是
p是質數,且gcd(a,p)=1,那么 a^(p-1)≡1(mod p)
費馬小定理降冪例子:
p=7
求(2^32)%p
由費馬小定理(2^6)%p=1
所以((2^6)^k)%p=1 (k=1.2.。。。n)
(2^32)%p=(2^[(6*5)+2])%p
=[2^(6*5)*2^2]%p
=[(2^(6*5)%p)(2^2%p)]%p //(a*b)%m=[(a%m)(b%m)]%m;
=[1*(2^2%p)]%p //2^(6*5)%p為對費馬小定理的應用
=2^2%p;

所以我們推廣一下
求(a^N)mod(p)
a^(p-1)≡1(mod p)
k=N%(p-1)
(a^N)mod(p)=(a^k)mod(p)

對于這一題
(2^(N-1))mod(10^9+7)
k=(N-1)mod(10^9+7-1)=Nmod(10^9+7-1)-1
最后求(2^k)mod(p)
由于N很大,我們只能用字符數組從一位開始%mod

參考代碼:
https://blog.csdn.net/valieli/article/details/52806037

#include <iostream> #include <stdio.h> using namespace std;char N[100005]; const int mod=1e9+7; //費馬小定理 long long PhiMa(char *N,int mod) {long long ans=0;for(int i=0;N[i]!='\0';i++)ans=(ans*10+N[i]-'0')%mod;return ans; } //快速冪取余運算 long long FastPow(long long a,long long b) {long long ans=1;while(b){if(b&1)ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans; } int main() {while(scanf("%s",N)!=EOF){long long n=PhiMa(N,mod-1)-1;long long pp=FastPow(2,n);printf("%lld\n",pp);}return 0; }

總結

以上是生活随笔為你收集整理的(小费马定理降幂)Sum的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。