Eratosthenes筛法
对于不超过n的每个非负整数p,删除2p, 3p,4p,...,当处理完所有数之后,还没有被删除的就是素数。用shu[i]表示i已经被删除。
vector<int> shu(n+1,1);
for (int i=2;i<=n;i++ )
{
if(shu[i])//i为素数
{
for(int j=i*i;j<=n;j+=i)
{
shu[j]=0;
}
}
}
给定外层循环变量i,内层循环的次数n/i。欧拉在1734年得到的结果:1+1/2 ...+1/= ln(n+1)+y,其中欧拉常数v=0.577218。可以在很短的时间内得到10^6以内的所有素数。