简单遍历:
bool isprime(int n) { int s = sqrt( double(n) )+1; //对n开根号 for(int i=2;i<=s;i++) //n除以每个比n开根号小比1大的自然数 if(n%i==0) //如果有能被整除的,则不是质数 return 0; return 1; }
筛法求素数:
筛法求质数,效率最高,但会比较浪费内存
首先建立一个boolean类型的数组,用来存储你要判断某个范围内自然数中的质数,例如,你要输出小于200的质数,你需要建立一个大小为201(建立201个存储位置是为了让数组位置与其大小相同)的boolean数组,初始化为true。
其次用第二种方法求的第一个质数(在此是2),然后将是2的倍数的数全置为false(2除外),即2、4、6、8……位置上置为false。然后是3的倍数的全置为false(3除外),一直到14(14是200的开平方),这样的话把不是质数的位置上置为false了,剩下的全是质数了,挑着是true的打印出来就行了。
#define N 1000 bool isprm[N]; void isprime() { int i,j,k=0; int s,e=sqrt( double(N) )+1; //sqrt是对于double数开平方 memset(isprm,1,sizeof(isprm)); //prm[k++]=2; isprm[0] = isprm[1] = 0; for(i=4 ;i < N; i=2+i) isprm[i]=0; for(i=3;i<e;i=2+i) if(isprm[i]) for(s=i*2,j=i*i;j<N;j=j+s) isprm[j]=0; //因为j是奇数,所以+奇数后是偶数,不必处理 }
今天再次看筛法求素数的竟然有点不懂了,首先不知道为什么
for(i=3;i<e;i=2+i) if(isprm[i]) for(s=i*2,j=i*i;j<N;j=j+s) isprm[j]=0;这里的i<e,其实这和下面的一个不懂得可以一并解决(也就是j的初始值为什么是i*i),可以举例来说明比如35就可以用5*5+2*5来看,这个5*5+2*5可以看做是在遍历以5为因数的时候做出的,所以说在遍历一个数的时候i的倍数只要时候小于i*i的时候就可以不管了,因为在以前遍历小因子的时候肯定是做过了,只要从i*i的开始就好了(因为i*i可能通过其他的乘积的不到,比如25),先这么理解着吧。。。