A022089 -id:A022089 - cp's OEIS frontend

A035513 Wythoff array read by falling antidiagonals.

1, 2, 4, 3, 7, 6, 5, 11, 10, 9, 8, 18, 16, 15, 12, 13, 29, 26, 24, 20, 14, 21, 47, 42, 39, 32, 23, 17, 34, 76, 68, 63, 52, 37, 28, 19, 55, 123, 110, 102, 84, 60, 45, 31, 22, 89, 199, 178, 165, 136, 97, 73, 50, 36, 25, 144, 322, 288, 267, 220, 157, 118, 81, 58, 41, 27, 233, 521

Offset: 1

N. J. A. Sloane

T(0,0)=1, T(0,1)=2,...; y^2-x^2-xy

Inverse of sequence A064274 considered as a permutation of the nonnegative integers. - Howard A. Landman, Sep 25 2001

The Wythoff array W consists of all the Wythoff pairs (x(n),y(n)), where x=A000201 and y=A001950, so that W contains every positive integer exactly once. The differences T(i,2j+1)-T(i,2j) form the Wythoff difference array, A080164, which also contains every positive integer exactly once. - Clark Kimberling, Feb 08 2003

For n>2 the determinant of any n X n contiguous subarray of A035513 (as a square array) is 0. - Gerald McGarvey, Sep 18 2004

From Clark Kimberling, Nov 14 2007: (Start)

Except for initial terms in some cases:

(Row 1) = A000045

(Row 2) = A000032

(Row 3) = A006355

(Row 4) = A022086

(Row 5) = A022087

(Row 6) = A000285

(Row 7) = A022095

(Row 8) = A013655 (sum of Fibonacci and Lucas numbers)

(Row 9) = A022112

(Row 10-19) = A022113, A022120, A022121, A022379, A022130, A022382, A022088, A022136, A022137, A022089

(Row 20-28) = A022388, A022096, A022090, A022389, A022097, A022091, A022390, A022098, A022092

(Column 1) = A003622 = AA Wythoff sequence

(Column 2) = A035336 = BA Wythoff sequence

(Column 3) = A035337 = ABA Wythoff sequence

(Column 4) = A035338 = BBA Wythoff sequence

(Column 5) = A035339 = ABBA Wythoff sequence

(Column 6) = A035340 = BBBA Wythoff sequence

Main diagonal = A020941. (End)

The Wythoff array is the dispersion of the sequence given by floor(n*x+x-1), where x=(golden ratio). See A191426 for a discussion of dispersions. - Clark Kimberling, Jun 03 2011

If u and v are finite sets of numbers in a row of the Wythoff array such that (product of all the numbers in u) = (product of all the numbers in v), then u = v. See A160009 (row 1 products), A274286 (row 2), A274287 (row 3), A274288 (row 4). - Clark Kimberling, Jun 17 2016

All columns of the Wythoff array are compound Wythoff sequences. This follows from the main theorem in the 1972 paper by Carlitz, Scoville and Hoggatt. For an explicit expression see Theorem 10 in Kimberling's paper from 2008 in JIS. - Michel Dekking, Aug 31 2017

The Wythoff array can be viewed as an infinite graph over the set of nonnegative integers, built as follows: start with an empty graph; for all n = 0, 1, ..., create an edge between n and the sum of the degrees of all i < n. Finally, remove vertex 0. In the resulting graph, the connected components are chains and correspond to the rows of the Wythoff array. - Luc Rousseau, Sep 28 2017

Suppose that h < k are consecutive terms in a row of the Wythoff array. If k is in an even numbered column, then h = floor(k/tau); otherwise, h = -1 + floor(k/tau). - Clark Kimberling, Mar 05 2020

From Clark Kimberling, May 26 2020: (Start)

For k > = 0, column k shows the numbers m having F(k+1) as least term in the Zeckendorf representation of m. For n >= 1, let r(n,k) be the number of terms in column k that are <= n. Then n/r(n,k) = n/(F(k+1)*tau + F(k)*(n-1)), by Bottomley's formula, so that the limiting ratio is 1/(F(k+1)*tau + F(k)). Summing over all k gives Sum_{k>=0} 1/(F(k+1)*tau + F(k)) = 1. Thus, in the limiting sense:

38.19...% of the numbers m have least term 1;

23.60...% have least term 2;

14.58...% have least term 3;

