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.

Previous Showing 41-50 of 296 results. Next

A157725 a(n) = Fibonacci(n) + 2.

Original entry on oeis.org

2, 3, 3, 4, 5, 7, 10, 15, 23, 36, 57, 91, 146, 235, 379, 612, 989, 1599, 2586, 4183, 6767, 10948, 17713, 28659, 46370, 75027, 121395, 196420, 317813, 514231, 832042, 1346271, 2178311, 3524580, 5702889, 9227467, 14930354, 24157819, 39088171, 63245988, 102334157
Offset: 0

Views

Author

N. J. A. Sloane, Jun 26 2010

Keywords

Comments

a(n) = A226649(2*n+1) - A226649(2*n). - Reinhard Zumkeller, Jul 30 2013

Crossrefs

Programs

  • Haskell
    a157725 = (+ 2) . a000045
    a157725_list = 2 : 3 : map (subtract 2)
                           (zipWith (+) a157725_list $ tail a157725_list)
    -- Reinhard Zumkeller, Jul 30 2013
  • Magma
    [ Fibonacci(n) + 2: n in [0..40] ]; // Vincenzo Librandi, Apr 24 2011
    
  • Mathematica
    Fibonacci[Range[0, 50]] + 2 (* or *)
    LinearRecurrence[{2, 0, -1}, {2, 3, 3}, 50] (* Paolo Xausa, Jul 28 2024 *)
  • PARI
    a(n)=fibonacci(n)+2 \\ Charles R Greathouse IV, Jul 02 2013
    

Formula

G.f.: -(1+x)*(3*x-2) / ( (x-1)*(x^2+x-1) ). - R. J. Mathar, Aug 08 2012
a(0) = 2, a(1) = 3, a(n) = a(n - 2) + a(n - 1) - 2. - Reinhard Zumkeller, Jul 30 2013
E.g.f.: 2*(exp(x) + exp(x/2)*sinh(sqrt(5)*x/2)/sqrt(5)). - Stefano Spezia, Apr 09 2025

A056588 Coefficient triangle of certain polynomials.

Original entry on oeis.org

1, 1, -1, 1, -2, -1, 1, -4, -4, 1, 1, -7, -16, 7, 1, 1, -12, -53, 53, 12, -1, 1, -20, -166, 318, 166, -20, -1, 1, -33, -492, 1784, 1784, -492, -33, 1, 1, -54, -1413, 9288, 17840, -9288, -1413, 54, 1, 1, -88, -3960, 46233, 163504, -163504, -46233, 3960, 88, -1
Offset: 0

Views

Author

Wolfdieter Lang, Jul 10 2000

Keywords

Comments

