OFFSET
1,3
COMMENTS
The definition of the sequence is unclear. - Michael Somos, Mar 07 2012
When n is a square, the trivial solution (x,y) = (1,0) is taken; otherwise we take the least nontrivial solution that satisfies one of the four equations with +1, -1, +4 or -4. - Ray Chandler, Aug 22 2015
REFERENCES
A. Cayley, Report of a committee appointed for the purpose of carrying on the tables connected with the Pellian equation ..., Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 13, pp. 430-443.
C. F. Degen, Canon Pellianus. Hafniae, Copenhagen, 1817.
D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, p. 55.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Ray Chandler, Table of n, a(n) for n = 1..10000
A. Cayley, Report of a committee appointed for the purpose of carrying on the tables connected with the Pellian equation ..., Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 13, pp. 430-443. (Annotated scanned copy)
MATHEMATICA
r[x_, n_] := Reduce[lhs = x^2 - n*y^2; y > 0 && (lhs == -1 || lhs == 1 || lhs == -4 || lhs == 4), y, Integers]; a[n_ /; IntegerQ[Sqrt[n]]] = 1; a[n_] := (x = 1; While[r[x, n] === False, x++]; x); yy[1, _] = 1; yy[x_, n_] := r[x, n][[2]]; A006704 = Table[x = a[n]; Print[{n, x, yy[x, n]}]; x, {n, 1, 55}] (* Jean-François Alcover, Mar 07 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
5 terms corrected by Jean-François Alcover, Mar 09 2012
Corrected a(47)=48 and extended by Ray Chandler, Aug 22 2015
STATUS
approved