bzoj2982: combination

Lucas定理模版题

T行，每行一个数，为C(n, m) mod 10007的答案。

 

#include<bits/stdc++.h>
#define rep(i,l,r) for(int i=l;i<=r;++i)
using namespace std;
typedef long long ll;
const ll p=10007;
ll n,m,T;
ll Pow(int x,int y){
    if(!y)return 1;
    ll t=Pow(x,y/2)%p; (t*=t)%=p;
    if(y&1)return t*x;else return t;
}
ll c(ll n,ll m){
    if(n<m)return 0;
    ll ans=1;
    rep(i,1,m){
        ans=ans*(((n-m+i)%p)*Pow(i,p-2)%p)%p;  
    }
    return ans;
}
ll lucas(ll n,ll m){
    if(!m)return 1;
    return lucas(n/p,m/p)*c(n%p,m%p)%p;
}
int main(){
    cin>>T;
    while(T--){
        cin>>n>>m;
        cout<<lucas(n,m)<<endl;
    }
}