Longge的问题 [Bzoj 2705，SDOI2012]

题目地址请点击——

---

Longge的问题

---

问题描述

$Longge$ 的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。
现在问题来了：给定一个整数 $N$，你需要求出 $\sum_{i=1}^Ngcd(i, N)$ 。

---

输入描述

一个整数，为 $N$。

---

输出描述

一个整数，为所求的答案。

---

样例输入

6

---

样例输出

15

---

HINT

对于 $60\%$ 的数据，$0<N<=2^{16}$。

对于 $100\%$ 的数据，$0<N<=2^{32}$。

---

Solution

设 $x_1$，$x_2$，$x_3$，……，$x_k$ 为 $n$ 的 $k$ 个因数，$f(i)$ 表示满足 $gcd(j,n)=i\ (1<=j<=n)$ 的 $j$ 的个数。

$\therefore\ ans=\sum_{i=1}^kf(x_i)·x_i$ 。

显然，$f(1)=\varphi(n)$。

那么，如何计算 $f(k)\ (k>1)$ 呢？

∵ gcd(j,n)=i (1<=j<=n)⇒gcd(ji,ni)=1⇒f(i)=φ(ni)

$$\begin{align*} &\because \ gcd(j,n)=i\ (1<=j<=n)\\ &\Rightarrow gcd(\frac{j}{i},\frac{n}{i})=1\\ &\Rightarrow f(i)=\varphi(\frac{n}{i})\\ \end{align*}$$

$\therefore\ ans=\sum_{i=1}^k \varphi(\frac{n}{x_i}) ·x_i$ 。

---

Code

[cpp]

1. #include <iostream>  
2. #include <cstdio>  
3. #include <cmath>  
4.   
5. #define LL long long  
6.   
7. using namespace std;  
8.   
9. LL n,sqr,ans;  
10.   
11. LL phi(LL x){  
12.     LL tmp=x,tmp2=x,sqr=(LL)sqrt((double)x);  
13.     for(LL i=2;i<=sqr;i++){  
14.         if(x%i==0){  
15.             tmp=tmp*(i-1)/i;  
16.             while(x%i==0)x/=i;  
17.         }  
18.     }  
19.     if(x!=1)tmp=tmp*(x-1)/x;  
20.     return tmp;  
21. }  
22.   
23. int main(){  
24.     scanf(”%lld”,&n);  
25.     sqr=(LL)sqrt((double)n);  
26.     for(LL i=1;i<=sqr;i++)  
27.         if(!(n%i)){  
28.             ans+=i*phi(n/i);  
29.             if(i!=n/i)ans+=(n/i)*phi(i);  
30.         }  
31.     printf(”%lld\n”,ans);  
32.     return 0;  
33. }