题意：
求 x 在[1, n]范围内，y 在[1, m]范围内的满足 gcd(x, y) = k 的x,y对数。
题解：
- 与GCD那题稍有不同。首先这道题(x,y)(y,x)认为不同。
- 这道题要用到一个分块的技巧。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int N = 1e6 + 5;
int T, k, n, m;
int pri[N], cnt;
int mu[N];
int vis[N];
int sum[N];
void prime(){
mu[1] = 1;
sum[1] = 1;
for(int i = 2; i <= 4e5; ++ i){
if(!vis[i]){
pri[++ cnt] = i;
mu[i] = -1;
}
for(int j = 1; j <= cnt && pri[j] * i <= 4e5; ++ j){
vis[pri[j] * i] = 1;
if(i % pri[j] == 0){
mu[pri[j] * i] = 0;
break;
}
mu[pri[j] * i] = -mu[i];
}
sum[i] = sum[i - 1] + mu[i];
}
}
int main()
{
prime();
scanf("%d",&T);
for(int Case = 1; Case <= T; ++ Case){
scanf("%d%d%d", &n, &m, &k);
if(!k){
printf("0
", Case);
continue;
}
n /= k, m /= k;
if(n > m) swap(n, m);
ll ans = 0;
for(int i = 1, last = 1; i <= n; i = last + 1){
last = min(n / (n / i), m / (m / i));
ans += (ll)(sum[last] - sum[i - 1]) * (n / i) * (m / i);
}
printf("%lld
", ans);
}
return 0;
}