題目

傳送門：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 = 7000;
ll a[maxn], top=0, p[maxn];
int fac(ll x) {
    int n=x, flag=0;
    ll end = sqrt(x);
    for(ll i=2;i<=end;i++) {
        if(x%i==0) {
            flag=1;
            a[++top] = i;
            if(i*i!=x) a[++top] = x/i;
            while(n%i==0
 
) n/=i;
        }
    }
    if(x==2) flag=0;
    return flag;
}
ll phi(ll x) {
    ll ans = x, end = sqrt(x)+1;
    for(ll i=2;i<=end;i++) {
        if(x%i==0) {
            ans = ans / i * (i-1);
            while(x%i==0) x/=i;
        }
    }
    if(x>1) ans = ans / x * (x-1);
    return
 
 ans;
}
void Ph() {
    for(int i=2;i<maxn-10;i++) p[i]=i;
    for(ll i=2;i<maxn-10;i++) {
        if(p[i]==i) {
            for(ll j=i;j<maxn-10;j+=i) {
                p[j] = p[j] / i * (i-1);
            }
        }
    }
}
int main() {
    int n; Ph();

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

【POJ】2480 Longge's problem（歐拉函數）