A022130 -id:A022130 - cp's OEIS frontend

A028387 a(n) = n + (n+1)^2.

1, 5, 11, 19, 29, 41, 55, 71, 89, 109, 131, 155, 181, 209, 239, 271, 305, 341, 379, 419, 461, 505, 551, 599, 649, 701, 755, 811, 869, 929, 991, 1055, 1121, 1189, 1259, 1331, 1405, 1481, 1559, 1639, 1721, 1805, 1891, 1979, 2069, 2161, 2255, 2351, 2449, 2549, 2651

Offset: 0

Patrick De Geest

a(n+1) is the least k > a(n) + 1 such that A000217(a(n)) + A000217(k) is a square. - David Wasserman, Jun 30 2005
Values of Fibonacci polynomial n^2 - n - 1 for n = 2, 3, 4, 5, ... - Artur Jasinski, Nov 19 2006
A127701 * [1, 2, 3, ...]. - Gary W. Adamson, Jan 24 2007
Row sums of triangle A135223. - Gary W. Adamson, Nov 23 2007
Equals row sums of triangle A143596. - Gary W. Adamson, Aug 26 2008
a(n-1) gives the number of n X k rectangles on an n X n chessboard (for k = 1, 2, 3, ..., n). - Aaron Dunigan AtLee, Feb 13 2009
sqrt(a(0) + sqrt(a(1) + sqrt(a(2) + sqrt(a(3) + ...)))) = sqrt(1 + sqrt(5 + sqrt(11 + sqrt(19 + ...)))) = 2. - Miklos Kristof, Dec 24 2009
When n + 1 is prime, a(n) gives the number of irreducible representations of any nonabelian group of order (n+1)^3. - Andrew Rupinski, Mar 17 2010
a(n) = A176271(n+1, n+1). - Reinhard Zumkeller, Apr 13 2010
The product of any 4 consecutive integers plus 1 is a square (see A062938); the terms of this sequence are the square roots. - Harvey P. Dale, Oct 19 2011
Or numbers not expressed in the form m + floor(sqrt(m)) with integer m. - Vladimir Shevelev, Apr 09 2012
Left edge of the triangle in A214604: a(n) = A214604(n+1,1). - Reinhard Zumkeller, Jul 25 2012
Another expression involving phi = (1 + sqrt(5))/2 is a(n) = (n + phi)(n + 1 - phi). Therefore the numbers in this sequence, even if they are prime in Z, are not prime in Z[phi]. - Alonso del Arte, Aug 03 2013
a(n-1) = n*(n+1) - 1, n>=0, with a(-1) = -1, gives the values for a*c of indefinite binary quadratic forms [a, b, c] of discriminant D = 5 for b = 2*n+1. In general D = b^2 - 4ac > 0 and the form [a, b, c] is a*x^2 + b*x*y + c*y^2. - Wolfdieter Lang, Aug 15 2013
a(n) has prime factors given by A038872. - Richard R. Forberg, Dec 10 2014
A253607(a(n)) = -1. - Reinhard Zumkeller, Jan 05 2015
An example of a quadratic sequence for which the continued square root map (see A257574) produces the number 2. There are infinitely many sequences with this property - another example is A028387. See Popular Computing link. - N. J. A. Sloane, May 03 2015
Left edge of the triangle in A260910: a(n) = A260910(n+2,1). - Reinhard Zumkeller, Aug 04 2015
Numbers m such that 4m+5 is a square. - Bruce J. Nicholson, Jul 19 2017
The numbers represented as 131 in base n: 131_4 = 29, 131_5 = 41, ... . If 'digits' larger than the base are allowed then 131_2 = 11 and 131_1 = 5 also. - Ron Knott, Nov 14 2017
From Klaus Purath, Mar 18 2019: (Start)
Let m be a(n) or a prime factor of a(n). Then, except for 1 and 5, there are, if m is a prime, exactly two squares y^2 such that the difference y^2 - m contains exactly one pair of factors {x,z} such that the following applies: x*z = y^2 - m, x + y = z with
x < y, where {x,y,z} are relatively prime numbers. {x,y,z} are the initial values of a sequence of the Fibonacci type. Thus each a(n) > 5, if it is a prime, and each prime factor p > 5 of an a(n) can be assigned to exactly two sequences of the Fibonacci type. a(0) = 1 belongs to the original Fibonacci sequence and a(1) = 5 to the Lucas sequence.
But also the reverse assignment applies. From any sequence (f(i)) of the Fibonacci type we get from its 3 initial values by f(i)^2 - f(i-1)*f(i+1) with f(i-1) < f(i) a term a(n) or a prime factor p of a(n). This relation is also valid for any i. In this case we get the absolute value |a(n)| or |p|. (End)
a(n-1) = 2*T(n) - 1, for n>=1, with T = A000217, is a proper subsequence of A089270, and the terms are 0,-1,+1 (mod 5). - Wolfdieter Lang, Jul 05 2019
a(n+1) is the number of wedged n-dimensional spheres in the homotopy of the neighborhood complex of Kneser graph KG_{2,n}. Here, KG_{2,n} is a graph whose vertex set is the collection of subsets of cardinality 2 of set {1,2,...,n+3,n+4} and two vertices are adjacent if and only if they are disjoint. - Anurag Singh, Mar 22 2021
Also the number of squares between (n+2)^2 and (n+2)^4. - Karl-Heinz Hofmann, Dec 07 2021
(x, y, z) = (A001105(n+1), -a(n-1), -a(n)) are solutions of the Diophantine equation x^3 + 4*y^3 + 4*z^3 = 8. - XU Pingya, Apr 25 2022
The least significant digit of terms of this sequence cycles through 1, 5, 1, 9, 9. - Torlach Rush, Jun 05 2024

			From _Ilya Gutkovskiy_, Apr 13 2016: (Start)
