Lucas定理学习小记

（1）Lucas定理：p为素数，则有:

（2）证明： n=(ak...a2,a1,a0)p = (ak...a2,a1)p*p + a0 = [n/p]*p+a0，m=[m/p]*p+b0其次，我们知道，对任意质数p有(1+x)^p=1+(x^p)(mod p) 。我们只要证明这个式子：C(n,m)=C([n/p],[m/p]) * C(a0,b0)(mod p)，那么就可以用归纳法证明整个定理。对于模p而言，我们有下面的式子成立：

上式左右两边的x的某项x^m(m<=n)的系数对模p同余。其中左边的x^m的系数是 C(n,m)。而由于a0和b0都小于p，因此右边的x^m 一定是由 x^([m/p]*p) 和 x^b0 (即i=[m/p] , j=b0 ) 相乘而得因此有：C(n,m)=C([n/p],[m/p]) * C(a0,b0) (mod p)。

（3）拓展应用：上面的p是素数，那么不是素数怎么办呢？若不是素数，将p分解质因数，将C(n,m)分别按照（1）中的方法求对p的质因数的模，然后用中国剩余定理合并。比如计算C(10,3)%14。C(10,3)=120,14有两个质因数2和7，120%2=0,120%7=1,这样用(2,0)(7,1)找到最小的正整数8即是答案，即C(10,3)%14=8。注意，这里只适用于p分解完质因数后每个质因数只出现一次，例如12=2*2*3就不行，因为2出现了两次。若p分解完质因数后，含有某个质因数出现多次，比如C(10,3)%98,其中98=2*7*7,此时就要把7*7看做一个数,即:120%2=0,120%49=22,用(2,0)(49,22)和中国剩余定理得到答案22，即C(10,3)%98=22。此时，你又会有疑问，C(10,3)%49不也是模一个非素数吗？此时不同的是这个非素数不是一般的非素数，而是某个素数的某次方。下面（4）介绍如何计算C(n,m)%p^t(t>=2,p为素数)。

（4）计算C(n,m)%p^t。我们知道，C(n,m)=n!/m!/(n-m)!，若我们可以计算出n!%p^t，我们就能计算出m!%p^t以及(n-m)!%p^t。我们不妨设x=n!%p^t,y=m!%p^t,z=(n-m)!%p^t,那么答案就是x*reverse(y,p^t)*reverse(z,p^t)(reverse(a,b)计算a对b的乘法逆元)。那么下面问题就转化成如何计算n!%p^t。比如p=3,t=2,n=19,

n!=1*2*3*4*5*6*7*8* ……*19

=[1*2*4*5*7*8*… *16*17*19]*(3*6*9*12*15*18)

=[1*2*4*5*7*8*… *16*17*19]*3^6(1*2*3*4*5*6)

然后发现后面的是(n/p)!，于是递归即可。前半部分是以p^t为周期的[1*2*4*5*7*8]=[10*11*13*14*16*17](mod 9)。下面是孤立的19，可以知道孤立出来的长度不超过 p^t,于是暴力即可。那么最后剩下的3^6啊这些数怎么办呢？我们只要计算出n!,m!,(n-m)!里含有多少个p(不妨设a,b,c)，那么a-b-c就是C(n,m)中p的个数，直接算一下就行。

至此整个计算C(n,m)%p(p为任意数)的问题完美解决。下面给出代码：

