传送门
Description
Longge的数学成绩非常好,并且他非常乐于挑战高难度的数学问题。现在问题来了:给定一个整数N,你需要求出∑gcd(i, N)(1<=i <=N)。
Input
一个整数,为N。
Output
一个整数,为所求的答案。
Sample Input
6
Sample Output
15
HINT
【数据范围】
对于60%的数据,0
题解
这道题一开始自己yy了好多种做法,然而在对拍的时候都出错了QAQ
于是只好去膜了一发黄学长的博客
题解实际上很简单,但又非常巧妙
设f(k)为所有i(1≤i≤N)且gcd(i,N)=k的i的个数,那么题目所求的就可以表示为
∑i=1Ni×f(i)
若gcd(x,N)=k,则
gcd(xk,Nk)=1
即求 Nk 的欧拉函数
枚举i的时候只需枚举到 N−−√ 就可以了
CODE:
#include<cmath>
#include<cstdio>
typedef long long ll;
ll N,n,ans;
inline ll phi(ll N)//注意是ll N 不要传进来int,否则会得95分
{
ll ans=N;
for(int i=2;i<=n;i++)
if(N%i==0)
{
ans=ans/i*(i-1);
while(N%i==0) N/=i;
}
if(N>1) ans=ans/N*(N-1);
return ans;
}
int main()
{
scanf("%lld",&N);
n=sqrt(N);
for(int i=1;i<=n;i++)
if(N%i==0)
{
ans+=1ll*i*phi(N/i);
if(i*i<N) ans+=1ll*(N/i)*phi(i);
}
printf("%lld",ans);
return 0;
}
总结
遇到需要仔细研究的题目的时候思路还是不准确,有时候正解其实很简单,但是还需要巧妙的转化。
可以将题目中出现有规律的东西提出来观察,比如说这道题中设出的f(k);可以其中不变的因素提出来,然后就可以发现一些美妙的规律,比如说这道题中最后的
gcd(xk,Nk)=1
,于是发现可以用欧拉函数求解。