9.01...% have least term 5, etc. (End)

Named after the Dutch mathematician Willem Abraham Wythoff (1865-1939). - Amiram Eldar, Jun 11 2021

From Clark Kimberling, Jun 04 2025: (Start)

Let u(n) = (T(n,1),T(n,2)) mod 2. The positive integers (A000027) are partitioned into 4 sets (sequences):

{n : u(n) = (0,0)} = (3, 5, 9, 15, 19, 25, 29,...) = 1 + 2*A190429

{n: u(n) = (0,1)} = (2, 6, 12, 16, 18, 22, 28,...) = A191331

{n : u(n) = (1,0)} = (1, 7, 11, 13, 17, 21, 23,...) = A086843

{n: u(n) = (1,1)} = (4, 8, 10, 14, 20, 24, 26,...) = A191330.

Let v(n) = (T(n,1),T(n,2)) mod 3. The positive integers are partitioned into 9 sets (sequences):

{n : v(n) = (0,0)} = (4, 13, 19, 28, 43, 52,...) = 1 + 3*A190434

{n: v(n) = (0,1)} = (3, 12, 27, 36, 42, 51,...) = 3*A140399

{n : v(n) = (0,2)} = (5, 11, 20, 35, 44, 50,...) = 2 + 3*A190439

{n: v(n) = (1,0)} = (9, 18, 24, 33, 48, 57,...) = 3*A140400

{n: v(n) = (1,1)} = (2, 8, 17, 26, 32, 41,...) = A384601

{n : v(n) = (1,2)} = (1, 10, 16, 25, 34, 40,...) = A384602

{n: v(n) = (2,0)} = (14, 23, 29, 38, 47, 53,...) = 2 + 3*A190438

{n: v(n) = (2,1)} = (7, 22, 31, 37, 46, 61,...) = 1 + 3*A190433

{n : v(n) = (2,2)} = (6, 15, 21, 30, 39, 45,...) = 3*A140398.

Conjecture: If m >= 2, then {(T(n,1), T(n,2)) mod m} has cardinality m^2. (End)

			The Wythoff array begins:
   1    2    3    5    8   13   21   34   55   89  144 ...
   4    7   11   18   29   47   76  123  199  322  521 ...
   6   10   16   26   42   68  110  178  288  466  754 ...
   9   15   24   39   63  102  165  267  432  699 1131 ...
  12   20   32   52   84  136  220  356  576  932 1508 ...
  14   23   37   60   97  157  254  411  665 1076 1741 ...
  17   28   45   73  118  191  309  500  809 1309 2118 ...
  19   31   50   81  131  212  343  555  898 1453 2351 ...
  22   36   58   94  152  246  398  644 1042 1686 2728 ...
  25   41   66  107  173  280  453  733 1186 1919 3105 ...
  27   44   71  115  186  301  487  788 1275 2063 3338 ...
  ...
The extended Wythoff array has two extra columns, giving the row number n and A000201(n), separated from the main array by a vertical bar:
0     1  |   1    2    3    5    8   13   21   34   55   89  144   ...
1     3  |   4    7   11   18   29   47   76  123  199  322  521   ...
2     4  |   6   10   16   26   42   68  110  178  288  466  754   ...
3     6  |   9   15   24   39   63  102  165  267  432  699 1131   ...
4     8  |  12   20   32   52   84  136  220  356  576  932 1508   ...
5     9  |  14   23   37   60   97  157  254  411  665 1076 1741   ...
6    11  |  17   28   45   73  118  191  309  500  809 1309 2118   ...
7    12  |  19   31   50   81  131  212  343  555  898 1453 2351   ...
8    14  |  22   36   58   94  152  246  398  644 1042 1686 2728   ...
9    16  |  25   41   66  107  173  280  453  733 1186 1919 3105   ...
10   17  |  27   44   71  115  186  301  487  788 1275 2063 3338   ...
11   19  |  30   49   79   ...
12   21  |  33   54   87   ...
13   22  |  35   57   92   ...
14   24  |  38   62   ...
15   25  |  40   65   ...
16   27  |  43   70   ...
17   29  |  46   75   ...
18   30  |  48   78   ...
19   32  |  51   83   ...
20   33  |  53   86   ...
21   35  |  56   91   ...
22   37  |  59   96   ...
23   38  |  61   99   ...
24   40  |  64   ...
25   42  |  67   ...
26   43  |  69   ...
27   45  |  72   ...
28   46  |  74   ...
29   48  |  77   ...
30   50  |  80   ...
31   51  |  82   ...
32   53  |  85   ...
33   55  |  88   ...
34   56  |  90   ...
35   58  |  93   ...
36   59  |  95   ...
37   61  |  98   ...
38   63  |     ...
   ...
