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万以内也很快:


  • 相关阅读:
    supervisor 安装与使用
    CF39C-Moon Craters【dp】
    NWERC2020J-Joint Excavation【构造,贪心】
    CF25E-Test【AC自动机,bfs】
    CF19E-Fairy【树形结构,差分】
    CF11D-A Simple Task【状压dp】
    CF5E-Bindian Signalizing【单调栈】
    P6628-[省选联考 2020 B 卷] 丁香之路【欧拉回路,最小生成树】
    CF666E-Forensic Examination【广义SAM,线段树合并】
    CF235D-Graph Game【LCA,数学期望】
  • 原文地址:https://www.cnblogs.com/linruier/p/9485191.html
Copyright © 2011-2022 走看看