如果一个数只有1和他本身两个因数,那这个数就是质数。
例:7 = 1 x 7
,5 = 1 x 5
。
如果一个数除了1和他本身,还有其他因数,那这个数就是合数。
例:8 = 1 x 8 = 2 x 4
,12 = 1 x 12 = 2 x 6 = 3 x 4
。
1既不是质数,也不是合数
就是从2开始,一直到根号 n \sqrt{n} n ,判断每一个数是不是n的因数。如果是的话n就不是质数了。
bool prime(int n){ for(int i = 2; i * i <= n; i++){//这样能避免精度误差 if(n % i == 0){ return false;//找到了因数,就直接返回合数。 } } return true;//没有找到其他因数,返回质数。 }
假如一个数是质数,那他的倍数(除了他本身) 都是 合数。
那我们可以用一个数组
p
r
i
m
e
prime
prime,用来记录一个数是不是质数。
int prime[MAXN]; int get_prime(int n){//返回值是质数的个数 memset(n, true, sizeof n); for(int i = 2; i <= n; i++){ if(prime[i]){//如果i是质数 cnt++; for(int j = 2; j * i <= n; j++){//枚举i的倍数 prime[i * j] = false; } } } return cnt; }
它的复杂度是 O ( n l o g n ) O(n log n) O(nlogn),详情见文章:Eratosthenes筛法(埃式筛法)时间复杂度分析。
欧拉欧拉欧拉! 欢迎欧拉再次登上c++的舞台!
回到正题。
看刚才的算法,我们发现重复的地方就在于一个合数可能被计算多次。
例如:24 = 1 × 24 = 2 × 12 = 3 × 8 = 4 × 6
。被计算了8次!这明显浪费时间了。
所以,我们统一把合数写成如下形式:合数 = 质数 × i
,也就是说把i的枚举倍数设定为质数,这样每个数都会只被判定一次了。
int prime[MAXN]; int get_prime(int n){//返回值是质数的个数 memset(n, true, sizeof n); for(int i = 2; i <= n; i++){ if(prime[i]){//如果i是质数 prime[++cnt] = i;//记录下质数 for(int j = 1; j <= cnt; j++){//枚举i的倍数 prime[i * prime[j]] = false; } } } return cnt; }
值得一提的是,prime[++cnt] = i;
这个语句必须在 j循环 的前面,不然所有质数的平方都逃走了。
时间复杂度:
O
(
n
)
O(n)
O(n)