PHP求最大公约数 && 最小公倍数 && 质数

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