POJ 2478 Farey Sequence（欧拉函数）

Description
给出n，求小于等于n的数中互素的整数对数
Input
多组输入，每组用例输入一整数n，以n=0结束输入
Output
对于每组输入，输出小于等于n的数中互素的整数对数
Sample Input
2
3
4
5
0
Sample Output
1
3
5
9
Solution
问题转化为前n项欧拉函数求和，注意此时euler[1]=0
Code

#include<cstdio>#include<iostream>#include<algorithm>#include<cstring>#include<cmath>#include<queue>#include<stack>#include<map>#include<vector>#include<string>using namespace std;#define maxn 1000001#define INF 0x3f3f3f3ftypedef long long ll;int euler[maxn],prime[maxn],cnt;void Get_Euler(int n)//求n以内所有数的欧拉函数值 {    memset(euler,0,sizeof(euler));    euler[1]=1;    cnt=0;     for(int i=2;i<=n;i++)    {        if(!euler[i])        {            euler[i]=i-1;            prime[cnt++]=i;        }         for(int j=0;j<cnt&&prime[j]*i<=n;j++)        {            if(i%prime[j])                 euler[prime[j]*i]=euler[i]*(prime[j]-1);            else             {                euler[prime[j]*i]=euler[i]*prime[j];                break;            }        }    }}ll dp[maxn]={0ll};int main(){    Get_Euler(maxn);    euler[1]=0;//注意此处第一项欧拉函数值为0     for(int i=1;i<maxn;i++)        dp[i]=dp[i-1]+euler[i];    int n;    while(scanf("%d",&n),n)        printf("%lld\n",dp[n]);    return 0;}