bzoj2982 -- Lucas定理

Lucas定理裸题。。

Lucas定理:C(n,m)=C(n%p,m%p)*C(n/p,m/p)%p

预处理出阶乘、逆元的阶乘就可以了。

 

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define p 10007
int i,j,k,n,m,inv[p],fac[p],fac2[p],t;
inline void Init(){
    for(fac[0]=fac[1]=1,i=2;i<p;i++)fac[i]=fac[i-1]*i%p;
    for(inv[0]=inv[1]=fac2[0]=fac2[1]=1,i=2;i<p;i++)inv[i]=inv[p%i]*(p-p/i)%p,fac2[i]=fac2[i-1]*inv[i]%p;
}
inline int C(int n,int m){
    if(n<m)return 0;
    return fac[n]*fac2[m]%p*fac2[n-m]%p;
}
inline int Get(int n,int m){
    if(m==0)return 1;
    return C(n%p,m%p)*Get(n/p,m/p)%p;
}
int main(){
    scanf("%d",&t);
    Init();
    while(t--){
        scanf("%d%d",&n,&m);
        printf("%d\n",Get(n,m));
    }
    return 0;
}