快速幂运算

quick_pow.cpp

#include <type_traits>

template <typename T, typename P>
T qpow(
  T a, 
  typename std::enable_if<std::is_integral<P>::value, P>::type n,
  T init) {
    T ans = init;
    while (n) {
        if (n & 1) {
            ans = ans * a;
        }
        n >>= 1;
        a = a * a;
    }
    return ans;
}

// usage
#include <iostream>

int main() {
    int res = qpow<int, int>(2, 11, 1);
    std::cout << res << std::endl;
    return 0;
}