POJ 2478 Farey Sequence（n级法雷数列的元素个数）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <climits>
#include <cmath>
#include <ctime>
#include <cassert>
#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
using namespace std;
typedef long long ll;
const int MAX_N = 1000010;

int n;
ll euler[MAX_N], f[MAX_N];

void GetEuler()
{
memset(euler, 0, sizeof(euler));
euler[1] = 1;
for(int i = 2; i < MAX_N; i++){
if(euler[i] == 0){
for(int j = i; j < MAX_N; j += i){
if(euler[j] == 0) euler[j] = j;
euler[j] = euler[j] / i * (i - 1);
}
}
}
}

void init()
{
GetEuler();
f[2] = 1;
for(int i = 3; i < MAX_N; i++){
f[i] = f[i - 1] + euler[i];
}
}

int main()
{
init();
while(~scanf("%d", &n) && n){
printf("%lld\n", f
);
}
return 0;
}