Koike's Trace Formula states that Tr((Uκp)n)=−∑0≤u<√pn (u,p)=1H(u2−4pn)γ(u)κγ(u)2−pn−1, where κ∈Cp, H(D) is the Hurwitz class number of D, and γ(u) is the unique p-adic unit root of the equation x2−ux+pn=0. Is there an implementation of this formula in sage ?