素数筛法求素数

素数筛类似于打表标记,预先处理掉非素数的数,即素数的倍数(任意非素数都可以由几个素数相乘得到),于是效率比暴力求解快得多。

埃氏筛法的效率为O(n loglog n),简单易懂,但是会重复标记,比如当i为2时,6会被标记掉,然而当i为3时,6又会被重复标记,这样的重复访问加大了时间复杂度,于是有了欧拉筛。

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<ctime> 
 5 using namespace std;
 6 bool v[100000055];
 7 void ss(int n)
 8 {
 9     memset(v,true,sizeof(v));//全定为素数 
10     for(int i=2;i<=sqrt(n);i++)
11     {
12         if(v[i])
13         {
14             for(int j=i*i;j<=n;j+=i)//筛掉i的倍数 
15             v[j]=0;
16         }
17     }
18 }
19 int main()
20 {
21      clock_t start,finish;
22    double totaltime;
23    
24     int n;
25     cin>>n;
26     start=clock();
27     ss(n);
28     int ans=0;
29     for(int i=2;i<=n;i++)
30     {
31         if(v[i])
32         ans++;
33     }
34     cout<<ans<<endl;
35       finish=clock();
36    totaltime=(double)(finish-start)/CLOCKS_PER_SEC;
37    cout<<"\n此程序的运行时间为"<<totaltime<<"秒!"<<endl;
38    return 0;
39 }

欧拉筛的时间复杂度为o(n),就是线性。注释为本人拙见,水平有限见谅。

 1 #include<iostream>
 2 #include<cstring>
 3 using namespace std;
 4 int v[10000000],f[10000000];
 5 int ss(int n)
 6 {
 7     memset(v,0,sizeof(v));//用0表示素数,1表示非素数,先全部标记为素数 
 8     int cnt=0;
 9     for(int i=2;i<n;i++)
10     {
11         if(!v[i])//当这个i为素数时 
12 13             f[cnt++]=i;//存素数 
14             for(int j=0;j<cnt&&i*f[j]<n;j++)
15             {
16                 v[i*f[j]]=1;//筛掉素数i的倍数,任意非质数都可以由质数相乘得到 
17                 if(i%f[j]==0)//当素数i是前几个素数中某个数的倍数时,f[j+1]*i已经被f[j]乘以某个数筛掉了
18                 break;
19             }
20 21     }
22     return cnt;//n以内素数的个数 
23 }
24 int main()
25 {
26     int n;
27     cin>>n;
28     int cnt=ss(n); 
29     for(int i=0;i<cnt;i++)
30     cout<<f[i]<<' ';
31     cout<<endl;
32 }

 

posted @ 2019-03-14 21:21  魔法少女郭德纲  阅读(163)  评论(0编辑  收藏  举报