BZOJ 2226: [Spoj 5971] LCMSum（数论+欧拉函数）

题目描述

传送门

题目大意：问你1到n中的数和n的lcm之和，T组询问。T <= 300000，n <= 1000000。

---

思路

一看就是比较套路的题。随便化简一下，大概是这样

$$\sum_{i=1}^n[n,i]=\sum_{i=1}^n{ni\over(n,i)}=n\sum_{i=1}^n\sum_{d=1}^n[(i,n)=d]{i\over d}$$

然后先枚举$d$（只能是$n$的约数），令$i=id$（将$i$看成$id$），变成

$$n\sum_{d|n}\sum_{i=1}^{n\over d}[(i,{n \over d})=1]i$$
接下来还要不要 反演呢？其实到这里已经大功告成了。根据我不知道从哪里偷来的公式 $\sum_{i=1}^n[(i,n)=1]i={{\varphi(n)*n+[n=1]}\over2}$，答案就是
$$n\sum_{d|n}{\varphi(d)*d+[d=1] \over 2}$$

至于这个公式怎么证？显然。（证它干嘛呢）

于是总时间就变成了$O(n+T\sqrt n)$，理论跑不过，但政治书上说，理论要联系实际。

---

代码

#include <bits/stdc++.h>
#define maxn 1000100
#define temp (i * prime[j])
using namespace std;

typedef long long LL;
int n, T, cnt;
int prime[maxn], phi[maxn];
bool Vis[maxn];
LL ans;

void Da(){
    phi[1] = 1;
    for(int i = 2; i < maxn; i++){
        if(!Vis[i]){
            prime[++cnt] = i;
            phi[i] = i - 1;
        }
        for(int j = 1; j <= cnt && temp < maxn; j++){
            Vis[temp] = true;
            if(i % prime[j] == 0){
                phi[temp] = phi[i] * prime[j];
                break;
            }
            else  phi[temp] = phi[i] * (prime[j] - 1);
        }
    }
}

int main(){

    scanf("%d", &T);

    Da();

    while(T --){
        scanf("%d", &n);
        ans = 0LL;
        for(int i = 1; i * i <= n; i++)
            if(n % i == 0){
                ans += (1LL * phi[i] * i + ((i == 1) ? 1 : 0)) >> 1LL;
                if(n / i != i)  ans += (1LL * phi[n/i] * (n/i)) >> 1LL;
            }
        printf("%lld\n", n * ans);  
    }

    return 0;
}

---

明敕星驰封宝剑，辞君一夜取楼兰。