排列组合

　　Lucas定理模板题……

　　感觉我做题顺序有点问题啊……应该是BZOJ 2982-->HDOJ 3037-->BZOJ 1272

　　好吧这个现在来看就有些水了……

　　预处理一下fact和inv即可

 1 /**************************************************************
 2     Problem: 2982
 3     User: Tunix
 4     Language: C++
 5     Result: Accepted
 6     Time:4 ms
 7     Memory:1352 kb
 8 ****************************************************************/
 9
10 //BZOJ 2982
11 #include<cstdio>
12 #include<cstdlib>
13 #include<cstring>
14 #include<iostream>
15 #include<algorithm>
16 #define rep(i,n) for(int i=0;i<n;++i)
17 #define F(i,j,n) for(int i=j;i<=n;++i)
18 #define D(i,j,n) for(int i=j;i>=n;--i)
19 using namespace std;
20
21 int getint(){
22     int v=0,sign=1; char ch=getchar();
23     while(ch<‘0‘||ch>‘9‘) {if (ch==‘-‘) sign=-1; ch=getchar();}
24     while(ch>=‘0‘&&ch<=‘9‘) {v=v*10+ch-‘0‘; ch=getchar();}
25     return v*=sign;
26 }
27 /*******************tamplate********************/
28 const int N=10086,P=10007;
29 int fac[N],inv[N],n,m;
30 int pow(int a,int b){
31     int r=1;
32     for(;b;b>>=1,a=a*a%P)
33         if (b&1) r=r*a%P;
34     return r;
35 }
36 void calc(){
37     fac[0]=1;
38     F(i,1,P-1) fac[i]=fac[i-1]*i%P;
39     inv[P-1]=pow(fac[P-1],P-2); inv[0]=1;
40     D(i,P-2,1) inv[i]=inv[i+1]*(i+1)%P;
41 }
42 inline int c(int n,int m){
43     if (n<m) return 0;
44     if (n<P && m<P) return fac[n]*inv[m]%P*inv[n-m]%P;
45     return c(n%P,m%P)*c(n/P,m/P)%P;
46 }
47 int main(){
48     int T=getint();
49     calc();
50     while(T--){
51         n=getint(); m=getint();
52         printf("%d\n",c(n,m));
53     }
54     return 0;
55 }