beats

3509. 【NOIP2013模拟11.5B组】倒霉的小C(beats) 
(File IO): input:beats.in output:beats.out

Time Limits:  1000 ms  Memory Limits:  262144 KB  Detailed Limits  

Goto ProblemSet

Description

小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也吸光。

Input

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

Output

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

Sample Input

4

Sample Output

9

Data Constraint

对于20%的数据，满足：n<=100；

对于40%的数据，满足：n<=10^5；

对于60%的数据，满足：n<=10^7；

对于100%的数据，满足：n<=10^10。

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
long long a,b,d,ans;
long long fai(long long);
int main()
{
	freopen("beats.in","r",stdin);
	freopen("beats.out","w",stdout);
	scanf("%lld",&a);
	for (b=1;b<=trunc(sqrt(a));b++)
	{
		if (a%b==0)
		{
			d=b;
			ans=ans+d*fai(a/d);
			if (d!=a/d) 
			{
				d=a/b;
				ans=ans+d*fai(a/d);
			}
		}
	}
	ans++;
	printf("%lld",ans);
	return 0;
}
long long fai(long long x)
{
	int i;
	long long t;
	t=x;
	i=2;
	while (x>1) 
	{
		if (x%i==0) 
        {
			t=(t-t/i);
			while (x%i==0) 
    		{
				x=x/i;
			}
		}
		i++;
	}
	return t;
}