常见算法题（1）求最大公因数

最大公因（约）数，是指两个或多个整数共有因数中最大的一个。两个数a，b的最大公因数记作(a,b)。

求最大公因数的算法有：

(1)辗转相除法
有两整数a和b：
    ① a%b得余数c
    ② 若c=0，则b即为两数的最大公约数
    ③ 若c≠0，则a=b，b=c，再回去执行①
例如求27和15的最大公约数过程为：
    27÷15 余12 15÷12余3 12÷3余0。因此，3即为最大公约数

⑵ 更相减损法
有两整数a和b：
    ① 若a>b，则a=a-b
    ② 若a<b，则b=b-a
    ③ 若a=b，则a（或b）即为两数的最大公约数
    ④ 若a≠b，则再回去执行①
例如求27和15的最大公约数过程为：
    27－15＝12( 15>12 ) 15－12＝3( 12>3 )
    12－3＝9( 9>3 ) 9－3＝6( 6>3 )
    6－3＝3( 3==3 )
    因此，3即为最大公约数
    
⑶穷举法
有两整数a和b：
    ① i=1
    ② 若a，b能同时被i整除，则t＝i
    ③ i++
    ④ 若 i <= a(或b)，则再回去执行②
    ⑤ 若 i > a(或b)，则t即为最大公约数，结束
改进：
    ① i= a(或b)
    ② 若a，b能同时被i整除，则i即为最大公约数，
    结束
    ③ i--，再回去执行②

下面以一个实际例子来解释第一种算法。

【例1】实际问题，求a和b^b次方的最大公因数。

描述

    As we know, gcd(A, B) means the Greatest Common Divisor (GCD) of A and B.

    But zy thinks it is so easy to just let you calculate it.
    So now your task is to calculate gcd(A, B^B),it is easy right?

    Notice: B ^ B means B multiply by himself B times.

输入

    Multiply test case, each case per line.
    Each line 2 integers A, B, 1 <= A, B <= 1000

输出

    For each test case output one line( one number ), indicating the value of gcd (A, B^B)

样例输入

100 1
100 2

样例输出1

1
4

解决方案分析：

    ！！！注意！！！
    这里有个大坑！！！
    这里是求a与b^b次方的最大公因数，想当然的肯定是先求出b^b的值，再求a与b^b的最大公因数。
    ！！！但是！！！
    当b取值很大时，b^b的值会非常之大！！！而且pow函数的返回的可表示最大数的数值类型为long double（12字节），当求幂的数过大时，pow无法给出正确计算结果。

    因此，这里采取一个比较巧妙的方法：

    记a与b的最大公因数为(a,b)，可知a与b^b的最大公因数(a,b^b)。

算法如下：

1. 求出(a,b)，并使用result = result * (a,b)保存连乘的结果。转到2。
2. 如果(a,b) == 1,结束。(a,b^b) = result。
3. 如果(a,b) != 1,令a = a/(a,b)，转到1。

具体代码如下：

#include<iostream>
using namespace std;

int gcd(int a, int b){
    int r;
    while(b!=0){
        r = a % b;
        a = b;
        b = r;
    }
    return a;
}

int main(){
    int a,b,result,flag;
    while(cin >> a >> b){
        result = 1;
        if(a<1 || b>1000)
            exit(0);
        for(int i=0; i<b; i++){
            flag = gcd(a,b);
            result = result*flag;
            if(flag == 1)
                break;
            else
                a = a/flag;
        }    
        cout << result << endl;
    }
    return 0;
}

也可以写成如下的简洁形式：

#include<iostream>
using namespace std;

int gcd(int a, int b){
    int r;
    while(b!=0){
        r = a % b;
        a = b;
        b = r;
    }
    return a;
}

int main(){
    int a,b,result,flag;
    while(cin >> a >> b){
        result = 1;
        if(a<1 || b>1000)
            exit(0);
         for(; gcd(a,b)!=1; a=a/gcd(a,b)){
                result = result*gcd(a,b);
        }
        cout << result << endl;
    }
    return 0;
}