bzoj2705 [SDOI2012]Longge的问题

传送门
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的个数，那么题目所求的就可以表示为

$$\sum_{i=1}^N i×f(i)$$
若gcd(x,N)=k，则
$$gcd(\frac{x}{k},\frac{N}{k})=1$$
即求 $\frac{N}{k}$的欧拉函数
枚举i的时候只需枚举到 $\sqrt{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(\frac{x}{k},\frac{N}{k})=1$，于是发现可以用欧拉函数求解。