poj 2262（筛法求素数）

题意：就是一百万内的哥德巴赫猜想，每个整数都能拆成两个素数之和，输入一个数，把这个等式打印出来，输出整数=最小+最大就可以了

这里必须要用到筛法：

用筛法求素数的基本思想是：把从1开始的、某一范围内的正整数从小到大顺序排列， 1不是素数，首先把它筛掉。剩下的数中选择最小的数是素数，然后去掉它的倍数。依次类推，直到筛子为空时结束。如有：

1 2 3 4 5 6 7 8 9 10

11 12 13 14 15 16 17 18 19 20

21 22 23 24 25 26 27 28 29 30

1不是素数，去掉。剩下的数中2最小，是素数，去掉2的倍数，余下的数是：

3 5 7 9 11 13 15 17 19 21 23 25 27 29

剩下的数中3最小，是素数，去掉3的倍数，如此下去直到所有的数都被筛完，求出的素数为：

2 3 5 7 11 13 17 19 23 29

 

#include<iostream>
#include<string>
#define M 1000002           //最大一百万，多个2保险
using namespace std;

int prime[M];

void GetPrime()
{
	int i,j;
	memset(prime,1,sizeof(prime));  //将每个数字都先放进素数表中
		
	prime[0]=prime[1]=0;             //1和0不是素数
	for( i=2;i<=M;i++)               
		if(prime[i])                   //如果是素数
		{                           
			for( j=i*2;j<=M;j+=i)     //从i*2开始将i的倍数都筛除数组  
				prime[j]=0;
		}
}

int main()
{
	int n,i;
	GetPrime(); //构建素数表过程
	while(scanf("%d",&n)!=EOF)
	{
		for(i=2;i<=n/2;i++)
			if(prime[i]&&prime[n-i])           //如果i是素数并且n-i也是素数，则满足定理
			{
				printf("%d = %d + %d\n",n,i,n-i);
				break;
			}
	}
}