for (i=2;i<=n;i++)
    if (!p[i])
        for (j=i*i;j<=n;j+=i)
            p[j] = 1;