luogu P2480 [SDOI2010]古代猪文(Lucas定理&中国剩余定理)

luogu P2480 [SDOI2010]古代猪文(Lucas定理&中国剩余定理)

题目大意

iPig打算研究古时某个朝代的猪文文字。根据相关文献记载，那个朝代流传的猪文文字恰好为远古时期的k分之一，其中k是N的一个正约数（可以是1和N）。不过具体是哪k分之一，以及k是多少，由于历史过于久远，已经无从考证了。

iPig觉得只要符合文献，每一种能整除N的k都是有可能的。他打算考虑到所有可能的k。显然当k等于某个定值时，该朝的猪文文字个数为N / k。然而从N个文字中保留下N / k个的情况也是相当多的。iPig预计，如果所有可能的k的所有情况数加起来为P的话，那么他研究古代文字的代价将会是G的P次方。

现在他想知道猪王国研究古代文字的代价是多少。由于iPig觉得这个数字可能是天文数字，所以你只需要告诉他答案除以999911659的余数就可以了。

解题思路

根据欧拉公式显然有模数有循环节phi(mod)，据此对组合数做Lucas即可，由于模数不是质数因此需要对每个质数分别求组合数最后用CRT合并

AC代码

#include<bits/stdc++.h>
using namespace std;
const int phi=999911658;
const int mod=999911659;
const int size=1e5+5;
int a[10],b[10],k;
typedef long long LL;
int fac[size];
int quick_pow(int a,int b,int p)
{
	int ans=1;
	while(b)
	{
		if(b&1) ans=(LL)ans*a%p;
		a=(LL)a*a%p;
		b>>=1;
	}
	return ans;
}
#define int long long 
void exgcd(int a,int b,int &x,int &y)
{
    if(b==0){ x=1; y=0; return;}
    exgcd(b,a%b,x,y);
    int tp=x;
    x=y; y=tp-a/b*y;
}

int CRT()
{
    int ans=0,lcm=1,x,y;
    for(int i=1;i<=k;++i) lcm*=b[i];
    for(int i=1;i<=k;++i)
    {
        int tp=lcm/b[i];
        exgcd(tp,b[i],x,y);
        x=(x%b[i]+b[i])%b[i];//x要为最小非负整数解
        ans=(ans+tp*x*a[i])%lcm;
    }
    return (ans+lcm)%lcm;
}
#undef int
void init(int p)
{
	fac[0]=1;
	for(int i=1;i<=p;i++) fac[i]=(LL)fac[i-1]*(i)%p;
}
LL combi(int n,int m,int p)
{
	if(n>m) return 0;
	return (LL)fac[m]*quick_pow(fac[n],p-2,p)%p*quick_pow(fac[m-n],p-2,p)%p;
} 
LL Lucas(LL n,LL m,int p)
{
	if(n>m) return 0;if(n==0||n==m) return 1;
	return Lucas(n/p,m/p,p)*combi(n%p,m%p,p)%p;
}
int main()
{
	k=4;
	b[1]=2,b[2]=3,b[3]=4679,b[4]=35617;
	LL n,g;
	scanf("%lld%lld",&n,&g);
	if(g%mod==0) 
	{
		puts("0");
		return 0;
	}
	int up=sqrt(n);
	for(int i=1;i<=k;i++)
	{
		init(b[i]);
		for(int j=1;j<=up;j++)
		{
			if(n%j==0)
			{
				a[i]=(a[i]+Lucas(j,n,b[i]))%b[i];
				if(j*j!=n) a[i]=(a[i]+Lucas(n/j,n,b[i]))%b[i];
			}
		}
	}
	
	printf("%d\n",quick_pow(g,CRT(),mod));
	return 0;
}