题意为给定整数n,求解n以内的整数中有多少个素数(质数)
1.最朴素的做法是:先从2到n-1判断是不是素数,每个数字i都用2到i-1的数字去除,有一个可以整除它,就不是质数。
bool isPrime(int x){
if(x<=1)return false;
for(int i=2;i<x;i++)
if(x%i==0)return false;
return true;}
int countPrimes(int n) {
int k=0;
for(int i=0;i<n;i++)
if(isPrime(i))k++;
return k;
}
这种方法时间复杂度O(n^2), Time Limit Exceeded
2.基于上面的改善:
1. 不用2到i-1的除。从2到sqrt(i)即可,很容易理解吧。
2. 如果2不能整除,也不用尝试4、6、8……所有的偶数了。
3. 那么除了2以外,只要看质数。但是如果3不能整除,好像6、9、12……也肯定不能整除。5同理。
事实上,这是一个定理。将该数N用小于等于sqrt{N}的所有素数去试除,若均无法整除,则N为素数。原理就是,一个合数的因数,肯定能分解成n个质数的乘积。
也就是,只要用质数去整除他就可以了,当然这里包含了2。
2. 如果2不能整除,也不用尝试4、6、8……所有的偶数了。
3. 那么除了2以外,只要看质数。但是如果3不能整除,好像6、9、12……也肯定不能整除。5同理。
事实上,这是一个定理。将该数N用小于等于sqrt{N}的所有素数去试除,若均无法整除,则N为素数。原理就是,一个合数的因数,肯定能分解成n个质数的乘积。
也就是,只要用质数去整除他就可以了,当然这里包含了2。
public class Solution {
public int countPrimes(int n) {
int count = 0;
List<Integer> primes = new ArrayList<Integer>();
for(int i = 2; i < n; i++) {
boolean isPrime = true;
// 将该数N用小于等于sqrt{N}的所有素数去试除,若均无法整除,则N为素数
for(Integer prime : primes) {
if(prime > Math.sqrt(i)) {
break;
}
if(i % prime == 0) {
isPrime = false;
break;
}
}
if(isPrime) {
count++;
primes.add(i);
}
}
return count;
}
}
3.上面的解法还会重复判断。对后面的每个数字,还是要从前面已经得出的所有prime里面,循环再判断一遍。故有埃拉托斯特尼筛法。埃氏筛,是一种公元前250年由古希腊数学家埃拉托斯特尼所提出的一种简单检定素数的算法。
给出要筛数值的范围n,找出\sqrt{n}以内的素数p_{1},p_{2},\dots,p_{k}。先用2去筛,即把2留下,把2的倍数剔除掉;再用下一个素数,也就是3筛,把3留下,把3的倍数剔除掉;接下去用下一个素数5筛,把5留下,把5的倍数剔除掉;不断重复下去......。
详细列出算法如下:
列出2以后的所有序列:
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
标出序列中的第一个素数,也就是2,序列变成:
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
将剩下序列中,划摽2的倍数(用红色标出),序列变成:
本例中,因为25大于2的平方,我们返回第二步:
剩下的序列中第一个素数是3,将主序列中3的倍数划出(绿色),主序列变成:
25仍然大于3的平方,所以我们还要返回第二步:
现在序列中第一个素数是5,同样将序列中5的倍数(蓝色)划出,主序列成了:
因为25等于5的平方,跳出循环.
结论:去掉红色的数字,2到25之间的素数是:2 3 5 7 11 13 17 19 23。
列出2以后的所有序列:
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
标出序列中的第一个素数,也就是2,序列变成:
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
将剩下序列中,划摽2的倍数(用红色标出),序列变成:
- 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
本例中,因为25大于2的平方,我们返回第二步:
剩下的序列中第一个素数是3,将主序列中3的倍数划出(绿色),主序列变成:
- 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
25仍然大于3的平方,所以我们还要返回第二步:
现在序列中第一个素数是5,同样将序列中5的倍数(蓝色)划出,主序列成了:
- 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
因为25等于5的平方,跳出循环.
结论:去掉红色的数字,2到25之间的素数是:2 3 5 7 11 13 17 19 23。
int countPrimes(int n) {
if(n<2)return 0;
int *prime=(int *)malloc(n*sizeof(int));
for(int i=0;i<n;i++)
prime[i]=1;
prime[0]=prime[1]=0;
int k=0;
int j;
for(int i=2;i*i<n;i++){
if(prime[i]!=1)continue;//是素数?
for(j=i*i;j<n;j+=i)//小于n?,小于的话将该素数的平方及其平方后+该素数的定位合数
prime[j]=0;
}
for(j=0;j<n;j++)
if(prime[j]==1)k++;
return k;
}