【BZOJ】【中国剩余定理】1951: [Sdoi2010]古代猪文

最新推荐文章于 2019-07-05 14:08:00 发布

Sdywolf

最新推荐文章于 2019-07-05 14:08:00 发布

阅读量333

点赞数

分类专栏： BZOJ 中国剩余定理

本文链接：https://blog.csdn.net/FYOIER/article/details/82118601

版权

BZOJ 同时被 2 个专栏收录

46 篇文章 0 订阅

订阅专栏

中国剩余定理

1 篇文章 0 订阅

订阅专栏

题意

求

G \sum i | n (n i) (mod 999911659)

$G^{\sum_{i|n}{n\choose i}}\pmod{999911659}$

题解

由费马小定理

a p - 1 \equiv 1 (mod p)

$a^{p-1}\equiv 1\pmod p$
得

G \sum i | n (n i) (mod t t) = G \sum i | n (n i) (mod t t - 1) (mod t t)

$G^{\sum{i|n}{n\choose i}}\pmod{tt}=G^{\sum{i|n}{n\choose i}\pmod {tt-1}}\pmod{tt}$
但是，由于

n n $n$ 很大，所以无法直接求组合数，考虑将

t t - 1

$tt-1$ 质因数分解后用

Lucas L u c a s $Lucas$ 算出每个答案对每个质因数取模的余数，再用中国剩余定理合并。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int maxn=16,maxs=35617,tt=999911659;
LL n,m,G,M[maxn],a[maxn],pr[maxs+6],f[maxn][maxs+6],inv[maxn][maxs+6];
bool vis[maxs+6];
void make_pr(){
    for(int i=2;i<=maxs;i++){
        if(!vis[i])pr[++pr[0]]=i;
        for(int j=1;j<=pr[0]&&pr[j]<=maxs/i;j++){
            vis[pr[j]*i]=1;
            if(!(i%(pr[j])))break;
        }
    }
}
LL C(LL x,LL y,int t){if(x<y)return 0;return f[t][x]*inv[t][y]%M[t]*inv[t][x-y]%M[t];}
LL lucas(LL x,LL y,int t){
    if(!y)return 1;
    return C(x%M[t],y%M[t],t)*lucas(x/M[t],y/M[t],t)%M[t];
}
LL power(LL x,LL y){
    LL w=x,sum=1;if(!(x%tt))sum=0;
    while(y){
        if(y&1)(sum*=w)%=tt;
        y>>=1;(w*=w)%=tt;
    }return sum;
}
LL gcd(LL a,LL b){return !b?a:gcd(b,a%b);}
LL exgcd(LL a,LL b,LL&x,LL&y){if(!b)return x=1,y=0,a;LL r=exgcd(b,a%b,y,x);return y-=a/b*x,r;}
LL lcm(LL a,LL b){return a/gcd(a,b)*b;}
LL solve(LL a,LL b,LL c){
    b%=c;if(b<0)b+=c;;LL x,y,r=exgcd(a,c,x,y);if(b%r)return -1;
    c/=r;return (x*(b/r)%c+c)%c;
}
LL CRT(int len,LL*a,LL*m){
    LL A=a[1],M=m[1];
    for(int i=2;i<=len;i++){
        LL x=solve(m[i],A-a[i],M);if(x==-1)return -1;
        M=lcm(M,m[i]);A=m[i]*x%M+a[i];if(A>=M)A-=M;
    }return A;
}
int main(){
    freopen("pig.in","r",stdin);freopen("pig.out","w",stdout);
    scanf("%lld%lld",&n,&G);make_pr();LL X=tt-1;
    for(int i=1;i<=pr[0]&&pr[i]<=sqrt(tt-1);i++) if(!(X%pr[i])){
        M[++m]=pr[i];
        while(X>1&&!(X%pr[i]))X/=pr[i];
    }
    if(X>1)M[++m]=X;
    for(int t=1;t<=m;t++){
        f[t][0]=inv[t][1]=inv[t][0]=1;
        for(int i=1;i<=maxs;i++)f[t][i]=f[t][i-1]*i%M[t];
        for(int i=2;i<=maxs;i++)inv[t][i]=(M[t]-(LL)M[t]/i*inv[t][M[t]%i]%M[t])%M[t];
        for(int i=2;i<=maxs;i++)(inv[t][i]*=inv[t][i-1])%=M[t];
    }
    for(int i=1;i<=sqrt(n);i++) if(!(n%i)){
        for(int j=1;j<=m;j++) (a[j]+=lucas(n,i,j))%=M[j];
        if(i*i<n)for(int j=1;j<=m;j++)(a[j]+=lucas(n,n/i,j))%=M[j];
    }
    printf("%lld",power(G,CRT(m,a,M)));
    return 0;
}