用 Binary powering

int powmod(int a, int b) {
  int bin=1, i, res=1, tmp=a%10007;

  for(i=0; i<32; i++) {
    if(b & bin)
      res = res * tmp % 10007;
    bin <<= 1;
    tmp = tmp * tmp % 10007;
  }
  return res;
}