lucas定理小结

最新推荐文章于 2020-10-22 13:22:58 发布

whai362

最新推荐文章于 2020-10-22 13:22:58 发布

阅读量502

点赞数

分类专栏：组合数学文章标签： algorithm 组合数学

本文链接：https://blog.csdn.net/whai362/article/details/46536137

版权

组合数学专栏收录该内容

31 篇文章 0 订阅

订阅专栏

lucas定理小结
如果：
n = a[k]*p^k + a[k-1]*p^(k-1) + ... + a[1]*p + a[0]
m = b[k]*p^k + b[k-1]*p^(k-1) + ... + b[1]*p + b[0]
则：
C(n,m) = pe(i=0~k,C(a[i],b[i]))%p 其中pe表示连乘符号。
代码实现：

/*lucas定理
  如果：
  n = a[k]*p^k + a[k-1]*p^(k-1) + ... + a[1]*p + a[0]
  m = b[k]*p^k + b[k-1]*p^(k-1) + ... + b[1]*p + b[0]
  则：
  C(n,m) = pe(i=0~k,C(a[i],b[i]))%p 其中pe表示连乘符号。
 */
LL lucas(int n,int m,int MOD){
	if(n<m) return 0;	//注意
	LL ret=1;
	while(n && m){
		LL a=n%MOD,b=m%MOD;
		if(a<b) return 0;
		ret=ret*C(a,b,MOD)%MOD;
		n/=MOD;
		m/=MOD;
	}
	return ret;
}