51nod 1040最大公约数之和 欧拉函数

51nod 1040最大公约数之和 欧拉函数

我是真的菜啊，为啥一直局限于想$\phi (n)$，没有想到其他的数应该怎么解决呢。。。找规律也找不到（废话，欧拉函数怎么肉眼看出来 ）
算了不说了自闭了，欧拉函数白学了………………

题意：给一个数n，计算下式:

$$ans=\sum _{i=1}^{n}gcd(i,n)$$

做法：分解定理分解，计算每一个质因子的贡献和

假设$gcd(a,n)=x且a<=n$，即问题转化为当gcd为x时，求此时的a的数量记作该因子的贡献
若$gcd(a/x,n/x)=1$,即求不大于n/x的且与n/x互质的个数，即$\phi (\frac{n}{x})$，
那么原式即求$\phi (\frac{n}{x})\ast x$

时间复杂度是$O(2logn\sqrt{n})$

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<stack>
#include<bitset>
#include<cstdlib>
#include<cmath>
#include<set>
#include<list>
#include<deque>
#include<queue>
#include<map>
#define ll long long
#define pb push_back
#define rep(x,a,b) for (int x=a;x<=b;x++)
#define repp(x,a,b) for (int x=a;x<b;x++)
#define W(x) printf("%d\n",x)
#define WW(x) printf("%lld\n",x)
#define pi 3.14159265358979323846
#define mem(a,x) memset(a,x,sizeof a)
#define lson rt<<1,l,mid
#define rson rt<<1|1,mid+1,r
using namespace std;
const int maxn=2e6+7;
const int INF=1e9;
const ll INFF=1e18;
ll phi(ll n)
{
    ll res=n;
    for (int i=2;i*i<=n;i++)
    {
        if (n%i==0)
        {
            res=res/i*(i-1);
            while(n%i==0)n/=i;
        }
    }
    if (n>1)res=res/n*(n-1);
    return res;
}

ll solve(ll n)
{
    ll ans=0;
    for (int i=1;i*i<=n;i++)
    {
        if (i*i==n)ans+=phi(i)*i;
        else if (n%i==0)
        {
            ans+=phi(n/i)*(i);
            ans+=phi(i)*(n/i);
        }
    }
    return ans;
}
int main()
{
    ll n;
    scanf("%lld",&n);
    WW(solve(n));
}