洛谷-2480 [SDOI2010]古代猪文

题目描述
猪王国的文明源远流长，博大精深。
iPig在大肥猪学校图书馆中查阅资料，得知远古时期猪文文字总个数为N。当然，一种语言如果字数很多，字典也相应会很大。当时的猪王国国王考虑到如果修一本字典，规模有可能远远超过康熙字典，花费的猪力、物力将难以估量。故考虑再三没有进行这一项劳猪伤财之举。当然，猪王国的文字后来随着历史变迁逐渐进行了简化，去掉了一些不常用的字。
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的余数就可以了。
输入格式
输入文件ancient.in有且仅有一行：两个数N、G，用一个空格分开。
输出格式
输出文件ancient.out有且仅有一行：一个数，表示答案除以999911659的余数。

输入输出样例
输入 #1
4 2

输出 #1
2048

说明/提示
数据规模
10%的数据中，1 <= N <= 50；
20%的数据中，1 <= N <= 1000；
40%的数据中，1 <= N <= 100000；
100%的数据中，1 <= G <= 1000000000，1 <= N <= 1000000000。

解释：求解$G^{{\sum }_{d|n}{C}_n^d}\%mod=G^{{\sum }_{d|n}{C}_n^d\%\varphi (mod)}\%mod$
=$G^{{\sum }_{d|n}{C}_n^d\%(mod-1)}\%mod$,又$ $mod-1$可以分解为$2,3,4679,35617$
$x=a_1\mathrm{m}\mathrm{o}\mathrm{d}2$
$x=a_2\mathrm{m}\mathrm{o}\mathrm{d}3$
$x=a_3\mathrm{m}\mathrm{o}\mathrm{d}4679$
$x=a_4\mathrm{m}\mathrm{o}\mathrm{d}35617$
通过中国剩余定理求出$x$，在计算$G^x$就是答案

#include<bits/stdc++.h>
#define N 55555
#define Mod 999911659LL
using namespace std;
long long a[N];
int d[N]={0};
long long ai[10],bi[10];
int t=0;
int p[10]={2,3,4679,35617};
long long pow(long long y,int z,int p){
    y%=p;long long ans=1;
    for(int i=z;i;i>>=1,y=y*y%p)if(i&1)ans=ans*y%p;
    return ans;
}
long long C(long long n,long long m,long long p){
    if(m>n)return 0;
    return ((a[n]*pow(a[m],p-2,p))%p*pow(a[n-m],p-2,p)%p);
}
long long Lucas(long long n,long long m,long long p){
    if(!m)return 1;
    return C(n%p,m%p,p)*Lucas(n/p,m/p,p)%p;
}
inline int read(){
    int f=1,x=0;char ch;
    do{ch=getchar();if(ch=='-')f=-1;}while(ch<'0'||ch>'9');
    do{x=x*10+ch-'0';ch=getchar();}while(ch>='0'&&ch<='9');
    return f*x;
}
long long n,g;
void init(long long p){
    a[0]=1;for(int i=1;i<=p;i++)a[i]=(a[i-1]*i)%p;
}
long long mul(long long a,long long b,long long mod){
    long long res=0;
    while(b>0){
        if(b&1) res=(res+a)%mod;
        a=(a+a)%mod;
        b>>=1;
    }
    return res;
}
long long exgcd(long long a,long long b,long long &x,long long &y){
    if(b==0){x=1;y=0;return a;}
    long long gcd=exgcd(b,a%b,x,y);
    long long tp=x;
    x=y; y=tp-a/b*y;
    return gcd;
}
long long excrt(int n){
    long long x,y;
    long long M=bi[1],ans=ai[1];
    for(int i=2;i<=n;i++){
        long long a=M,b=bi[i],c=(ai[i]-ans%b+b)%b;
        long long gcd=exgcd(a,b,x,y),bg=b/gcd;
        if(c%gcd!=0) return -1;
        x=mul(x,c/gcd,bg);
        ans+=x*M;
        M*=bg;
        ans=(ans%M+M)%M;
    }
    return (ans%M+M)%M;
}
int main(){
    n=read();g=read();
    for(int i=1;i*i<=n;i++){
       if(n%i==0){
            d[t++]=i;
            if(i*i!=n){
                d[t++]=n/i;
            }
       }
    }
    if(g%Mod==0){
        cout<<0<<endl;
        return 0;
    }
    for(int i=0;i<4;i++){
        long long sum=0;
        init(p[i]);
        for(int j=0;j<t;j++){
            sum+=Lucas(n,d[j],p[i]);
            sum%=p[i];
        }
        bi[i+1]=p[i];
        ai[i+1]=sum;
    }
    long long ret=excrt(4);
    ret=pow(g,ret,Mod);
    cout<<ret<<endl;
    return 0;
}