异或【数论】

>Description
SarvaTathagata是个神仙，一天他在研究数论时，书上有这么一个问题：求不超过n两两的数的gcd。
SarvaTathagata这么神仙的人当然觉得这个是sb题啦。学习之余，他还发现gcd的某一个特别好的性质：如果有两个数i,j满足gcd(i,j)=i^j(这里的为c++中的异或)的话，那么这两个数组成的数对(i,j)就是一个nb的数对(这里认为(i,j)和(j,i)为相同的，并不需要计算2次)。
当然，SarvaTathagata并不会只满足于判断一个数对是否nb，他还想知道满足两个数都是不超过n并且nb的数对有多少个。
由于SarvaTathagata实在是太神仙了，他认为这种题实在是太简单了。于是他找到了你，看看你是否能解决这个问题。

---

>Input
共一行一个整数n，含义如题所述。

>Output
一行一个整数，表示nb的数对的个数。

---

>Sample Input
123456

>Sample Output
214394

---

>解题思路
一道数论题，直接暴力是肯定不行的→→

题目要求$gcd(a,b)=axorb$，异或的意思是两数在二进制的状态下相对位置相同为0，不同为1

默认情况a<b
∵
a xor b ≥ b - a
gcd(a,b) ≤ b - a
又∵
要求a xor b = gcd(a,b)
∴
a xor b = gcd(a,b) = b - a

所以只用枚举gcd(a,b)和b就行了（这里gcd(a,b)就是a）

---

>代码

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;

int n, ans;

int main()
{
	scanf ("%d", &n);
	for (int i = 1; i <= n; i++) //枚举a(gcd)
	  for (int j = 2; i * j <= n; j++) //b一定是a的倍数，因为要使gcd(a,b)=b-a
	    if (i * j - i == (i ^ (i * j))) ans++; //gcd(a,b)已等于b-a，所以判断a xor b是否等于b-a
	printf ("%d", ans);
	return 0;
}