mpir-2.6.0/mpz/n_pow_ui.cの153行目の関数
void mpz_n_pow_ui (mpz_ptr r, mp_srcptr bp, mp_size_t bsize, mpir_ui e)
がべき乗計算のメインだろ。
ただし単独では動かないだろうから、他ソースも参照の必要あり。


MPIR Multiple Precision Integers and Rationals
http://www.mpir.org/