【BZOJ】【2982】Combination

排列组合

　　Lucas定理模板题……

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

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

　　预处理一下fact和inv即可

/**************************************************************
    Problem: 2982
    User: Tunix
    Language: C++
    Result: Accepted
    Time:4 ms
    Memory:1352 kb
****************************************************************/
 
//BZOJ 2982
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
using namespace std;
 
int getint(){
    int v=0,sign=1; char ch=getchar();
    while(ch<'0'||ch>'9') {if (ch=='-') sign=-1; ch=getchar();}
    while(ch>='0'&&ch<='9') {v=v*10+ch-'0'; ch=getchar();}
    return v*=sign;
}
/*******************tamplate********************/
const int N=10086,P=10007;
int fac[N],inv[N],n,m;
int pow(int a,int b){
    int r=1;
    for(;b;b>>=1,a=a*a%P)
        if (b&1) r=r*a%P;
    return r;
}
void calc(){
    fac[0]=1;
    F(i,1,P-1) fac[i]=fac[i-1]*i%P;
    inv[P-1]=pow(fac[P-1],P-2); inv[0]=1;
    D(i,P-2,1) inv[i]=inv[i+1]*(i+1)%P;
}
inline int c(int n,int m){
    if (n<m) return 0;
    if (n<P && m<P) return fac[n]*inv[m]%P*inv[n-m]%P;
    return c(n%P,m%P)*c(n/P,m/P)%P;
}
int main(){
    int T=getint();
    calc();
    while(T--){
        n=getint(); m=getint();
        printf("%d
",c(n,m));
    }
    return 0;
}