In this article we present the post-Newtonian (pN) coefficients of the energy
flux (and angular momentum flux) at infinity and event horizon for a particle
in circular, equatorial orbits about a Kerr black hole (of mass M and
spin-parameter a) up to 20-pN order. When a pN term is not a polynomial in
a/M and includes irrational functions (like polygamma functions), it is
written as a power series of a/M. This is achieved by calculating the fluxes
numerically with an accuracy greater than 1 part in 10600. Such high
accuracy allows us to extract analytical values of pN coefficients that are
linear combinations of transcendentals like the Euler constant, logarithms of
prime numbers and powers of π. We also present the 22-pN expansion
(spin-independent pN expansion) of the ingoing energy flux at the event horizon
for a particle in circular orbit about a Schwarzschild black hole