【POJ2478】Farey Sequence

【题目大意】

求在[1,n]范围内的(a,b)且a≠b，满足gcd(a,b)=1的个数

 

【题解】

欧拉函数裸题，无需多言。

 1 /*************
 2   POJ 2478
 3   by chty
 4   2016.11.3
 5 *************/
 6 #include<iostream>
 7 #include<cstdio>
 8 #include<cstring>
 9 #include<cstdlib>
10 #include<ctime>
11 #include<cmath>
12 #include<algorithm>
13 using namespace std;
14 #define MAXN 1000010
15 int n,phi[MAXN],check[MAXN],prime[MAXN];
16 inline int read()
17 {
18     int x=0,f=1;  char ch=getchar();
19     while(!isdigit(ch))  {if(ch=='-')  f=-1;  ch=getchar();}
20     while(isdigit(ch))  {x=x*10+ch-'0';  ch=getchar();}
21     return x*f;
22 }
23 void get()
24 {
25     phi[1]=1;  int cnt=0;
26     for(int i=2;i<=MAXN;i++)
27     {
28         if(!check[i])  {prime[++cnt]=i;  phi[i]=i-1;}
29         for(int j=1;j<=cnt&&prime[j]*i<=MAXN;j++)
30         {
31             check[i*prime[j]]=1;
32             if(i%prime[j])  phi[i*prime[j]]=phi[i]*(prime[j]-1);
33             else {phi[i*prime[j]]=phi[i]*prime[j];  break;}
34         }
35     }
36 }
37 int main()
38 {
39     freopen("cin.in","r",stdin);
40     freopen("cout.out","w",stdout);
41     get();
42     while(scanf("%d",&n)&&n)
43     {
44         long long ans=0;
45         for(int i=2;i<=n;i++)  ans+=phi[i];
46         printf("%lld\n",ans);
47     }
48     return 0;
49 }

 

posted @ 2016-11-03 14:41  chty  阅读(174)  评论(0)  编辑  收藏  举报