验证哥德巴赫猜想

//验证2到一百万（任何一个大于2的偶数都可以表示成两个素数之和） 
#include <iostream>
#include <cstring>
#include <cmath>
using namespace std;

const int N = 1000000;

bool vis[N] = {0};
void init_prim()
{
    int i,j,k;
    int m = (int)(sqrt(N)+0.5); 
    memset(vis,0,sizeof(vis));
    vis[0] = 1;
    vis[1] = 1;
    for(i=2; i<m; i++)
    if(!vis[i])
    {
        for(j=2*i; j<N; j+=i)
            vis[j] = 1;
    }  
}

int main()
{
    int i,j,k;
    init_prim();
    for(i=4; i<N; i+=2)
    {
        for(j=2; j<i; j++)
        {
            int temp = i-j;
            if(!vis[j]&&!vis[temp])//为加快速度，可以把素数单独存在数组里 
            {
                cout<<i<<" = "<<j<<" + "<<temp<<endl;
                break;
            }
        }
        if(i==j)
        {
            cout<<"在"<<i<<"处定理不成立！"<<endl;
            exit(1);//在此异常退出时并未打印----------- 
        }
    } 
    cout<<"--------------"<<endl;
    while(1); 
    return 0;
}