数论-卢卡斯定理（lucas)与拓展卢卡斯定理 (exlucas)

最新推荐文章于 2021-11-19 12:57:47 发布

2码8改

最新推荐文章于 2021-11-19 12:57:47 发布

阅读量703

点赞数 2

分类专栏：数论文章标签：卢卡斯定理拓展卢卡斯定理大数取模

喜欢随便转，， https://blog.csdn.net/qq_40744093

本文链接：https://blog.csdn.net/qq_40744093/article/details/88294340

版权

数论专栏收录该内容

8 篇文章 0 订阅

订阅专栏

卢卡斯定理（lucas)

[用途]

求解 $C_{n}^{m}\% p$ ,其中m,n较大，p较小且为素数

[结论]

$C_{n}^{m}\equiv C_{n/p}^{m/p}C_{n\%p}^{m\%p}(mod \quad p)$

[证明]

假设： $\left\{\begin{array}{cc} n=sp+q\\ m=tp+r \end{array}\right.$
即证明 $C_{n}^{m}\equiv C_{s}^{t}C_{q}^{r}(mod \quad p)$
使用构造法：
$(1+x)^{n}\equiv(1+x)^{sp+q}\equiv[(1+x)^{p}]^{s}(1+x)^{q}(mod \quad p)$
根据二项式定理
$(1+x)^{p}=\sum\limits_{i=0}^p{C_{p}^{i}x^{i}}$
除了第一项 $C_{p}^{0}x^{0}（即1）$ 与最后一项 $C_{p}^{p}x^{p}（即x^{p})$ 不能整除p，中间项均能整除p,所以在 $\quad p)$ 意义下 $1+x)^{p}=1+x^{p}$
则 $(1+x)^{n}\equiv[(1+x)^{p}]^{s}(1+x)^{q}\equiv(1+x^{p})^{s}(1+x)^{q}(mod \quad p)$
即 $\sum\limits_{k=0}^n{C_{n}^{k}x^{k}}\equiv\sum\limits_{i=0}^s{C_{s}^{i}x^{pi}}\sum\limits_{j=0}^q{C_{q}^{j}x^{j}}(mod \quad p)$
我们在等式两边找 $x^m$ 项的系数
左:含 $x^m$ 项显然是 $C_{n}^{m}$
右:含 $x^m$ 项只能由 $x^{pi}x^{j}=x^{pi+j}$ 得到，即 $m = p i + j$ ，显然 $i = t, j = r$ ，此时对应系数恰好是 $C_{s}^{t}C_{q}^{r}$
所以 $C_{n}^{m}\equiv C_{s}^{t}C_{q}^{r}(mod \quad p)$ 成立
即 $C_{n}^{m}\equiv C_{n/p}^{m/p}C_{n\%p}^{m\%p}(mod \quad p)$

[代码]

#include<iostream>
using namespace std;
typedef long long LL;
const LL N=1e5+2;
LL a[N];
void init(LL p)
{
	a[1]=1;
	for(int i=2;i<=p;++i)a[i]=a[i-1]*i%p;
}
void exgcd(LL a,LL b,LL &x,LL &y)
{
	if(!b){
		x=1;
		y=0;
		return;
	}
	exgcd(b,a%b,y,x);
	y-=a/b*x;
}
LL ksm(LL x,LL n,LL mod)
{
	LL ans=1;
	while(n){
		if(n&1)ans=ans*x%mod;
		n>>=1;
		x=x*x%mod;
	}
	return ans;
}
LL C(LL n,LL m,LL p)
{
	if(n==m||m==0)return 1;
	if(n<m)return 0;
	if(m*2>n)m=n-m;						  /*C(n,m)=c(n,n-m)*/
	return a[n]*ksm(a[m]*a[n-m],p-2,p)%p; /*求(a[m]*a[n-m])在(mod p)意义下的乘法逆元*/
										  /*拓展欧几里得与费马小定理均可*/ 
	/*LL x,y;
	exgcd(a[m]*a[n-m],p,x,y);
	return (a[n]*x%p+p)%p;*/ 
}
LL lucas(LL n,LL m,LL p)
{
	if(!m)return 1;
	return lucas(n/p,m/p,p)*C(n%p,m%p,p)%p;
}
int main()
{
	ios::sync_with_stdio(false);
	LL T,n,m,p;
	cin>>T;
	while(T--){
		cin>>n>>m>>p;
		init(p);
		cout<<lucas(n+m,m,p)<<endl;
	}
	return 0;
}

拓展卢卡斯定理（exlucas)

[用途]

求解 $C_{n}^{m}\% P$ ,其中m,n较大，P较小且不一定为素数

[解析]

