https://www.luogu.org/problem/show?pid=2398
很有意思的题目;
这种题特别好;
代码不长,思路不浅;
题解细细看能看懂;
AC之后有所收获;
我们设f[i]表示gcd为i的对数;
那么对答案的贡献就是f[i]*i;
怎么求f[i]呢;
我们看
包含i的因子的数是不是有(n/i)个;
那显然包含i因子的两个数 的gcd一定包含i;
换句话说,那个gcd一定是i的倍数;
所以
f[i]=(n/i)*(n/i)-f[i*2]-f[i*3]-f[i*4]….
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#define Ll long long
using namespace std;
Ll n;
Ll f[100001],ans;
int main()
{
scanf("%lld",&n);
for(int i=n;i;i--){
f[i]=(n/i)*(n/i);
for(int j=2;i*j<=n;j++)f[i]-=f[j*i];
ans+=f[i]*i;
}
printf("%lld",ans);
}