【NOIP2016提高A组模拟8.19】(雅礼联考day2)公约数

题目

给定一个正整数，在[1,n]的范围内，求出有多少个无序数对(a,b)满足gcd(a,b)=a xor b。

分析

显然a=b是一定不满足，

我们设a>b，

易得gcd(a,b)<=a-b、a xor b>=a-b

那么gcd(a,b)=a xor b=a-b

gcd(a, a xor c)=c,而c是a的约数

设a-b=c，我们枚举它

a=i*c。

那么就只用判断a xor c=a-c即可。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=1000000007;
const int N=50005;
using namespace std;
int ans,n;
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
for(int j=2;j<=n/i;j++)
{
int a=i*j;
if((a^i)==(a-i))
ans++;
}
printf("%d",ans);
}