欧拉函数——龙哥的问题（暴力）

传送门：221. 龙哥的问题 - AcWing题库

思路：

在原问题中，设d=gcd(i,N),能发现，d一定是N的约数，同时能发现1=gcd(i/d,N/d),i/d与N/d互质。

其中i/d的数目就是N/d的欧拉值，而d的数目就是N的约数的数目。

转化一下得 ∑（1≤i≤N）gcd(i，N)  =  Σ（d|N）φ(N/d)∗(d)

至此，解题步骤如下

第一步：枚举n的约数d，这里因为约数总是成对出现，所以可以一起处理（n/d）

第二步：求n/d或n/(n/d)的欧拉值

第三步：答案累加φ(n/d)∗(d)和φ(d)∗(n/d)

注：求一个数的约数只需要枚举1~√N的范围，题目的2^31，实际只需要枚举

1~√21亿==14491.3767即可，而20亿以内的的数字约数最多的数的月数个数为1536，实际要算的就更少了。

代码：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<unordered_map>
using namespace std;
typedef pair<int,int>PII;
typedef long long ll;
const int N=1e7+10,M=1e7+10;

int get(int n)
{
     int sum=n;
    for(int i=2;(ll)i*i<=n;i++)
    {
        if(n%i==0)
        {
            sum=sum/i*(i-1);
            while(n%i==0) n/=i;
        }
    }
    if(n>1) sum=sum/n*(n-1);
    return sum;
}
int main()
{
    int n;
    scanf("%d",&n);
    ll ans=0;
    for(int i=1;(ll)i*i<=n;i++)
    {
        if(n%i==0)
        {
            ans+=get(i)*(n/i);
            if(i!=n/i) ans+=get(n/i)*i;
        }
    }
    cout<<ans<<endl;
    return 0;
}