51nod1040最大公约数之和（欧拉函数）

题目：51nod1040.
题目大意：给定$n$，求$1$~$n$与$n$的最大公约数之和.
$1\le n\le 10^9$.

先将式子写出来：
$$\sum _{i=1}^{n}gcd(i,n)$$

考虑从枚举$1$~$n$转变为枚举$n$的因数：
$$\sum _{i=1}^{n}gcd(i,n)=\sum _{d|n}d\sum _{i=1}^n[gcd(i,n)=d]$$

我们设$m=\frac{n}{d}$，那么：
$$\begin{gathered} \sum _{i=1}^{n}gcd(i,n) \\ =\sum _{d|n}d\sum _{i=1}^m[gcd(i,m)=1] \\ =\sum _{d|n}d\varphi (m) \\ =\sum _{d|n}d\varphi (\frac{n}{d}) \end{gathered}$$

由于$n$很大，我们没有办法直接线性筛筛出$\varphi$值，所以直接根号求解就可以了.

那么时间复杂度为$O({\sqrt{n}}^2)=O(n)$，但是这个上界很松，事实上复杂度明显跑不满，所以可以过.

不过一个更正经的做法是，大力分解$n$的所有质因数，然后计算欧拉函数时直接用$n$的质因数计算，时间复杂度$O(\sqrt{n}\log n)$.

代码如下：

#include<bits/stdc++.h>
using namespace std;

#define Abigail inline void
typedef long long LL;

const int N=50;

int n;
int pr[N+9],cp;

void Get_pr(int n){
  int now=n;
  for (int i=2;i*i<=n;++i)
    if (now%i==0)
      for (pr[++cp]=i;now%i==0;now/=i);
  if (now>1) pr[++cp]=now;
}

int Get_phi(int n){
  int res=n;
  for (int i=1;i<=cp;++i)
    if (n%pr[i]==0) res=res/pr[i]*(pr[i]-1);
  return res;
}

LL ans;

Abigail into(){
  scanf("%d",&n);
}

Abigail work(){
  Get_pr(n);
  for (int i=1;i*i<=n;++i)
    if (n%i==0){
      ans+=(LL)i*Get_phi(n/i);
      if (i*i^n) ans+=(LL)n/i*Get_phi(i);
    }
}

Abigail outo(){
  printf("%lld\n",ans);
}

int main(){
  into();
  work();
  outo();
  return 0;
}