素数表学习总结

素数表的打法多种多样，每种都有其特点，学习了几种，总结一下。

思路：

第四种：从小到大判断素数，再用素数的倍数是合数筛去合数，保留素数 ，重复性高；

第三种是根据一个合数的最小质因子<=sqrt(N) 来筛选；

第一种将不必要的访问去除，通过一个合数的最小素因子将其筛去，大大优化了时间复杂度；

如果需要将素数存起来，这时标记和保存数组就可以用一个了，素数0，1位置开始存都可以 prime[0]=2,prime[1]=3,都小于指针位置。
 

 第一种：自己写的，优点：适用于10^6的素数表，在时间上与第二种差别不大，空间上将vis数组省略

void primetable()//10^6:15ms 10^7: 200ms 10^8:1570ms
{
	int i=0,j=0,cnt=0;//memset(prime,0,sizeof(prime));
	for(i=2;i<MAXN;++i)//从2开始 
	{
		if(!prime[i])//1.找素数 
			prime[cnt++]=i;//找到覆盖前面的值 从0开始 满足2,3 
		for(j=0;j<cnt&&i*prime[j]<MAXN;++j)//2.筛素数 
		{                        //利用已知素数 筛选(i,i*i]的素数 
			prime[i*prime[j]]=1;//标记 
			if(i%prime[j]==0)  break;//利用最小的素因子筛	
		} 
	}
} 

第二种：优点：耗时最少，适用于10^7 时间比空间要求更紧时 可用此种方法替代第一种； 

 

const int MAXN=1e7+10;//167ms 10^6 13ms  10^8: 1287ms
int prime[MAXN];
bool vis[MAXN];//相对于第一种多开了一个标记数组 

void primetable()
{
	int i=0,j=0,cnt=0;
	memset(vis,false,sizeof(vis));// vis[0]=vis[1]=true;
	for(i=2;i<MAXN;++i)
	{
		if(!vis[i])
		{
			prime[++cnt]=i;
		}
		for(j=1;j<=cnt&&i*prime[j]<MAXN;++j)
		{
			vis[i*prime[j]]=true;
			if(i%prime[j]==0) break;
		 } 
	}
}

第三种：不如用1,2；

void primetable()//思路：算数基本定理：每个合数都存在最小质因子且小于或等于sqrt(n)
{
	int i=0,j=0,cnt=0;
	//memset(prime,0,sizeof(prime)); prime[0]=prime[1]=1;
	for(i=2;i*i<MAXN;++i)//先打出一个sqrt(MAXN)的素数表，用他来筛选MAXN范围的素数
	{
		if(!prime[i])
		{
			for(j=i*i;j<MAXN;j+=i)
			{
				prime[j]=1;
			}
		}
	}
	for(i=2;i<MAXN;++i)
	if(!prime[i]) prime[cnt++]=i;//printf("%d\n",i); 
}

第四种：前面的都是根据这一种的优化

void primetable()//思路：将所有素数的倍数筛去（重复性较高） 
{
	int i=0,j=0,cnt=0; 
    //memset(prime,0,sizeof(prime));	prime[0]=prime[1]=1;
	for(i=2;i<MAXN;++i)
	{
		if(!prime[i])
		{
			prime[cnt++]=i;//printf("%d\n",prime[cnt-1]);
			for(j=2*i;j<MAXN;j+=i)
		    prime[j]=1;
		}
	 } 
}

优化版详细讲解：http://blog.csdn.net/morewindows/article/details/7347459

题目：ZZULI 1437

测试代码：

#include<cstdio> //注：b版都是对a版空间的优化，其他都一样！
#include<cstring>
#include<time.h>
#include<algorithm>
using namespace std;

const int MAXN=1e5+10;//1000ms左右 10^6 120ms 
int prime[MAXN+333];
bool vis[MAXN];//相对于第一种多开了一个标记数组 

void primetable1_a()
{
	int i=0,j=0,cnt=0;//memset(prime,0,sizeof(prime));
	for(i=2;i<MAXN;++i)//从2开始 
	{
		if(!prime[i])//1.找素数 
			prime[cnt++]=i;//找到覆盖前面的值 从0开始 满足2,3 
		for(j=0;j<cnt&&i*prime[j]<MAXN;++j)//2.筛素数 
		{                        //利用已知素数 筛选(i,i*i]的素数 
			prime[i*prime[j]]=1;//标记 
			if(i%prime[j]==0)  break;//利用最小的素因子筛	
		} 
	}
} 

void primetable1_b()
{
	int i=0,j=0,cnt=0;
	memset(vis,false,sizeof(vis));// vis[0]=vis[1]=true;
	for(i=2;i<MAXN;++i)
	{
		if(!vis[i])
		{
			prime[++cnt]=i;
		}
		for(j=1;j<=cnt&&i*prime[j]<MAXN;++j)
		{
			vis[i*prime[j]]=true;
			if(i%prime[j]==0) break;
		 } 
	}
}
 
void primetable2_a()//10^6:16ms 10^7:183ms 10^8:3000ms 
{
	int i=0,j=0,cnt=0;
	memset(vis,false,sizeof(vis)); vis[0]=vis[1]=true;
	for(i=2;i*i<MAXN;++i)
	{
		if(!vis[i])
		{
			for(j=i*i;j<MAXN;j+=i)
		    vis[j]=true;
		}
	}
	for(i=2;i<MAXN;++i)
	if(!vis[i]) prime[++cnt]=i;
}

void primetable2_b()
{
	int i=0,j=0,cnt=0;
	//memset(prime,0,sizeof(prime)); prime[0]=prime[1]=1;
	for(i=2;i*i<MAXN;++i)
	{
		if(!prime[i])
		{
			for(j=i*i;j<MAXN;j+=i)
			{
				prime[j]=1;
			}
		}
	}
	for(i=2;i<MAXN;++i)
	if(!prime[i]) prime[cnt++]=i;//printf("%d\n",i); 
}

void primetable3_a()
{
	int i=0,j=0,cnt=0; memset(vis,false,sizeof(vis));
	for(i=2;i<MAXN;++i)
	{
		if(!vis[i])
		{
			prime[++cnt]=i;
			for(j=2*i;j<MAXN;++j)
		    vis[j]=true;
		}
	}
}

void primetable3_b()
{
	int i=0,j=0,cnt=0; 
    //memset(prime,0,sizeof(prime));	prime[0]=prime[1]=1;
	for(i=2;i<MAXN;++i)
	{
		if(!prime[i])
		{
			prime[cnt++]=i;//printf("%d\n",prime[cnt-1]);
			for(j=2*i;j<MAXN;j+=i)
		    prime[j]=1;
		}
	 } 
}


int main()
{
	int N,T,ans;
//	clock_t b1,b2,e1,e2;
//	b1=clock();
	primetable2_b();
//	e1=clock();
//	
//	b2=clock();
//	primetable1() ;
//	e2=clock();
//	printf("%d\n\n",e1-b1);
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&N);
		int pos=lower_bound(prime,prime+1000+1,N)-prime;// printf("%d\n",pos);
	    if(prime[pos]==N) 
		ans=N;
	    else
	    {
	    	int L=pos-1;
	    	if(prime[pos]-N==N-prime[L]) ans=prime[pos];
	    	else ans=prime[pos]-N<N-prime[L]?prime[pos]:prime[L];
		}
		printf("%d\n",ans);
	}
	return 0;
}