This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A006704 M0119 #31 Jul 08 2025 16:53:27
%S A006704 1,1,2,1,1,5,8,2,1,3,10,4,3,15,4,1,4,17,170,4,5,197,24,5,1,5,26,16,5,
%T A006704 11,1520,6,23,35,6,1,6,37,25,6,32,13,3482,20,7,24335,48,7,1,7,50,36,7,
%U A006704 485,89,15,151,99,530,8,39,63,8,1,8,65,48842,8,25,251,3480,17,1068,43
%N A006704 Solution to Pellian: x such that x^2 - n y^2 = +- 1, +- 4.
%C A006704 The definition of the sequence is unclear. - _Michael Somos_, Mar 07 2012
%C A006704 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
%D A006704 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.
%D A006704 C. F. Degen, Canon Pellianus. Hafniae, Copenhagen, 1817.
%D A006704 D. H. Lehmer, Guide to Tables in the Theory of Numbers. Bulletin No. 105, National Research Council, Washington, DC, 1941, p. 55.
%D A006704 N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H A006704 Ray Chandler, <a href="/A006704/b006704.txt">Table of n, a(n) for n = 1..10000</a>
%H A006704 A. Cayley, <a href="/A002349/a002349.pdf">Report of a committee appointed for the purpose of carrying on the tables connected with the Pellian equation ...</a>, Collected Mathematical Papers. Vols. 1-13, Cambridge Univ. Press, London, 1889-1897, Vol. 13, pp. 430-443. (Annotated scanned copy)
%t A006704 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 *)
%Y A006704 Cf. A006705.
%K A006704 nonn
%O A006704 1,3
%A A006704 _N. J. A. Sloane_
%E A006704 5 terms corrected by _Jean-François Alcover_, Mar 09 2012
%E A006704 Corrected a(47)=48 and extended by _Ray Chandler_, Aug 22 2015