快速线性筛法求素数&&一般筛法求素数代码

来源:互联网 发布:淘宝下载2016官方 编辑:程序博客网 时间:2024/06/06 02:24

一般筛法求素数代码

思路:素数的倍数不是素数。

#include<iostream>using namespace std;const int M = 2000000;int a[M] = {1,1};int prime[M]={0};int k=0;void init(){    for(int i=2;i<=M;i++)    {        if(!a[i])        {            prime[k++] = i;            for(int j=i+i;j<=M;j+=i)                    a[j]= 1;        }    }}int main(){    int n;    while(cin>>n)    {        for(int i=0;i<n;i++)        cout<<prime[i]<<" ";        cout<<endl;    }}

快速线性筛法

原理:1.一个合数是由n个素数的乘积所组成的

2.素数的倍数不是素数。

建议无论手模拟一个两个过程 就明白谁原理和优势了。

#include<iostream>using namespace std;const int Max=2000000;long long prime[Max] = {0};int k=0;int a[Max]={1,1};void init(){    for(long long i=2;i<Max;i++)    {        if(!a[i])//        prime[k++] = i;        for(long long j=0;j<k&&i*prime[j]<Max;j++)        {            a[i*prime[j]] = 1;            if(i%prime[j]==0)//                break;        }    }}int main(){    int n;    init();    while(cin>>n)    {        for(int i=0;i<n;i++)       cout<<prime[i]<<" ";        cout<<endl;    }}
推荐:http://blog.csdn.net/dinosoft/article/details/5829550

0 0