G.f. for column m: see column sequences: A000012, A000071, A056589-91, for m=0..4.
The row polynomials p(n,x) := sum(a(n,m)*x^m) occur as numerators of the g.f. for the (n+1)-th power of Fibonacci numbers A000045. The corresponding denominator polynomials are the row polynomials q(n+2,x) = Sum_{m=0..n+2} A055870(n+2, m)*x^m (signed Fibonomial triangle).
The row polynomials p(n,x) and the companion denominator polynomials q(n,x) can be deduced from Riordan's recursion result.
The explicit formula is found from the recursion relation for powers of Fibonacci numbers (see Knuth's exercise with solution). - Roger L. Bagula, Apr 03 2010

Examples

			Row polynomial for n=4: p(4,x) = 1 - 7*x - 16*x^2 + 7*x^3 + x^4. x*p(4,x) is the numerator of the g.f. for A056572(n), n >= 0 (fifth power of Fibonacci numbers) {0,1,1,32,243,...}. The denominator polynomial is Sum_{m=0..6} A055870(6,m)*x^m (n=6 row polynomial of signed fibonomial triangle).
From _Roger L. Bagula_, Apr 03 2010: (Start)
1;
1,  -1;
1,  -2,    -1;
1,  -4,    -4,     1;
1,  -7,   -16,     7,      1;
1, -12,   -53,    53,     12,      -1;
1, -20,  -166,   318,    166,     -20,     -1;
1, -33,  -492,  1784,   1784,    -492,    -33,    1;
1, -54, -1413,  9288,  17840,   -9288,  -1413,   54,  1;
1, -88, -3960, 46233, 163504, -163504, -46233, 3960, 88, -1; (End)
		

References

  • D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, 1969, Vol. 1, p. 84, (exercise 1.2.8. Nr. 30) and p. 492 (solution).

Crossrefs

Programs

  • Maple
    A056588 := proc(n,k)
        if k = 0 then
            1;
        elif k >n then
            0;
        else
            combinat[fibonacci](k+1)^(n+1)+add( A055870(n+2, j)*(combinat[fibonacci](k+1-j)^(n+1)), j=1..k) ;
        end if;
    end proc: # R. J. Mathar, Jun 14 2015
  • Mathematica
    p[x_, n_] = Sum[(((1 + Sqrt[5])^k - (1 - Sqrt[5])^k)/(2^k*Sqrt[5]))^n*x^k, {k, 0, Infinity}];
    a = Table[CoefficientList[FullSimplify[Numerator[p[ x, n]]/x], x]/2^(1 + Floor[n/2]), {n, 1, 10}];
    Table[a[[n]]/a[[n]][[1]], {n, 1, 10}];
    Flatten[%] (* Roger L. Bagula, Apr 03 2010 *)
  • PARI
    S(n, k) = (-1)^floor((k+1)/2)*(prod(j=0, k-1, fibonacci(n-j))/prod(j=1, k, fibonacci(j)));
    T(n, k) = sum(j=0, k, fibonacci(k+1-j)^(n+1) * S(n+2, j));
    tabl(m) = for (n=0, m, for (k=0, n, print1(T(n, k), ", ")); print);
    tabl(9); \\ Tony Foster III, Aug 20 2018

Formula

a(n, m)=0 if nA000045(n) (Fibonacci) and sfibonomial(n, m) := A055870(n, m).
From Roger L. Bagula, Apr 03 2010: (Start)
p(x,n) = Sum_{k>=0} (((1 + sqrt(5))^k - (1 - sqrt(5))^k)/(2^k*sqrt(5)))^n*x^k;
t(n,m) = Numerator_coefficients(p(x,n)/x)/2^(1 + floor(n/2));
out(n,m) = t(n,m)/t(n,1). (End)
T(n, k) = Sum_{j=0..k} Fibonacci(k+1-j)^(n+1) * A055870(n+2, j). - Tony Foster III, Aug 20 2018
Sum_{j=0..n-1} a(n-1, n-1-j)*A010048(k+j, n) = Fibonacci(k)^n. - Tony Foster III, Jul 24 2018

A105809 Riordan array (1/(1 - x - x^2), x/(1 - x)).

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 3, 4, 3, 1, 5, 7, 7, 4, 1, 8, 12, 14, 11, 5, 1, 13, 20, 26, 25, 16, 6, 1, 21, 33, 46, 51, 41, 22, 7, 1, 34, 54, 79, 97, 92, 63, 29, 8, 1, 55, 88, 133, 176, 189, 155, 92, 37, 9, 1, 89, 143, 221, 309, 365, 344, 247, 129, 46, 10, 1, 144, 232, 364, 530, 674, 709, 591
Offset: 0

Views

Author

Paul Barry, May 04 2005

Keywords

Comments

Previous name was: A Fibonacci-Pascal matrix.
From Wolfdieter Lang, Oct 04 2014: (Start)
In the column k of this triangle (without leading zeros) is the k-fold iterated partial sums of the Fibonacci numbers, starting with 1. A000045(n+1), A000071(n+3), A001924(n+1), A014162(n+1), A014166(n+1), ..., n >= 0. See the Riordan property.
For a combinatorial interpretation of these iterated partial sums see the H. Belbachir and A. Belkhir link. There table 1 shows in the rows these columns. In their notation (with r = k) f^(k)(n) = T(k, n+k).
The A-sequence of this Riordan triangle is [1, 1] (see the recurrence for T(n, k), k >= 1, given in the formula section). The Z-sequence is A165326 = [1, repeat(1, -1)]. See the W. Lang link under A006232 for Riordan A- and Z-sequences. (End)

Examples

			The triangle T(n,k) begins:
n\k   0   1   2    3    4    5    6    7    8   9  10 11 12 13 ...
0:    1
1:    1   1
2:    2   2   1
3:    3   4   3    1
4:    5   7   7    4    1
5:    8  12  14   11    5    1
6:   13  20  26   25   16    6    1
7:   21  33  46   51   41   22    7    1
8:   34  54  79   97   92   63   29    8    1
9:   55  88 133  176  189  155   92   37    9   1
10:  89 143 221  309  365  344  247  129   46  10   1
11: 144 232 364  530  674  709  591  376  175  56  11  1
12: 233 376 596  894 1204 1383 1300  967  551 231  67 12  1
13: 377 609 972 1490 2098 2587 2683 2267 1518 782 298 79 13  1
... reformatted and extended - _Wolfdieter Lang_, Oct 03 2014
------------------------------------------------------------------
Recurrence from Z-sequence (see a comment above): 8 = T(0,5) = (+1)*5 + (+1)*7 + (-1)*7 + (+1)*4 + (-1)*1 = 8. - _Wolfdieter Lang_, Oct 04 2014
		

Crossrefs

Cf. A165326 (Z-sequence), A027934 (row sums), A010049(n+1) (antidiagonal sums), A212804 (alternating row sums), inverse is A105810.
Some other Fibonacci-Pascal triangles: A027926, A036355, A037027, A074829, A109906, A111006, A114197, A162741, A228074.

Programs

  • Haskell
    a105809 n k = a105809_tabl !! n !! k
    a105809_row n = a105809_tabl !! n
    a105809_tabl = map fst $ iterate
       (\(u:_, vs) -> (vs, zipWith (+) ([u] ++ vs) (vs ++ [0]))) ([1], [1,1])
    -- Reinhard Zumkeller, Aug 15 2013
  • Maple
    T := (n,k) -> `if`(n=0,1,binomial(n,k)*hypergeom([1,k/2-n/2,k/2-n/2+1/2], [k+1,-n], -4)); for n from 0 to 13 do seq(simplify(T(n,k)),k=0..n) od; # Peter Luschny, Oct 10 2014
  • Mathematica
    T[n_, k_] := Sum[Binomial[n-j, k+j], {j, 0, n}]; Table[T[n, k], {n, 0, 11}, {k, 0, n}] (* Jean-François Alcover, Jun 11 2019 *)

Formula

Riordan array (1/(1-x-x^2), x/(1-x)).
Triangle T(n, k) = Sum_{j=0..n} binomial(n-j, k+j); T(n, 0) = A000045(n+1);
T(n, m) = T(n-1, m-1) + T(n-1, m).
T(n, k) = Sum_{j=0..n} binomial(j, n+k-j). - Paul Barry, Oct 23 2006
G.f. of row polynomials Sum_{k=0..n} T(n, k)*x^k is (1 - z)/((1 - z - z^2)*(1 - (1 + x)*z)) (Riordan property). - Wolfdieter Lang, Oct 04 2014
T(n, k) = binomial(n, k)*hypergeom([1, k/2 - n/2, k/2 - n/2 + 1/2],[k + 1, -n], -4) for n > 0. - Peter Luschny, Oct 10 2014
From Wolfdieter Lang, Feb 13 2025: (Start)
Array A(k, n) = Sum_{j=0..n} F(j+1)*binomial(k-1+n-j, k-1), k >= 0, n >= 0, with F = A000045, (from Riordan triangle k-th convolution in columns without leading 0s).
A(k, n) = F(n+1+2*k) - Sum_{j=0..k-1} F(2*(k-j)-1) * binomial(n+1+j, j), (from iteration of partial sums).
Triangle T(n, k) = A(k, n-k) = Sum_{j=k..n} F(n-j+1) * binomial(j-1, k-1), 0 <= k <= n.
T(n, k) = F(n+1+k) - Sum_{j=0..k-1} F(2*(k-j)-1) * binomial(n - (k-1-j), j). (End)
T(n, k) = A027926(n, n+k), for 0 <= k <= n. - Wolfdieter Lang, Mar 08 2025

Extensions

Use first formula as a more descriptive name, Joerg Arndt, Jun 08 2021

A005968 Sum of cubes of first n Fibonacci numbers.

Original entry on oeis.org

0, 1, 2, 10, 37, 162, 674, 2871, 12132, 51436, 217811, 922780, 3908764, 16558101, 70140734, 297121734, 1258626537, 5331629710, 22585142414, 95672204155, 405273951280, 1716768021816, 7272346018247, 30806152127640, 130496954475672, 552793970116297, 2341672834801754
Offset: 0

Views

Author

Keywords

Comments

From Alexander Adamchuk, Aug 07 2006: (Start)
The only two prime terms are a(2) = 2 and a(4) = 37.
The prime p divides a(p-1) iff p is in A045468.
The prime p divides a((p-1)/2) iff p is in A047650.
3^4 divides a(p) iff p is in A003628.
3^5 divides a(p) for p = {37,53,109,181,197,269,397,431,541,...}.
3^6 divides a(p) for p = {109,541,...}.
3^7 divides a(p) for p = {557,...}. (End)

References

  • Art Benjamin, Timothy A. Carnes, and Benoit Cloitre, Recounting the Sums of Cubes of Fibonacci Numbers, Congressus Numerantium, Proceedings of the Eleventh International Conference on Fibonacci Numbers and their Applications, (William Webb, ed.), Vol 194, pp. 45-51, 2009.
  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, p. 14.
  • A. Brousseau, Fibonacci and Related Number Theoretic Tables. Fibonacci Association, San Jose, CA, 1972, p. 18.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Partial sums of A056570. Cf. A119284 (alternating sum).
Sums of other powers: A000071, A001654, A005969, A098531, A098532, A098533, A128697.

Programs

  • Magma
    [(1/10)*( Fibonacci(3*n+2)-(-1)^(n)*6*Fibonacci(n-1)+5 ): n in [0..30]]; // G. C. Greubel, Jan 17 2018
  • Maple
    with(combinat): l[0] := 0: for i from 1 to 50 do l[i] := l[i-1]+fibonacci(i)^3; printf(`%d,`,l[i]) od: # James Sellers, May 29 2000
    A005968:=(-1+2*z+z**2)/(z-1)/(z**2+4*z-1)/(z**2-z-1); # conjectured by Simon Plouffe in his 1992 dissertation
  • Mathematica
    f[n_]:=(Fibonacci[n]*Fibonacci[n+1]^2+(-1)^(n-1)*Fibonacci[n-1]+1)/2;Table[f[n],{n,0,5!}] (* Vladimir Joseph Stephan Orlovsky, Nov 22 2010 *)
    Accumulate[Fibonacci[Range[0,20]]^3]
    CoefficientList[Series[x*(1-2*x-x^2)/((1-x)*(1+x-x^2)*(1-4*x-x^2)), {x, 0, 50}], x] (* Vincenzo Librandi, Jun 09 2013 *)
  • PARI
    a(n)=(fibonacci(n)*fibonacci(n+1)^2+(-1)^(n-1)*fibonacci(n-1)+1)/2
    
  • PARI
    a(n)=(fibonacci(3*n+2)-(-1)^(n)*6*fibonacci(n-1)+5)/10
    
  • PARI
    a(n)=sum(i=1,n,fibonacci(i)^3)
    

Formula

G.f.: x*(1-2*x-x^2)/((1-x)*(1+x-x^2)*(1-4*x-x^2)). - Ralf Stephan, Apr 23 2004
a(n) = (1/2)*(F(n)*F(n+1)^2 + (-1)^(n-1)*F(n-1) + 1). - Benoit Cloitre, Aug 06 2004
a(n) = Sum_{i=1..n} A000045(i)^3.
a(n) = (1/10)*(F(3*n+2) - (-1)^(n)*6*F(n-1) + 5). - Art Benjamin and Timothy A. Carnes
a(n+5) = 4*a(n+4) + 3*a(n+3) - 9*a(n+2) + 2*a(n+1) + a(n). - Benoit Cloitre, Sep 12 2004

Extensions

More terms from James Sellers, May 29 2000

A032858 Numbers whose base-3 representation Sum_{i=0..m} d(i)*3^i has d(m) > d(m-1) < d(m-2) > ...

Original entry on oeis.org

0, 1, 2, 3, 6, 7, 10, 11, 19, 20, 23, 30, 33, 34, 57, 60, 61, 69, 70, 91, 92, 100, 101, 104, 172, 173, 181, 182, 185, 208, 209, 212, 273, 276, 277, 300, 303, 304, 312, 313, 516, 519, 520, 543, 546, 547, 555, 556, 624, 627, 628, 636, 637
Offset: 1

Views

Author

Keywords

Comments

Every other base-3 digit must be strictly less than its neighbors. - M. F. Hasler, Oct 05 2018
The terms can be generated in the following way: if A(n) are the terms with n digits in base 3, the terms with n+2 digits are obtained by prefixing them with '10' and with '20', and prefixing '21' to those starting with a digit '2'. It is easy to prove that #A(n) = A000045(n+2), since from the above we have #A(n+2) = 2*#A(n) + #A(n-1) = #A(n) + #A(n+1). (The #A(n-1) numbers starting with '2' are #A(n-2) numbers prefixed with '20' and #A(n-3) prefixed with '21'.) - M. F. Hasler, Oct 05 2018

Examples

			The base-3 representation of the initial terms is 0, 1, 2, 10, 20, 21, 101, 102, 201, 202, 212, 1010, 1020, 1021, 2010, 2020, 2021, 2120, 2121, 10101, 10102, ...
		

Crossrefs

Cf. A032859 .. A032865 for base-4 .. 10 variants.
Cf. A000975 (or A056830 in binary) for the base-2 analog.
Cf. A306105 for these terms written in base 3.

Programs

  • Mathematica
    sdQ[n_]:=Module[{s=Sign[Differences[IntegerDigits[n, 3]]]}, s==PadRight[{}, Length[s], {-1, 1}]]; Select[Range[0, 700], sdQ] (* Vincenzo Librandi, Oct 06 2018 *)
  • PARI
    is(n,b=3)=!for(i=2,#n=digits(n,b),(n[i-1]-n[i])*(-1)^i>0||return) \\ M. F. Hasler, Oct 05 2018

Formula

a(A000071(n+3)) = floor(3^(n+1)/8) = A033113(n). - M. F. Hasler, Oct 05 2018

Extensions

Definition edited, cross-references and a(1) = 0 inserted by M. F. Hasler, Oct 05 2018

A119282 Alternating sum of the first n Fibonacci numbers.

Original entry on oeis.org

0, -1, 0, -2, 1, -4, 4, -9, 12, -22, 33, -56, 88, -145, 232, -378, 609, -988, 1596, -2585, 4180, -6766, 10945, -17712, 28656, -46369, 75024, -121394, 196417, -317812, 514228, -832041, 1346268, -2178310, 3524577, -5702888, 9227464, -14930353, 24157816, -39088170, 63245985, -102334156, 165580140, -267914297, 433494436, -701408734, 1134903169, -1836311904, 2971215072, -4807526977, 7778742048
Offset: 0

Views

Author

Stuart Clary, May 13 2006

Keywords

Comments

Apart from signs, same as A008346.
Natural bilateral extension (brackets mark index 0): ..., 88, 54, 33, 20, 12, 7, 4, 2, 1, 0, [0], -1, 0, -2, 1, -4, 4, -9, 12, -22, 3, ... This is A000071-reversed followed by A119282.
Alternating sums of rows of the triangle in A141169. - Reinhard Zumkeller, Mar 22 2011

Crossrefs

Programs

  • Magma
    [0] cat [(&+[(-1)^k*Fibonacci(k):k in [1..n]]): n in [1..30]]; // G. C. Greubel, Jan 17 2018
  • Mathematica
    FoldList[#1 - Fibonacci@ #2 &, -Range@ 50] (* Michael De Vlieger, Jan 27 2016 *)
    Accumulate[Table[(-1)^n Fibonacci[n], {n, 0, 49}]] (* Alonso del Arte, Apr 25 2017 *)
  • PARI
    a(n) = sum(k=1, n, (-1)^k*fibonacci(k)); \\ Michel Marcus, Jan 27 2016
    

Formula

Let F(n) be the Fibonacci number A000045(n).
a(n) = Sum_{k = 1..n} (-1)^k F(k).
Closed form: a(n) = (-1)^n F(n-1) - 1 = (-1)^n A008346(n-1).
Recurrence: a(n) - 2 a(n-2) + a(n-3)= 0.
G.f.: A(x) = -x/(1 - 2 x^2 + x^3) = -x/((1 - x)(1 + x - x^2)).
Another recurrence: a(n) = a(n-2) - a(n-1) - 1. - Rick L. Shepherd, Aug 12 2009

A167616 a(n) = Fibonacci(n) - 5.

Original entry on oeis.org

0, 3, 8, 16, 29, 50, 84, 139, 228, 372, 605, 982, 1592, 2579, 4176, 6760, 10941, 17706, 28652, 46363, 75020, 121388, 196413, 317806, 514224, 832035, 1346264, 2178304, 3524573, 5702882, 9227460, 14930347, 24157812, 39088164, 63245981, 102334150
Offset: 5

Views

Author

N. J. A. Sloane, Jun 26 2010

Keywords

Crossrefs

Programs

Formula

a(n) = a(n-1) + a(n-2) + 5. - Zak Seidov, Jun 27 2010
a(5)=0, a(6)=3, a(7)=8, a(n) = 2*a(n-1) - a(n-3). - Harvey P. Dale, Dec 31 2011
G.f.: x^6*(3+2*x)/((1-x)*(1-x-x^2)). - L. Edson Jeffery, Mar 17 2013

A011794 Triangle defined by T(n+1, k) = T(n, k-1) + T(n-1, k), T(n,1) = 1, T(1,k) = 1, T(2,k) = min(2,k).

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 3, 4, 5, 1, 3, 6, 7, 8, 1, 4, 7, 11, 12, 13, 1, 4, 10, 14, 19, 20, 21, 1, 5, 11, 21, 26, 32, 33, 34, 1, 5, 15, 25, 40, 46, 53, 54, 55, 1, 6, 16, 36, 51, 72, 79, 87, 88, 89, 1, 6, 21, 41, 76, 97, 125, 133, 142, 143, 144, 1, 7, 22, 57, 92, 148, 176, 212, 221, 231, 232, 233
Offset: 1

Views

Author

Keywords

Examples

			matrix(10,10,n,k,a(n-1,k-1))
  [ 0 0 0 0 0 0 0 0 0 0 ]
  [ 0 1 1 1 1 1 1 1 1 1 ]
  [ 0 1 2 2 2 2 2 2 2 2 ]
  [ 0 1 2 3 3 3 3 3 3 3 ]
  [ 0 1 3 4 5 5 5 5 5 5 ]
  [ 0 1 3 6 7 8 8 8 8 8 ]
Triangle begins as:
  1;
  1, 2;
  1, 2,  3;
  1, 3,  4,  5;
  1, 3,  6,  7,  8;
  1, 4,  7, 11, 12, 13;
  1, 4, 10, 14, 19, 20, 21;
  1, 5, 11, 21, 26, 32, 33, 34;
  1, 5, 15, 25, 40, 46, 53, 54, 55;
  1, 6, 16, 36, 51, 72, 79, 87, 88, 89;
		

Crossrefs

Columns include A008619 and (essentially) A055802, A055803, A055804, A055805, A055806.
Essentially a reflected version of A055801.
Sums include: A039834 (signed row), A131913 (row).

Programs

  • Magma
    function T(n,k) // T = A011794(n,k)
      if k eq 1 or n eq 1 then return 1;
      elif n eq 2 then return Min(2, k);
      else return T(n-1,k-1) + T(n-2,k);
      end if;
    end function;
    [T(n,k): k in [1..n], n in [1..15]]; // G. C. Greubel, Oct 21 2024
    
  • Mathematica
    T[n_, k_]:= T[n, k]= T[n-1, k-1] + T[n-2, k]; T[n_, 1] = 1; T[1, k_] = 1; T[2, k_] := Min[2, k]; Table[T[n, k], {n,15}, {k,n}]//Flatten (* Jean-François Alcover, Feb 26 2013 *)
  • PARI
    T(n,k)=if(n<=0 || k<=0,0, if(n<=2 || k==1, min(n,k), T(n-1,k-1)+T(n-2,k)))
    
  • SageMath
    def T(n, k): # T = A011794
        if (k==1 or n==1): return 1
        elif (n==2): return min(2,k)
        else: return T(n-1, k-1) + T(n-2, k)
    flatten([[T(n, k) for k in range(1,n+1)] for n in range(1,16)]) # G. C. Greubel, Oct 21 2024

Formula

T(n,n) = Fibonacci(n+1). - Jean-François Alcover, Feb 26 2013
From G. C. Greubel, Oct 21 2024: (Start)
Sum_{k=1..n} T(n, k) = A131913(n-1).
Sum_{k=1..n} (-1)^(k-1)*T(n, k) = A039834(n).
Sum_{k=1..floor((n+1)/2)} T(n-k+1,k) = (1/2)*((1-(-1)^n)*A074878((n+3)/2) + (1+(-1)^n)*A008466((n+6)/2)) (diagonal row sums).
Sum_{k=1..floor((n+1)/2)} (-1)^(k-1)*T(n-k+1,k) = (-1)^floor((n-1)/2)*A103609(n) + [n=1] (signed diagonal row sums). (End)

Extensions

Entry improved by comments from Michael Somos
More terms added by G. C. Greubel, Oct 21 2024

A019274 Number of recursive calls needed to compute the n-th Fibonacci number F(n), starting with F(1) = F(2) = 1.

Original entry on oeis.org

0, 0, 2, 4, 8, 14, 24, 40, 66, 108, 176, 286, 464, 752, 1218, 1972, 3192, 5166, 8360, 13528, 21890, 35420, 57312, 92734, 150048, 242784, 392834, 635620, 1028456, 1664078, 2692536, 4356616, 7049154, 11405772, 18454928, 29860702, 48315632
Offset: 1

Views

Author

Kim Trammell (kim(AT)coc.com) and others

Keywords

Comments

Let g = F(2) + F(3) + ... + F(n) = F(n+2) - 2. Some numbers in the range [0,g] have unique representations of the form Sum_{i=1..n} a(i)*F(i) where each a(i) is 1 or -1. These numbers have the form g-k for k in the sequence. - Louis ten Bosch (louis_ten_bosch(AT)hotmail.com), Jan 01 2003

Crossrefs

Cf. A000045.
Antidiagonal sums of array A017125.

Programs

Formula

a(n) = a(n-1) + a(n-2) + 2. a(n) = 2*F(n) - 2 = 2*A000071(n).
a(n+2) = Sum_{k=0..n} Fibonacci(n-k) + k*Fibonacci(n-k).
a(n) = Sum_{k=0..n} (2 - 2*0^(n-k))*F(k). - Paul Barry, Oct 24 2007
a(n) = F(n) + F(n+3) - 2, n>=-1 (where F(n) is the n-th Fibonacci number). - Zerinvary Lajos, Jan 31 2008
G.f.: 2*x^3 / ( (x-1)*(x^2+x-1) ). - R. J. Mathar, Jul 01 2012
a(1)=0, a(2)=0, a(3)=2, a(n) = 2*a(n-1) - a(n-3). - Harvey P. Dale, Oct 16 2012

A053538 Triangle: a(n,m) = ways to place p balls in n slots with m in the rightmost p slots, 0<=p<=n, 0<=m<=n, summed over p, a(n,m)= Sum_{k=0..n} binomial(k,m)*binomial(n-k,k-m), (see program line).

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 5, 5, 4, 1, 1, 8, 10, 7, 5, 1, 1, 13, 18, 16, 9, 6, 1, 1, 21, 33, 31, 23, 11, 7, 1, 1, 34, 59, 62, 47, 31, 13, 8, 1, 1, 55, 105, 119, 101, 66, 40, 15, 9, 1, 1, 89, 185, 227, 205, 151, 88, 50, 17, 10, 1, 1, 144, 324, 426, 414, 321, 213, 113, 61, 19, 11, 1, 1
Offset: 0

Views

Author

Wouter Meeussen, May 23 2001

Keywords

Comments

Riordan array (1/(1-x-x^2), x(1-x)/(1-x-x^2)). Row sums are A000079. Diagonal sums are A006053(n+2). - Paul Barry, Nov 01 2006
Subtriangle of the triangle given by (0, 1, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 05 2012
Mirror image of triangle in A208342. - Philippe Deléham, Mar 05 2012
A053538 is jointly generated with A076791 as an array of coefficients of polynomials u(n,x): initially, u(1,x)=v(1,x)=1, for n>1, u(n,x) = x*u(n-1,x) + v(n-1,x) and v(n,x) = u(n-1,x) + v(n-1,x). See the Mathematica section at A076791. - Clark Kimberling, Mar 08 2012
The matrix inverse starts
1;
-1, 1;
-1, -1, 1;
1, -2, -1, 1;
3, 1, -3, -1, 1;
1, 6, 1, -4, -1, 1;
-7, 4, 10, 1, -5, -1, 1;
-13, -13, 8, 15, 1, -6, -1, 1;
3, -31, -23, 13, 21, 1, -7, -1, 1; - R. J. Mathar, Mar 15 2013
Also appears to be the number of subsets of {1..n} containing n with k maximal anti-runs of consecutive elements increasing by more than 1. For example, the subset {1,3,6,7,11,12} has maximal anti-runs ((1,3,6),(7,11),(12)) so is counted under a(12,3). For runs instead of anti-runs we get A202064. - Gus Wiseman, Jun 26 2025

Examples

			n=4; Table[binomial[k, j]binomial[n-k, k-j], {k, 0, n}, {j, 0, n}] splits {1, 4, 6, 4, 1} into {{1, 0, 0, 0, 0}, {3, 1, 0, 0, 0}, {1, 4, 1, 0, 0}, {0, 0, 3, 1, 0}, {0, 0, 0, 0, 1}} and this gives summed by columns {5, 5, 4, 1, 1}
Triangle begins :
   1;
   1,  1;
   2,  1,  1;
   3,  3,  1, 1;
   5,  5,  4, 1, 1;
   8, 10,  7, 5, 1, 1;
  13, 18, 16, 9, 6, 1, 1;
...
(0, 1, 1, -1, 0, 0, 0, ...) DELTA (1, 0, -1, 1, 0, 0, 0, ...) begins :
  1;
  0,  1;
  0,  1,  1;
  0,  2,  1,  1;
  0,  3,  3,  1, 1;
  0,  5,  5,  4, 1, 1;
  0,  8, 10,  7, 5, 1, 1;
  0, 13, 18, 16, 9, 6, 1, 1;
		

Crossrefs

Column k = 1 is A000045.
Row sums are A000079.
Column k = 2 is A010049.
For runs instead of anti-runs we have A202064.
For integer partitions see A268193, strict A384905, runs A116674.
A034839 counts subsets by number of maximal runs.
A384175 counts subsets with all distinct lengths of maximal runs, complement A384176.
A384877 gives lengths of maximal anti-runs in binary indices, firsts A384878.
A384893 counts subsets by number of maximal anti-runs.

Programs

  • GAP
    Flat(List([0..12], n-> List([0..n], k-> Sum([0..n], j->  Binomial(j,k)*Binomial(n-j,j-k)) ))); # G. C. Greubel, May 16 2019
  • Magma
    [[(&+[Binomial(j,k)*Binomial(n-j,j-k): j in [0..n]]): k in [0..n]]: n in [0..12]]; // G. C. Greubel, May 16 2019
    
  • Maple
    a:= (n, m)-> add(binomial(k, m)*binomial(n-k, k-m), k=0..n):
    seq(seq(a(n,m), m=0..n), n=0..12);  # Alois P. Heinz, Sep 19 2013
  • Mathematica
    Table[Sum[Binomial[k, m]*Binomial[n-k, k-m], {k,0,n}], {n,0,12}, {m,0,n}]
  • PARI
    {T(n,k) = sum(j=0,n, binomial(j,k)*binomial(n-j,j-k))}; \\ G. C. Greubel, May 16 2019
    
  • Sage
    [[sum(binomial(j,k)*binomial(n-j,j-k) for j in (0..n)) for k in (0..n)] for n in (0..12)] # G. C. Greubel, May 16 2019
    

Formula

From Philippe Deléham, Mar 05 2012: (Start)
T(n,k) = T(n-1,k) + T(n-1,k-1) + T(n-2,k) - T(n-2,k-1), T(0,0) = T(1,0) = T(1,1) = 1 and T(n,k) = 0 if k<0 or if k>n.
G.f.: 1/(1-(1+y)*x-(1-y)*x^2).
Sum_{k, 0<=k<=n} T(n,k)*x^k = A077957(n), A000045(n+1), A000079(n), A001906(n+1), A007070(n), A116415(n), A084326(n+1), A190974(n+1), A190978(n+1), A190984(n+1), A190990(n+1), A190872(n+1) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 respectively. (End)
Previous Showing 41-50 of 296 results. Next