A152163 -id:A152163 - cp's OEIS frontend

A029907 a(n+1) = a(n) + a(n-1) + Fibonacci(n), with a(0) = 0 and a(1) = 1.

0, 1, 2, 4, 8, 15, 28, 51, 92, 164, 290, 509, 888, 1541, 2662, 4580, 7852, 13419, 22868, 38871, 65920, 111556, 188422, 317689, 534768, 898825, 1508618, 2528836, 4233872, 7080519, 11828620, 19741179, 32916068, 54835556, 91276202, 151814645, 252318312

Offset: 0

N. J. A. Sloane

Number of matchings of the fan graph on n vertices, n>0 (a fan is the join of the path graph with one extra vertex).
a(n+1) gives row sums of A054450. - Paul Barry, Oct 23 2004
Number of parts in all compositions of n into odd parts. Example: a(5)=15 because the compositions 5, 311, 131, 113, and 11111 have a total of 1+3+3+3+5=15 parts.
a(n-1) is the number of compositions of n that contain one even part; for example, a(5-1)=a(4)=8 counts the compositions 1112, 1121, 1211, 14, 2111, 23, 32, 41. - Joerg Arndt, May 21 2013

			a(4)=8 because matchings of fan graph with edges {OA,OB,OC,AB,AC} are: {},{OA},{OB},{OC},{AB},{AC},{OA,BC},{OC,AB}.

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Jean-Luc Baril and Nathanaël Hassler, Intervals in a family of Fibonacci lattices, Univ. de Bourgogne (France, 2024). See p. 5.
David Broadhurst, Zig-zag resistance and OEIS sequence A029907, Apr 11 2024
Jia Huang, Compositions with restricted parts, arXiv:1812.11010 [math.CO], 2018.
Mengmeng Liu and Andrew Yezhou Wang, The Number of Designated Parts in Compositions with Restricted Parts, J. Int. Seq., Vol. 23 (2020), Article 20.1.8.
Index entries for linear recurrences with constant coefficients, signature (2,1,-2,-1).

Cf. A000032, A000045, A001629, A010049, A054450, A152163, A240847.

a029907 n = a029907_list !! n
a029907_list = 0 : 1 : zipWith (+) (tail a000045_list)
                      (zipWith (+) (tail a029907_list) a029907_list)
-- Reinhard Zumkeller, Nov 01 2013

[((n+4)*Fibonacci(n)+2*n*Fibonacci(n-1))/5: n in [0..40]]; // Vincenzo Librandi, Feb 25 2018

with(combinat); A029907 := proc(n) options remember; if n <= 1 then n else procname(n-1)+procname(n-2)+fibonacci(n-1); fi; end;

CoefficientList[Series[x(1-x^2)/(1-x-x^2)^2, {x, 0, 37}], x] (* or *)
a[n_]:= a[n]= a[n-1] +a[n-2] +Fibonacci[n-1]; a[0]=0; a[1]=1; Array[a, 37] (* or *)
LinearRecurrence[{2,1,-2,-1}, {0,1,2,4}, 38] (* Robert G. Wilson v, Jun 22 2014 *)

alias(F,fibonacci); a(n)=((n+4)*F(n)+2*n*F(n-1))/5;

def A029907(n): return (1/5)*(n*lucas_number2(n, 1, -1) + 4*fibonacci(n))
[A029907(n) for n in (0..40)] # G. C. Greubel, Apr 06 2022

