一个数（3000内）等于两个素数之和

 

2. #include "stdafx.h"
3. #include "iostream"
5. using namespace std;
7. unsigned int a[500]; // 用于存放前500个素数
8. unsigned int f(unsigned int x);  
9. bool is(unsigned int x);
11. void init() //  初始化数组全零
12. {
13.     for(int i=0; i<500; i++)
14.     {
15.         a[i]=0;
16.     }
17. }
19. void  init2()  //  用素数序列来初始化数组
20. {
21.     int j=0;
22.     for(unsigned int i; i<=4294967294; i++)
23.     {
24.         if(is(i))
25.         {
26.             a[j]=i;
27.             j++;
28.         }
30.         if(j>500)
31.         {
32.             break;
33.         }
34.     }
36. }
38. bool is(unsigned int x) //  判断一个数x是不是素数，是返回true
39. {
40.     unsigned int s;
42.     if(x<2) return false;
43.     if(x == 2 ) return true;
45.     s = f(x);
47.     for(unsigned int i=2 ; i<=s ; i++)
48.     {
49.         if(x%i==0) return false;
50.     }
52.     return true;
53. }
55. unsigned int f(unsigned int x)  //返回一个数，这个数大于他的平方根
56. {
57.     unsigned int i=0;
58.     unsigned int s=1;
60.     while(x>=s)
61.     {
62.         x-=s;
63.         s+=2;
64.         i++;
65.     }
67.     i++;
68.     return  i;
69. }
72.     
74. int main()
75. {
76.     init();
77.     init2();
78.     int i;
79.     int x;
80.     int y;
82.     cout << a[499] << endl;  //打印数组中最大的素数
83. loop:
84.     cout << "input a number : " ;
85.     cin >> i ;
87.     if(i<=2) 
88.     {
89.         cout << "The number must bigger than 3";
90.         goto loop;
91.     }
93.     for(x=0; x<=499; x++)
94.     {
96.         for(y=x+1; y<=499; y++)
97.         {
98.             if(a[x]+a[y]==i) goto end;
99.         }
100.     }
102.     cout << "error " <<endl;
104. end:
105.     cout << i << "=" << a[x] << " + " <<  a[y] << endl;
106.   
107.     cout << endl;
109. }

一个数（3000内）等于两个素数之和