Hi, David. This is roughly Robin's comment, just with longer discussions, and, I am afraid, links to my own answers. That's life. I think you would probably enjoy J.H. Conway, "The Sensual Quadratic Form," especially PSL_2(Z) on pages 27-33. I always like "The Markoff and Lagrange Spectra" by Thomas W. Cusick and Mary E. Flahive and "Binary Quadratic Forms" by Duncan A. Buell. On particular issues I think you are raising, answers by me and by Dror in: Reasons for switching from simple continued fractions for $\sqrt D$ to reduced forms: https://mathoverflow.net/questions/22811/ Numbers (here primes) occurring as "diagonal" coefficients: https://mathoverflow.net/questions/29850/ Other: https://mathoverflow.net/questions/22552/ In my language, I think you are asking mostly about small numbers occurring as diagonal coefficients, just not necessarily prime. For me, such numbers need also to be *primitively* represented by the forms in the cycle (for you I guess it is by $ x^2 - D y^2 ,$ using $a_0 = \lfloor \sqrt D \rfloor$ your quadratic form is equivalent to the reduced form $\langle 1,\; 2 a_0, \; a_0^2 - D \rangle.$ ) phoebus:~/Cplusplus> ./Pell Input n for Pell 67 0 form 1 16 -3 delta -5 1 form -3 14 6 delta 2 2 form 6 10 -7 delta -1 3 form -7 4 9 delta 1 4 form 9 14 -2 delta -7 5 form -2 14 9 delta 1 6 form 9 4 -7 delta -1 7 form -7 10 6 delta 2 8 form 6 14 -3 delta -5 9 form -3 16 1 delta 16 10 form 1 16 -3 disc 268 Automorph, written on right of Gram matrix: -1106 -17901 -5967 -96578 Pell automorph -48842 -399789 -5967 -48842 Pell unit -48842^2 - 67 * -5967^2 = 1 ========================================= phoebus:~/Cplusplus> Anyway, let me know if you want to see any more worked examples. Or the computer program. It is C++ so the numbers are bounded. Easy enough in other languages, of course.