Illustration of initial terms:
                                        o               o
                        o           o   o o           o o
            o       o   o o       o o   o o o       o o o
    o   o   o o   o o   o o o   o o o   o o o o   o o o o
o   o o o   o o o o o   o o o o o o o   o o o o o o o o o
n=0  n=1       n=2           n=3               n=4
(End)
From _Klaus Purath_, Mar 18 2019: (Start)
Examples:
a(0) = 1: 1^1-0*1 = 1, 0+1 = 1 (Fibonacci A000045).
a(1) = 5: 3^2-1*4 = 5, 1+3 = 4 (Lucas A000032).
a(2) = 11: 4^2-1*5 = 11, 1+4 = 5 (A000285); 5^2-2*7 = 11, 2+5 = 7 (A001060).
a(3) = 19: 5^2-1*6 = 19, 1+5 = 6 (A022095); 7^2-3*10 = 19, 3+7 = 10 (A022120).
a(4) = 29: 6^2-1*7 = 29, 1+6 = 7 (A022096); 9^2-4*13 = 29, 4+9 = 13 (A022130).
a(11)/5 = 31: 7^2-2*9 = 31, 2+7 = 9 (A022113); 8^2-3*11 = 31, 3+8 = 11 (A022121).
a(24)/11 = 59: 9^2-2*11 = 59, 2+9 = 11 (A022114); 12^2-5*17 = 59, 5+12 = 17 (A022137).
(End)

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Patrick De Geest, World!Of Numbers, Palindromes of the form n+(n+1)^2.
Adalbert Kerber, A matrix of combinatorial numbers related to the symmetric groups, Discrete Math., 21 (1978), 319-321. [Annotated scanned copy]
Clark Kimberling, Complementary Equations, Journal of Integer Sequences, Vol. 10 (2007), Article 07.1.4.
Nandini Nilakantan and Anurag Singh, Homotopy type of neighborhood complexes of Kneser graphs, KG_{2,k}, Proceeding-Mathematical Sciences, 128, Article number: 53(2018).
Yanni Pei and Jiang Zeng, Counting signed derangements with right-to-left minima and excedances, arXiv:2206.11236 [math.CO], 2022.
Popular Computing (Calabasas, CA), The CSR Function, Vol. 4 (No. 34, Jan 1976), pages PC34-10 to PC34-11. Annotated and scanned copy.
Zdzislaw Skupień and Andrzej Żak, Pair-sums packing and rainbow cliques, in Topics In Graph Theory, A tribute to A. A. and T. E. Zykovs on the occasion of A. A. Zykov's 90th birthday, ed. R. Tyshkevich, Univ. Illinois, 2013, pages 131-144, (in English and Russian).
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Complement of A028392. Third column of array A094954.
Cf. A000217, A002522, A062392, A062786, A127701, A135223, A143596, A052905, A162997, A062938 (squares of this sequence).
A110331 and A165900 are signed versions.
Cf. A002327 (primes), A094210.
Frobenius number for k successive numbers: this sequence (k=2), A079326 (k=3), A138984 (k=4), A138985 (k=5), A138986 (k=6), A138987 (k=7), A138988 (k=8).
Cf. also A135223, A176271, A214604, A105728, A005408, A002378, A084990, A253607, A260910, A089270.

