素数筛选模板

1.普通素数筛选模板

最普通的方法。。。

//判断是否是一个素数  Mark 标记数组 index 素数个数
int Prime(){
    int index = 0;
    memset(Mark,0,sizeof(Mark));
    prime[index++]=2;
    for(int i = 3;i < MAXSIZE;i+=2){
        if(Mark[i] != 1)
        {
            prime[index++] = i;
            for(int j = i+i;j < MAXSIZE;j += i){
                Mark[j] = 1;
            }
        }
    }
    return index;
}

 2.快速素数筛选

这种方法经证明只需要2n的复杂度,但是因为有*,/,%等运算,所以只比一般素数筛选快3倍左右。

关于复杂度的证明,用整体的思路就很好证明,算法由一层循环O(n)再加上所有产生的合数,可以证明每个合数都会被产生一次且仅一次。

所以这个算法每一次运算都产生了一个合数(真是牛逼),没有丝毫重复的操作。最终复杂度为O(n)

int Mark[MAXSIZE];
int prime[MAXSIZE];

//判断是否是一个素数  Mark 标记数组 index 素数个数
int Prime(){
    int index = 0;
    memset(Mark,0,sizeof(Mark));
    for(int i = 2; i < MAXSIZE; i++)
    {
        //如果未标记则得到一个素数
        if(Mark[i] == 0){
            prime[index++] = i;
        }
        //标记目前得到的素数的i倍为非素数
        for(int j = 0; j < index && (long long)prime[j] * i < MAXSIZE; j++)
        {
            Mark[i * prime[j]] = 1;
            if(i % prime[j] == 0){
                break;
            }
        }
    }
    return index;
}

 

posted @ 2016-03-01 17:20  chenhuan001  阅读(504)  评论(0编辑  收藏  举报