组合数取模

原文链接：http://blog.csdn.net/wang907553141/article/details/52143128

原文参考博客：http://blog.csdn.net/acdreamers/article/details/8037918

下面附上Lucas定理的代码：

const LL mod=1e9+7;typedef long long LL;LL n,m;LL quickPow(LL a,LL k){//快速幂求a^k%mod    LL ans=1;    while(k)    {        if(k&1)//相与运算,如果k是奇数那么结果为真,执行下面操作        ans=(ans*a)%mod;        a=(a*a)%mod;        k>>=1;    }    return ans;}LL C(LL a,LL b,LL p){//求从n个数中挑出m个数的组合数    if(a<b)      return 0;    if(b>a-b)      b=a-b;    LL up=1,down=1;    for(LL i=0;i<b;i++)    {        up=up*(a-i)%p;        down=down*(i+1)%p;    }    return up*quickPow(down,p-2)%p;//乘法逆元(费马小定理)}LL lucas(LL a,LL b,LL p){    if(b==0)        return 1;    return C(a%p,b%p,p)*lucas(a/p,b/p,p)%p;}