1 //验证2到一百万（任何一个大于2的偶数都可以表示成两个素数之和） 2 #include <iostream> 3 #include <cstring> 4 #include <cmath> 5 using namespace std; 6 7 const int N = 1000000; 8 9 bool vis[N] = {0}; 10 void init_prim() 11 { 12 int i,j,k; 13 int m = (int)(sqrt(N)+0.5); 14 memset(vis,0,sizeof(vis)); 15 vis[0] = 1; 16 vis[1] = 1; 17 for(i=2; i<m; i++) 18 if(!vis[i]) 19 { 20 for(j=2*i; j<N; j+=i) 21 vis[j] = 1; 22 } 23 } 24 25 int main() 26 { 27 int i,j,k; 28 init_prim(); 29 for(i=4; i<N; i+=2) 30 { 31 for(j=2; j<i; j++) 32 { 33 int temp = i-j; 34 if(!vis[j]&&!vis[temp])//为加快速度，可以把素数单独存在数组里 35 { 36 cout<<i<<" = "<<j<<" + "<<temp<<endl; 37 break; 38 } 39 } 40 if(i==j) 41 { 42 cout<<"在"<<i<<"处定理不成立！"<<endl; 43 exit(1);//在此异常退出时并未打印----------- 44 } 45 } 46 cout<<"--------------"<<endl; 47 while(1); 48 return 0; 49 }