Fast GCD

template<typename T>
T gcd(T a, T b) {
    while((a %= b) && (b %= a));
    return a | b;
}