G.f.: x*(1-x^2)/(1-x-x^2)^2.
a(n) = ((n+4)*Fibonacci(n) + 2*n*Fibonacci(n-1))/5.
a(n+1) = Sum_{k=0..n} Sum_{j=0..floor(k/2)} binomial(n-j, j). - Paul Barry, Oct 23 2004
a(n) = A010049(n+1) + A152163(n+1). - R. J. Mathar, Dec 10 2011
a(n) = F(n) + Sum_{k=1..n-1} F(k)*F(n-k), where F=Fibonacci. - Reinhard Zumkeller, Nov 01 2013
a(n) = (1/5)*(n*A000032(n) + 4*A000045(n)). - G. C. Greubel, Apr 06 2022
a(n) = A001629(n+1) - A001629(n-1), where A001629 is the first convolution of the Fibonacci numbers. - Gregory L. Simay, Aug 30 2022
E.g.f.: exp(x/2)*(5*x*cosh(sqrt(5)*x/2) + sqrt(5)*(5*x + 8)*sinh(sqrt(5)*x/2))/25. - Stefano Spezia, Dec 04 2023

Additional formula from Wolfdieter Lang, May 02 2000

Additional comments from Michael Somos, Jul 23 2002

A147703 Triangle [1,1,1,0,0,0,...] DELTA [1,0,0,0,...] with Deléham DELTA defined in A084938.

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 5, 9, 5, 1, 13, 27, 20, 7, 1, 34, 80, 73, 35, 9, 1, 89, 234, 252, 151, 54, 11, 1, 233, 677, 837, 597, 269, 77, 13, 1, 610, 1941, 2702, 2225, 1199, 435, 104, 15, 1, 1597, 5523, 8533, 7943, 4956, 2158, 657, 135, 17, 1

Offset: 0

Paul Barry, Nov 10 2008

Equal to A062110*A007318 when A062110 is regarded as a triangle read by rows.

			Triangle begins
   1;
   1,   1;
   2,   3,   1;
   5,   9,   5,   1;
  13,  27,  20,   7,  1;
  34,  80,  73,  35,  9,  1;
  89, 234, 252, 151, 54, 11, 1;

Indranil Ghosh, Rows 0..100, flattened

Row sums are A006012. Diagonal sums are A147704.

Cf. A062110, A007318, A206800, A001519, A321620.

# The function RiordanSquare is defined in A321620:
RiordanSquare(1 / (1 - x / (1 - x / (1 - x))), 10); # Peter Luschny, Jan 26 2020

nmax=9; Flatten[CoefficientList[Series[CoefficientList[Series[(1 - 2*x)/(1 - (3 + y)*x + (1 + y)*x^2), {x, 0, nmax}], x], {y, 0, nmax}], y]] (* Indranil Ghosh, Mar 11 2017 *)

Riordan array ((1-2x)/(1-3x+x^2), x(1-x)/(1-3x+x^2)).
Sum_{k=0..n} T(n,k)*x^k = A152239(n), A152223(n), A152185(n), A152174(n), A152167(n), A152166(n), A152163(n), A000007(n), A001519(n), A006012(n), A081704(n), A082761(n), A147837(n), A147838(n), A147839(n), A147840(n), A147841(n), for x = -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively. - Philippe Deléham, Dec 01 2008
G.f.: (1-2*x)/(1-(3+y)*x+(1+y)*x^2). - Philippe Deléham, Nov 26 2011
T(n,k) = 3*T(n-1,k) + T(n-1,k-1) - T(n-2,k) - T(n-2,k-1), for n > 1. - Philippe Deléham, Feb 12 2012
The Riordan square of the odd indexed Fibonacci numbers A001519. - Peter Luschny, Jan 26 2020

A103631 Triangle read by rows: T(n,k) = abs(qStirling2(n,k,q)) for q = -1, with 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 2, 1, 0, 1, 1, 3, 2, 1, 0, 1, 1, 4, 3, 3, 1, 0, 1, 1, 5, 4, 6, 3, 1, 0, 1, 1, 6, 5, 10, 6, 4, 1, 0, 1, 1, 7, 6, 15, 10, 10, 4, 1, 0, 1, 1, 8, 7, 21, 15, 20, 10, 5, 1, 0, 1, 1, 9, 8, 28, 21, 35, 20, 15, 5, 1, 0, 1, 1, 10, 9, 36, 28, 56, 35, 35, 15, 6, 1

