1.最大公约数
greatest common divisor
(1)更相减损法
大减小,一直减,减到出现相同的数,即为最小公倍数
BUG:不能求负数的最大公约数
function gcd($a,$b){
while ($a != $b) {
if($a > $b) $a -= $b;
else $b -= $a;
}
return $a;
}
echo gcd(10,15);
(2)辗转相除法
用b除a,得a=bq…r1(0≤r)。若r1=0,则(a,b)=b;若r1≠0,则再用r1除b,得b=r1q…r2(0≤r2).若r2=0,则(a,b)=r1,若r2≠0,则继续用r2除r1,……如此下去,直到能整除为止。其最后一个非零余数即为(a,b)
该算法可以求所有,包括负数
//递归解法
function gcd($a, $b){
if($b==0)
return $a;
else
return gcd($b, $a%$b);
}
//非递归解法
function gcd($a, $b){
while ($b != 0) {
$tmp = $a % $b;
$a = $b;
$b = $tmp;
}
return $a;
}
echo gcd1(-10,-15);
2.最小公倍数
分解因数法未解
(1)穷举
function gbs($a, $b ){
if(!is_int($a) || !is_int($b)){
return false; //不是整数,返回false
}
if($a <= 0 || $b <= 0){
return false; //不是正整数,返回false
}
$max = $a >= $b ? $a : $b;
$min = $a < $b ? $a : $b;
//比较大的数,一直翻倍,直到能整除小的数
for($i = $max; $i <= $max * $min; $i+= $max){
if( $i % $min == 0 ){
return $i; //这就是要求的最小公倍数
}
}
}
echo gbs(10,15); //30
(2)公式法
两个数的乘积等于这两个数的最大公约数与最小公倍数的乘积
function GongBeiShu($a,$b){
$gys = gcd($a, $b);
return $a*$b/$gys;
}
//最大公约数-递归解法
function gcd($a, $b){
if($b==0)
return $a;
else
return gcd($b, $a%$b);
}
echo GongBeiShu(10,15); //30
3.质数
除了1和其本身外,不被其他数整除的数
(1)是否是质数
求一个数是否是质数,只需遍历开根后的范围即可,不必从头到尾遍历,因为一个数由两个数相乘,到某一个因数的平方后,两个乘数就会翻转。
function isZhiShu($num){
$len = floor(sqrt($num));
for($i = 2;$i<=$len;++$i){
if($num%$i == 0)
return '否';
}
return '是';
}
echo isZhiShu(87);//否
(2)求范围内的质数(如2-200)
利用上述证明是否是质数的函数去求
PS:2是最小的质数,也是唯一一个偶质数
function isZhiShu($num){
$len = floor(sqrt($num));
for($i = 2;$i<=$len;++$i){
if($num%$i == 0)
return false;
}
return true;
}
function zhiShu($low, $high){
for($i = $low; $i<=$high;++$i){
if(isZhiShu($i)){
echo $i,' ';
}
}
}
echo zhiShu(2,200);
//2 3 5 7 11 13 17 19 23 29 31 37 41 43 47
//53 59 61 67 71 73 79 83 89 97 101 103 107
//109 113 127 131 137 139 149 151 157 163
//167 173 179 181 191 193 197 199