一分钟学会使用卢卡斯定理(仅限于应用层面,无数学证明)

在网上找了一大堆卢卡斯的描述,很难懂。
其实公式非常简短
C m n m o d p = C m / p n / p ∗ C m   m o d   p n   m o d   p m o d   p ( p 为 质 数 ) C_m^nmodp=C_{m/p}^{n/p}∗C_{m\ mod\ p}^{n \ mod \ p}mod \ p(p为质数) Cmnmodp=Cm/pn/pCm mod pn mod pmod pp其中我们对 C m / p n / p C_{m/p}^{n/p} Cm/pn/p执行递归操作
C m   m o d   p n   m o d   p C_{m\ mod\ p}^{n \ mod \ p} Cm mod pn mod p进行暴力求解
递归操作代码:

int Lucas(int n, int m, int p) {
     if (!m) return 1;
     return (LL)C(n%p, m%p)*Lucas(n/p, m/p, p)%p;
}

暴力求解——直接运用乘法逆元求解

ll inv(ll x)//ll->long long
{
	return qpow(x,mod-2);
}
void init()
{
	f[0] = 1;
	for(ll i=1;i<=200000;i++)
	{
 		f[i] = (f[i-1]*i)%p;
	}
}
ll C(int n,int m)
{
	return f[m]*inv(f[m-n+1])%p*inv(f[n-1])%p;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值