a028387 n = n + (n + 1) ^ 2  -- Reinhard Zumkeller, Jul 17 2014

[n + (n+1)^2: n in [0..60]]; // Vincenzo Librandi, Apr 26 2011

FoldList[## + 2 &, 1, 2 Range@ 45] (* Robert G. Wilson v, Feb 02 2011 *)
Table[n + (n + 1)^2, {n, 0, 100}] (* Vincenzo Librandi, Oct 17 2012 *)
Table[ FrobeniusNumber[{n, n + 1}], {n, 2, 30}] (* Zak Seidov, Jan 14 2015 *)

a(n)=n^2+3*n+1 \\ Charles R Greathouse IV, Jun 10 2011

def a(n): return (n**2+3*n+1) # Torlach Rush, May 07 2024

[n+(n+1)^2 for n in range(0,48)] # Zerinvary Lajos, Jul 03 2008

a(n) = sqrt(A062938(n)). - Floor van Lamoen, Oct 08 2001
a(0) = 1, a(1) = 5, a(n) = (n+1)*a(n-1) - (n+2)*a(n-2) for n > 1. - Gerald McGarvey, Sep 24 2004
a(n) = A105728(n+2, n+1). - Reinhard Zumkeller, Apr 18 2005
a(n) = A109128(n+2, 2). - Reinhard Zumkeller, Jun 20 2005
a(n) = 2*T(n+1) - 1, where T(n) = A000217(n). - Gary W. Adamson, Aug 15 2007
a(n) = A005408(n) + A002378(n); A084990(n+1) = Sum_{k=0..n} a(k). - Reinhard Zumkeller, Aug 20 2007
Binomial transform of [1, 4, 2, 0, 0, 0, ...] = (1, 5, 11, 19, ...). - Gary W. Adamson, Sep 20 2007
G.f.: (1+2*x-x^2)/(1-x)^3. a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - R. J. Mathar, Jul 11 2009
a(n) = (n + 2 + 1/phi) * (n + 2 - phi); where phi = 1.618033989... Example: a(3) = 19 = (5 + .6180339...) * (3.381966...). Cf. next to leftmost column in A162997 array. - Gary W. Adamson, Jul 23 2009
a(n) = a(n-1) + 2*(n+1), with n > 0, a(0) = 1. - Vincenzo Librandi, Nov 18 2010
For k < n, a(n) = (k+1)*a(n-k) - k*a(n-k-1) + k*(k+1); e.g., a(5) = 41 = 4*11 - 3*5 + 3*4. - Charlie Marion, Jan 13 2011
a(n) = lower right term in M^2, M = the 2 X 2 matrix [1, n; 1, (n+1)]. - Gary W. Adamson, Jun 29 2011
G.f.: (x^2-2*x-1)/(x-1)^3 = G(0) where G(k) = 1 + x*(k+1)*(k+4)/(1 - 1/(1 + (k+1)*(k+4)/G(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 16 2012
Sum_{n>0} 1/a(n) = 1 + Pi*tan(sqrt(5)*Pi/2)/sqrt(5). - Enrique Pérez Herrero, Oct 11 2013
E.g.f.: exp(x) (1+4*x+x^2). - Tom Copeland, Dec 02 2013
a(n) = A005408(A000217(n)). - Tony Foster III, May 31 2016
From Amiram Eldar, Jan 29 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2).
Product_{n>=1} (1 - 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2)/6. (End)
a(5*n+1)/5 = A062786(n+1). - Torlach Rush, Jun 05 2024

Minor edits by N. J. A. Sloane, Jul 04 2010, following suggestions from the Sequence Fans Mailing List

A035513 Wythoff array read by falling antidiagonals.

Original entry on oeis.org

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

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 *)

A127830 a(n) = Sum_{k=0..n} (binomial(floor(k/2),n-k) mod 2).

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 2, 3, 3, 3, 2, 2, 3, 2, 3, 5, 5, 4, 4, 5, 4, 3, 3, 3, 4, 4, 3, 4, 5, 3, 5, 8, 8, 7, 6, 7, 7, 5, 6, 8, 7, 6, 5, 5, 5, 4, 4, 5, 6, 5, 5, 7, 6, 4, 5, 6, 7, 7, 5, 6, 8, 5, 8, 13, 13, 11, 10, 12, 11, 8, 9, 11, 11, 10, 8, 9, 10, 7, 9, 13, 12

Offset: 0

Paul Barry, Feb 01 2007

Row sums of number triangle A127829.
From Johannes W. Meijer, Jun 05 2011: (Start)
The Ze3 and Ze4 triangle sums, see A180662 for their definitions, of Sierpinski's triangle A047999 equal this sequence.
The sequences A127830(2^n-p), p>=0, are apparently all Fibonacci like sequences, i.e., the next term is the sum of the two nonzero terms that precede it; see the crossrefs. (End)

Cf.: A000045 (p=0), A000204 (p=7), A001060 (p=13), A000285 (p=14), A022095 (p=16), A022120 (p=24), A022121 (p=25), A022113 (p=28), A022096 (p=30), A022097 (p=31), A022098 (p=32), A022130 (p=44), A022137 (p=48), A022138 (p=49), A022122 (p=52), A022114 (p=53), A022123 (p=56), A022115 (p=60), A022100 (p=62), A022101 (p=63), A022103 (p=64), A022136 (p=79), A022388 (p=80), A022389 (p=88). - Johannes W. Meijer, Jun 05 2011

A127830 := proc(n) local k: option remember: add(binomial(floor(k/2), n-k) mod 2, k=0..n) end: seq(A127830(n), n=0..80); # Johannes W. Meijer, Jun 05 2011

Table[Sum[Mod[Binomial[Floor[k/2],n-k],2],{k,0,n}],{n,0,80}] (* James C. McMahon, Jan 04 2025 *)

def A127830(n): return sum(not ~(k>>1)&n-k for k in range(n+1)) # Chai Wah Wu, Jul 29 2025

a(2^n) = F(n); a(2^(n+1)+1) = L(n).

a(n) mod 2 = A000931(n+5) mod 2 = A011656(n+4).

A199535 Clark Kimberling's even first column Stolarsky array read by antidiagonals.

Original entry on oeis.org

1, 2, 4, 3, 7, 6, 5, 11, 9, 10, 8, 18, 15, 17, 12, 13, 29, 24, 27, 19, 14, 21, 47, 39, 44, 31, 23, 16, 34, 76, 63, 71, 50, 37, 25, 20, 55, 123, 102, 115, 81, 60, 41, 33, 22, 89, 199, 165, 186, 131, 97, 66, 53, 35, 26, 144, 322, 267, 301, 212, 157, 107, 86, 57, 43, 28

Offset: 1

Casey Mongoven, Nov 07 2011

The rows of the array can be seen to have the form A(n, k) = p(n)*Fibonacci(k) + q(n)*Fibonacci(k+1) where p(n) is the sequence {0, 1, 3, 3, 3, 5, 7, 7, 9, 9, 11, 11, 13, 13, 15, 15, 17, ...}{n >= 1} and q(n) is the sequence {1, 3, 3, 7, 2, 9, 9, 13, 13, 17, 17, 19, 19, 23, 23, 25, ...}{n >= 1}. - G. C. Greubel, Jun 23 2022

			The even first column stolarsky array (EFC array), northwest corner:
  1......2.....3.....5.....8....13....21....34....55....89...144 ... A000045;
  4......7....11....18....29....47....76...123...199...322...521 ... A000032;
  6......9....15....24....39....63...102...165...267...432...699 ... A022086;
  10....17....27....44....71...115...186...301...487...788..1275 ... A022120;
  12....19....31....50....81...131...212...343...555...898..1453 ... A013655;
  14....23....37....60....97...157...254...411...665..1076..1741 ... A000285;
  16....25....41....66...107...173...280...453...733..1186..1919 ... A022113;
  20....33....53....86...139...225...364...589...953..1542..2495 ... A022096;
  22....35....57....92...149...241...390...631..1021..1652..2673 ... A022130;
Antidiagonal rows (T(n, k)):
   1;
   2,   4;
   3,   7,   6;
   5,  11,   9,  10;
   8,  18,  15,  17, 12;
  13,  29,  24,  27, 19, 14;
  21,  47,  39,  44, 31, 23, 16;
  34,  76,  63,  71, 50, 37, 25, 20;
  55, 123, 102, 115, 81, 60, 41, 33, 22;

Clark Kimberling, The first column of an interspersion, Fibonacci Quarterly 32 (1994), pp. 301-314.

Cf. A000032, A000045, A000285, A013655, A022086, A022088, A022095.
Cf. A022096, A022097, A022113, A022120, A022121, A022122, A022130.
Cf. A022388, A022390, A035506, A035513, A097135, A199536, A199537.

From G. C. Greubel, Jun 23 2022: (Start)
T(n, 1) = A000045(n+1).
T(n, 2) = A000032(n+1), n >= 2.
T(n, 3) = A022086(n) = A097135(n), n >= 3.
T(n, 4) = A022120(n-2), n >= 4.
T(n, 5) = A013655(n-1), n >= 5.
T(n, 6) = A000285(n-2), n >= 6.
T(n, 7) = A022113(n-4), n >= 7.
T(n, 8) = A022096(n-4), n >= 8.
T(n, 9) = A022130(n-6), n >= 9.
T(n, 10) = A022098(n-5), n >= 10.
T(n, 11) = A022095(n-7), n >= 11.
T(n, 12) = A022121(n-8), n >= 12.
T(n, 13) = A022388(n-10), n >= 13.
T(n, 14) = A022122(n-10), n >= 14.
T(n, 15) = A022097(n-10), n >= 15.
T(n, 16) = A022088(n-10), n >= 16.
T(n, 17) = A022390(n-14), n >= 17.
T(n, n) = A199536(n).
T(n, n-1) = A199537(n-1), n >= 2. (End)

More terms added by G. C. Greubel, Jun 23 2022

A280931 a(n) = 2F(n-1) + 9F(n-4) + 9*F(n-7) where n >= 7 and F = A000045.

Original entry on oeis.org

34, 62, 96, 158, 254, 412, 666, 1078, 1744, 2822, 4566, 7388, 11954, 19342, 31296, 50638, 81934, 132572, 214506, 347078, 561584, 908662, 1470246, 2378908, 3849154, 6228062, 10077216, 16305278, 26382494, 42687772, 69070266, 111758038, 180828304, 292586342

Offset: 7

Vincenzo Librandi, Jan 24 2017

Vincenzo Librandi, Table of n, a(n) for n = 7..1100
H. Zhao and X. Li, On the Fibonacci numbers of trees, The Fibonacci Quarterly, Vol. 44, Number 1 (2006), page 37.
Index entries for linear recurrences with constant coefficients, signature (1,1).

Cf. A000045, A022125, A022130, A101156, A280932.

[2*Fibonacci(n-1)+9*Fibonacci(n-4)+9*Fibonacci(n-7): n in [7..40]];

a0:=34; a1:=62; [GeneralizedFibonacciNumber(a0, a1, n): n in [0..40]];

Mathematica
```
LinearRecurrence[{1, 1}, {34, 62}, 35]
```

G.f.: 2*x^7*(17 + 14*x)/(1 - x - x^2).
a(n) = a(n-1) + a(n-2).
From the g.f.: a(n) = 34*F(n-6) + 28*F(n-7) = 28*F(n-5) + 6*F(n-6) = 6*F(n-4) + 22*F(n-5) = 22*F(n-3) - 16*F(n-4) = -16*F(n-2) + 38*F(n-3) = 38*F(n-1) - 54*F(n-2) = -54*F(n) + 92*F(n-1), and so on.
a(n) = 2*A022125(n-5).
a(n) = F(n+2) + F(n-3) + F(n-11). - Greg Dresden, Jul 07 2022

Corrected and extended by Bruno Berselli, Jan 24 2017

A280932 a(n) = 2F(n-1) + 2F(n-3) + 10F(n-5) + 9F(n-8) where n >= 8 and F = A000045.

Original entry on oeis.org

56, 97, 153, 250, 403, 653, 1056, 1709, 2765, 4474, 7239, 11713, 18952, 30665, 49617, 80282, 129899, 210181, 340080, 550261, 890341, 1440602, 2330943, 3771545, 6102488, 9874033, 15976521, 25850554, 41827075, 67677629, 109504704, 177182333, 286687037

Offset: 8

Vincenzo Librandi, Jan 24 2017

Vincenzo Librandi, Table of n, a(n) for n = 8..1100
H. Zhao and X. Li, On the Fibonacci numbers of trees, The Fibonacci Quarterly, Vol. 44, Number 1 (2006), page 37.
Index entries for linear recurrences with constant coefficients, signature (1,1).

Cf. A000045, A022130, A101156, A280931.

[2*Fibonacci(n-1)+2*Fibonacci(n-3)+10*Fibonacci(n-5)+9*Fibonacci(n-8): n in [8..40]];

a0:=56; a1:=97; [GeneralizedFibonacciNumber(a0, a1, n): n in [0..40]];

Mathematica
```
LinearRecurrence[{1, 1}, {56, 97}, 35]
```

G.f.: x^8*(56 + 41*x)/(1 - x - x^2).
a(n) = a(n-1) + a(n-2).
From the g.f.: a(n) = 56*F(n-7) + 41*F(n-8) = 41*F(n-6) + 15*F(n-7) = 15*F(n-5) + 26*F(n-6) = 26*F(n-4) - 11*F(n-5) = -11*F(n-3) + 37*F(n-4) = 37*F(n-2) - 48*F(n-3) = -48*F(n-1) + 85*F(n-2) = 85*F(n) - 133*F(n-1), and so on.

Corrected and extended by Bruno Berselli, Jan 24 2017

A187892 a(0)=0, a(1)=6, a(n)=a(n-1)+a(n-2)-1.

Original entry on oeis.org

0, 6, 5, 10, 14, 23, 36, 58, 93, 150, 242, 391, 632, 1022, 1653, 2674, 4326, 6999, 11324, 18322, 29645, 47966, 77610, 125575, 203184, 328758, 531941, 860698, 1392638, 2253335, 3645972, 5899306, 9545277, 15444582, 24989858, 40434439, 65424296, 105858734

Offset: 0

Vladimir Joseph Stephan Orlovsky, Mar 15 2011

Index entries for linear recurrences with constant coefficients, signature (2,0,-1).

Cf. A000071, A001611, A001612, A187890.

Join[{a=0,b=6},Table[c=a+b-1;a=b;b=c,{n,100}]]
LinearRecurrence[{2,0,-1},{0,6,5},40] (* Harvey P. Dale, Aug 17 2019 *)

G.f.: -x^2*(-6+7*x) / ( (x-1)*(x^2+x-1) ). - R. J. Mathar, Mar 15 2011

a(n) = 1+A022130(n-3), n>2. - R. J. Mathar, Mar 15 2011

Showing 1-8 of 8 results.

cp's OEIS Frontend

A028387 a(n) = n + (n+1)^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

Python

Sage

Formula

Extensions

A035513 Wythoff array read by falling antidiagonals.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Python

Formula

Extensions

A126714 Dual Wythoff array read along antidiagonals.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

A127830 a(n) = Sum_{k=0..n} (binomial(floor(k/2),n-k) mod 2).

Views

Author

Keywords

Comments

Crossrefs

Programs

Maple

Mathematica

Python

Formula

A199535 Clark Kimberling's even first column Stolarsky array read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A280931 a(n) = 2*F(n-1) + 9*F(n-4) + 9*F(n-7) where n >= 7 and F = A000045.

Views

Author

Keywords

Links

Crossrefs

Programs

Magma

Magma

Mathematica

Formula

Extensions

A280931 a(n) = 2F(n-1) + 9F(n-4) + 9*F(n-7) where n >= 7 and F = A000045.

A280932 a(n) = 2F(n-1) + 2F(n-3) + 10F(n-5) + 9F(n-8) where n >= 8 and F = A000045.