普通素数筛选与快速素数筛选模板

1.普通素数筛选模板

最普通的方法 

      说明：解决这个问题的诀窍是如何安排删除的次序，使得每一个非质数都只被删除一次。 中学时学过一个因式分解定理，他说任何一个非质（合）数都可以分解成质数的连乘积。例如，16=2^4，18=2 * 3^2，691488=2^5 * 3^2 * 7^4等。如果把因式分解中最小质数写在最左边，有16=4^2，18=2*9,691488=2^5 * 21609,；换句话说，把合数N写成N=p^k * q，此时q当然是大于p的，因为p是因式分解中最小的质数。由于因式分解的唯一性，任何一个合数N，写成N=p^k * q;的方式也是唯一的。 由于q>=p的关系，因此在删除非质数时，如果已知p是质数，可以先删除p^2,p^3,p^4,... ，再删除pq,p^2*q,p^3*q,...,（q是比p大而没有被删除的数），一直到pq>N为止。

      因为每个非质数都只被删除一次，可想而知，这个程序的速度一定相当快。依据Gries与Misra的文章，线性的时间，也就是与N成正比的时间就足够了（此时要找出2N的质数）。

#include <iostream>
#include<cstring>
using namespace std;

const int nmax=100;
int prime[nmax];//素数数组
int mark[nmax];//全体数的标记数组 

int Prime(int n){
	int index=0;//index记录素数个数
	memset(mark,0,sizeof(mark));
	//prime[index++]=2; //素数数组中第一个素数为2，更新index
	for(int i=2;i<n;i++){
		if(mark[i]!=1){//如果此时序列周未标记 
		 	prime[index++]=i;//此时未标记序列中最小的那个是素数
			for(int j=2*i;j<n;j+=i){//将该数的整数倍都标记 
				mark[j]=1;
			}
		}
	} 
	return index;
}
int main(int argc, char** argv) {
	int ans=Prime(100);
    printf("%d\n",ans);
	int cnt=0;
	/*for(int i=2;i<nmax;i++){
		if(mark[i]!=1){
			printf("%d ",i);
			cnt++;
		}
	}
	printf("\n%d\n",cnt);*/
	return 0;
}

2.快速素数筛选

      这种方法经证明只需要2n的复杂度，但是因为有*,/,%等运算，所以只比一般素数筛选快3倍左右。
关于复杂度的证明，用整体的思路就很好证明，算法由一层循环O(n)再加上所有产生的合数，可以证明每个合数都会被产生一次且仅一次。所以这个算法每一次运算都产生了一个合数（真是牛逼），没有丝毫重复的操作。最终复杂度为O(n)

#include <iostream>
#include<cstring>

using namespace std;
const int nmax=100;
int Mark[nmax];
int prime[nmax];

//判断是否是一个素数  Mark 标记数组 index 素数个数
int Prime(int n){
    int index = 0;
    memset(Mark,0,sizeof(Mark));
    for(int i = 2; i <= n; i++)
    {
        //如果未标记则得到一个素数
        if(Mark[i] == 0){
            prime[index++] = i;
        }
        
        //标记目前得到的素数的i倍为非素数
        for(int j = 0; j < index && (long long)prime[j] * i <=n; j++)
        {
            Mark[i * prime[j]] = 1;
            if(i % prime[j] == 0){//如果倍数i能整除已有的素数，则已经标记过了 
                break;
            }
        }
    }
    return index;
}
int main(int argc, char** argv) {
	int ans=Prime(100);
	printf("%d\n",ans);
	int cnt=0;
	for(int i=2;i<=100;i++){
		if(Mark[i]!=1){
			printf("%d ",i);
			cnt++;
		}
	}
	printf("\n%d\n",cnt);
	return 0;
}