POJ - 2480 Longge‘s problem（积性函数+欧拉函数）

传送门

---

$gcd(x,n)$ 当$n$固定时，该函数是积性函数，因为对$n$质因数分解得到$p_1^{a_1}{p}_2^{a_2}...p_n^{a_n}$，那么$gcd(x,n)=gcd(x,p_1^{a_1})\ast gcd(x,p_2^{a_2})...gcd(x,p_n^{a_n})$成立，故为积性函数

设$gcd(x,n)=g,1\le x\le n$，那么$g\ast gcd(x/g,n/g)=1$，而满足$gcd(x/g,n/g)=1$的所有数，实际上就是小于$n/g$且与$n/g$的互质个数，即欧拉函数$\phi (n/g)$

那么问题就变成了这样一个公式${\sum }_{d|n}d\ast \phi (n/d)$，只需要枚举$n$的所有因数，即可求出答案

因为$2^{10}\ast 3^{10}=6e7$，对某个小于$10$的质数再次求幂，这个新的幂次肯定也为个位数，大胆猜想小于$int$的数最多只有几百个因数，那么总的时间复杂度为$O(k\ast \sqrt{n})$

#include <cstdio>
#include <math.h>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn = 2e5 + 10;

ll euler_phi(ll n) {
    int m = sqrt(n + 0.5);
    ll ans = n;
    for (int i = 2; i <= m; i++) {
        if (n % i == 0) {
            ans = ans / i * (i - 1);
            while (n % i == 0) n /= i;
        }
    }
    if (n > 1) ans = ans / n * (n - 1);
    return ans;
}

int main() {
    ll n;
    vector<int> fac;
    //cout<<euler_phi(7)<<endl;
    while (cin >> n) {
        fac.clear();
        for (ll i = 1; i * i <= n; i++)
            if (n % i == 0) {
                if (i * i == n) fac.push_back(i);
                else {
                    fac.push_back(i);
                    fac.push_back(n / i);
                }
            }
        ll ans = 0;
        for (int i = 0; i < fac.size(); i++) {
            ans += 1LL * fac[i] * euler_phi(n / fac[i]);
        }
        cout << ans << "\n";
    }
    return 0;
}