bzoj 2705: [SDOI2012]Longge的问题

2705: [SDOI2012]Longge的问题
Time Limit: 3 Sec Memory Limit: 128 MB
Submit: 2295 Solved: 1398
[Submit][Status][Discuss]
Description
Longge的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。现在问题来了：给定一个整数N，你需要求出∑gcd(i, N)(1<=i <=N)。
Input
一个整数，为N。
Output
一个整数，为所求的答案。
Sample Input
6
Sample Output
15
HINT

【数据范围】

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

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

---

【分析】

极其鬼畜的一道题，理解起来巨困难

考虑转化问题，从1到n的范围内，有许多个i的gcd(i,n)等于同一个n的因数。我们可以枚举n的每一个因数k，累计“以该数k为解的gcd(i,n)的个数s(k)”乘以该数，就能得到答案。

若有gcd(n,m)=k，那么n和m同除公约数k后，可以得到gcd(n/k,m/k)=1。由前式可知(m/k)与(n/k)互质。满足条件的(m/k)个数，也就是s(k)就等于phi(n/k) 即欧拉函数

---

【代码】

//2705: [SDOI2012]Longge的问题
//需要用到公式法求欧拉函数
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define LL long long
#define fo(i,j,k) for(long long i=j;i<=k;i++)
using namespace std;
LL n,m,ans;
LL fai(LL x)
{
    LL a=x;
    fo(i,2,m)
      if(x%i==0)
      {
        a=a/i*(i-1);
        while(x%i==0) x=x/i;
      }
    if(x>1) a=a/x*(x-1);
    return a;
}
int main()
{
    scanf("%lld",&n);
    m=sqrt(n);
    fo(i,1,m)
      if(n%i==0)
      {
        ans+=fai(n/i)*i;
        ans+=fai(i)*(n/i);
      }
    printf("%lld\n",ans);
    return 0;
}

（需要用欧拉函数的公式，还要进行细节处理）