cp's OEIS Frontend

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.

Showing 1-6 of 6 results.

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

Views

Author

Keywords

Comments

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

Examples

			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)
		

References

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

Crossrefs

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.

Programs

  • 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
  • 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
    
  • Python
    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

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

Extensions

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

Keywords

Comments

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

References

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

Crossrefs

Essentially the same as A097135.
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).

Programs

  • Magma
    [3*Fibonacci(n): n in [0..40]]; // Vincenzo Librandi, Dec 31 2016
    
  • Maple
    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
  • Mathematica
    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 *)
  • PARI
    a(n)=3*fibonacci(n) \\ Charles R Greathouse IV, Nov 06 2014
    
  • SageMath
    def A022086(n): return 3*fibonacci(n)
    print([A022086(n) for n in range(41)]) # G. C. Greubel, Apr 10 2025

Formula

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

Author

Reinhard Zumkeller, May 20 2006

Keywords

Examples

			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;
		

Crossrefs

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

Programs

  • Magma
    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
    
  • Mathematica
    (* 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 *)
  • SageMath
    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

Formula

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)

A280154 a(n) = 5*Lucas(n).

Original entry on oeis.org

10, 5, 15, 20, 35, 55, 90, 145, 235, 380, 615, 995, 1610, 2605, 4215, 6820, 11035, 17855, 28890, 46745, 75635, 122380, 198015, 320395, 518410, 838805, 1357215, 2196020, 3553235, 5749255, 9302490, 15051745, 24354235, 39405980, 63760215, 103166195, 166926410, 270092605, 437019015
Offset: 0

Author

Bruno Berselli, Dec 27 2016

Keywords

Comments

Fibonacci sequence beginning 10, 5.
After 5, the sequence provides the 3rd column of the rectangular array in A213590.
After 5, all terms belong to A191921 because a(n) = Lucas(n+4) - 3*Lucas(n-1).
From G. C. Greubel, Dec 27 2016: (Start)
{a(n) mod 3} yields (1,2,0,2,2,1,0,1), repeated, and is given as A082115.
{a(n) mod 6} yields (4,5,3,2,5,1,0,1,1,2,3,5,2,1,3,4,1,5,0,5,5,4,3,1) and is given as A082117. (End)

Crossrefs

Subsequence of A084176.
Cf. A022088: 5*Fibonacci(n).
Cf. A022359: Lucas(n+5) + Lucas(n-5).
Cf. sequences with formula Fibonacci(n+k) + Fibonacci(n-k): A006355 (k=0, without the initial 1), A000032 (k=1), A022086 (k=2), A022112 (k=3, with an initial 4), A022090 (k=4), this sequence (k=5), A022352 (k=6).

Programs

  • Magma
    [5*Lucas(n): n in [0..40]];
    
  • Maple
    F := n -> combinat:-fibonacci(n):
    seq(F(n+5) + F(n-5), n=0..38); # Peter Luschny, Dec 29 2016
  • Mathematica
    Table[5 LucasL[n], {n, 0, 40}]
  • PARI
    vector(40, n, n--; fibonacci(n+5)+fibonacci(n-5))
    
  • Sage
    def A280154():
        x, y = 10, 5
        while True:
            yield x
            x, y = y, x + y
    a = A280154(); print([next(a) for  in range(39)]) # _Peter Luschny, Dec 29 2016

Formula

G.f.: 5*(2 - x)/(1 - x - x^2).
a(n) = a(n-1) + a(n-2) for n>1.
a(n) = Fibonacci(n+5) + Fibonacci(n-5), with Fibonacci(-k) = -(-1)^k*Fibonacci(k) for the negative indices.

A294157 Fibonacci sequence beginning 2, 8.

Original entry on oeis.org

2, 8, 10, 18, 28, 46, 74, 120, 194, 314, 508, 822, 1330, 2152, 3482, 5634, 9116, 14750, 23866, 38616, 62482, 101098, 163580, 264678, 428258, 692936, 1121194, 1814130, 2935324, 4749454, 7684778, 12434232, 20119010, 32553242, 52672252, 85225494, 137897746, 223123240
Offset: 0

Author

Bruno Berselli, Oct 24 2017

Keywords

References

  • Steven Vajda, Fibonacci and Lucas Numbers, and the Golden Section: Theory and Applications, Dover Publications (2008), page 24 (formula 8).

Crossrefs

Similar sequences listed in A294116.

