质数的判定 —— 试除法(O(sqrt(n))
质数在大于1的整数中,如果只包含1和本身这两个约数,就被称为质数,或者叫素数。
有3中形式
- 由于sqrt()计算较慢。(不推荐)。
bool is_prime(int n){
if (n < 2) return false;
for (int i = 2; i <= sqrt(n); ++i)
if (n % i == 0) return false;
return true;
}
- 由于
i*i
有可能会大于n,而当n较大时,可能会爆int,从而成为负值。(不推荐)
bool is_prime(int n){
if (n < 2) return false;
for (int i = 2; i * i <= n; ++i)
if (n % i == 0) return false;
return true;
}
- 推荐 —— 试除法
bool is_prime(int n){
if (n < 2) return false;
for (int i = 2; i <= n / i; ++i)
if (n % i == 0) return false;
return true;
}
分解质因数 —— 试除法(O(sqrt(n))
质因数是指能整除给定正整数的质数。
正整数的因数分解可将正整数表示为一连串的质因子相乘,质因子如重复可以用指数表示。
n中最多只包含一个大于 n \sqrt{n} n的质因子
void divide(int x){
for(int i = 2; i <= n / i; ++i) # i 一定是质数
if (n % i == 0)
{
int s = 0;
while(n % i == 0)
{
n /= i;
s ++ ;
}
printf("%d %d\n", i, s);
}
if (n > 1) printf("%d %d\n", n, 1);
printf("\n");
}
筛质数
要求出从1到n中所有的质数的个数
思路是将从2到n中质数的倍数全部删去。剩下的就是质数。
埃氏筛法
O
(
n
l
o
g
l
o
g
n
)
O(nlog{log{n}})
O(nloglogn)
void get_prime(int x){
for(int i = 2; i <= n; ++i){
if(!st[i]){
prime[idx++] = i;
for(int j = i + i; j <= n; j += i) st[j] = true;
}
}
cout << idx << endl;
}
线性筛法
原理:
n只会被最小质因子筛掉
由于当i模primes[j]等于0时就结束,则primes[j]一定是i的最小质因子
- i % pj == 0
pj 一定是i的最小质因子,pj一定是pj * i 的最小质因子。 - i % pj != 0
pj一定小于i的所有质因子,pj也一定是pj*i的最小质因子。
任何一个合数x,假设pj是x的最小质因子,当i枚举到x/pj的时候就会被筛掉
void get_prime(){
for(int i = 2; i <= n; ++i){
if(!st[i]) prime[cnt++] = i;
for(int j = 0; prime[j] <= n / i; ++j){
st[primes[j] * i] = true;
if(i % primes[j] == 0) break; // 关键
}
}
}