lucas定理概述

继续学习数论：lucas定理。
我也不会证，记个结论吧：$C(a,b)=C(a$%$p$,$b$%$p)*C(a$/$p$,$b$/$p)$
然后我们递归就好了。
当然要求是$p$它是一个质数。如果$p$不是一个质数怎么办咧？
那我们把$p$分解质因数一下就可以使用lucas定理了。

哦，等一下，那我们怎么合并咧？我们知道了答案模各种质数的值，是不是很像中国剩余定理？我们得到了一系列同余方程，然后把它们解出来即可。

还是做一道题吧：1951: [Sdoi2010]古代猪文
首先我们使用欧拉定理，只用求出答案模上999911658就行。对于$n$暴力分解。然而这玩意儿不是一个质数，我们用上面的方法，可以求答案。注意特判 G G <script type="math/tex" id="MathJax-Element-1618">G</script>等于999911659的情况，这样不满足欧拉定理（真坑）。
code:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define md 999911659
using namespace std;

long long n,g,s,x,y;
long long ans;
long long mod[5]={0,2,3,4679,35617};
long long inv[5][35618],rv[5][35618];

long long qui(long long a,long long b,long long p){
    long long ret=1;
    while(b!=0){
        if((b&1)==1) ret=ret*a%p;
        a=a*a%p;
        b=b>>1;
    }
    return ret;
}

long long exgcd(long long a,long long b){
    if(b==0){
        x=1,y=0;
        return a;
    }
    long long ret=exgcd(b,a%b);
    long long tep=y;
    y=x-a/b*y;
    x=tep;
    return ret;
}

long long C(long long a,long long b,int ty){
    if(a==0) return 1;
    return rv[ty][b]*inv[ty][b-a]%mod[ty]*inv[ty][a]%mod[ty];
}

long long lucas(long long a,long long b,int ty){
    if(a<mod[ty]&&b<mod[ty]){
        if(a>b) return 0;
        return C(a,b,ty);
    }
    return lucas(a%mod[ty],b%mod[ty],ty)*lucas(a/mod[ty],b/mod[ty],ty)%mod[ty];
}

long long add(long long u){
    long long m1=mod[1],r1=lucas(u,n,1),m2,r2;
    for(int i=2;i<=4;i++){
        m2=mod[i],r2=lucas(u,n,i);
        long long d=exgcd(m1,-m2);
        long long t=abs((-m2)/d);
        x=x*(r2-r1)/d;
        x=(x%t+t)%t;
        t=x*m1+r1;
        m1=m1*m2/exgcd(m1,m2);r1=t;
    }
    return (r1%m1+m1)%m1;
}

int main()
{
    for(int i=1;i<=4;i++){
      rv[i][1]=1;
      for(int j=2;j<=mod[i]-1;j++){
        rv[i][j]=rv[i][j-1]*j%mod[i];
      }
      inv[i][mod[i]-1]=qui(mod[i]-1,mod[i]-2,mod[i]);
      for(int j=mod[i]-2;j>=1;j--)
        inv[i][j]=(inv[i][j+1]*(j+1))%mod[i];
      inv[i][0]=1;
    }
    scanf("%lld%lld",&n,&g);
    if(g==md){
      printf("0");
      return 0;
    }
    s=sqrt(n);
    for(int i=1;i<=s;i++)
      if(n%i==0){
        ans=(ans+add(i))%(md-1);
        if(i<n/i) ans=(ans+add(n/i))%(md-1);
      }
    printf("%lld",qui(g,ans,md));
}