判断一个数n是否是素数,只需要判断它是否能被2到n之间的数整除就行了,若不能被整除,则说明是素数。考虑到某数大于n/2时,n不可能被该数整除,故只需遍历2~n/2即可。更进一步,由合数定理可知,若一个数是合数,则它的最小质因数必小于等于该数的平方根,由此可得更为高效的代码。
/*代码如下*/
#include<iostream>
#include<cmath>
using namespace std;
int main()
{
int n;
cin>>n;
int i;
for(i = 2;i <= sqrt(n); i++){
if(n % i == 0) //如果可以整除,则说明不是素数,立即结束循环
break;
}
if(i <= sqrt(n))
cout<<"n is not a prime number!"<<endl;
else
cout<<"n is a prime number!"<<endl;
}
当所求的数值较大时,上述代码依然会更为费时,而求素数最高效的算法应为筛选法(埃拉托色尼筛选法):
(1)先把1删除(现今数学界1既不是质数也不是合数)
(2)读取队列中当前最小的数2,然后把2的倍数删去
(3)读取队列中当前最小的数3,然后把3的倍数删去
(4)读取队列中当前最小的数5,然后把5的倍数删去
(5)如上所述直到需求的范围内所有的数均删除或读取
/*筛选法求素数的具体代码,此处将不符合要求的数置零*/
#include<iostream>
#include<cmath>
using namespace std;
int main()
{
int i,n = 1000;
int a[1001];
for(i = 0;i <= 1001; i++)
a[i] = i;
a[1] = 0;
for(i = 2;i < sqrt(n); i++){
for(int j = i + 1;j <= n ; j++){
if(a[j] != 0 && a[j] % i == 0)
a[j] = 0;
}
}
for(i = 1,n = 0; i<= 1000; i++){
if(a[i] != 0){
cout<<a[i]<<"\t";
if(++n % 10 == 0)
cout<<endl;
}
}
}