UVA12716 GCD XOR

UVA12716 GCD XOR

---

题目

洛谷原题
说白了就是求 ${\sum }_{i=1}^n{\sum }_{j=i}^n[\gcd (i,j)==ixorj]$

---

题解

当 $i=j$ 时，条件 $\gcd (i,j)=ixorj$ 不可能成立，钦定 $i>j$。
分析：
$\because$ $\gcd (i,j)=\gcd (i-j,j),\gcd (i-j,j)\le i-j$
$\therefore$ $\gcd (i,j)\le i-j$

$\because$ 异或相当于不退位的减法
$\therefore$ $ixorj\ge i-j$

$\therefore$ $\gcd (i,j)\le i-j\le ixorj$

$\because$ $\gcd (i,j)=ixorj,\gcd (i,j)\le i-j\le ixorj$
$\therefore$ $ixorj=i-j$

$\therefore$ 问题变成了${\sum }_{i=1}^n{\sum }_{j=i}^n[ixorj==i-j]$
然而直接暴力枚举 $i$ 和 $j$ 任会超时。

---

$\because$ $\gcd (i,j)=ixorj,ixorj=i-j$
$\therefore$ $\gcd (i,j)=i-j$
$\because$ $\gcd (i,j)=\gcd (i,i-j)$
$\therefore$ $\gcd (i,i-j)=i-j$
$\therefore$ $i$ 为 $i-j$ 的倍数
$\therefore$ 我们只需枚举 $j$，再枚举 $j$ 的倍数 $i$，再判断 $ixorj=i-j$ 就行。
然而还是会超时······

---

于是就需要预处理出$ans$：

for(int j=1; j<=Max/2; j++)
	for(int i=j<<1;i<=Max;i+=j)
		if((i^j)==i-j)
			ans[i]++;

再求前缀和：

for(int i=2;i<=Max;i++)
	ans[i]+=ans[i-1];

$ans[i]$表示$i$为$1$~$i$的答案。

---

代码

贴上十分简短的代码：

#include<cstdio>
#include<iostream>
#define Max 30000001
#define Max2 15000000
using namespace std;
int t,n,ans[Max+5];
int main()
{
	for(int j=1; j<=Max2; j++)
		for(int i=j<<1;i<=Max;i+=j)
			if((i^j)==i-j)
				ans[i]++;
	for(int i=2;i<=Max;i++)
		ans[i]+=ans[i-1];
	scanf("%d",&t);
	for(int i=1; i<=t; i++)
	{
		scanf("%d",&n);
		printf("Case %d: %d\n",i,ans[n]);
	}
}