Offset: 0

Paul Barry, Feb 11 2005

Previous name: An invertible triangle whose row sums are F(n+1).
Triangle inverse has general term (-1)^(n-k)*binomial(floor(n/2),n-k). Diagonal sums are A103632.
Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 08 2005
Row sums are Fibonacci numbers (A000045).
Another version of triangle in A065941. - Philippe Deléham, Jan 01 2009
From Johannes W. Meijer, Aug 11 2011: (Start)
The T(n,k) coefficients appear in appendix 2 of Parks's remarkable article "A new proof of the Routh-Hurwitz stability criterion using the second method of Liapunov" if we assume that the b(r) coefficients are all equal to 1; see the second Maple program.
The T(n,k) triangle is related to a linear (n+1)-th order differential equation with coefficients a(n,k), see triangle A194005.
Parks's triangle appears to be an appropriate name for the triangle given above. (End)

			From _Paul Barry_, Oct 02 2009: (Start)
Triangle begins:
  1,
  0, 1,
  0, 1, 1,
  0, 1, 1, 1,
  0, 1, 1, 2, 1,
  0, 1, 1, 3, 2,  1,
  0, 1, 1, 4, 3,  3,  1,
  0, 1, 1, 5, 4,  6,  3,  1,
  0, 1, 1, 6, 5, 10,  6,  4, 1,
  0, 1, 1, 7, 6, 15, 10, 10, 4, 1
Production matrix is:
  0, 1,
  0, 1, 1,
  0, 0, 0, 1,
  0, 0, 0, 1, 1,
  0, 0, 0, 0, 0, 1,
  0, 0, 0, 0, 0, 1, 1,
  0, 0, 0, 0, 0, 0, 0, 1,
  0, 0, 0, 0, 0, 0, 0, 1, 1,
  0, 0, 0, 0, 0, 0, 0, 0, 0, 1 (End)

Reinhard Zumkeller, Rows n = 0..150 of triangle, flattened
Henry W. Gould, A Variant of Pascal's Triangle, The Fibonacci Quarterly, Vol. 3, Nr. 4, Dec. 1965, pp. 257-271, with corrections.
P. C. Parks, A new proof of the Routh-Hurwitz stability criterion using the second method of Liapunov , Math. Proc. of the Cambridge Philosophical Society, Vol. 58, Issue 04 (1962) pp. 694-702.

Cf. A103633 (signed version).

a103631 n k = a103631_tabl !! n !! k
a103631_row n = a103631_tabl !! n
a103631_tabl = [1] : [0,1] : f [1] [0,1] where
   f xs ys = zs : f ys zs where
     zs = zipWith (+)  ([0,0] ++ xs)  (ys ++ [0])
-- Reinhard Zumkeller, May 07 2012

/* As triangle: */ [[Binomial(Floor((2*n-k-1)/2), n-k): k in [0..n]]: n in [0..15]]; // Vincenzo Librandi, Aug 28 2016

From Johannes W. Meijer, Aug 11 2011: (Start)
A103631 := proc(n,k): binomial(floor((2*n-k-1)/2),n-k) end: seq(seq(A103631(n,k), k=0..n), n=0..12);
nmax:=12: for n from 0 to nmax+1 do b(n):=1 od: A103631 := proc(n,k) option remember: local j: if k=0 and n=0 then b(1) elif k=0 and n>=1 then 0 elif k=1 then b(n+1) elif k=2 then b(1)*b(n+1) elif k>=3 then expand(b(n+1)*add(procname(j,k-2), j=k-2..n-2)) fi: end: for n from 0 to nmax do seq(A103631(n,k), k=0..n) od: seq(seq(A103631(n,k),k=0..n), n=0..nmax); # (End)

