2017.2.25【初中部 提高组】模拟赛B组 倒霉的小C(beats) 题解

原题：

http://172.16.0.132/senior/#contest/show/1929/1

题目描述：

小G最近迷上了岛国动漫《Angel Beats》，她为了画出一个更霸气的Angel Beats的logo，想了如下办法：
从（0，0）开始，画到（n，1），再从（n，1），画到（2*n，-1），再到（3*n，2），再到（4*n，-2），依此类推，即每次画出一个（n，(-1)^(i+1)*i）的向量，一共画出n个这样的向量。现在小G想让小C求出这个图形穿过了多少格点（坐标都是整数）。
由于小C想要认真地听他的数学课并且想自己在接力赛中因RP暴光而发生接力棒传错这类的糗事，所以这个问题就交给你啦。小G说，如果连你也解决不好，就把你的RP也吸光。

输入：

输入文件中仅一行为一个整数n。

输出：

输出文件中仅一行为一个数，表示穿过的格点数。

样例输出：

4

样例输出：

9

数据范围限制：

对于20%的数据，满足：n<=100；
对于40%的数据，满足：n<=10^5；
对于60%的数据，满足：n<=10^7；
对于100%的数据，满足：n<=10^10。

分析：

通过简单观察可以发现，每次画出向量（n，i）经过的格点个数为gcd(i，n)，那么答案就等于Ans=1+
直接求解的时间复杂度是O(n)的。
那么Ans=1+，其中d为n的约数。fai(n)表示1～n中与n互质的数的个数。通过这样的变形，我们就可以得到时间复杂度为O(C*sqrt(n))的算法，C为n的约数个数。

实现：

#include<cstdio>    
#include<cmath>   

long long i,n,ans,k;
long long dg(long long x)  
{  
    long long k=x,i=2;  
    for(i=2;x>1;i++)
        if(x%i==0)
        {
            k-=k/i;
            while(x%i==0) x/=i;
        }
    return k;  
}  
int main()  
{
    freopen("beats.in","r",stdin);freopen("beats.out","w",stdout);
    scanf("%lld",&n);
    for(i=1;i<=trunc(sqrt(n));i++)  
        if(n%i==0)   
        {  
              k=i;  
              ans+=k*dg(n/k);  
              if(k!=n/k)  
              {  
                    k=n/i;  
                    ans+=k*dg(n/k);  
              }  
        }  
    printf("%lld\n",ans+1);
}