poj 2478 Farey Sequence 欧拉函数

题目链接：http://poj.org/problem?id=2478

/*
* 题目：
*       求0<a<b<=n中，所有的最简分数a/b的个数
*
* 分析：
*      先化成求sigma( phi(n) )，然后转化为快速求欧拉函数。
*
* */

#include <cstdio>
#include <cstring>
#include <cmath>
#include <iostream>

using namespace std;

const int X = 1000005;

typedef long long ll;

ll ans[X];
bool use[X];

void init(){
	memset(use,false,sizeof(use));
	for(int i=1;i<X;i++)
		ans[i] = i;
	for(int i=2;i<X;i++)
		if(!use[i]){
			for(int j=i;j<X;j+=i){
				use[j] = true;
				ans[j] = ans[j]/i*(i-1);
			}
		}
	ans[2] = 1;
	for(int i=3;i<X;i++)
		ans[i] += ans[i-1];
}

int main(){
#ifndef ONLINE_JUDGE
	freopen("sum.txt","r",stdin);
#endif
	int n;
	init();
	while(scanf("%d",&n),n)
		printf("%lld\n",ans[n]);

	return 0;
}