[UVA-12716] GCD XOR

问题：

$\forall a < b \le N$, 求有多少对 $(a, b)$ 使得 $a \,\mathrm{xor}\, b = \mathrm{gcd}(a, b)$。

推导：

满足条件的 $(a, b)$ 一定有 $a \,\mathrm{xor}\, b = \mathrm{gcd}(a, b) = b - a$。

证明：

• $\mathrm{gcd}(a, b) \le b - a$：
  

  $a, b$ 均为 $\mathrm{gcd}(a, b)$ d的倍数，则 $b-a$ 也为其倍数，所以 $\mathrm{gcd}(a,b) \le b - a$。

• $a \,\mathrm{xor}\, b \ge b - a$：
  二进制下一共有4种减法情况：

b	a	c
0	0	0
0	1	1(带借位)
1	0	1
1	1	0

可以看出结果总是 a, b 的异或，并且可能带有借位，减小了实际值，所以有$a \,\mathrm{xor}\, b \ge b - a$

• 设 $g = a \,\mathrm{xor}\, b = \mathrm{gcd}(a, b)$, 则有 $g = b - a$

算法：

枚举 $\hat{g}$，对其倍数 $\hat{b}$，得到唯一的 $\hat{a} = \hat{g} \,\mathrm{xor}\, \hat{b}$，检验是否有 $\hat{g} = \hat{b} - \hat{a}$。

代码：

#include <bits/stdc++.h>

using namespace std;
const int N_MAX = 3e7 + 10;
int N;

int ans[N_MAX];
inline int est() {
  int cnt = 0;
  for (int c = 1; c <= N; c++) {
    for (int a = c + c; a <= N; a += c) {
      // divisor of a is c
      if ((a - c) == (a ^ c)) {
        cnt++;
        ans[a]++;
      }
    }
  }
  return cnt;
}
int main() {
  N = N_MAX - 1;
  est();
  for (int i = 1; i < N_MAX; i++) ans[i] += ans[i-1];
  int T;
  scanf("%d", &T);
  int ks = 0;
  while (T--) {
    int x; scanf("%d", &x);
    printf("Case %d: %d\n", ++ks, ans[x]);
  }
}