RosettaCodeData/Task/Exponentiation-operator/PARI-GP/exponentiation-operator-1.parigp

ex(a, b)={
  my(c = 1);
  while(b > 1,
    if(b % 2, c *= a);
    a = a^2;
    b >>= 1
  );
  a * c
};