整数分块、欧几里得、逆元、欧拉定理(2023-03-26一周小结)

一、整数分块——

       求从1到n的n/i的和,即求\sum_{i=1}^{n}\left \lfloor \frac{n}{i} \right \rfloor

代码段:

ll ans=0;
for(int l=1;l<=n;l=r+1)
{
    r=n/(n/l);
    ans+=(r-l+1)*(n/l);
}

二、欧几里得——

       用于求两个数的最大公约数,也称辗转相除法。

代码段:

ll gcd(ll a,ll b)
{
    if(!b) return a;
    return gcd(b,a%b);
}(有模取模)

       扩展欧几里得:用于求逆元,求同余方程,求ax+by=gcd(x,y)。

代码段:

ll exgcd(ll a,ll b,ll &x, ll &y)
{//a的逆元对b取模
    if(!b)
    {
        x=1;
        y=0;//特解情况
        return a;
    }
    ll other,r=exgcd(b,a%b,x,y);
    other =x;
    x=y;
    y=other-a/b*y;//通解
    return r;
}

三、逆元——

       求逆元可用扩展欧几里得,也可以用费马小定理,以下是费马小定理求解:

ll Inverse(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a;//奇数
        b>>=1;
         a=a*a;//偶数
    }
    return ans;
}(有模取模)

四、欧拉定理——a^{\varphi (x)}\equiv 1(mod p)

五、二分乘与二分幂

ll midmul(ll a,ll b,ll mod)
{//二分乘
    if(b==0) return 0;
    if(b&1) return (a+midmul((a<<1)%mod,b>>1)%mod)%mod;//奇数
    else return midmul((a<<1)%mod,b>>1)%mod;//偶数
}

ll midpow(ll a,ll b,ll mod)
{//二分幂
    if(b==0) return 1;
    if(b&1) return a*midpow(a*a%mod,b>>1,mod)%mod;//这里的乘法用二分乘更好,直接乘大数据容易出错
    else return return midpow(a*a%mod,b>>1,mod)%mod;
}

六、Lucas定理

       用于求组合数,根据Lucas定理有:C_{a}^{b} %mod=C_{a/mod}^{b/mod}   *  C_{a%mod}^{b%mod}  %mod.

代码段:

ll C(ll a,ll b, ll mod)
{//Ca_b=a!/(b!*(a-b)!)
    ll ans=1;
    ll sum=1;
    for(int i=0;i<b;i++)
    {
        sum=sum*(a-i)%mod;//a*(a-1)*(a-2)*(a-3)*……*(a-b+1)=a!/(a-b)!
        ans=ans*(i+1)%p;//b!
    }
    ans=Inverse(ans,mod-2,mod);
    return sum*ans%mod;
}

ll Lucas(ll a,ll b,ll mod)
{
    if(!b) return 1;
    return Lucas(a/mod,b/mod,mod)*C(a%mod,b%mod,mod)%mod;
}

不对的地方欢迎指正

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值