欧拉函数线性筛法

2018-06-17 22:12:28来源:未知 阅读 ()

新老客户大回馈,云服务器低至5折

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cmath>
 5 #define lli long long int 
 6 using namespace std;
 7 const int MAXN=10000001;
 8 void read(int &n)
 9 {
10     char c='+';int x=0;bool flag=0;
11     while(c<'0'||c>'9'){c=getchar();if(c=='-')flag=1;}
12     while(c>='0'&&c<='9')
13     x=(x<<1)+(x<<3)+c-48,c=getchar();
14     flag==1?n=-x:n=x;
15 }
16 int n,m;
17 bool check[MAXN];
18 int prime[MAXN];
19 int phi[MAXN];
20 int tot=0;
21 int  main()
22 {
23            cin>>n;
24            phi[1]=1;
25         for(int i=2;i<=n;i++)
26         {
27             if(!check[i])
28                 prime[++tot]=i,phi[i]=i-1;// 只有i与它互质 
29             for(int j=1;j<=tot;j++)
30             {
31                 if(i*prime[j]>n)
32                     break;
33                 check[i*prime[j]]=1;
34                 if(i%prime[j]==0)
35                 {
36                     phi[i*prime[j]]=phi[i]*prime[j]; 
37                     break;
38                 }
39                 else
40                     phi[i*prime[j]]=phi[i]*(prime[j]-1);
41             }
42         }
43         printf("%d\n",phi[n]);
44     return 0;
45 }

注意 if(i%prime[j]==0)

不要写成if(!i%prime[j])

标签:

版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有

上一篇:04:最长公共子上升序列

下一篇:Common Subsequence