51nod1040最大公约数之和

题目大意：

给出一个n，求1-n这n个数，同n的最大公约数的和。比如：n = 6

1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6，加在一起 = 15

Input

1个数N(N <= 10^9)

Output

公约数之和

思路：将n分解质因子，每个质因子贡献的是多少，累加求和。每个质因子贡献的次数可以用欧拉函数来求。若n%i==0 ， 且n/i！=i，则ans += i*（Euler(n/i)) + (n/i) * Euler(i);

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <set>
#include <iomanip>

using namespace std;
//#pragma comment(linker, "/STACK:102400000,102400000")
#define maxn 1000005
#define MOD 1000000007
#define mem(a , b) memset(a , b , sizeof(a))
#define LL long long
#define INF 100000000

int n;
int Euler(int n)
{
    int res = n ;
    for(int i = 2 ; i * i <= n ; i ++ )
    {
        if(n % i == 0)
        {
            n /= i;
            res -= res / i ;
            while(n % i == 0)
                n /= i;
        }
    }
    if(n > 1)
        res -= res / n  ;
    return res;
}

int main()
{
    while(scanf("%d" , &n) != EOF)
    {
        LL ans = 0;
        for(int i = 1 ; i * i <= n ; i ++)
        {
            if(i * i == n) ans += i*Euler(i);
            else if(n % i == 0)
            {
                ans = ans + i * Euler(n/i) + (n/i) * Euler(i);
            }
        }
        cout << ans << endl;
    }
    return 0;
}