【数论-Lucas定理】

旧城冷巷雨未停 提交于 2019-12-07 12:41:56

1.写在前面:我始终觉得,对于一个问题要知其然,更要知其所以然。Lucas定理在刚刚接触数论的时候就知道了,因为这是一个很常用的定理,常常和中国剩余定理放在一起考。最近在组合数学上出现了很多问题,但是都是找个结论就过去了。浑浑噩噩并不懂其中原理,感觉自己的数学直觉一直在下降,以前我甚至能够从数据中看出来数学表达式,后来学了很多的工具,就变懒了,甚至连一个稍微复杂一点的积分都懒得动手,直接自适应Simpson积过去。今天又拿起了笔,准备把Lucas定理自己证明一遍。证完之后,我似乎又看到之前的那个影子:不畏惧任何数学表达式,即使推导了十几页也毫不畏惧。废话到此结束,正题开始。

2.Lucass定理:

      

Lucas定理是同余理论中的一个很重要的定理,用于组合数取模。常常使用在问题的结果需要对某个数据取模,n,m很大,达到1e15以上,但是p在1e9以内。一般来说最好的效果实在1e5以内。看到这个式子,一个十分典型的递归,并且注意到p是素数。所以我们就会有很多疑问:这个定理怎么来的?p为什么一定要是素数?p不是素数可以吗?我们来解决这些问题。

3.在解决这个问题之前,我想说一个问题:记得是ZOJ的某次月赛吧,问题就是求:

.  其实说白了,就是统计杨辉三角形每一行中有多少个奇数。这个问题当时我打了个表,看出了规律。结论和i和n的二进制有关,但是为什么呢?

4.首先证明Lucas定理。

                  Lucas定理在数论中表述是这样的:

                     

解释一下:其实就是把a,b在p进制下表示出来,然后发现C(a,b)是和后边的式子同余的。这个式子写成递归就是上边的递归形式下的Lucas定理。递归形式下的更加易于编程,但是由于递归的存在,p不能过大。问题来了,这个结论怎么来的?继续证明。

----多项式同余-----

首先介绍一个辅助工具:多项式下的同余,其实也很容易:

         

.   其中对于任意的

都有:ai和bi对同余,那么就说这两个多项式是同余的。

知道了这些,我们还需要明白这样一个结论:

%n=0,对于所有的

成立的条件是:n是质数。这个很容易证明,就不再证了。其实这里也说明了,为什么Lucas定理要求p一定是质数,而且,从这里我们也可以找出来p不是质数该如何处理。这个我们后边再讲。这样我们来看一个性质:

     

。这个证明可以用上边的结论证;

  

.  我们知道中间的所有项mod p是0,所以只有第一项和最后一项。那么上式得证。

那么我们对于一般的式子:

  

我们把a拆分成p进制:

                 

对于这个等式,在等式的左边,我们考虑x的b次方的系数:C(a,b)那么对于右边,其实就是在每个乘积式子中拿出一些凑出b。前边公式我们知道:

.也就是b在p进制下的拆分,所以对于每一个乘积式子,拿出bi。由乘法原理:C(ai,bi)的乘积即使左边x的b次方的系数。根据多项式相等的原则,我们推出:

                         

到了这里,Lucas定理得证。我们回答了第一个问题,Lucas定理是怎么来的。其实也回答了第二个问题,为什么p要是素数?因为只有p是素数才满足:对于所有的j<=p-1是成立的。那么p不是素数行吗?这个问题我们最后回答,先回答第四个问题:二项式系数:C(n,i)里面对于每一个n到底有多少个奇数?

5.二项式系数奇偶以及更加一般的C(a,b)%p=0的问题。在二项式系数中,判断是不是奇偶,其实是就是在判断:C(n,i)%2==0?.那么根据Lucas定理,这里的素数就是2.我们考虑这个式子:把n,i的二进制数位分别写出来,如果存在ai<bi,这个组合数就是0.显然是个偶数。所以我们猜想充分必要条件就是对于每一个数位,ai必须要大于bi.现在我们正面考虑:对于每一个ai>=bi,相等的时候是显然成立的。大于的时候,这个地方就是C(1,0)=1,乘起来必然是一个奇数。所以猜想得证:C(n,i)%2=1的充分必要条件是:对于i和n的二进制数位,n的每一位必须大于等于i的每一位。然后考虑每一层二项式奇数有多少?假设我用f(x)表示x的二进制数位中1的个数,对于所有的i属于[0,n],C(n,i)%2是奇数的个数就是:pow(2,f(n-1))。为什么?上边我已经证明了80%,接下来可以自己证明。那么这个结论是不是可以推广呢?答案是可以的,因为我们上来其实就说明了一个很重要的观点:a,b需要在p进制下观察,才能更好的认识这个问题。在我们证明二进制下的结论的时候,并没有很特殊的假设,只是建立在2是个质数这个结论上,所以对于所有的质数,都是成立的。那么我们可以得出下边的推论:

对于C(a,b)%p=0,p是质数。这个式子成立的条件是,对于a,b在p进制下的数位:ai,bi.只要存在一对数位使得ai<bi即可。

6.p不是质数的处理:p不是质数,就显得很棘手了。但是我们没有办法修改Lucas定理,所以只能把P变成质数,通过质因数分解

                                

那么列出同余方程:

                                  

得到了这个,我们就可以用中国剩余定理来合并答案。问题转化为如何求:这个式子的值。

我们把组合数展开:

                                  

尝试着把n!中所有pi的倍数分离,然后神奇的性质就出现了。发现n!%pi是有循环节的。这样,我们一边递归n!,一边快速幂即可求出答案。

最后放几个链接:

HDU3037:这是个裸题,直接就是结论了,C(n+m,m)%p.采用Lucas定理即可。

二项式系数奇偶问题

洛谷Lucas定理模板

最后,给出HDU3037我AC的代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll  n, m, p;
ll Ext_gcd(ll a, ll b, ll &x, ll &y) 
{
	if (b == 0) { x = 1, y = 0; return a; }
	ll ret = Ext_gcd(b, a%b, y, x);
	y -= a / b * x;
	return ret;
}
ll Inv(ll a, int m) 
{   
	ll d, x, y, t = (ll)m;
	d = Ext_gcd(a, t, x, y);
	if (d == 1) return (x%t + t) % t;
	return -1;
}

ll Cm(ll n, ll m, ll p)  
{
	ll a = 1, b = 1;
	if (m > n) return 0;
	while (m)
	{
		a = (a*n) % p;
		b = (b*m) % p;
		m--;
		n--;
	}
	return (ll)a*Inv(b, p) % p;  
}

int Lucas(ll n, ll m, ll p)  
{
	if (m == 0) return 1;
	return (ll)Cm(n%p, m%p, p)*(ll)Lucas(n / p, m / p, p) % p;
}

int main()
{
	int  T;
	cin >> T;
	while (T--)
	{
		scanf("%lld%lld%lld", &n, &m, &p);
		printf("%d\n", Lucas(n + m, m, p));
	}
	return 0;
}

参考书籍:冯志刚《初等数论》

 

 

 

 

 

 

 

易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!