p[x, -1] = 0; p[x, 0] = 1; p[x, 1] = x; p[x, 2] = x + x^2; p[x_, n_] := p[x, n] = p[x, n - 1] + x^2*p[x, n - 2]; (* with *) Table[ExpandAll[p[x, n]], {n, 0, 10}]; (* or *) a = Table[CoefficientList[p[x, n], x], {n, 0, 10}]; Flatten[a] (* Roger L. Bagula, Apr 27 2008 *)
Table[Binomial[Floor[(2*n - k - 1)/2], n - k], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, Aug 27 2016 *)
qStirling2[n_, k_, q_] /; 1 <= k <= n := q^(k - 1) qStirling2[n - 1, k - 1, q] + Sum[q^j, {j, 0, k - 1}] qStirling2[n - 1, k, q];
qStirling2[n_, 0, _] := KroneckerDelta[n, 0];
qStirling2[0, k_, _] := KroneckerDelta[0, k];
qStirling2[, , _] = 0;
Table[Abs[qStirling2[n, k, -1]], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 10 2020 *)

from sage.combinat.q_analogues import q_stirling_number2
for n in (0..9):
    print([abs(q_stirling_number2(n,k).substitute(q=-1)) for k in [0..n]])
# Peter Luschny, Mar 09 2020