Programs

  • Magma
    a0:=2; a1:=8; [GeneralizedFibonacciNumber(a0, a1, n): n in [0..40]];
    
  • Maple
    f:= gfun:-rectoproc({a(n)=a(n-1)+a(n-2),a(0)=2,a(1)=8},a(n),remember):
    map(f, [$0..100]); # Robert Israel, Oct 24 2017
  • Mathematica
    LinearRecurrence[{1, 1}, {2, 8}, 40]
  • PARI
    Vec(2*(1 + 3*x)/(1 - x - x^2) + O(x^40)) \\ Colin Barker, Oct 25 2017
    
  • Sage
    a = BinaryRecurrenceSequence(1, 1, 2, 8)
    print([a(n) for n in range(38)]) # Peter Luschny, Oct 25 2017

Formula

G.f.: 2*(1 + 3*x)/(1 - x - x^2).
a(n) = a(n-1) + a(n-2).
a(n) = 2*A000285(n).
Let g(r,s;n) be the n-th generalized Fibonacci number with initial values r, s. We have:
a(n) = Lucas(n) + g(0,7;n), see A022090;
a(n) = Fibonacci(n) + g(2,7;n), see A022113;
a(n) = 2*g(1,8;n) - g(0,8;n);
a(n) = g(1,k;n) + g(1,8-k;n) for all k in Z.
a(h+k) = a(h)*Fibonacci(k-1) + a(h+1)*Fibonacci(k) for all h, k in Z (see S. Vajda in References section). For h=0 and k=n:
a(n) = 2*Fibonacci(n-1) + 8*Fibonacci(n).
Sum_{j=0..n} a(j) = a(n+2) - 8.
a(n) = (2^(-n)*((1-sqrt(5))^n*(-7+sqrt(5)) + (1+sqrt(5))^n*(7+sqrt(5)))) / sqrt(5). - Colin Barker, Oct 25 2017

A168622 Triangle read by rows: T(n, k) = [x^k]( 7*(1+x)^n - 6*(1+x^n) ) with T(0, 0) = 1.

Original entry on oeis.org

1, 1, 1, 1, 14, 1, 1, 21, 21, 1, 1, 28, 42, 28, 1, 1, 35, 70, 70, 35, 1, 1, 42, 105, 140, 105, 42, 1, 1, 49, 147, 245, 245, 147, 49, 1, 1, 56, 196, 392, 490, 392, 196, 56, 1, 1, 63, 252, 588, 882, 882, 588, 252, 63, 1, 1, 70, 315, 840, 1470, 1764, 1470, 840, 315, 70, 1
Offset: 0

Author

Roger L. Bagula and Gary W. Adamson, Dec 01 2009

Keywords

Examples

			Triangle begins as:
  1;
  1,  1;
  1, 14,   1;
  1, 21,  21,   1;
  1, 28,  42,  28,    1;
  1, 35,  70,  70,   35,    1;
  1, 42, 105, 140,  105,   42,    1;
  1, 49, 147, 245,  245,  147,   49,   1;
  1, 56, 196, 392,  490,  392,  196,  56,   1;
  1, 63, 252, 588,  882,  882,  588, 252,  63,  1;
  1, 70, 315, 840, 1470, 1764, 1470, 840, 315, 70, 1;
		

Crossrefs

Columns (essentially): A008589 (k=1), A024966 (k=2).

Programs

  • Magma
    A168622:= func< n,k | k eq 0 or k eq n select 1 else 7*Binomial(n,k) >;
    [A168622(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Apr 10 2025
    
  • Mathematica
    (* First program *)
    p[x_, n_]:= With[{m=3}, If[n==0, 1, (2*m+1)(1+x)^n - 2*m*(1+x^n)]];
    Table[CoefficientList[p[x,n], x], {n,0,12}]//Flatten
    (* Second program *)
    A168622[n_, k_]:= If[k==0 || k==n, 1, 7*Binomial[n,k]];
    Table[A168622[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Apr 10 2025 *)
  • SageMath
    def A168622(n,k):
        if k==0 or k==n: return 1
        else: return 7*binomial(n,k)
    print(flatten([[A168622(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Apr 10 2025

Formula

From G. C. Greubel, Apr 10 2025: (Start)
T(n, k) = 7*binomial(n, k), with T(n, 0) = T(n, n) = 1.
T(n, n-k) = T(n, k).
Sum_{k=0..n} T(n, k) = 2*A048489(n-1) + 6*[n=0].
Sum_{k=0..n} (-1)^k*T(n, k) = -6*(1 + (-1)^n) + 13*[n=0].
Sum_{k=0..floor(n/2)} T(n-k, k) = A022090(n+1) - 3*(3 + (-1)^n) + 6*[n=0].
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = (14/sqrt(3))*(-1)^n*cos((4*n+1)*Pi/6) - 6*(1 + (-1)^n*cos(n*Pi/2)) + 6*[n=0]. (End)
Showing 1-6 of 6 results.