zoukankan      html  css  js  c++  java
  • 素数筛法求素数

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

    埃氏筛法的效率为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<<"
    此程序的运行时间为"<<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 }
  • 相关阅读:
    rocketMQ配置事故
    微信网页授权问题记录
    记一次Spring配置事故
    文件转换
    Java对象空间分配流程
    mysql(六)索引的数据结构
    mysql(五)查询缓存
    mysql(四)log
    从项目中加载文件
    cookie
  • 原文地址:https://www.cnblogs.com/spongeb0b/p/10533574.html
Copyright © 2011-2022 走看看