题目链接
解体思路
直接进行莫比乌斯反演
过程如下
我们枚举d,把
μ
\mu
μ求一个前缀和,然后通过分块来做这道题,不分块我不知道能不能过。
AC代码
#include <iostream>
using namespace std;
const int MAX_N = 1e7+8;
int vis[MAX_N], prime[MAX_N], mu[MAX_N];
long long f[MAX_N];
int ans = 0;
void init(int n)
{
mu[1] = 1;
for (int i=2; i<=n; i++)
{
if (!vis[i])
{
prime[ans++] = i;
mu[i] = -1;
}
for (int j=0; j<ans && i*prime[j] <= n; j++)
{
vis[i*prime[j]] = 1;
if (i%prime[j])
{
mu[i*prime[j]] = -mu[i];
}
else
{
mu[i*prime[j]] = 0;
break;
}
}
mu[i] += mu[i-1];
}
}
int main()
{
long long n;
long long sum = 0;
scanf("%lld", &n);
init(n);
for (int i=0; i<ans; i++)
{
long long m=n/prime[i];
for (long long l=1, r; l<=m; l=r+1)
{
r = min(m/(m/l), m);
sum+=(m/l)*(m/l)*(mu[r]-mu[l-1]);
}
}
printf("%lld\n", sum);
return 0;
}