根据题目要求,我们要求出 ∑ gcd(i,N) 。初看这题,仿佛很难下手,因为N的范围是到Max_int,答案已经达到longlong型,不可能全部枚举。 我们不妨枚举N的约数,这样枚举的效率是O (log n)。然后我们要求出针对每一个约数k, gcd(i,N)=k的答案数。
那么我们通过如下过程:
∵gcd(i,N)=k ∴ gcd(i/k,N/k)=1
又∵i>0,i<=N
∴ 我们用S(k)表示gcd(i/k,N/k)=1这样的i的个数
则S(k)即 满足1<=i<=N/k的与N/k互素的i总数
那么最终的答案Ans = ∑(k*S(k))
我们用欧拉函数求出每一个约数k的gcd(i,N)=k的答案数,就能得到答案
#include "iostream"
using namespace std;
typedef long long ll;
ll sum=0,n;
inline ll Eular (int x){
ll ans=x;
for (int i=2;i*i<=x;i++)
{ if(x%i==0) ans=ans/i*(i-1);
while(x%i==0) x/=i;}
if(x>1) ans=ans/x*(x-1);
return ans;
}
int main(){
int i; cin>>n;
for (i=1;i*i<=n;i++)
if(n%i==0){
sum+=n/i*Eular(i);
if(i*i-n)
sum+=i*Eular(n/i);
}
cout<<sum<<endl;
return 0;
}