We limit input to m'n' ≤ 2000.
The paper shows that there is a quasi-period consisting of d=mn linear expressions of the form pt22qx+qt, t=1,...,d. A criteria is given for all the pt to be equal to 1. In the case of (m/1)e, the primes m = 7, 11, 13, 31, 41, ... give at least one progression with pt > 1.
In the first version of this program, a zero partial quotient was present, being located at the end of the factorizations of the first or the start of the factorization of the second matrix in Theorem 1 (b) of the paper. This is now removed, thanks to an elegant proof of Alan Offer on 30th June 2010. We also have improved the code for getting the "initial" matrix, by using a suggestion of Alan that we use the recurrence relations of his proof.
Last modified 9th January 2013
Return to main page