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 A228405 #30 Oct 29 2021 23:16:43 %S A228405 0,1,1,0,1,2,2,2,3,5,0,2,4,7,12,4,4,6,10,17,29,0,4,8,14,24,41,70,8,8, %T A228405 12,20,34,58,99,169,0,8,16,28,48,82,140,239,408,16,16,24,40,68,116, %U A228405 198,338,577,985,0,16,32,56,96,164,280,478,816,1393,2378 %N A228405 Pellian Array, A(n, k) with numbers m such that 2*m^2 +- 2^k is a square, and their corresponding square roots, read downward by diagonals. %C A228405 The left column, A(n,0), is A000129(n), Pell Numbers. %C A228405 The top row, A(0,k), is A077957(k) plus an initial 0, which is the inverse binomial transform of A000129. %C A228405 These may be considered initializing values, or results, depending the perspective taken, since there are several ways to generate the array. See Formula section for details. %C A228405 The columns of the array hold all values, in sequential order, of numbers m such that 2m^2 + 2^k or 2m^2 - 2^k are squares, and their corresponding square roots in the next column, which then form the "next round" of m values for k+1. %C A228405 For example A(n,0) are numbers such that 2m^2 +- 1 are squares, the integer square roots of each are in A(n,1), which are then numbers m such that 2m^2 +- 2 are squares, with those square roots in A(n,2), etc. %C A228405 A(n, k)/A(n,k-2) = 2; A(n,k)/A(n,k-1) converges to sqrt(2) for large n. %C A228405 A(n,k)/A(n-1,k) converges to 1 + sqrt(2) for large n. %C A228405 The other columns of this array hold current OEIS sequences as follows: %C A228405 A(n,1) = A001333(n); A(n,2) = A163271(n); A(n,3) = A002203(n); %C A228405 Bisections of these column-oriented sequences also appear in the OEIS, corresponding to the even and odd rows of the array, which in turn correspond to the two different recursive square root equations in the formula section below. %C A228405 Farey fraction denominators interleave columns 0 and 1, and the corresponding numerators interleave columns 1 and 2, for approximating sqrt(2). See A002965 and A119016, respectively. %C A228405 The other rows of this array hold current OEIS sequences as follows: %C A228405 A(1,k) = A016116(k); A(2,k) = A029744(k) less the initial 1; %C A228405 A(3,k) = A070875(k); A(4,k) = A091523(k) less the initial 8. %C A228405 The Pell Numbers (A000219) are the only initializing set of numbers where the two recursive square root equations (see below) produce exclusively integer values, for all iterations of k. For any other initial values only even iterations (at k = 2, 4, ...) produce integers. %C A228405 The numbers in this array, especially the first three columns, are also integer square roots of these expressions: floor(m^2/2), floor(m^2/2 + 1), floor (m^2/2 - 1). See A227972 for specific arrangements and relationships. Also: ceiling(m^2/2), ceiling(m^2/2 + 1), ceiling (m^2/2 -1), m^2+1, m^2-1, m^2*(m^2-1)/2, m^2*(m^2-1)/2, in various different arrangements. Many of these involve: A000129(2n)/2 = A001109(n). %C A228405 A001109 also appears when multiplying adjacent columns: A(n,k) * A(n,k+1) = (k+1) * A001109(n), for all k. %H A228405 Seiichi Manyama, <a href="/A228405/b228405.txt">Antidiagonals n = 0..139, flattened</a> %H A228405 MacTutor, <a href="https://mathshistory.st-andrews.ac.uk/Extras/Thompson_irrationals/">D'Arcy Thompson on Greek irrationals</a> %H A228405 D'Arcy Thompson, <a href="https://www.jstor.org/stable/2249223">Excess and Defect: Or the Little More and the Little Less</a>, Mind, New Series, Vol. 38, No. 149 (Jan., 1929), pp. 43-55 (13 pages). See page 50. %F A228405 If using the left column and top row to initialize: A(n,k) = A(n,k-1) + A(n-1,k-1). %F A228405 If using only the top row to initialize, then each column for k = i is the binomial transform of A(0,k) restricted to k=> i as input to the transform with an appropriate down shift of index. The inverse binomial transform with a similar condition can produce each row from A000129. %F A228405 If using only the first two rows to initialize then the Pell equation produces each column, as: A(n,k) = 2*A(n-1, k) + A(n-2, k). %F A228405 If using only the left column (A000219(n) = Pell Numbers) to initialize then the following two equations will produce each row: %F A228405 A(n,k) = sqrt(2*A(n,k-1) + (-2)^(k-1)) for even rows %F A228405 A(n,k) = sqrt(2*A(n,k-1) - (-2)^(k-1)) for odd rows. %F A228405 Interestingly, any portion of the array can also be filled "backwards" given the top row and any column k, using only: A(n,k-1) = A(n-1,k-1) + A(n-1, k), or if given any column and its column number by rearranging the sqrt recursions above. %e A228405 With row # as n. and column # as k, and n, k =>0, the array begins: %e A228405 0, 1, 0, 2, 0, 4, 0, 8, ... %e A228405 1, 1, 2, 2, 4, 4, 8, 8, ... %e A228405 2, 3, 4, 6, 8, 12, 16, 24, ... %e A228405 5, 7, 10, 14, 20, 28, 40, 56, ... %e A228405 12, 17, 24, 34, 48, 68, 96, 136, ... %e A228405 29, 41, 58, 82, 116, 164, 232, 328, ... %e A228405 70, 99, 140, 198, 280, 396, 560, 792, ... %e A228405 169, 239, 338, 478, 676, 956, 1352, 1912, ... %e A228405 408, 577, 816, 1154, 1632, 2308, 3264, 4616, ... %Y A228405 Cf. A000219, A077957, A001333, A163271, A002203, A001109, A227972. %Y A228405 Main diagonal is A007070. %K A228405 nonn,tabl,look %O A228405 0,6 %A A228405 _Richard R. Forberg_, Aug 21 2013