打素数表

常用方法

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