lucas定理

tanjunming2020

于 2023-02-14 19:50:11 发布

阅读量1.1k

点赞数 5

分类专栏：定理文章标签： c++

本文链接：https://blog.csdn.net/tanjunming2020/article/details/129032365

版权

定理专栏收录该内容

4 篇文章 0 订阅

订阅专栏

Lucas定理提供了一种在模质数p的情况下计算组合数C_n^m的方法，特别是当n和m非常大而质数p较小时。通过将n和m表示为p进制，可以将问题分解为更小的组合数计算，简化了大整数的处理。证明基于二项式定理和模运算性质，保证了计算的正确性。示例代码展示了如何实际应用该定理解决相关问题。

摘要由CSDN通过智能技术生成

介绍

当正整数 $n, m$ 很大，且质数 $p$ 较小的时候，要求 $C_n^m$ 对 $p$ 取模后的值，可以用lucas定理。

定理内容如下：

若 $p$ 为质数， $n, m$ 为正整数，且 $n\geq m$ ，设 $n = a p + b, m = c p + d$ ，其中 $0\leq b<p,0\leq d<p$ ，则有

$C_n^m\equiv C_a^b\times C_c^d \pmod p$

递归下去，我们可以发现，若 $n$ 的 $p$ 进制为 $(\overline{n_1,n_2\dots n_k})_p$ ， $m$ 的 $p$ 进制为 $(\overline{m_1,m_2\dots m_k})_p$ ，则有

$C_n^m\equiv C_{n_k}^{m_k}\times C_{n_{k-1}}^{m_{k-1}}\times \cdots\times C_{n_0}^{m_0} \pmod p$

这样就可以方便地解决大整数取模的问题了，不过前提是 $p$ 为较小的质数。

证明

根据二项式定理， $1+x)^n$ 中 $x^m$ 的系数为 $C_n^m$ 。

为了方便叙述，以下式子都是在模 $p$ 的意义下进行的。

$(1+x)^n\equiv (1+x)^{ap+b}\equiv ((1+x)^p)^a(1+x)^b\equiv (1+x^p)^a(1+x)^b$

最后一步成立是因为当 $0 < i < p$ 时， $C_p^i=\dfrac{p!}{i!(p-i)!}$ ， $p$ 为质数，原式的分子中可以提出一个 $p$ ，所以此时 $C_p^i$ 是 $p$ 的倍数，可以被模去。

继续推导

$(1+x^p)^a(1+x)^b\equiv \sum\limits_{i=1}^aC_a^ix^{pi}\sum\limits_{j=0}^bC_b^jx^j$

因为 $b < p$ ，所以 $j < p$ ， $m = c p + d$ ，考虑 $x^m$ 的系数， $x^m$ 的系数只能为 $C_a^c\times C_b^d$ ，得证 $C_n^m\equiv C_a^b\times C_c^d \pmod p$ 。

例题

洛谷P3807 【模板】卢卡斯定理/Lucas 定理

根据上述方法直接做即可。

code

#include<bits/stdc++.h>
using namespace std;
long long p;
long long jc[100005],ny[100005];
long long mi(long long t,long long v){
	if(v==0) return 1;
	long long re=mi(t,v/2);
	re=re*re%p;
	if(v&1) re=re*t%p;
	return re;
}
long long C(long long n,long long m){
	if(n<m) return 0;
	return jc[n]*ny[m]%p*ny[n-m]%p;
}
long long lucas(long long n,long long m){
	if(n<m) return 0;
	if(m==0) return 1;
	return C(n%p,m%p)*lucas(n/p,m/p)%p;
}
int main()
{
	int t;
	long long n,m,a,b;
	scanf("%d",&t);
	while(t--){
		scanf("%lld%lld%lld",&a,&b,&p);
		n=a+b;m=a;
		jc[0]=1;
		for(int i=1;i<p;i++) jc[i]=jc[i-1]*i%p;
		ny[p-1]=mi(jc[p-1],p-2);
		for(int i=p-2;i>=0;i--) ny[i]=ny[i+1]*(i+1)%p;
		printf("%lld\n",lucas(n,m));
	}
	return 0;
}