题目描述:
统计所有小于非负整数 n
的质数的数量。
暴击法(会超时,不能过):
class Solution {
public:
bool isprime(int x)
{
for (int i = 2; i * i <= x; i++)
{
if (x % i == 0) //说明不是质数
{
return false;
}
}
return true;
}
int countPrimes(int n) {
int num = 0;
for (int i = 2; i < n; i++)
{
num += isprime(i);
}
return num;
}
};
关于暴力法,虽然不能过,但自己在写暴力法时,寻找因子的区间是在[0, x/2]之间,看完题解之后才发现想错了:考虑到如果 y 是 x 的因数,那么 x/y 也必然是 x 的因数,因此我们只要校验 y 或者 x/y 即可。而如果我们每次选择校验两者中的较小数,则不难发现较小数一定落在[2, √x] 的区间中。
时间复杂度:O(n√n)。单个数检查的时间复杂度为O(√n) ,一共要检查 O(n) 个数,因此总时间复杂度为 O(n√n)。空间复杂度:O(1)。
方法一:
关于埃氏筛法:
我们设 isPrime[i] 表示数 i 是不是质数,如果是质数则为 1,否则为 0。从小到大遍历每个数,如果这个数为质数,则将其所有的倍数都标记为合数(除了该质数本身),即 0,这样在运行结束的时候我们即能知道质数的个数。
常规的埃氏筛:
class Solution {
public:
int countPrimes(int n) {
int num = 0;
vector<int> isprime(n, 1);
for (int i = 2; i < n; i++)
{
if (isprime[i]) //如果是质数
{
num++;
}
for (int j = 2 * i; j < n; j+=i) //将其倍数标记为合数
{
isprime[j] = 0;
}
}
return num;
}
};
注意到,对于一个质数 x,如果按上文说的我们从 2x 开始标记其实是冗余的,应该直接从 x·x 开始标记,因为 2x,3x,… 这些数一定在 x·x 之前就被其他数的倍数标记过了,例如 2 的所有倍数,3 的所有倍数等
于是改进的埃氏筛:
class Solution {
public:
int countPrimes(int n) {
int num = 0;
vector<int> isprime(n, 1);
for (int i = 2; i < n; i++)
{
if (isprime[i]) //如果是质数
{
num++;
}
if ((long long)i * i < n) //先转化成longlong类型判断,否则i*i可能越界超过int的最大范围
{
for (int j = i * i; j < n; j+=i) //将其倍数标记为合数
{
isprime[j] = 0;
}
}
}
return num;
}
};
时间复杂度:O(nloglogn)。空间复杂度:O(n)。