Description
Input
仅包含一行,为两个整数n和m。
Output
仅包含一个整数,表示总共产生的能量损失。
Sample Input
5 4
【样例输入2】
3 4
Sample Output
36
【样例输出2】
20
对于100%的数据:1 ≤ n, m ≤ 100,000。
题解
可以发现(x,y)的能量损耗就是gcd(x,y)-1,所以
$$\begin{aligned}
ans &= \sum_{i, j} 2(gcd(i,j)-1)+1\\
&= 2\left(\sum_{i, j}gcd(i, j)\right)-nm\\
&= 2\left(\sum_{d}d\sum_{i, j}\left[gcd(i, j)=d\right]\right)-nm\\
&= 2\left(\sum_{d}d\sum_{d|t}\mu\left(\frac td\right)\sum_{i, j}\left[t|gcd(i, j)\right]\right)-nm\\
&= 2\left(\sum_{d}d\sum_{d|t}\mu\left(\frac td\right)\left\lfloor\frac nt\right\rfloor\left\lfloor\frac mt\right\rfloor\right)-nm\\
&= 2\left(\sum_{t}\left\lfloor\frac nt\right\rfloor\left\lfloor\frac mt\right\rfloor\sum_{d|t}d\mu\left(\frac td\right)\right)-nm\\
&= 2\left(\sum_{t}\left\lfloor\frac nt\right\rfloor\left\lfloor\frac mt\right\rfloor\phi(t)\right)-nm
\end{aligned}$$
枚举t,计算即可。
附代码:
#include <algorithm>
#include <cstdio>
#include <cstring>
typedef long long LL;
const int N = 100050;
int phi[N];
void getPhi() {
for (int i = 1; i < N; ++i) phi[i] = i;
for (int i = 2; i < N; ++i) if (phi[i] == i)
for (int j = i; j < N; j += i)
phi[j] -= phi[j] / i;
}
int main() {
int n, m;
scanf("%d%d", &n, &m);
LL ans = 0;
getPhi();
for (int i = 1; i <= n && i <= m; ++i)
ans += (LL)phi[i] * (n / i) * (m / i);
printf("%lld\n", ans * 2 - (LL)n * m);
return 0;
}