【hdu4704】 Sum——费马小定理+快速幂取模

题目链接

题目的大意就是要求把一个数划分为m（1<=m<=n）个数相加之和的方案数（答案取模1e9+7）。

容易想到答案就是sum（C（n-1,i）（0<=i<=n-1）），也就是2^（n-1）（组合数第n行的所有数的和为2^（n-1）)。

既然答案是2^n-1，第一反应应该是快速幂取模。很可惜，这里的n巨大无比，这种做法显然是不现实的。于是我们就需要用到费马小定理啦。

我们知道 a^p-1%p=1，那么我们令n=t*(p-1)+k，即n-1=t*(p-1)+(k-1)。那么2^n-1=2^t*(p-1) * 2^k-1，重点来了！根据同余的性质，2^t*(p-1) %p=1。也就是说我们只需要求出2^k-1%p的值，这个值就是最后的答案，变形到这里，求出k的值就可以用快速幂了。

#include<cstdio>
#include<cstring>
#include<algorithm>
const int mod=1e9+7;
typedef long long LL;
using namespace std;
char ch[1000005];
int len;
LL kuai(int aa,LL b)
{
    LL res=1;LL a=aa;
    while(b){
        if(b&1)res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
LL fei(int p)
{
    LL ans=0;
    for(int i=1;i<=len;i++)
        ans=(ans*10+ch[i]-48)%p;
    return ans;
}
int main()
{
    while(~scanf("%s",ch+1)){
        len=strlen(ch+1);
        LL k=fei(mod-1);
        LL an=kuai(2,k-1);
        printf("%lld
",an);
    }
    return 0;
}