题意:求
∑i=1ngcd(i,n)
这种题目居然也是省选题qwq..
考虑到n<=1e9,于是不可能直接暴力搞。
直接求出每一种约数的贡献然后直接乘起来就可以了
问题就转化成求每一个约数出现的次数,即
∑gcd(i,j)=d
继续转化就是gcd(i/d,j/d)=1,那么贡献为phi(max(i,j)/d).
但是时间复杂度很鬼畜啊,感觉根号套根号不久退化成O(n)了,为啥能过了。。。
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e5+5;
typedef long long ll;
ll n,m,ans;
ll phi(ll x)
{
ll y=x;
fo(i,2,m)
if (x%i==0)
{
y=y/i*(i-1);
while(x%i==0)x/=i;
}
if (x>1)y=y/x*(x-1);
return y;
}
int main()
{
scanf("%lld",&n);
m=1ll*sqrt(n);
fo(i,1,m)
{
if (n%i==0)
{
ans+=i*phi(n/i);
if (i*i<n)ans+=(n/i)*phi(i);
}
}
printf("%lld\n",ans);
}