POJ2480：求所有最大公因数之和（附带详细数学推导）

Description

Longge is good at mathematics and he likes to think about hard mathematical problems which will be solved by some graceful algorithms. Now a problem comes: Given an integer N(1 < N < 2^31),you are to calculate ∑gcd(i, N) 1<=i <=N.
“Oh, I know, I know!” Longge shouts! But do you know? Please solve it.

Input
Input contain several test case.
A number N per line.

Output
For each N, output ,∑gcd(i, N) 1<=i <=N, a line

Sample Input
2
6

Sample Output
3
15

推导（其中所用欧拉函数与积性函数知识可自行查阅）：

$$设f(N)=\sum _{i=1}^{N}gc{d}_{(i,N)}$$可知f(n)是积性函数（此处不给出证明）
不妨令N的素因子分解式N=$p_1^{a_1}$·$p_2^{a_2}$·…·$p_n^{a_n}$（$p_i$是素数）

所以由积性函数的性质可知
$f(N)=f(p_1^{a_1})\cdot f(p_2^{a_2})\cdot \dots \cdot f(p_n^{a_n})$，

$$对f(p_1^{a_1})=\sum _{i=1}^{p_1^{a_1}}gc{d}_{(i,p_1^{a_1})}$$
和$p_1^{a_1}$互素的数的个数有$\phi$($p_1^{a_1}$)个，他们的最大公因数是1；
和$p_1^{a_1}$的最大公因数为$p_1$的数都可以写成$p_1\cdot k$且k与$p_1$互素，故k可以取[1,$p_1^{a_1-1}$]中与$p_1$互素的数，而与$p_1^{a_1-1}$不互素的数均为$p_1$的倍数，所以满足条件的k有$\phi$($p_1^{a_1-1}$)个，即和$p_1^{a_1}$的最大公因数为$p_1$的数有$\phi$($p_1^{a_1-1}$)个；
依此类推……
和$p_1^{a_1}$的最大公因数为$p_1^{a_1-1}$的数有$\phi (p_1)$个。
和$p_1^{a_1}$的最大公因数为$p_1^{a_1}$的数有$\phi$(1)=1个。
将它们求和，即可得
$f(p_1^{a_1})=1\cdot \phi (p_1^{a_1})+p_1\cdot \phi (p_1^{a_1-1})+\dots +p_1^{a_1-1}\cdot \phi (p_1)+p_1^{a_1}\cdot \phi (1)$

由欧拉函数公式$\phi (x)=x\prod _{i=1}^n(1-\frac{1}{p_i})$
(其中$p_i$是x的所有质因数，x是不为0的整数)
所以$f(p_1^{a_1})=a_1\cdot p_1^{a_1}\cdot (1-\frac{1}{p_1})+p_1^{a_1}$
$=p_1^{a_1}[a_1(1-\frac{1}{p_1})+1]$

所以$f(N)=f(p_1^{a_1})\cdot f(p_2^{a_2})\cdot \dots \cdot f(p_n^{a_n})$
=$N\cdot [1+a_1(1-\frac{1}{p_1})]\cdot \dots \cdot [1+a_n(1-\frac{1}{p_n})]$
=$N\sum _{i=1}^n(1+a_i(1-\frac{1}{p_i}))$

解题代码如下，代码中为最大程度减少循环次数，还运用了定理：对任意正整数n，最多只有一个素因子大于$\sqrt{n}$，且如果有素因子大于$\sqrt{n}$，那么这个素因子的次数一定是1

#include <iostream>
#include <cstring>
#include <cmath>
using namespace std;

int main()
{
    int n;
    long long ans;
    while(cin>>n) {
        ans=n;
        int b=sqrt(1.0*n);//sqrt的形参是double类型的。
        for(int i=2; i<=b; i++)
        {
            if(n%i==0)//说明i是其中一个素因子
            {
                int a=0;
                while(n%i==0)//计算对应这个素因子的次数
                {
                    n/=i;
                    a++;
                }
                ans=ans+ans*a*(i-1)/i;//技巧：像此题除法一定要保证是整除
            }
        }
        if(n!=1)
            ans=ans+ans*(n-1)/n;//此处是在判断是否存在大于根号n的质因子
        cout<<ans<<endl;
    }
    return 0;
}