最大公约数

编程之美2.7节。解法三即stein算法。下面是其迭代式实现。

template<class T>
T gcdStein(T a, T b) {
    for(T r = 1;;) {
        if (!a) return r *= b;
        if (!b) return r *= a;

        T o1 = 0;
        while (!(a&1)) a>>=1, ++o1;
        T o2 = 0;
        while (!(b&1)) b>>=1, ++o2;
        r <<= min(o1, o2);

        if (a == 1 || b == 1) return r;
        if (a > b) a -= b;
        else b -= a;
    }
}

int main() {
    int a, b;
    while (cin>>a>>b)
        cout<<gcdStein<int>(a, b)<<'\n';
    return 0;
}