【POJ】2480 Longge's problem（欧拉函数）

题目

传送门：QWQ

分析

题意就是求∑gcd(i, N) 1<=i <=N.。

显然$ gcd(i,n) = x $时，必然$x|n$。

所以我们枚举一下n的约数，对于每个约数x，显然$ gcd(i/x,n/x)=1$

所以我们计算一下n/x的欧拉函数就ok了。

联赛前刷水题qwq

代码

// #include <bits/stdc++.h>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn = ;
ll a[maxn], top=, p[maxn];
int fac(ll x) {
    int n=x, flag=;
    ll end = sqrt(x);
    for(ll i=;i<=end;i++) {
        if(x%i==) {
            flag=;
            a[++top] = i;
            if(i*i!=x) a[++top] = x/i;
            while(n%i==) n/=i;
        }
    }
    if(x==) flag=;
    return flag;
}
ll phi(ll x) {
    ll ans = x, end = sqrt(x)+;
    for(ll i=;i<=end;i++) {
        if(x%i==) {
            ans = ans / i * (i-);
            while(x%i==) x/=i;
        }
    }
    if(x>) ans = ans / x * (x-);
    return ans;
}
void Ph() {
    for(int i=;i<maxn-;i++) p[i]=i;
    for(ll i=;i<maxn-;i++) {
        if(p[i]==i) {
            for(ll j=i;j<maxn-;j+=i) {
                p[j] = p[j] / i * (i-);
            }
        }
    }
}
int main() {
    int n; Ph();

    while(scanf("%lld",&n)==) {
        if(n==){
            printf("1\n"); continue;
        }
        top=; int q=fac(n);
        if(q == ) {
            printf("%lld\n",phi(n)+n); continue;
        }
        sort(a+,a++top);
        int pp = top;
        ll ans=phi(n)+n;
        for(int i=;i<=pp;i++) {
            ll qwq;
            if(n/a[i] < maxn-) qwq = p[n/a[i]];
            else qwq = phi(n/a[i]);
            ans += a[i] * qwq;
        }
        printf("%lld\n",ans);
    }
}