组合数

int fac[N],invfac[N];
void init(){
    fac[0]=1;
    for(int i=1;i<N;i++)fac[i]=fac[i-1]*i%MOD;
    invfac[N-1]=inv(fac[N-1]);
    for(int i=N-2;i>=0;i--)invfac[i]=invfac[i+1]*(i+1)%MOD;
}