power.h: extract from euler.cc