寻找素数对
题目: http://acm.hdu.edu.cn/showproblem.php?pid=1262
Problem Description
哥德巴赫猜想大家都知道一点吧.我们现在不是想证明这个结论,而是想在程序语言内部能够表示的数集中,任意取出一个偶数,来寻找两个素数,使得其和等于该偶数.
做好了这件实事,就能说明这个猜想是成立的.
由于可以有不同的素数对来表示同一个偶数,所以专门要求所寻找的素数对是两个值最相近的.
Input
输入中是一些偶整数M(5<M<=10000).
Output
对于每个偶数,输出两个彼此最接近的素数,其和等于该偶数.
Sample Input
20 30 40
Sample Output
7 13
13 17
题目: http://acm.hdu.edu.cn/showproblem.php?pid=1262
Problem Description
哥德巴赫猜想大家都知道一点吧.我们现在不是想证明这个结论,而是想在程序语言内部能够表示的数集中,任意取出一个偶数,来寻找两个素数,使得其和等于该偶数.
做好了这件实事,就能说明这个猜想是成立的.
由于可以有不同的素数对来表示同一个偶数,所以专门要求所寻找的素数对是两个值最相近的.
Input
输入中是一些偶整数M(5<M<=10000).
Output
对于每个偶数,输出两个彼此最接近的素数,其和等于该偶数.
Sample Input
20 30 40
Sample Output
7 13
13 17
17 23
做了一通DFS、BFS休息会,刚好看到书上有关于素数的,做一下。
中文题目,不解释。
搜索从m/2开始搜索,小的放前面。
#include <iostream>
#include <string.h>
using namespace std;
bool prim[100001];
/*
// 素性测试O(n^0.5)
bool is_prim(int n)
{
int i;
for(i=2; i*i<=n; ++i)
if(n%i==0) return false;
return n!=1;
}
*/
// 埃氏筛法
void make_prim(int n)
{
memset(prim,0,sizeof(prim));
int i,j;
prim[0]=prim[1]=1;
for(i=2;i<=n;++i)
{
if(prim[i]) continue;
for(j=i+i;j<=n;j+=i)
prim[j]=1;
}
}
int main()
{
int m,a,b;
make_prim(10000);
while(cin>>m)
{
for(a=m/2;a<m;--a)
{
b=m-a;
if(!prim[a] && !prim[b]) break;
}
cout<<a<<" "<<b<<endl;
}
return 0;
}