如果要写一个素数算法,我想大家第一反应会是如下的代码:
bool IsPrime(unsigned uValue){
for(unsigned int i= 2; i< uValue; ++i){
if(uValue%i== 0){
return false;
}
}
return true;
}
这个代码的正确性肯定是不需要否定的,但是,我们根据数学中的唯一分解定理可以推出,如果一个数N不是素数,那么一定存在x使得x|N(注意'|'读作整除),那么知道x和N/x都是N的约数,自然的有min(x, N/x)<= 根号(N)。所以这个代码可以往下优化如下的代码:
bool IsPrime(unsigned uValue){
unsigned int uMax= sqrt(uValue); //在math.h里面声明了
for(unsigned int i= 2; i<= uMax; ++i){
if(uValue%i== 0){
return false;
}
}
return true;
}
这个代码基本上是大家都写的代码,其正确性也是勿用质疑的。然后,这个程序仍然不是最简的,事实上