组合数学lucas定理 BZOJ2982 combination

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 using namespace std;
 6 const int mod=10007;
 7 int t,n,m;
 8 int qpow(int a,int b){
 9     int ans;
10     for(ans=1;b;b>>=1,a=a*a%mod)
11         if(b&1) ans=ans*a%mod;
12     return ans;
13 }
14 int zuhe(int nn,int mm){
15     if(nn<mm) return 0;
16     if(mm>nn-mm) mm=nn-mm;
17     long long s1=1,s2=1;
18     for(int i=0;i<mm;i++){
19         s1=s1*(nn-i)%mod;
20         s2=s2*(i+1)%mod;
21     }
22     return s1*qpow(s2,mod-2)%mod;//这里在int范围内 
23 }
24 int lucas(int nn,int mm){
25     if(!mm) return 1;
26     return zuhe(nn%mod,mm%mod)*lucas(nn/mod,mm/mod)%mod;
27 }
28 int main(){
29     scanf("%d",&t);
30     while(t--){
31         scanf("%d%d",&n,&m);
32         printf("%d\n",lucas(n,m));
33     } 
34     return 0;
35 }