luogu UVA12716 GCD等于XOR GCD XOR

背景：

一年之前看到的题了。

题目传送门：

https://www.luogu.org/problem/UVA12716

题意：

求$1\le b\le a\le n$的数对满足$\gcd (a,b)=a\text{ xor }b$的组数。

思路：

设$\gcd (a,b)=a\text{ xor }b=c$
考虑两个不等式：
$[1]$：
考虑更相减损术，有：$\gcd (a,b)=\gcd (b,a-b)=c$
因此显然有$a-b\ge c$。
$[2]$：
结论：$a-b\le c$。
将$a,b$转成二进制，考虑当前这一位$i$，由于$a\ge b$，因此若$a_i=1$，则$b_i=0,1$显然成立；若$a_i=0,b_i=0$，显然成立；若$a_i=0,b_i=1$，则前面一定存在一个$j<i$，且$a_j=1,b_j=0$，则这一段的异或结果一定形如$100..001$，而这一段的减法结果一定形如$01...11$，成立。

综上$a-b=c$。
因此$\gcd (a,b)=a\text{ xor }b=a-b$时成立。
考虑枚举$\gcd (a,b)$作为$b$，枚举$\gcd (a,b)$的倍数作为$a$，此时已经满足了$\text{gcd}$的限制，直接判断$a\text{ xor }b,a-b$是否相等即可。
可以类似埃氏筛来做。
离线，前缀和一下即可。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
	int n;
	int f[30000010];
void init(int ma)
{
	for(int i=1;i<=ma;i++)
		for(int j=(i<<1);j<=ma;j+=i)
			if(j-i==(i^j)) f[j]++;
	for(int i=1;i<=ma;i++)
		f[i]+=f[i-1];
}
int main()
{
	int T;
	init(30000000);
	scanf("%d",&T);
	for(int i=1;i<=T;i++)
	{
		scanf("%d",&n);
		printf("Case %d: %d\n",i,f[n]);
	}
}