标记素数法(模板)+素数打表

#include <stdio.h>
#include <string.h>
#define N 3000000

int f[3000000];

int main()
{
	memset(f, 0, sizeof(f));
	int i, j;
	f[0]=1;
	f[1]=1;
	for(i=2; i<=N; i++)
	{
		if(f[i]==0)
		{
			for(j=i*2; j<=N; j+=i)
			{
				f[j]=1; //不是素数
			}
		}
	} // 打印所有发f[i]==0, 即为素数
	for(i=2; i<=N; i++)
	{
		if(f[i]==0)
			printf("%d  ", i );
	}
	return 0;
}

 

#include <stdio.h>
#include <string.h>
#include <math.h>
#define N 3000000

int c=0;
int f[3000000];
int prime[1000000];

int main()
{
	int i, j;

    memset(f, 0, sizeof(f));
	f[0]=1;
	f[1] =1;
	for(i=2; i<=N; i++)
	{
		if(f[i]==0)
		{
			prime[c++] = i;
			for(j=2*i; j<=N; j+=i )
			{
				f[j]=1;
			}
		}
	}
	for(i=0; i<=1000; i++)
	{
			printf("%d  ", prime[i] );
	}
	return 0;
}

 素数打表法二:

#include<stdio.h>  
#include<string.h>  

#define N 1000005

int prime[100005];  
int flag[1000005];  
int e;
  
void getP()  // 素数打表,找出素数存栈  
{  
	int i, j;  
	e = 0;  
	memset(flag, 0, sizeof(flag) ); //标记初始化
	
	for ( i=2; i<N; i++)  
	{  
		if ( flag[i]==0 ) 
		{
			prime[e++] = i; //进栈  
		}
		for ( j=0; j<e && i*prime[j]<N; j++ )  
		{  
			flag[ i * prime[j] ] = 1;  
		}  
	}  
}  

 

posted @ 2014-08-09 08:18  我喜欢旅行  阅读(816)  评论(0编辑  收藏  举报