JZOJ 3509. 【NOIP2013模拟11.5B组】倒霉的小C

目录：

---

题目：

单击查看题目

---

分析：

【算法分析与设计】
通过简单观察可以发现，每次画出向量（n，i）经过的格点个数为gcd(i，n)，那么答案就等于

直接求解的时间复杂度是O(n)的。
那么，其中d为n的约数。fai(n)表示1～n中与n互质的数的个数。通过这样的变形，我们就可以得到时间复杂度为$O(C*sqrt(n))$的算法，C为n的约数个数。

---

代码：

#pragma GCC optimize("3")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<map>
#define LL long long
using namespace std;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
LL phi(LL x)
{
    LL k=x;
    LL i=2;
    while(x>1)
    {
        if(x%i==0)
        {
            k-=k/i;
            while(x%i==0) x/=i;
        }
        i++;
    }
    return k;
}
LL ans=0;
int main()
{
//  freopen("beats.in","r",stdin);
//  freopen("beats.out","w",stdout);
    LL n=read();
    for(LL i=1;i<=trunc(sqrt(n));i++)
    {
        if(n%i==0)
        {
            LL k=i;
            ans+=k*phi(n/k);
            if(n/k!=k)
            {
                k=n/i;
                ans+=k*phi(n/k);
            }
        }
    }
    printf("%lld",ans+1);
    fclose(stdin);
    fclose(stdout); 
    return 0;
}