牛客小白赛 18 F.Forsaken的位运算魔法（类欧几里得算法（模板））

---

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

根据题解计算有多少个 $ij$ 第 p 位为 1，等价于计算 $\sum _{i=0}^n\sum _{j=0}^n(\lfloor \frac{ij}{2^p}\rfloor -2\lfloor \frac{ij}{2^{p+1}}\rfloor )$
当 $i$ 固定时，会发现这个是类欧几里得算法形式，套个类欧的模板即可。

类欧算法参考：https://www.cnblogs.com/encodetalker/p/11037506.html
https://blog.csdn.net/dreaming__ldx/article/details/86768953

计算有多少个 $ij$ 第 p 位为0，可以用总数减一下 为 1的个数得到
复杂度：$O(cn{\log }^{2}n)$，$c$ 是一个较大的常数

---

代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e5 + 10;
const int mod = 1e9 + 7;
typedef long long ll;
ll ans[100],n,k,inv2;
ll fpow(ll a,ll b) {
	ll r = 1;
	while(b) {
		if(b & 1) r = r * a % mod;
		b >>= 1;
		a = a * a % mod;
	}
	return r;
}
ll cal(ll x) {
	return x * (x + 1) % mod * inv2 % mod;
}
ll f(ll a,ll b,ll c,ll n) {
	if(a >= c || b >= c) {
		return (f(a % c,b % c,c,n) + cal(n) * (a / c) % mod + (b / c) * (n + 1) % mod) % mod;
	} else {
		ll m = (a * n + b) / c;
		if(a == 0) return 0;
		else return (n * m % mod + mod - f(c,c - b - 1,a,m - 1)) % mod;
	}
}
/*ll f(ll a,ll b,ll c,ll n) {		//另一种写法
	if(!a) return (n + 1) * (b / c) % mod;
	if(a >= c || b >= c) {
		return (f(a % c,b % c,c,n) + cal(n) * (a / c) % mod + (b / c) * (n + 1) % mod) % mod;
	} else {
		ll m = (a * n + b) / c;
		return (n * m % mod + mod - f(c,c - b - 1,a,m - 1)) % mod;
	}
} */
int main() {
	scanf("%lld%lld",&n,&k);
	inv2 = fpow(2,mod - 2);
	for(int i = 0; i <= 34; i++)
		for(int a = 0; a <= n; a++)
			ans[i] = (ans[i] + f(a,0,(1ll << i),n)) % mod;
		
	ll res = 0,tot = (n + 1) * (n + 1);
	for(int i = 0; i <= 34; i++) {
		if(k >> i & 1) {
			ll cnt = (ans[i] - 2 * ans[i + 1] + mod) % mod;
			res += (tot - cnt + mod) % mod * (1ll << i) % mod;
			res %= mod;
		} else {
			ll cnt = (ans[i] - 2 * ans[i + 1] + mod) % mod;
			res += cnt * (1ll << i) % mod;
			res %= mod;
		}
	} 
	printf("%lld\n",res);
	return 0;
}