标签:lucas int 定理 base mod ll jc
定义
证明
代码实现
//lucas定理
//1.qpow,阶乘(init)
//2.lucas
//3.组合排列C(nm)
#include<bits/stdc++.h>
#define ll long long
using namespace std ;
const int mod = 110119 ;
ll jc[mod+1],nc[mod+1] ;
ll qpow(ll x,int k)
{
ll base = 1 ;
while(k)
{
if(k&1) base = base*x%mod ;
k>>=1 ;
x = x*x%mod ;
}
return base ;
}
void init()
{
jc[0] = 1 , nc[0] = 1 ;
for(int i = 1 ; i < mod ; i++)
{
jc[i] = jc[i-1]*i%mod ;
nc[i] = qpow(jc[i],mod-2) ;//费马小定理
}
}
ll compose(ll n,ll m)
{
if(n<0 || m<0 || n<m) return 0 ;
if(n==m) return 1 ;
return jc[n]*nc[n-m]%mod*nc[m]%mod ;
}
ll lucas(ll n,ll m)
{
if(n<0 || m<0 || n<m) return 0 ;
ll ans = 1 ;
while(n||m)
{
ans = ans*compose(n%mod,m%mod)%mod ;
n/=mod ;
m/=mod ;
}
return ans ;
}
int main()
{
init() ;
ll n,m ;
while(~scanf("%lld%lld",&n,&m))
{
printf("%lld\n",lucas(n,m)) ;
}
return 0 ;
}
参考:(9条消息) 卢卡斯定理(十分钟带你看懂)_Combatting 的博客-CSDN博客_卢卡斯定理
标签:lucas,int,定理,base,mod,ll,jc 来源: https://blog.csdn.net/weixin_50553852/article/details/120751777
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。