Forsaken的位运算魔法类欧几里得

最新推荐文章于 2020-06-08 09:56:05 发布

henucm

最新推荐文章于 2020-06-08 09:56:05 发布

阅读量368

点赞数

分类专栏：牛客网

本文链接：https://blog.csdn.net/henucm/article/details/103150490

版权

牛客网专栏收录该内容

26 篇文章 0 订阅

订阅专栏

该博客探讨了如何从位运算的角度出发，利用类欧几里得算法解决特定问题。博主详细解释了如何遍历数字的每一位，通过判断某位是否为0来统计相应数量，并指出这种问题可以转化为类欧几里得算法的形式，从而提供了一种解决方案。

摘要由CSDN通过智能技术生成

传送

从贡献角度考虑。遍历 k 的每一位，若这一位为 0，计算有多少个 ij 这一位为1。反之计算这一位有多少个数字为 0。

计算有多少个 ij 第 p 位为 1，等价于计算我们固定i，会发现这个是类欧几里得算法形式，套个类欧的模板即可。原博主

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int mod=1e9+7;
const int inv2=500000004;
ll f(ll a, ll b, ll c, ll n) {
    if(a == 0)
        return (n+1) * (b/c)%mod;
    if(a < c && b < c) {
        ll m = (a*n+b)/c;
        if(m == 0)
            return 0;
        return (n*(m%mod)%mod - f(c, c-b-1, a, m-1)+mod)%mod;
    }
    return (f(a%c, b%c, c, n) + (n+1)*(b/c)%mod + (n+1)*n%mod*inv2%mod*(a/c)%mod)%mod;
}
int main()
{
    ll n,k;
    cin>>n>>k;
    ll ans=0;
    for(int j=0;j<=35;j++)
    {
        int flag=(k>>j)&1;
        ll s=1ll<<j;
        for(int i=0;i<=n;i++) //枚举i
        {
            ll x=f(i,0,s,n);
            ll y=f(i,0,s+s,n)*2%mod;
            x=((x-y)%mod+mod)%mod; //1的个数
            if(flag)
            {
                ans=(ans+(n-x+1)*s%mod)%mod;
            }
            else
            {
                ans=(ans+x*s%mod)%mod;
            }
        }
    }
    cout<<(ans+mod)%mod<<endl;
    return 0;
}