LINK
题意:求1~n, 和 1 ~ m 中互质的本质不同的数对的数量。
设n < m, 则可以转变成求 1 ~ n 与 1 ~ n 的数对 和 1 ~ n 与 n + 1 ~ m 的 数对, 前者直接时欧拉函数的前缀和即可(phi[i] :小于 i 且与 i 互质的数的数量。
所以主要要考虑后者,很容易想到1 ~ n 和 j 互质的数 = n - (1 ~ n 与 j 不互质的数), 那怎么求不互质的数呢,
可以将j分解为质因数, 显然质因数的倍数和j不互质,
比如 2x, 3x, 所以我们要减去这部分, n / 2, n / 3, 但有一个关键点, 这样就有可能产生重复, 比如 6 既是2的倍数, 又是3的倍数, 这样就必须要用容斥定理去重 减去
n / (2 * 3)
这个方法网上有三种, 这里借鉴的是 DFS
#include <bits/stdc++.h>
using namespace std;
const int N = 1e5 + 10;
#define ll long long
int num[N] ={0};// i的素因子个数
int p[N][20]; // i的素因子
ll phi[N] = {0};
void phi_table()
{
phi[1] = 1;
for(int i = 2; i < N; i ++)
{
if(phi[i])
{
phi[i] += phi[i-1];
continue;
}
for(int j = i; j < N; j += i)
{
if(!phi[j]) phi[j] = j;
phi[j] = phi[j] * (i - 1) / i;
p[j][num[j] ++] = i;
}
phi[i] += phi[i - 1];
}
}
int DFS(int k, int b, int cur)// 求[1, b]中与cur不互质的数的个数
{
int rel = 0;
for(int i = k; i < num[cur]; i ++)
{
rel += b / p[cur][i] - DFS(i + 1, b / p[cur][i], cur);
}
return rel;
}
int main()
{
phi_table();
int t, caseno = 1;
cin >> t;
while(t --)
{
int a, b, c, d, k;
cin >> a >> b >> c >> d >> k;
if(k == 0)
{
cout << "Case " << caseno ++ << ": " << 0 << endl;
continue;
}
if(b > d) swap(b, d);
b /= k, d /= k;
ll ans = phi[b];
for(int i = b + 1; i <= d; i ++)
{
ans += b - DFS(0, b, i);
}
cout << "Case " << caseno ++ << ": " << ans << endl;
}
}