bzoj2705 [SDOI2012]Longge的问题

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2705

题目大意：给定一个N，要求sigma(gcd(i,N))的值，1<=i<=N。

题目分析：我们枚举gcd(i,N)的值，即，枚举N的约数，时间复杂度为sqrt(N),然后对于每一个枚举值p，我们计算phi(N/p)，由于N的值比较大，我们不能预处理，所以我们就暴力在sqrt(N)的时间内求出phi(N/p)，那么gcd(i,N)为p的那一部分的和就等于p*phi(N/p)。我们把这些结果累加起来即可。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#define maxn 11000000
using namespace std;
typedef long long LL;
LL d;
int n;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
LL phi(LL x)
{
    LL t=x;
    for(LL i=2;i<=d;i++)
        if(x%i==0)
        {
            t=t/i*(i-1);
            while(x%i==0)x/=i;
        }
    if(x>1)t=t/x*(x-1);
    return t;
}
int main()
{
    scanf("%d",&n); 
    LL ans=0;
	d=sqrt(n);
	for (int i=1;i<=d;i++)
	{
	    if (n%i==0)
		{
		    ans+=phi(i)*n/i;
			if (i!=n/i) ans+=phi(n/i)*i;	
		}	
	}
	cout<<ans<<endl;
	return 0;
}