OFFSET
1,2
COMMENTS
T(0,0)=1, T(0,1)=2,...; y^2-x^2-xy<y if and only if there exist (i,j) with x=T(i,2j) and y=T(i,2j+1). - Claude Lenormand (claude.lenormand(AT)free.fr), Mar 17 2001
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
(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
REFERENCES
John H. Conway, Posting to Math Fun Mailing List, Nov 25 1996.
Clark Kimberling, "Stolarsky interspersions," Ars Combinatoria 39 (1995) 129-138.
LINKS
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.
Pedro Zanetti, C++ code snippet, for generating this sequence.
FORMULA
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
EXAMPLE
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)
MAPLE
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
MATHEMATICA
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
PROG
(PARI) 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
(Python)
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
CROSSREFS
KEYWORD
AUTHOR
EXTENSIONS
Comments about the extended Wythoff array added by N. J. A. Sloane, Mar 07 2016
STATUS
approved