题目大意:
给出一个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;
}