思路:
在原问题中,设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;
}