BZOJ 2982 combination

lucas定理裸题，表达式：C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p
由于p为素数
根据费马小定理：
已知(a, p) = 1，则 a^p-1 ≡ 1 (mod p),  所以 a*a^p-2 ≡ 1 (mod p)。
也就是 (m!(n-m)!)的逆元为 (m!(n-m)!)^p-2

#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;
const int maxn = 10005;
const ll p = 10007;
int T;
ll fac[maxn],n,m;
ll pow(ll a,ll b)
{
    ll ret=1;
    while(b){
        if(b&1) ret=ret*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ret;
}
void fact()
{
    fac[0]=1;
    for(int i=1;i<=p;i++) fac[i]=fac[i-1]*i%p;
}
ll lucas()
{
    ll ret=1;
    while(n&&m)
    {
        ll a=n%p,b=m%p;
        if(a<b) return 0;
        ret=(ret*fac[a]*pow(fac[b]*fac[a-b]%p,p-2))%p;
        n/=p;m/=p;
    }
    return ret;
}
int main()
{
    scanf("%d",&T);
    fact();
    while(T--)
    {
        scanf("%lld%lld",&n,&m);
        printf("%lld\n",lucas());
    }
    return 0;
}