埃拉托斯特尼筛法:
这种方法的思想是设置一个标志数组isPrimes[n],标志数组的每一位标示相应的数是不是素数,初始化为全true。
算法从i=2开始,依次将质数的倍数标记为非素数,即将标记数组的相应位改为false。
标记质数的倍数的时候从i*i也就是i的平方开始标记,不需要从i*j(且j<i
)开始,因为i*j,在遇到j时已经被标记了,因为j比i小,所以遇到j比遇到i要早。
因为从i*i开始标记,所以i终止条件也为sqrt(n),否则i*i将大于n。
此方法时间复杂度为O(n*lglgn),空间复杂度为O(n)。
void isprime(int n)
{
if (n < 2)
return ;
vector<int> vec(n+1,1);
int sq = pow(n+1, 0.5);
int j;
for (auto i = 2; i < sq; ++i)
{
if (vec[i] == 0)
continue;
for (j = i*i; j <= n; j += i)
{
vec[j] = 0;
}
}
stringstream str_stream;
for (auto i = 2; i < vec.size(); ++i)
{
if (vec[i] == 1)
str_stream <<i<<" ";
}
str_stream << endl;
cout << str_stream.str();
}
int main()
{
clock_t clockBegin, clockEnd;
clockBegin = clock();
isprime(50000);
clockEnd = clock();
printf("%d\n", clockEnd - clockBegin);
return 0;
}