Each row of the extended Wythoff array also satisfies the Fibonacci recurrence, and may be extended to the left using this recurrence backwards.
From _Peter Munn_, Jun 11 2021: (Start)
The Wythoff array appears to have the following relationship to the traditional Fibonacci rabbit breeding story, modified for simplicity to be a story of asexual reproduction.
Give each rabbit a number, 0 for the initial rabbit.
When a new round of rabbits is born, allocate consecutive numbers according to 2 rules (the opposite of many cultural rules for inheritance precedence): (1) newly born child of Rabbit 0 gets the next available number; (2) the descendants of a younger child of any given rabbit precede the descendants of an older child of the same rabbit.
Row n of the Wythoff array lists the children of Rabbit n (so Rabbit 0's children have the Fibonacci numbers: 1, 2, 3, 5, ...). The generation tree below shows rabbits 0 to 20. It is modified so that each round of births appears on a row.
                                                                 0
                                                                 :
                                       ,-------------------------:
                                       :                         :
                       ,---------------:                         1
                       :               :                         :
              ,--------:               2               ,---------:
              :        :               :               :         :
        ,-----:        3         ,-----:         ,-----:         4
        :     :        :         :     :         :     :         :
     ,--:     5     ,--:     ,---:     6     ,---:     7     ,---:
     :  :     :     :  :     :   :     :     :   :     :     :   :
  ,--:  8  ,--:  ,--:  9  ,--:  10  ,--:  ,--:  11  ,--:  ,--:  12
  :  :  :  :  :  :  :  :  :  :   :  :  :  :  :   :  :  :  :  :   :
  : 13  :  : 14  : 15  :  : 16   :  : 17  : 18   :  : 19  : 20   :
The extended array's nontrivial extra column (A000201) gives the number that would have been allocated to the first child of Rabbit n, if Rabbit n (and only Rabbit n) had started breeding one round early.
(End)

John H. Conway, Posting to Math Fun Mailing List, Nov 25 1996.
Clark Kimberling, "Stolarsky interspersions," Ars Combinatoria 39 (1995) 129-138.

Alois P. Heinz, Table of n, a(n) for n = 1..5151
Peter G. Anderson, More Properties of the Zeckendorf Array, Fib. Quart. 52-5 (2014), 15-21.
L. Carlitz, Richard Scoville, and V. E. Hoggatt, Jr., Fibonacci representations, Fib. Quart., Vol. 10, No. 1 (1972), pp. 1-28.
Code Golf Stack Exchange, New Order #5: where Fibonacci and Beatty meet at Wythoff.
J. H. Conway, Allan Wechsler, Marc LeBrun, Dan Hoey, and N. J. A. Sloane, On Kimberling Sums and Para-Fibonacci Sequences, Correspondence and Postings to Math-Fun Mailing List, Nov 1996 to Jan 1997.
John Conway and Alex Ryba, The extra Fibonacci series and the Empire State Building, Math. Intelligencer, Vol. 38, No. 1 (2016), pp. 41-48.
Larry Ericksen and Peter G. Anderson, Patterns in differences between rows in k-Zeckendorf arrays, The Fibonacci Quarterly, Vol. 50, No. 1 (February 2012), pp. 11-18. - _N. J. A. Sloane_, Jun 10 2012
Clark Kimberling, Interspersions.
Clark Kimberling, The Zeckendorf array equals the Wythoff array, Fibonacci Quarterly, Vol. 33, No. 1 (1995), pp. 3-8.
Clark Kimberling, Complementary equations and Wythoff Sequences, JIS, Vol. 11 (2008), Article 08.3.3.
Clark Kimberling, Lucas Representations of Positive Integers, J. Int. Seq., Vol. 23 (2020), Article 20.9.5.
Clark Kimberling and Kenneth B. Stolarsky, Slow Beatty sequences, devious convergence, and partitional divergence, Amer. Math. Monthly, 123, No. 2 (2016), pp. 267-273.
Stéphane Legendre, Labeled Fibonacci Trees, Fibonacci Quart. 53 (2015), no. 2, 152-167.
A. J. Macfarlane, On the fibbinary numbers and the Wythoff array, arXiv:2405.18128 [math.CO], 2024. See pages 1-2.
Casey Mongoven, Sonification of multiple Fibonacci-related sequences, Annales Mathematicae et Informaticae, Vol. 41 (2013), pp. 175-192.
N. J. A. Sloane, My favorite integer sequences, in Sequences and their Applications (Proceedings of SETA '98).
N. J. A. Sloane, Classic Sequences
Sam Vandervelde, On the divisibility of Fibonacci sequences by primes of index two, The Fibonacci Quarterly, Vol. 50, No. 3 (2012), pp. 207-216. See Figure 1.
Eric Weisstein's World of Mathematics, Wythoff Array.
Wikipedia, Wythoff array.
Pedro Zanetti, C++ code snippet, for generating this sequence.
Index entries for sequences that are permutations of the natural numbers

See comments above for more cross-references.
Cf. A003622, A064274 (inverse), A083412 (transpose), A000201, A001950, A080164, A003603, A265650, A019586 (row that contains n).
For two versions of the extended Wythoff array, see A287869, A287870.

W:= proc(n,k) Digits:= 100; (Matrix([n, floor((1+sqrt(5))/2* (n+1))]). Matrix([[0,1], [1,1]])^(k+1))[1,2] end: seq(seq(W(n, d-n), n=0..d), d=0..10); # Alois P. Heinz, Aug 18 2008
A035513 := proc(r, c)
    option remember;
    if c = 1 then
        A003622(r) ;
    else
        A022342(1+procname(r, c-1)) ;
    end if;
end proc:
seq(seq(A035513(r,d-r),r=1..d-1),d=2..15) ; # R. J. Mathar, Jan 25 2015

W[n_, k_] := Fibonacci[k + 1] Floor[n*GoldenRatio] + (n - 1) Fibonacci[k]; Table[ W[n - k + 1, k], {n, 12}, {k, n, 1, -1}] // Flatten

T(n,k)=(n+sqrtint(5*n^2))\2*fibonacci(k+1) + (n-1)*fibonacci(k)
for(k=0,9,for(n=1,k, print1(T(n,k+1-n)", "))) \\ Charles R Greathouse IV, Mar 09 2016

from sympy import fibonacci as F, sqrt
import math
tau = (sqrt(5) + 1)/2
def T(n, k): return F(k + 1)*int(math.floor(n*tau)) + F(k)*(n - 1)
for n in range(1, 11): print([T(k, n - k + 1) for k in range(1, n + 1)]) # Indranil Ghosh, Apr 23 2017

from math import isqrt, comb
from gmpy2 import fib2
def A035513(n):
    a = (m:=isqrt(k:=n<<1))+(k>m*(m+1))
    x = n-comb(a,2)
    b, c = fib2(a-x+2)
    return b*(x+isqrt(5*x*x)>>1)+c*(x-1) # Chai Wah Wu, Jun 26 2025

T(n, k) = Fib(k+1)*floor[n*tau]+Fib(k)*(n-1) where tau = (sqrt(5)+1)/2 = A001622 and Fib(n) = A000045(n). - Henry Bottomley, Dec 10 2001

T(n,-1) = n-1. T(n,0) = floor(n*tau). T(n,k) = T(n,k-1) + T(n,k-2) for k>=1. - R. J. Mathar, Sep 03 2016

Comments about the extended Wythoff array added by N. J. A. Sloane, Mar 07 2016

A022086 Fibonacci sequence beginning 0, 3.

Original entry on oeis.org

0, 3, 3, 6, 9, 15, 24, 39, 63, 102, 165, 267, 432, 699, 1131, 1830, 2961, 4791, 7752, 12543, 20295, 32838, 53133, 85971, 139104, 225075, 364179, 589254, 953433, 1542687, 2496120, 4038807, 6534927, 10573734, 17108661, 27682395, 44791056, 72473451, 117264507

Offset: 0

N. J. A. Sloane

First differences of A111314. - Ross La Haye, May 31 2006
Pisano period lengths: 1, 3, 1, 6, 20, 3, 16, 12, 8, 60, 10, 6, 28, 48, 20, 24, 36, 24, 18, 60, ... . - R. J. Mathar, Aug 10 2012
For n>=6, a(n) is the number of edge covers of the union of two cycles C_r and C_s, r+s=n, with a single common vertex. - Feryal Alayont, Oct 17 2024

A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 7,17.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (1,1).

Essentially the same as A097135.
Cf. A026390, A036999, A111314, A119457, A187893.
Sequences of the form Fibonacci(n+k) + Fibonacci(n-k) are listed in A280154.
Sequences of the form m*Fibonacci: A000045 (m=1), A006355 (m=2), this sequence (m=3), A022087 (m=4), A022088 (m=5), A022089 (m=6), A022090 (m=7), A022091 (m=8), A022092 (m=8), A022093 (m=10), A022345...A022366 (m=11...32).

[3*Fibonacci(n): n in [0..40]]; // Vincenzo Librandi, Dec 31 2016

BB := n->if n=0 then 3; > elif n=1 then 0; > else BB(n-2)+BB(n-1); > fi: > L:=[]: for k from 1 to 34 do L:=[op(L),BB(k)]: od: L; # Zerinvary Lajos, Mar 19 2007
with (combinat):seq(sum((fibonacci(n,1)),m=1..3),n=0..32); # Zerinvary Lajos, Jun 19 2008

LinearRecurrence[{1, 1}, {0, 3}, 40] (* Arkadiusz Wesolowski, Aug 17 2012 *)
Table[Fibonacci[n + 4] + Fibonacci[n - 4] - 4 Fibonacci[n], {n, 0, 40}] (* Bruno Berselli, Dec 30 2016 *)
Table[3 Fibonacci[n], {n, 0, 40}] (* Vincenzo Librandi, Dec 31 2016 *)

a(n)=3*fibonacci(n) \\ Charles R Greathouse IV, Nov 06 2014

def A022086(n): return 3*fibonacci(n)
print([A022086(n) for n in range(41)]) # G. C. Greubel, Apr 10 2025

a(n) = 3*Fibonacci(n).
a(n) = F(n-2) + F(n+2) for n>1, with F=A000045.
a(n) = round( ((6*phi-3)/5) * phi^n ) for n>2. - Thomas Baruchel, Sep 08 2004
a(n) = A119457(n+1,n-1) for n>1. - Reinhard Zumkeller, May 20 2006
G.f.: 3*x/(1-x-x^2). - Philippe Deléham, Nov 19 2008
a(n) = A187893(n) - 1. - Filip Zaludek, Oct 29 2016
E.g.f.: 6*sinh(sqrt(5)*x/2)*exp(x/2)/sqrt(5). - Ilya Gutkovskiy, Oct 29 2016
a(n) = F(n+4) + F(n-4) - 4*F(n), F = A000045. - Bruno Berselli, Dec 29 2016

A119457 Triangle read by rows: T(n, 1) = n, T(n, 2) = 2*(n-1) for n>1 and T(n, k) = T(n-1, k-1) + T(n-2, k-2) for 2 < k <= n.

Original entry on oeis.org

1, 2, 2, 3, 4, 3, 4, 6, 6, 5, 5, 8, 9, 10, 8, 6, 10, 12, 15, 16, 13, 7, 12, 15, 20, 24, 26, 21, 8, 14, 18, 25, 32, 39, 42, 34, 9, 16, 21, 30, 40, 52, 63, 68, 55, 10, 18, 24, 35, 48, 65, 84, 102, 110, 89, 11, 20, 27, 40, 56, 78, 105, 136, 165, 178, 144, 12, 22, 30, 45, 64, 91, 126, 170, 220, 267, 288, 233

Offset: 1

Reinhard Zumkeller, May 20 2006

			Triangle begins as:
   1;
   2,  2;
   3,  4,  3;
   4,  6,  6,  5;
   5,  8,  9, 10,  8;
   6, 10, 12, 15, 16, 13;
   7, 12, 15, 20, 24, 26,  21;
   8, 14, 18, 25, 32, 39,  42,  34;
   9, 16, 21, 30, 40, 52,  63,  68,  55;
  10, 18, 24, 35, 48, 65,  84, 102, 110,  89;
  11, 20, 27, 40, 56, 78, 105, 136, 165, 178, 144;
  12, 22, 30, 45, 64, 91, 126, 170, 220, 267, 288, 233;

G. C. Greubel, Rows n = 1..50 of the triangle, flattened
Eric Weisstein's World of Mathematics, Fibonacci Number

Main diagonal: A023607(n).
Sums: A001891 (row), A355020 (signed row).
Columns: A000027(n) (k=1), A005843(n-1) (k=2), A008585(n-2) (k=3), A008587(n-3) (k=4), A008590(n-4) (k=5), A008595(n-5) (k=6), A008603(n-6) (k=7).
Diagonals: A000045(n+1) (k=n), A006355(n+1) (k=n-1), A022086(n-1) (k=n-2), A022087(n-2) (k=n-3), A022088(n-3) (k=n-4), A022089(n-4) (k=n-5), A022090(n-5) (k=n-6).
Cf. A023652, A112469.

A119457:= func< n,k | (n-k+1)*Fibonacci(k+1) >;
[A119457(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Apr 16 2025

(* First program *)
T[n_, 1] := n;
T[n_ /; n > 1, 2] := 2 n - 2;
T[n_, k_] /; 2 < k <= n := T[n, k] = T[n - 1, k - 1] + T[n - 2, k - 2];
Table[T[n, k], {n, 1, 12}, {k, 1, n}] // Flatten (* Jean-François Alcover, Dec 01 2021 *)
(* Second program *)
A119457[n_,k_]:= (n-k+1)*Fibonacci[k+1];
Table[A119457[n,k], {n,13}, {k,n}]//Flatten (* G. C. Greubel, Apr 16 2025 *)

def A119457(n,k): return (n-k+1)*fibonacci(k+1)
print(flatten([[A119457(n,k) for k in range(1,n+1)] for n in range(1,13)])) # G. C. Greubel, Apr 16 2025

T(n, k) = (n-k+1)*T(k,k) for 1 <= k < n, with T(n, n) = A000045(n+1).
From G. C. Greubel, Apr 15 2025: (Start)
T(n, k) = (n-k+1)*Fibonacci(k+1).
Sum_{k=1..floor((n+1)/2)} T(n-k+1, k) = (1/2)*(1-(-1)^n)*A023652(floor((n+1)/2)) + (1+(-1)^n)*A001891(floor(n/2)).
Sum_{k=1..floor((n+1)/2)} (-1)^(k-1)*T(n-k+1, k) = (1/2)*(1-(-1)^n)*A112469(floor((n-1)/2)) + (1+(-1)^n)*A355020(floor((n-2)/2)). (End)

A126714 Dual Wythoff array read along antidiagonals.

Original entry on oeis.org

1, 2, 4, 3, 6, 7, 5, 10, 11, 9, 8, 16, 18, 14, 12, 13, 26, 29, 23, 19, 15, 21, 42, 47, 37, 31, 24, 17, 34, 68, 76, 60, 50, 39, 27, 20, 55, 110, 123, 97, 81, 63, 44, 32, 22, 89, 178, 199, 157, 131, 102, 71, 52, 35, 25, 144, 288, 322, 254, 212, 165, 115, 84, 57, 40, 28

Offset: 1

R. J. Mathar, Feb 12 2007

The dual Wythoff array is the dispersion of the sequence w given by w(n)=2+floor(n*x), where x=(golden ratio), so that w=2+A000201(n). For a discussion of dispersions, see A191426. - Clark Kimberling, Jun 03 2011

			Array starts
1 2 3 5 8 13 21 34 55 89 144
4 6 10 16 26 42 68 110 178 288 466
7 11 18 29 47 76 123 199 322 521 843
9 14 23 37 60 97 157 254 411 665 1076
12 19 31 50 81 131 212 343 555 898 1453
15 24 39 63 102 165 267 432 699 1131 1830
17 27 44 71 115 186 301 487 788 1275 2063
20 32 52 84 136 220 356 576 932 1508 2440
22 35 57 92 149 241 390 631 1021 1652 2673
25 40 65 105 170 275 445 720 1165 1885 3050
28 45 73 118 191 309 500 809 1309 2118 3427

Clark Kimberling, "Stolarsky Interspersions," Ars Combinatoria 39 (1995) 129-138. See page 135 for the dual Wythoff array and other dual arrays. - Clark Kimberling, Oct 29 2009

P. Hegarty, U. Larsson, Permutations of the natural numbers with prescribed difference multisets, Electr. J. Combin. Numb. Theory 6 (2006) #A03.

First three rows identical to A035506. First column is A007066. First row is A000045. 2nd row is essentially A006355. 3rd row is essentially A000032. 4th row essentially A000285. 5th row essentially A013655 or A001060. 6th row essentially A022086 or A097135. 7th row essentially A022120. 8th row essentially A022087. 9th row essentially A022130. 10th row essentially A022088. 11th row essentially A022095. 12th row essentially A022089 etc.

Cf. A035513 (Wythoff array).

Tn1 := proc(T,nmax,row) local n,r,c,fnd; n := 1; while true do fnd := false; for r from 1 to row do for c from 1 to nmax do if T[r,c] = n then fnd := true; fi; od; if T[r,nmax] < n then RETURN(-1); fi; od; if fnd then n := n+1; else RETURN(n); fi; od; end; Tn2 := proc(T,nmax,row,ai1) local n,r,c,fnd; for r from 1 to row do for c from 1 to nmax do if T[r,c]+1 = ai1 then RETURN(T[r,c+1]+1); fi; od; od; RETURN(-1); end; T := proc(nmax) local a,col,row; a := array(1..nmax,1..nmax); for col from 1 to nmax do a[1,col] := combinat[fibonacci](col+1); od; for row from 2 to nmax do a[row,1] := Tn1(a,nmax,row-1); a[row,2] := Tn2(a,nmax,row-1,a[row,1]); for col from 3 to nmax do a[row,col] := a[row,col-2]+a[row,col-1]; od; od; RETURN(a); end; nmax := 12; a := T(nmax); for d from 1 to nmax do for row from 1 to d do printf("%d, ",a[row,d-row+1]); od; od;

(* program generates the dispersion array T of the complement of increasing sequence f[n] *)
r = 40; r1 = 12;  (* r=# rows of T, r1=# rows to show *)
c = 40; c1 = 12;   (* c=# cols of T, c1=# cols to show *)
x = GoldenRatio; f[n_] := Floor[n*x + 2]
(* f(n) is complement of column 1 *)
mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
rows = {NestList[f, 1, c]};
Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
t[i_, j_] := rows[[i, j]];  (* the array T *)
TableForm[Table[t[i, j], {i, 1, 10}, {j, 1,10}]]
(* Dual Wythoff array, A126714 *)
Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* array as a sequence *)
(* Program by Peter J. C. Moses, Jun 01 2011; added here by Clark Kimberling, Jun 03 2011 *)

A192958 Constant term of the reduction by x^2 -> x+1 of the polynomial p(n,x) defined at Comments.

Original entry on oeis.org

1, -1, 3, 7, 17, 33, 61, 107, 183, 307, 509, 837, 1369, 2231, 3627, 5887, 9545, 15465, 25045, 40547, 65631, 106219, 171893, 278157, 450097, 728303, 1178451, 1906807, 3085313, 4992177, 8077549, 13069787, 21147399, 34217251, 55364717, 89582037

Offset: 0

Clark Kimberling, Jul 13 2011

sign

The titular polynomials are defined recursively: p(n,x) = x*p(n-1,x) - 2 + n^2, with p(0,x)=1. For an introduction to reductions of polynomials by substitutions such as x^2 -> x+1, see A192232 and A192744.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (3,-2,-1,1).

Cf. A000045, A192232, A192744, A192951, A192959.

F:=Fibonacci;; List([0..40], n-> 6*F(n+1)-(2*n+5)); # G. C. Greubel, Jul 12 2019

F:=Fibonacci; [6*F(n+1)-(2*n+5): n in [0..40]]; // G. C. Greubel, Jul 12 2019

(* First program *)
q = x^2; s = x + 1; z = 40;
p[0, x]:= 1;
p[n_, x_]:= x*p[n-1, x] + n^2 - 2;
Table[Expand[p[n, x]], {n, 0, 7}]
reduce[{p1_, q_, s_, x_}]:= FixedPoint[(s PolynomialQuotient @@ #1 + PolynomialRemainder @@ #1 &)[{#1, q, x}] &, p1]
t = Table[reduce[{p[n, x], q, s, x}], {n, 0, z}];
u1 = Table[Coefficient[Part[t, n], x, 0], {n, 1, z}] (* A192958 *)
u2 = Table[Coefficient[Part[t, n], x, 1], {n, 1, z}] (* A192959 *)
(* Second program *)
With[{F=Fibonacci}, Table[6*F[n+1]-(2*n+5), {n,0,40}]] (* G. C. Greubel, Jul 12 2019 *)

vector(40, n, n--; f=fibonacci; 6*f(n+1)-(2*n+5)) \\ G. C. Greubel, Jul 12 2019

f=fibonacci; [6*f(n+1)-(2*n+5) for n in (0..40)] # G. C. Greubel, Jul 12 2019

a(n) = 3*a(n-1) - 2*a(n-2) - a(n-3) + a(n-4).
From R. J. Mathar, May 09 2014: (Start)
G.f.: (1 -4*x +8*x^2 -3*x^3)/((1-x-x^2)*(1-x)^2).
a(n) - 2*a(n-1) +a(n-2) = A022089(n-3). (End)
a(n) = 6*Fibonacci(n+1) - (2*n+5). - G. C. Greubel, Jul 12 2019

A180251 Decimal expansion of 6*(phi+1)/5, where phi is (1 + sqrt(5))/2.

Original entry on oeis.org

3, 1, 4, 1, 6, 4, 0, 7, 8, 6, 4, 9, 9, 8, 7, 3, 8, 1, 7, 8, 4, 5, 5, 0, 4, 2, 0, 1, 2, 3, 8, 7, 6, 5, 7, 4, 1, 2, 6, 4, 3, 7, 1, 0, 1, 5, 7, 6, 6, 9, 1, 5, 4, 3, 4, 5, 6, 2, 5, 3, 8, 3, 4, 7, 2, 4, 6, 3, 1, 2, 5, 5, 5, 3, 8, 2, 6, 8, 2, 9, 3, 9, 6, 4, 8, 6, 4, 8, 6, 4, 5, 0, 2, 7, 2, 6, 9, 3, 6, 4, 9, 8, 1, 7, 0, 4, 9, 0, 5, 6, 9, 0, 4, 6

Offset: 1

Grant Garcia, Jan 16 2011

This is an approximation to Pi.

6*(phi+1)/5 is not equal to Pi, although some have claimed this (see Dudley). - Kellen Myers, Oct 04 2013

			3.141640786499873817845504201238765741264371015766915434562538347246312555382...

Underwood Dudley, Mathematical Cranks, MAA 1992, pp. 247, 292.
Alfred S. Posamentier and Ingmar Lehmann, The (Fabulous) Fibonacci Numbers, New York, Prometheus Books, 2007, p. 119.

G. C. Greubel, Table of n, a(n) for n = 1..10000
Hung Viet Chu, Square the Circle in One Minute, arXiv:1908.01202 [math.GM], 2019.
Futility Closet, A Surprise Visitor

Cf. A022088, A022089, A001622, A000796.

(3/10)*(1 + Sqrt(5))^2; // G. C. Greubel, Jan 17 2018

RealDigits[(6/5)GoldenRatio^2, 10, 100][[1]] (* Alonso del Arte, Apr 09 2012 *)

3*(3+sqrt(5))/5 \\ Charles R Greathouse IV, Sep 13 2013

Limit of A022089(n+2)/A022088(n) as n approaches infinity.
6*(phi + 1)/5 = 6*phi^2/5 = 3(3 + sqrt(5))/5 = 9/5 + sqrt(9/5). - Charles R Greathouse IV, Sep 13 2013
Equals 24/(5-sqrt(5))^2. - Joost Gielen, Sep 20 2013

Showing 1-6 of 6 results.

cp's OEIS Frontend

A035513 Wythoff array read by falling antidiagonals.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Python

Formula

Extensions

A022086 Fibonacci sequence beginning 0, 3.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

SageMath

Formula

A119457 Triangle read by rows: T(n, 1) = n, T(n, 2) = 2*(n-1) for n>1 and T(n, k) = T(n-1, k-1) + T(n-2, k-2) for 2 < k <= n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Magma

Mathematica

SageMath

Formula

A126714 Dual Wythoff array read along antidiagonals.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

A192958 Constant term of the reduction by x^2 -> x+1 of the polynomial p(n,x) defined at Comments.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Mathematica

PARI

Sage

Formula

A180251 Decimal expansion of 6*(phi+1)/5, where phi is (1 + sqrt(5))/2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs