【BZOJ】2982 combination

【算法】组合数取模——lucas定理

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int p=;
int fac[p+],inv[p+];
void gcd(int a,int b,int &x,int &y)
{
    if(b==){x=;y=;}
    else{gcd(b,a%b,y,x);y-=x*(a/b);}
}
void pre_inv()
{
    fac[]=fac[]=;
    for(int i=;i<p;i++)fac[i]=(fac[i-]*i)%p;
    int xx,yy;
    gcd(fac[p-],p,xx,yy);
    inv[p-]=((xx%p)+p)%p;//扩欧解不一定是最小非负解！
    for(int i=p-;i>=;i--)inv[i]=(inv[i+]*(i+))%p;
}
int C(int n,int m)
{
    if(n<m)return ;
    if(n<p&&m<p)return (1ll*fac[n]*inv[m]*inv[n-m])%p;
    return (C(n%p,m%p)*C(n/p,m/p))%p;
}
int main()
{
    pre_inv();
    int T,n,m;
    scanf("%d",&T);
    //for(int i=0;i<10;i++)printf("inv_%d=%d\n",i,inv[i]);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        printf("%d\n",C(n,m));
    }
    return ;
}