UVa--12716 GCD XOR(math)

UVa 12716

题意

给定整数n，求有多少对整数(a, b)满足: 1 <= b <= a <= n，且 gcd(a, b) = x ^ b.

题解

设 $c = gcd(a, b) = a \text{ Xor } b$,
则:  c 是 a 的约数，$b = a \text{ Xor } c$
因此，可以枚举 a 和 c，为降低枚举的复杂度，利用 c 是 a 的约数的关系，枚举 a 和 a 的约数，但约数并不好枚举，所以反过来，枚举 c 和 c 的倍数，这就类似于素数筛法。
此外，满足此等式的 $a, b, c$，满足: $c = a - b$
最后算法为：枚举 c 和 a, 计算 b = a - c，判断 $a \text{^} b == c$。

#include <bits/stdc++.h>
using namespace std;

const int maxn = 30000000;
int  sum[maxn];

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

int main()
{
    #ifdef LOCAL
    freopen("data.in", "r", stdin);
    #endif // LOCAL
    int t, n, kcase = 1;

    init();
    for(cin >> t; t--; )
    {
        cin >> n;
        printf("Case %d: %lld\n", kcase++, sum[n]);
    }
}