码迷,mamicode.com
首页 > 其他好文 > 详细

# 数论-组合数+lacus定理

时间:2019-07-26 10:42:43      阅读:87      评论:0      收藏:0      [点我收藏+]

标签:质数   math   lis   turn   费马小定理   ==   lin   合数   范围   

数论-组合数+lacus定理

组合数计算

  • 为避免爆long long,\(20!\)就达到了long long 的范围,采用边乘边除的思想
ll C(ll n,ll m){
    if(n<m)return 0;
    ll ans=1;
    for(ll i=1;i<=m;i++)
        ans=ans*(n-m+i)/i;   
    return ans;
}

lacus定理-大组合数取模

  • 卢卡斯定理:C(n,m)%p=C(n/p,m/p)*C(n%p,m%p)%p,适合p的范围在10^5且p为素数的情况

  • 费马小定理:当p为质数时候, a^(p-1)≡1(mod p), a*a^(p-2)≡1(mod p)

ll qpow(ll a,ll b){    
    ll ans=1;
    a%=p;
    while(b){
        if(b&1)ans=ans*a%p;
        b>>=1;
        a=a*a%p
    }
    return ans;
}
ll C(ll n,ll m){
    if(n<m)return 0;
    ll ans=1;
    for(ll i=1;i<=m;i++)
        ans=ans*((n-m+i)%p*qpow(i,p-2)%p)%p;//费马小定理求逆元 
    return ans;
}
 
ll lacus(ll n,ll m){
    if(m==0)
        return 1;
    return (lacus(n/p,m/p)*C(n%p,m%p))%p;
}

# 数论-组合数+lacus定理

标签:质数   math   lis   turn   费马小定理   ==   lin   合数   范围   

原文地址:https://www.cnblogs.com/sstealer/p/11248668.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!