【BZOJ1951】古代猪文（CRT，卢卡斯定理）
题面
题解
要求什么很显然吧。。。
[Ans=G^{sum_{k|N}{C_N^k}}
]
给定的模数是一个质数，要求解的东西相当于是上面那坨东西的结果对于(varphi)的取值。
但是(varphi)不是质数，不好直接(Lucas)定理，把(varphi)分解质因数之后，
直接(CRT)合并结果就好了，所以这个就是(ex\_Lucas)
#include<iostream>
#include<cstdio>
using namespace std;
#define ll long long
#define RG register
#define MAX 50000
#define MOD (999911659)
#define phi (999911658)
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
int fpow(int a,int b,int P)
{
int s=1;if(!a)return 0;
while(b){if(b&1)s=1ll*s*a%P;a=1ll*a*a%P;b>>=1;}
return s;
}
int pri[5]={0,2,3,4679,35617},tot=4;
int jc[5][MAX],jv[5][MAX],N,G,ans;
void pre(int P)
{
jc[P][0]=1;
for(int i=1;i<=pri[P];++i)jc[P][i]=1ll*jc[P][i-1]*i%pri[P];
jv[P][pri[P]-1]=fpow(jc[P][pri[P]-1],pri[P]-2,pri[P]);
for(int i=pri[P]-2;~i;--i)jv[P][i]=1ll*jv[P][i+1]*(i+1)%pri[P];
}
int C(int n,int m,int P){return 1ll*jc[P][n]*jv[P][m]%pri[P]*jv[P][n-m]%pri[P];}
int Lucas(int n,int m,int P)
{
if(m>n)return 0;if(!m)return 1;
if(n<pri[P]&&m<pri[P])return C(n,m,P);
return 1ll*Lucas(n/pri[P],m/pri[P],P)*Lucas(n%pri[P],m%pri[P],P)%pri[P];
}
int exgcd(int a,int b,int &x,int &y)
{
if(!b){x=1;y=0;return a;}
int d=exgcd(b,a%b,y,x);
y-=a/b*x;return d;
}
int CRT(int k)
{
int x,y,a,ret=0;
for(int i=1;i<=4;++i)
{
a=Lucas(N,k,i);
exgcd(phi/pri[i],pri[i],x,y);
x=(x%pri[i]+pri[i])%pri[i];
ret=(ret+1ll*a*x%phi*(phi/pri[i])%phi)%phi;
}
return (ret+phi)%phi;
}
int main()
{
pre(1);pre(2);pre(3);pre(4);
N=read();G=read()%MOD;
for(int i=1;i*i<=N;++i)
if(N%i==0)
{
ans=(ans+CRT(i))%phi;
if(i*i!=N)ans=(ans+CRT(N/i))%phi;
}
ans=fpow(G,ans,MOD);printf("%d
",ans);
return 0;
}