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

    素数筛法听起来很高大上的算法,实际上算法思想却十分简单:对前面取得的素数,依次"去掉"取出的这个素数后面的所有该素数的倍数。最后剩下的素数表里面的就是要求的素数,效率很高。

    #include<iostream>
    using namespace std;
    const int MAX_N=100000000;
    int prime[MAX_N];//素数表
    bool is_prime[MAX_N+1];
    //返回n以内的素数的个数
    int sieve(int n)
    {
        int p=0;
        for(int i=0;i<=n;i++)is_prime[i]=true;
        is_prime[0]=is_prime[1]=false;
        for(int i=2;i<=n;i++)
        {
            if(is_prime[i])
            {
                prime[p++]=i;//素数打表
                for(int j=2*i;j<=n;j+=i)is_prime[j]=false;//去掉已有素数的倍数
            }
        }
        return p;
    }
    int main()
    {
        int N;
        cout<<"求2到多少以内的素数个数?"<<endl;
        cin>>N;
        cout<<"2到"<<N<<"共有"<<sieve(N)<<"个素数"<<endl;
        return 0;

    }



    //或者让我们直接输出n以内的所有素数看一看:

    #include<iostream>
    #include<iomanip>
    using namespace std;
    const int MAX_N=100000000;
    int prime[MAX_N];
    bool is_prime[MAX_N+1];
    //返回以内的素数的个数
    int sieve(int n)
    {
        int p=0;
        for(int i=0;i<=n;i++)is_prime[i]=true;
        is_prime[0]=is_prime[1]=false;
        for(int i=2;i<=n;i++)
        {
            if(is_prime[i])
            {
                prime[p++]=i;//素数打表
                for(int j=2*i;j<=n;j+=i)is_prime[j]=false;//去掉已有素数的倍数
            }
        }
        return p;
    }
    int main()
    {
        int N;
        cout<<"求2到多少以内的素数?"<<endl;
        cin>>N;
        sieve(N);
        int cnt=0;
        for(int i=0;i<=N;i++)
        {
            if(is_prime[i])
            {cout<<setw(9)<<i;cnt++;if(cnt%8==0)cout<<endl;}
        }
        return 0;
    }


    甚至求解10万以内也很快:


  • 相关阅读:
    wait与sleep区别?
    oracle死锁查询
    atomic 原子操作的类
    买票问题
    0001.第一个多线程demo--分批处理数据
    01: JavaScript实例
    01: 运维工作梳理
    04: 使用BeautifulSoup封装的xss过滤模块
    04: 打开tornado源码剖析处理过程
    03: 自定义异步非阻塞tornado框架
  • 原文地址:https://www.cnblogs.com/linruier/p/9485191.html
Copyright © 2011-2022 走看看