UVA12716 - GCD XOR（埃氏筛思想+数论）

题目链接 https://cn.vjudge.net/status/#un=&OJId=UVA&probNum=&res=0&language=&onlyFollowee=false

【题意】
输入整数n（1<=n<=3e7）问有多少对整数（a,b）满足1<=b<=a<=n 且 gcd(a,b)=a xor b

【思路】
暴力打表可以发现当gcd(a,b)=a xor b时，gcd(a,b)=a-b（a>b）可以证明这个结论是对的
首先当a>b时，有如下结论
a-b<=a xor b
若a xor b=c，则a xor c=b
那么如果gcd(a,b)=a xor b=c成立，则说明a-b<=c，又因为c是a,b的最大公约数，那么假设a=k1*c,b=k2*c,k1>k2
则a-b=(k1-k2)*c>=c，因此a-b=c成立
特殊的，当a=b时gcd(a,b)=a,a^b=0，不符合题意

然后用埃氏筛的思想从小到大枚举最大公约数c以及所有c的倍数a，只要a^(a-c)=c说明找到了一组解
最后求再一下前缀和

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

const int maxn=30000005;

int ans[maxn];

void init(){
    for(int c=1;c<maxn;++c){
        for(int a=c*2;a<maxn;a+=c){
            int b=a-c;
            if((a^b)==c) ++ans[a];
        }
    }
    for(int i=2;i<maxn;++i) ans[i]+=ans[i-1];
}

int main(){
    init();
    int T;
    scanf("%d",&T);
    for(int kase=1;kase<=T;++kase){
        int n;
        scanf("%d",&n);
        printf("Case %d: %d\n",kase,ans[n]);
    }
    return 0;
}