素数打表

#define N 50000  //质数范围   
int prime[1000000];  //prime[0]=2,prime[1]=3,prime[2]=5,……

void init_prime()  
{  
    int i, j;  
    for(i = 2;i <= sqrt(N*1.0); ++i)  
    {  
        if(!prime[i])  
            for(j = i * i; j < N; j += i)  
                prime[j] = 1;  
    }  
    j = 0;  
    for(i = 2;i <= N; ++i)  
        if(!prime[i])   
            prime[j++] = i;  
}