[math-fun] How can I accurately calculate (p/q)^n mod 1, for large n