We run the algorithm for c = br, r = m,...,n (suggested by Alan Offer).
Here a, b, m, n are positive integers satisfying a > b > 1, and 1 ≤ m ≤ n.
We employ the criterion Ai,c > c + b√c in the main loop and usually correct partial quotients are returned. However the example (a,b,m,n)=(991,2,146,148) shows that there are exceptions.
If we use the stronger cutoff condition Ai,c > c + b2√c in the main loop, as in http://www.numbertheory.org/php/log3.html, this anomaly disappears. However fewer partial quotients will be returned.
Last modified 6th January 2019
Return to main page