poj 2478 Farey Sequence

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

题目大意：给定一个整数 n ( 2≤n≤10⁶ ), 求 phi(2)+phi(3)+···+phi(n).

解题思路：简单题，主要是求 1-10⁶ 的欧拉函数 phi(i), 方法与素数筛法类似。

///////////////////////////////////////////////////////////////////////////
//problem_id: poj 2478
//user_id: SCNU20102200088
///////////////////////////////////////////////////////////////////////////

#include <algorithm>
#include <iostream>
#include <iterator>
#include <iomanip>
#include <cstring>
#include <cstdlib>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cmath>
#include <queue>
#include <stack>
#include <list>
#include <set>
#include <map>
using namespace std;

///////////////////////////////////////////////////////////////////////////
typedef long long LL;
const double PI=acos(-1.0);

const int x4[]={-1,0,1,0};
const int y4[]={0,1,0,-1};
const int x8[]={-1,-1,0,1,1,1,0,-1};
const int y8[]={0,1,1,1,0,-1,-1,-1};

typedef int T;
T max(T a,T b){ return a>b? a:b; }
T min(T a,T b){ return a<b? a:b; }
///////////////////////////////////////////////////////////////////////////

///////////////////////////////////////////////////////////////////////////
//Add Code:
const int maxn=1000000;
LL phi[maxn+5],ans[maxn+5];
///////////////////////////////////////////////////////////////////////////

int main(){
    ///////////////////////////////////////////////////////////////////////
    //Add code:
    int n,i,j;
    for(i=1;i<=maxn;i++) phi[i]=i;
    for(i=2;i<=maxn;i++){
        if(phi[i]==i){
            for(j=i;j<=maxn;j+=i) phi[j]=phi[j]/i*(i-1);
        }
    }
    ans[2]=phi[2];
    for(i=3;i<=maxn;i++) ans[i]=ans[i-1]+phi[i];
    while(scanf("%d",&n)!=EOF){
        if(n==0) break;
        printf("%lld
",ans[n]);
    }
    ///////////////////////////////////////////////////////////////////////
    return 0;
}

///////////////////////////////////////////////////////////////////////////
/*
Testcase:
Input:
2
3
4
5
0
Output:
1
3
5
9
*/
///////////////////////////////////////////////////////////////////////////