社区讨论

Lucas板子,p过大也没问题

P3807【模板】卢卡斯定理 / Lucas 定理参与者 1已保存回复 0

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
0 条
当前快照
1 份
快照标识符
@lyx4trqi
此快照首次捕获于
2024/07/22 23:18
2 年前
此快照最后确认于
2024/07/23 09:19
2 年前
查看原帖
CPP
/*卢卡斯定理*/
//背景:求C(n,m)%p,p是质数
//公式:结论1 -> C(p,x)%p=0;
//     结论2 -> [(1+x)^p]%p=1+x^p
//     结论3 -> 卢卡斯:C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p
//注意:这里p是质数,如果不是质数,移步扩展卢卡斯

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll ksm(ll a,ll b,ll p)
{
	ll ans=1,base=a;
	while(b)
	{
		if(b&1)ans*=base,ans%=p;
		base*=base,base%=p;b>>=1;
	}
	return ans;
}
ll getC(ll a, ll b,ll p)   
{
	ll ans = 1;
	for (ll i = 1, j = a; i <= b; i ++, j -- )
	{
		ans = ans * j % p;
		ans = ans * ksm(i, p - 2,p) % p;
	}
	return ans;
}
ll lucas(ll n,ll m,ll p)
{
	if(m==0)return 1;
	return lucas(n/p,m/p,p)*getC(n%p,m%p,p)%p;
}
void solve()
{
	ll n,m,p;
	cin>>n>>m>>p;
	cout<<lucas(n,m,p);
}


回复

0 条回复,欢迎继续交流。

正在加载回复...