BZOJ 1951 古代猪文 鲁卡斯定理+费马小定理+中国剩余定理

先给这道题的题目描述跪了（此处吐槽10分钟）
——————————————————————————
吐槽完毕，这道题首先看得出要算的结果是一个数的巨大无比次方，简单的快速幂显然不行了，但是最后的结果要mod一个质数，于是我们可以引入费马小定理：
a^(p-1)≡1（mod p）（当p为质数并且a不为p的倍数时成立），这能说明什么呢，就是说在mod p的意义下乘以一个数的p-1次方等于没乘！所以我们就可以将一个数巨大无比的次方mod（p-1）就行了。然而我们又发现了一个问题，这个次方是用组合数算出来的，由于组合数的递推公式无法使用（n，m过大），所以只能使用组合数的阶乘公式，然而就算我们用阶乘公式也要将阶乘从1一直乘到p-1，复杂度爆炸，所以我们又想到了另一个牛*的定理：中国剩余定理。
我们发现p-1恰好可以分解为35617*4679*3*2，那么我们只需算出这个次方分别在取余这4个质数之下的余数用中国剩余定理一合就可以了，中间还要穿插一下阶乘求逆（其实很简单）和鲁卡斯定理，鲁卡斯定理的内容就是从n个数中选出m个数的方案mod p==从n/p个数中选出m/p方案数*从n mod p个数中选出m mod p个数的方案，这样这道题就解决了。
鲁卡斯定理有毒…

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<iostream>
#include<iomanip>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
const long long maxn=999911659;
long long ni[4][40000];
long long jc[4][40000];
long long modd[4];
long long x[4];
long long Mi[4];
long long M;
void jisuan(long long num)
{
    long long mod=modd[num];
    long long mid=1;
    jc[num][0]=1;
    jc[num][1]=1;
    for(int i=2;i<mod;i++)
    {
        mid=(mid*i)%mod;
        jc[num][i]=(jc[num][i-1]*i)%mod;
    }
    ni[num][mod-1]=mod-1;
    for(int i=mod-2;i>=0;i--) ni[num][i]=(ni[num][i+1]*(i+1))%mod;
}
long long zuhe(long long p,long long q,long long o)
{
    long long sum=1;
    while(p && q)
    {
        long long x=p%modd[o],y=q%modd[o];
        if(y>x) return 0;
        sum=(sum*(jc[o][x]*(ni[o][x-y]*ni[o][y])%modd[o])%modd[o])%modd[o];
        p/=modd[o];
        q/=modd[o];
    }
    return sum;
}
long long qiuni(long long yuan,long long p)
{
    long long t=p-2;
    long long mid=yuan;
    long long ans=1;
    while(t)
    {
        if(t&1) ans=(ans*mid)%p;
        mid=(mid*mid)%p;
        t>>=1;
    }
    return ans;
}
long long china()
{
    long long t=0;
    for(int i=0;i<=3;i++)
    {
        long long mid=(x[i]*Mi[i])%M;
        mid=(mid*qiuni(Mi[i],modd[i]))%M;
        t=(t+mid)%M;
    }
    return t;
}
int main()
{
    long long n,g;
    scanf("%lld%lld",&n,&g);
    if(g%maxn==0) {cout<<0;return 0;}
    modd[0]=35617;
    modd[1]=4679;
    modd[2]=3;
    modd[3]=2;
    jisuan(0);jisuan(1);jisuan(2);jisuan(3);
    long long gen=sqrt(n);
    for(int i=1;i<=gen;i++)
        if(n%i==0)
        {
            for(int j=0;j<=3;j++)
            {
                x[j]=(x[j]+zuhe(n,i,j))%modd[j];
                if(i*i<n) x[j]=(x[j]+zuhe(n,n/i,j))%modd[j];
            }
        }
    M=modd[0]*modd[1]*modd[2]*modd[3];
    for(int i=0;i<=3;i++) Mi[i]=M/modd[i];
    long long goal=china();
    long long t=goal;
    long long ans=1;
    long long mid=g;
    while(t)
    {
        if(t&1) ans=(ans*mid)%maxn;
        mid=(mid*mid)%maxn;
        t>>=1;
    }
    cout<<ans<<endl;
    return 0;
}