NOIP 提高 公约数

Description

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

Solution

我们可以发现对于一个a、b，gcd（a,b）=gcd(b,a-b),又gcd(b,a-b)<=min(b,a-b),所以gcd(a,b)<=a-b。

而考虑异或操作，对于a,b两个值的异或，一定是相同位为0，否则为1，那么对于a-b，转换成二进制上的加减，一定也是相同位为0，否则为1，不同的是，它可能还要向上一位借一个1，所以一定有a-b<=a xor b。

中和上面推导，gcd(a,b)<=a-b<=a xor b,因为我们要求gcd(a，b)=a xor b，所以只能是gcd(a,b)=a-b=a xor b，所以我们枚举a-b，设a-b=k，gcd(a,b)=k，所以设a=kx,b=ky,gcd(x,y)=1,a-b=kx-ky=k(x-y)=k,所以x-y=1，但是当gcd（x,y）=z（z!=1）时，gcd（a,b）=kz,而a-b=k* z* x1-k* z* y1=k* z(x1-y1),所以gcd（x,y）=z（z!=1）仍然是一个可能正确的方案，只需要判断a xor b是否与gcd(a,b)相等即可。

代码

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
int n,i,t,j,k,l,ans;
int main(){
    scanf("%d",&n);
    i=2;
    while(i<n){
        t=i+1;
        k=n/t;ans++;
        for (j=2;j<=k;j++)
            if (((i*j)^(t*j))==j) ans++;
        i+=2;
    }
    printf("%d\n",ans);
}