SP5971 LCMSUM-LCM Sum (función de Euler, derivación de fórmulas)

SP5971 LCMSUM - Suma de LCM

Significado: preguntar

Arribista

$$\begin{gathered} \sum _{i=1}^{n}lcm(yo,n) \\ \sum _{i=1}^n\frac{in}{gcd(yo,n)} \\ norte\sum _{i=1}^n\frac{yo}{gcd(yo,n)} \\ norte\underset{d|n}{\Sigma }\sum _{i=1}^n\frac{yo}{re}\left[gcd\right(yo,n)=d] \\ norte\underset{d|n}{\Sigma }\sum _{i=1}^{\frac{n}{re}}yo\left[gcd\right(yo,\frac{n}{re})=1] \\ norte\underset{d|n}{\Sigma }\frac{\left(\varphi \right(d)+[d=1\left]\right)\times d}{2} \end{gathered}$$

La primera fórmula es el problema; la
segunda fórmula es poner $lcm\; está$ desensamblado; la
tercera fórmula es sacar la constante n; la
cuarta fórmula es enumerar todos los factores de n, que se pueden transformar en la cuarta fórmula; la
quinta fórmula es ad de$\left[\right[gcd(i,n)=d]$ Propóngalo dentro, y luego seleccione el rango de valores, número enumerado$\frac{yo}{re}$Transformación $yo$和$\left[gcd\right(yo,n)=d]\; se$ convierte en$\left[gcd\right(yo,\frac{n}{re})=1]$ ; La
sexta fórmula es usar una propiedad de la función de Euler:$ComosiGCD(yo,n)=1,luegogcd(n-i,n)=1(yo\ne 1)$ , y luego se puede transformar en la sexta fórmula;
final

$Code$

#include<bits/stdc++.h>
using namespace std;
const int N = 1e6+10;
typedef long long ll;
ll st[N], p[N], phi[N], sum[N], ans[N], tot;
//phi：欧拉函数，sum函数第六个公式中的欧拉函数phi*d，ans是存放答案的数组；
void init() {
    
    
    phi[1] = 1, st[1] = 1, sum[1] = 1;//sum[1]=1,把d=1的情况考虑了进去
    for(int i=2; i<N; i++) {
    
    
        if(!st[i]) p[tot++] = i, phi[i] = i-1, sum[i] = 1ll*phi[i]*i/2;
        for(int j=0; j<tot&&i*p[j]<N; j++) {
    
    
            st[i*p[j]] = 1;
            if(i % p[j] == 0) {
    
    
                phi[i*p[j]] = phi[i]*p[j];
                sum[i*p[j]] = 1ll*phi[i*p[j]]*(i*p[j])/2; 
                break;
            }
            phi[i*p[j]] = phi[i]*(p[j]-1);
            sum[i*p[j]] = 1ll*phi[i*p[j]]*(i*p[j])/2; 
        }
    }
    //类埃氏筛，把n的全部因子d加起来ans[n]+=sum[d];
    for(int i=1; i<N; i++) {
    
    
        for(int j=i; j<N; j+=i)
        ans[j] += sum[i];
    }
}
inline ll read() {
    
    //快读
    ll x = 0, f = 1; char c = getchar();
    while(c < '0' || c > '9') {
    
    
        if(c == '-') f = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
    
    
        x = (x << 1) + (x << 3) + (c ^ 48);
        c = getchar();
    }
    return x * f;
}

int main() {
    
    
    init();//初始化
    ll t;
    t = read();
    while(t--) {
    
    
        int n;
        n = read();
        printf("%lld\n", ans[n]*n);
    }
    return 0;
}

Fin