产生MAXN以内的所有素数代码,C/C++代码库,德仔网

代码语言

代码分类

【C/C++】产生MAXN以内的所有素数

作者:首君 / 发布于2016/2/23/ 626

//产生MAXN以内的所有素数
//note:2863311530就是10101010101010101010101010101010
//给所有2的倍数赋初值
#include <cmath>
#include <iostream>
using namespace std;
#define MAXN 100000000
unsigned int plist[6000000],pcount;
unsigned int isprime[(MAXN>>5)+1];
#define setbitzero(a) (isprime[(a)>>5]&=(~(1<<((a)&31))))
#define setbitone(a) (isprime[(a)>>5]|=(1<<((a)&31)))
#define ISPRIME(a) (isprime[(a)>>5]&(1<<((a)&31)))
void initprime(){
    int i,j,m;
    int t=(MAXN>>5)+1;
    for(i=0;i<t;++i)isprime[i]=2863311530;
    plist[0]=2;setbitone(2);setbitzero(1);
    m=(int)sqrt(MAXN);
    for(pcount=1,i=3;i<=m;i+=2)
        if(ISPRIME(i))
            for(plist[pcount++]=i,j=i<<1;j<=MAXN;j+=i)
                setbitzero(j);
    if(!(i&1))++i;
    for(;i<=MAXN;i+=2)if(ISPRIME(i))plist[pcount++]=i;
}

试试其它关键字

　素数　

同语言下

可能有用的

首君贡献的其它代码(39)