T(n,k) = binomial(floor((2*n-k-1)/2), n-k).
A polynomial recursion which produces this triangle: p(x, n) = p(x, n - 1) + x^2*p(x, n - 2). - Roger L. Bagula, Apr 27 2008
Sum_{k=0..n} T(n,k)*x^k = A152163(n), A000007(n), A000045(n+1), A026597(n), A122994(n+1), A158608(n), A122995(n+1), A158797(n), A122996(n+1), A158798(n), A158609(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Jun 12 2009
G.f.: (1+(y-1)*x)/(1-x-y^2*x^2). - Philippe Deléham, Mar 09 2012
T(n,k) = T(n-1,k) + T(n-2,k-2), T(0,0) = 1, T(1,0) = 0, T(1,1) = 1, T(n,k) = 0 if k < 0 or if k > n. - Philippe Deléham, Mar 09 2012

New name from Peter Luschny, Mar 09 2020

A062110 A(n,k) is the coefficient of x^k in (1-x)^n/(1-2*x)^n for n, k >= 0; Table A read by descending antidiagonals.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 4, 5, 3, 1, 0, 8, 12, 9, 4, 1, 0, 16, 28, 25, 14, 5, 1, 0, 32, 64, 66, 44, 20, 6, 1, 0, 64, 144, 168, 129, 70, 27, 7, 1, 0, 128, 320, 416, 360, 225, 104, 35, 8, 1, 0, 256, 704, 1008, 968, 681, 363, 147, 44, 9, 1, 0, 512, 1536, 2400, 2528, 1970

Offset: 0

Henry Bottomley, May 30 2001

The triangular version of this square array is defined by T(n,k) = A(k,n-k) for 0 <= k <= n. Conversely, A(n,k) = T(n+k,n) for n,k >= 0. We have [o.g.f of T](x,y) = [o.g.f. of A](x*y, x) and [o.g.f. of A](x,y) = [o.g.f. of T](y,x/y). - Petros Hadjicostas, Feb 11 2021
From Paul Barry, Nov 10 2008: (Start)
As number triangle, Riordan array (1, x(1-x)/(1-2x)). A062110*A007318 is A147703.
[0,1,1,0,0,0,....] DELTA [1,0,0,0,.....]. (Philippe Deléham's DELTA is defined in A084938.) (End)
Modulo 2, this triangle T becomes triangle A106344. - Philippe Deléham, Dec 18 2008

			Table A(n,k) (with rows n >= 0 and columns k >= 0) begins:
  1, 0,  0,   0,   0,    0,    0,     0,     0,     0, ...
  1, 1,  2,   4,   8,   16,   32,    64,   128,   256, ...
  1, 2,  5,  12,  28,   64,  144,   320,   704,  1536, ...
  1, 3,  9,  25,  66,  168,  416,  1008,  2400,  5632, ...
  1, 4, 14,  44, 129,  360,  968,  2528,  6448, 16128, ...
  1, 5, 20,  70, 225,  681, 1970,  5500, 14920, 39520, ...
  1, 6, 27, 104, 363, 1182, 3653, 10836, 31092, 86784, ...
  ... - _Petros Hadjicostas_, Feb 15 2021
Triangle T(n,k) (with rows n >= 0 and columns k = 0..n) begins:
  1;
  0,   1;
  0,   1,   1;
  0,   2,   2,   1;
  0,   4,   5,   3,   1;
  0,   8,  12,   9,   4,   1;
  0,  16,  28,  25,  14,   5,   1;
  0,  32,  64,  66,  44,  20,   6,   1;
  0,  64, 144, 168, 129,  70,  27,   7,   1;
  0, 128, 320, 416, 360, 225, 104,  35,   8,   1;
  ... - _Philippe Deléham_, Nov 30 2008

Michael De Vlieger, Table of n, a(n) for n = 0..11475 (rows 0 <= n <= 150, flattened.)
Eunice Y. S. Chan, Robert M. Corless, Laureano Gonzalez-Vega, J. Rafael Sendra, and Juana Sendra, Upper Hessenberg and Toeplitz Bohemians, arXiv:1907.10677 [cs.SC], 2019.
Milan Janjić, Words and Linear Recurrences, J. Int. Seq., 21 (2018), #18.1.4.

Rows of A include A000007, A011782, A045623, A058396, A062109, A169792, A169793, A169794, A169795, A169796, A169797.
Columns of A include A000012, A001477, A000096, A000297.
Main diagonal of A is A002002.
Table A(n, k) is a multiple of 2^(k-n); dividing by this gives a table similar to A050143 except at the edges.
Essentially the same array as A105306, A160232.

t[n_, n_] = 1; t[n_, k_] := 2^(n-2*k)*k*Hypergeometric2F1[1-k, n-k+1, 2, -1]; Table[t[n, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Oct 30 2013, after Philippe Deléham + symbolic sum *)

a(i,j)=if(i<0 || j<0,0,polcoeff(((1-x)/(1-2*x)+x*O(x^j))^i,j))

Formulas for the square array (A(n,k): n,k >= 0):
A(n, k) = A(n-1, k) + Sum_{0 <= j < k} A(n, j) for n >= 1 and k >= 0 with A(0, k) = 0^k for k >= 0.
G.f.: 1/(1-x*(1-y)/(1-2*y)) = Sum_{i, j >= 0} A(i, j) x^i*y^j.
From Petros Hadjicostas, Feb 15 2021: (Start)
A(n,k) = 2^(k-n)*n*hypergeom([1-n, k+1], [2], -1) for n >= 0 and k >= 1.
A(n,k) = 2*A(n,k-1) + A(n-1,k) - A(n-1,k-1) for n,k >= 1 with A(n,0) = 1 for n >= 0 and A(0,k) = 0 for k >= 1. (End)
Formulas for the triangle (T(n,k): 0 <= k <= n):
From Philippe Deléham, Aug 01 2006: (Start)
T(n,k) = A121462(n+1,k+1)*2^(n-2*k) for 0 <= k < n.
T(n,k) = 2^(n-2*k)*k*hypergeom([1-k, n-k+1], [2], -1) for 0 <= k < n. (End)
Sum_{k=0..n} T(n,k)*x^k = A152239(n), A152223(n), A152185(n), A152174(n), A152167(n), A152166(n), A152163(n), A000007(n), A001519(n), A006012(n), A081704(n), A082761(n), A147837(n), A147838(n), A147839(n), A147840(n), A147841(n), for x = -7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Dec 09 2008
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k-1) for 1 <= k <= n-1 with T(0,0) = T(1,1) = T(2,1) = T(2,2) = 1, T(1,0) = T(2,0) = 0, and T(n,k) = 0 if k > n or if k < 0. - Philippe Deléham, Oct 30 2013
G.f.: Sum_{n.k>=0} T(n,k)*x^n*y^k = (1 - 2*x)/(x^2*y - x*y - 2*x + 1). - Petros Hadjicostas, Feb 15 2021

Various sections edited by Petros Hadjicostas, Feb 15 2021

A152185 a(n) = -3a(n-1) + 5a(n-2), n > 1; a(0)=1, a(1)=-5.

Original entry on oeis.org

1, -5, 20, -85, 355, -1490, 6245, -26185, 109780, -460265, 1929695, -8090410, 33919705, -142211165, 596232020, -2499751885, 10480415755, -43940006690, 184222098845, -772366329985, 3238209484180, -13576460102465

Offset: 0

Philippe Deléham, Nov 28 2008

Paolo Xausa, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (-3,5).

Cf. A152163, A152166, A152167, A152174.

LinearRecurrence[{-3, 5}, {1, -5}, 25] (* Paolo Xausa, Jan 19 2024 *)

Vec((1-2*x)/(1+3*x-5*x^2)+O(x^99)) \\ Charles R Greathouse IV, Jan 11 2012

G.f.: (1-2x)/(1+3x-5x^2).
a(n) = Sum_{k=0..n} A147703(n,k)*(-6)^k.
a(n) = (-1)^n*A152187(n). - Philippe Deléham, Nov 29 2008

A153764 Triangle T(n,k), 0 <= k <= n, read by rows, given by [1,0,-1,0,0,0,0,0,0,0,0,...] DELTA [0,1,0,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 1, 0, 1, 2, 3, 1, 1, 0, 1, 3, 3, 4, 1, 1, 0, 1, 3, 6, 4, 5, 1, 1, 0, 1, 4, 6, 10, 5, 6, 1, 1, 0, 1, 4, 10, 10, 15, 6, 7, 1, 1, 0, 1, 5, 10, 20, 15, 21, 7, 8, 1, 1, 0, 1, 5, 15, 20, 35, 21, 28, 8, 9, 1, 1, 0, 1, 6, 15, 35, 35, 56, 28, 36, 9, 10, 1, 1, 0

Offset: 0

Philippe Deléham, Jan 01 2009

A130595*A153342 as infinite lower triangular matrices. Reflected version of A103631. Another version of A046854. Row sums are Fibonacci numbers (A000045).

A055830*A130595 as infinite lower triangular matrices.

			Triangle begins:
  1;
  1, 0;
  1, 1, 0;
  1, 1, 1, 0;
  1, 2, 1, 1, 0;
  1, 2, 3, 1, 1, 0;
  1, 3, 3, 4, 1, 1, 0;
  ...

G. C. Greubel, Table of n, a(n) for the first 45 rows

Cf. A046854, A103631.

/* As triangle */ [[Binomial(Floor((n+k-1)/2),k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Aug 28 2016

Table[Binomial[Floor[(n + k - 1)/2], k], {n, 0, 45}, {k, 0, n}] // Flatten (* G. C. Greubel, Aug 27 2016 *)

T(n,k) = binomial(floor((n+k-1)/2),k).
Sum_{k=0..n} T(n,k)*x^k = A122335(n-1), A039834(n-2), A000012(n), A000045(n+1), A001333(n), A003688(n), A015448(n), A015449(n), A015451(n), A015453(n), A015454(n), A015455(n), A015456(n), A015457(n) for x = -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 respectively. - Philippe Deléham, Dec 17 2011
Sum_{k=0..n} T(n,k)*x^(n-k) = A152163(n), A000007(n), A000045(n+1), A026597(n), A122994(n+1), A158608(n), A122995(n+1), A158797(n), A122996(n+1), A158798(n), A158609(n) for x = -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Dec 17 2011
G.f.: (1+(1-y)*x)/(1-y*x-x^2). - Philippe Deléham, Dec 17 2011
T(n,k) = T(n-1,k-1) + T(n-2,k), T(0,0) = T(1,0) = T(2,0) = T(2,1) = 1, T(1,1) = T(2,2) = 0, T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Nov 09 2013

A225799 a(n) = Sum_{k=0..n} binomial(n,k) * 10^(n-k) * Fibonacci(n+k).

Original entry on oeis.org

0, 11, 143, 3058, 55341, 1052755, 19717984, 371084087, 6973353387, 131101759514, 2464418392865, 46327530894271, 870879506447808, 16371134451297043, 307750614069672631, 5785211638097121890, 108752568228856901349, 2044371455527726003547, 38430858858805840293152

Offset: 0

John Molokach, Jul 27 2013

This sequence is part of a family of Fibonacci-like sequences, where:
Sum_{k=0..n} binomial(n,k)*m^(n-k)*Fibonacci(n+k) produces a sequence whose terms are divisible by (m+1); m>=1.
A recurrence relation for a(n) (m not equal to zero) is:
a(n) = (m+3)*a(n-1) + (m^2+m-1)*a(n-2); a(0)=0, a(1)=m+1.
Notable values of m include:
  m = 1: Fibonacci(3n),
  m = 0: Fibonacci(2n) (using recurrence relation only - the sum above is undefined for m=0),
  m = -1: the zero sequence,
  m = -2: (-1)*Fibonacci(n), or A152163(n+2).
For any value of m, the sequence gives a(n*k) divisible by a(n); n>=1, k>=1, m not equal to -1 (zero is not divisible by zero).
Equivalent sequences are given by: Sum_{k=0..n} binomial(n,k) * (m+1)^k * Fibonacci(k).
When these sequences are divided by m+1, we obtain the family of sequences A057088, A015553, A087567, A087579, A087584, A087603, and so on.
Another interesting value of m, m = -3, gives a(2n-1)= -2 * 5^(n-1); a(2n)=0.

Index entries for linear recurrences with constant coefficients, signature (13,109).

Cf. A000045, A027941, A152163, A014445, A057088, A015553, A087567, A087579, A087584, A087603.

Table[Sum[Binomial[n, k]*10^(n - k)*Fibonacci[n + k], {k, 0, n}], {n, 0, 25}]
FullSimplify[Table[((13 + 11 Sqrt[5])^n - (13 - 11 Sqrt[5])^n)/(2^n Sqrt[5]), {n, 0, 25}]]
LinearRecurrence[{13,109},{0,11},30] (* Harvey P. Dale, Jul 31 2018 *)

a(n) = ((13 + 11*sqrt(5))^n - (13 - 11*sqrt(5))^n)/(2^n*sqrt(5)).
a(n) = 13*a(n-1) + 109*a(n-2); a(0)=0, a(1)=11.
G.f.: 11*x*/(1 - 13*x - 109*x^2). - Corrected by Georg Fischer, May 10 2019

cp's OEIS Frontend

A029907 a(n+1) = a(n) + a(n-1) + Fibonacci(n), with a(0) = 0 and a(1) = 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

SageMath

Formula

Extensions

A147703 Triangle [1,1,1,0,0,0,...] DELTA [1,0,0,0,...] with Deléham DELTA defined in A084938.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A103631 Triangle read by rows: T(n,k) = abs(qStirling2(n,k,q)) for q = -1, with 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Sage

Formula

Extensions

A062110 A(n,k) is the coefficient of x^k in (1-x)^n/(1-2*x)^n for n, k >= 0; Table A read by descending antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A152185 a(n) = -3*a(n-1) + 5*a(n-2), n > 1; a(0)=1, a(1)=-5.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A153764 Triangle T(n,k), 0 <= k <= n, read by rows, given by [1,0,-1,0,0,0,0,0,0,0,0,...] DELTA [0,1,0,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

A152185 a(n) = -3a(n-1) + 5a(n-2), n > 1; a(0)=1, a(1)=-5.