转化为中国剩余定理(crt）

根据算术基本定理

$P=M_{1}M_{2}...M_{k}=\prod\limits_{i=0}^{k}M_{i}$
其中 $M_{i}=p_{i}^{t_{i}},M_{i}M_{j}(i!=j)$ 两两互素
设 $X=C_{n}^{m}$
则 $A_{i}=C_{n}^{m}\%M_{i}$
所以转化为中国剩余定理
$X\ \equiv\left\{\begin{array}{cc} A_{1} \quad (mod \quad M_{1})\\ A_{2} \quad (mod \quad M_{2})\\ ...\\ A_{k} \quad (mod \quad M_{k}) \end{array}\right.$
$M_{i}$ 根据算术基本定理循环分解即可
重点的是 $A_{i}$
$\begin{aligned}A_{i} &=C_{n}^{m}\%M_{i}\\ &=\frac{n!}{m!(n-m)!}\%M_{i}\\ \end{aligned}$
可惜的是 $m!(n-m)!在模M_{i}(M_i=p_{i}^{t_{i}})$ 意义下的乘法逆元不一定存在
我们可以先把阶层中含 $p_{i}$ 的项约去，这样就可以求逆元了，然后再乘回来即可
所以
$\begin{aligned}A_{i} &=\frac{n!}{m!(n-m)!}\%M_{i}\\ &=\frac{\frac{n!}{p_i}}{\frac{m!}{p_i}\frac{(n-m)!}{p_i}}p_i^k\%M_{i}\\ &=\frac{n!}{p_i}inv(\frac{m!}{p_i},M_i)inv(\frac{(n-m)!}{p_i},M_i)p_i^k\%M_{i}\\ &=(\frac{n!}{p_i}\%M_{i})inv(\frac{m!}{p_i}\%M_{i},M_i)inv(\frac{(n-m)!}{p_i}\%M_{i},M_i)p_i^k\%M_{i} \end{aligned}$

令
$nn=\frac{n!}{p_i}\%M_{i}$
$mm=\frac{m!}{p_i}\%M_{i}$
$nm=\frac{(n-m)!}{p_i}\%M_{i}$
则 $A_i=nn*inv(mm,M_i)*inv(nm,M_i)*p_i^k\%M_i$
现在的关键又转换为求 $n n, m m, n m$ 了
就nn来看
当n=19时, $p_i=3,t_i=2$
$n!=1*...*19=3^6*6!*[(1*2*4*5*7*8)*(10*11*13*14*16*17)*19)]$
$\frac{n!}{p_i}\%M_{i}=[(1*2*4*5*7*8)*(10*11*13*14*16*17)*19)]*6！\%M_i$
而中括号这一堆是有循环节的 $1*2*4*5*7*8\equiv10*11*13*14*16*17(mod\quad M_i)$
循环节为 $n/M_i$ ，循环一组后然后快速幂就行了
剩下的个数就是 $n\%M_i$ 这一类直接循环相乘
而 $6!$ 恰好是 $n/p_i)!$ ,直接递归即可
即 $\frac{n!}{p_i}\%M_{i}=[...]*\frac{(n/p_i)!}{p_i}\%M_{i}$
$[. . .]$ 就是上面的循环节那块，分类进行求解,over！

[代码]

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL N=1e5+9;
LL A[N],M[N];
LL ksm(LL x,LL n,LL mod)
{
	LL ans=1;
	while(n){
		if(n&1)ans=ans*x%mod;
		n>>=1,x=x*x%mod;
	}
	return ans;
}
void exgcd(LL a,LL b,LL &x,LL &y)
{
	if(!b)x=1,y=0;
	else exgcd(b,a%b,y,x),y-=a/b*x;
}
LL inv(LL a,LL p)
{
	LL x,y;
	exgcd(a,p,x,y);
	return (x+p)%p?x:x+p;
}
LL get(LL n,LL pi,LL p)									/*求(与pi互素后的n!)%M[i]*/ 
{
	if(!n)return 1;
	LL ans=1;
	if(n/p){											/*判断有无循环节 */ 
		for(LL i=2;i<=p;++i)if(i%pi)ans=ans*i%p;
		ans=ksm(ans,n/p,p);
	}
	for(LL i=2;i<=n%p;++i)if(i%pi)ans=ans*i%p;			/*循环节剩余部分*/ 
	return ans*get(n/pi,pi,p)%p;
}
LL exlucas(LL n,LL m,LL pi,LL p)						/*求A[i]*/ 
{
	LL nn=get(n,pi,p);									/*求(与pi互素后的n)%M[i]*/ 
	LL mm=get(m,pi,p);									/*求(m!与pi互素后的m!)%M[i]*/ 
	LL nm=get(n-m,pi,p);								/*求(与pi互素后的(n-m)!)%M[i]*/ 
	LL k=0;												/*含质因数pi的数量*/ 
	for(LL i=n;i;i/=pi)k+=i/pi;
	for(LL i=m;i;i/=pi)k-=i/pi;
	for(LL i=n-m;i;i/=pi)k-=i/pi;
	return nn*inv(mm,p)*inv(nm,p)*ksm(pi,k,p)%p;
}
LL crt(LL len,LL Lcm)
{
	LL ans=0;
	for(LL i=1;i<=len;++i){
		LL Mi=Lcm/M[i];
		ans=((ans+A[i]*inv(Mi,M[i])*Mi)%Lcm+Lcm)%Lcm;
	}
	return ans;
}
int main()
{
	ios::sync_with_stdio(false);
	LL n,m,P,num;
	while(cin>>n>>m>>P){
		if(n<m){
			cout<<0<<endl;
			continue;
		}
		num=0;
		memset(A,0,sizeof(A));
		memset(M,0,sizeof(M));
		for(LL x=P,i=2;i<=P;++i)
			if(x%i==0){
				M[++num]=1;
				while(x%i==0){
					M[num]*=i;
					x/=i;
				}
				A[num]=exlucas(n,m,i,M[num])%P;
			} 
		cout<<crt(num,P)<<endl;
	}
	return 0;
}

2码8改

关注

2
点赞
踩
4

收藏

觉得还不错? 一键收藏
0
评论
数论-卢卡斯定理（lucas)与拓展卢卡斯定理 (exlucas)

卢卡斯定理（lukas)[用途]求解Cnm%pC_{n}^{m}\% pCnm%p,其中m,n较大，p较小且为素数[结论]Cnm≡Cn/pm/pCn%pm%p(modp)C_{n}^{m}\equiv C_{n/p}^{m/p}C_{n\%p}^{m\%p}(mod \quad p)Cnm≡Cn/pm/pCn%pm%p(modp)[证明]假设：{n=sp+qm=tp+r\lef...
复制链接

扫一扫