i64 POW(i64 a,i64 b,i64 mod)
{i64 ans=1;while(b){if(b&1) ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans;
}i64 POW(i64 a,i64 b)
{i64 ans=1;while(b){if(b&1) ans=ans*a;a=a*a;b>>=1;}return ans;
}i64 exGcd(i64 a,i64 b,i64 &x,i64 &y)
{i64 t,d;if(!b){x=1;y=0;return a;}d=exGcd(b,a%b,x,y);t=x;x=y;y=t-a/b*y;return d;
}bool modular(i64 a[],i64 m[],i64 k)
{i64 d,t,c,x,y,i;for(i=2;i<=k;i++){d=exGcd(m[1],m[i],x,y);c=a[i]-a[1];if(c%d) return false;t=m[i]/d;x=(c/d*x%t+t)%t;a[1]=m[1]*x+a[1];m[1]=m[1]*m[i]/d;}return true;
}i64 reverse(i64 a,i64 b)
{i64 x,y;exGcd(a,b,x,y);return (x%b+b)%b;
}i64 C(i64 n,i64 m,i64 mod)
{if(m>n) return 0;i64 ans=1,i,a,b;for(i=1;i<=m;i++){a=(n+1-i)%mod;b=reverse(i%mod,mod);ans=ans*a%mod*b%mod;}return ans;
}i64 C1(i64 n,i64 m,i64 mod)
{if(m==0) return 1;return C(n%mod,m%mod,mod)*C1(n/mod,m/mod,mod)%mod;
}i64 cal(i64 n,i64 p,i64 t)
{if(!n) return 1;i64 x=POW(p,t),i,y=n/x,temp=1;for(i=1;i<=x;i++) if(i%p) temp=temp*i%x;i64 ans=POW(temp,y,x);for(i=y*x+1;i<=n;i++) if(i%p) ans=ans*i%x;return ans*cal(n/p,p,t)%x;
}i64 C2(i64 n,i64 m,i64 p,i64 t)
{i64 x=POW(p,t);i64 a,b,c,ap=0,bp=0,cp=0,temp;for(temp=n;temp;temp/=p) ap+=temp/p;for(temp=m;temp;temp/=p) bp+=temp/p;for(temp=n-m;temp;temp/=p) cp+=temp/p;ap=ap-bp-cp;i64 ans=POW(p,ap,x);a=cal(n,p,t);b=cal(m,p,t);c=cal(n-m,p,t);ans=ans*a%x*reverse(b,x)%x*reverse(c,x)%x;return ans;
}//计算C(n,m)%mod
i64 Lucas(i64 n,i64 m,i64 mod)
{i64 i,t,cnt=0;i64 A[205],M[205];for(i=2;i*i<=mod;i++) if(mod%i==0){t=0;while(mod%i==0){t++;mod/=i;}M[++cnt]=POW(i,t);if(t==1) A[cnt]=C1(n,m,i);else A[cnt]=C2(n,m,i,t);}if(mod>1){M[++cnt]=mod;A[cnt]=C1(n,m,mod);}modular(A,M,cnt);return A[1];
}

ps：这篇文章的代码好像是有问题的，在做这道题时：http://blog.csdn.net/yuyanggo/article/details/47378767，就出现问题了，但他讲的算法是正确的，代码的问题暂时我也不知道。如果要看代码的话，可以参考我在刚刚那篇文章中所写的代码。

lucas定理及扩展lucas定理相关推荐

孙子定理及扩展孙子定理
孙子定理及扩展孙子定理中国剩余定理 (Chinese Remainder Theorem) CRT是解决下列同余方程问题的,mi,mjm_i,m_jmi,mj两两互质. 该定理只需要用到两个取余 ...
Lucas定理与扩展Lucas
之前看了乘法逆元(详见除法取模与逆元),发现不能处理不互质的情况,于是去找方法,最后找到了Lucas定理... 虽然与期待中的不一样,但是还是非常有用的. (1)Lucas定理: 若p为素数,则有: ...
专题·Lucas定理【including Lucas定理，扩展Lucas
初见安~这里是数论专题(6)[详见数论专栏本篇有前置知识点需要掌握,建议先了解下:费马小定理,中国剩余定理,乘法逆元一.Lucas定理 Lucas定理用于求解的组合数取模的问题.其中p为质数. 组 ...
组合数学 —— 组合数取模 —— 卢卡斯定理与扩展卢卡斯定理
[卢卡斯定理] 1.要求:p 是质数,m.n 很大但 p 很小或者 n.m 不大但大于 p 2.定理内容其中, 3.推论当将 n 写成 p 进制:,将 m 写成 p 进制: 时,有: 4.实现 ...
2015 ICL, Finals, Div. 1 Ceizenpok’s formula（组合数取模，扩展lucas定理）
J. Ceizenpok's formula time limit per test 2 seconds memory limit per test 256 megabytes input stand ...
[学习笔记]扩展LUCAS定理
可以先做这个题[SDOI2010]古代猪文此算法和LUCAS定理没有半毛钱关系. [模板]扩展卢卡斯不保证P是质数. $C_n^m=\frac{n!}{m!(n-m)!}$ 麻烦的是分母. 如果互 ...
Lucas定理扩展Lucas
Lucas定理求解 C n m m o d p C_n^m\mod p Cnmmodp (p为素数) 蒟蒻不会证明,记着递推公式就好=_= L u c a s ( n , m ) = C ( ...
BZOJ 4338 糖果（扩展Lucas定理+CRT）
题目链接:BZOJ 4338 题目大意:用数字1~k填一个n*m的表格,每种数字可用任意次,要求每行数字1~m列单调不减,任意两行不完全相同,求方案数对P取模的值. 题解:扩展Lucas+CRT模板题 ...
扩展Lucas定理扩展Lucas板子
题意概述:多组询问,给出N,K,M,要求回答C(N,K)%M,1<=N<=10^18,1<=K<=N,2<=M<=10^6 分析: 模数不为质数只能用扩展Lucas ...

lucas定理及扩展lucas定理

Lucas定理学习小记

lucas定理及扩展lucas定理相关推荐

最新文章

热门文章