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-2 of 2 results.

A213579 Rectangular array: (row n) = b**c, where b(h) = F(h), c(h) = n-1+h, where F=A000045 (Fibonacci numbers), n>=1, h>=1, and ** = convolution.

Original entry on oeis.org

1, 3, 2, 7, 5, 3, 14, 11, 7, 4, 26, 21, 15, 9, 5, 46, 38, 28, 19, 11, 6, 79, 66, 50, 35, 23, 13, 7, 133, 112, 86, 62, 42, 27, 15, 8, 221, 187, 145, 106, 74, 49, 31, 17, 9, 364, 309, 241, 178, 126, 86, 56, 35, 19, 10, 596, 507, 397, 295, 211, 146, 98, 63, 39, 21
Offset: 1

Views

Author

Clark Kimberling, Jun 18 2012

Keywords

Comments

Principal diagonal: A213580.
Antidiagonal sums: A053808.
Row 1, (1,1,2,3,5,...)**(1,2,3,4,...): A001924.
Row 2, (1,1,2,3,5,...)**(2,3,4,5,...): A023548.
Row 3, (1,1,2,3,5,...)**(3,4,5,6,...): A023552.
Row 4, (1,1,2,3,5,...)**(4,5,6,7,...): A210730.
Row 5, (1,1,2,3,5,...)**(5,6,7,8,...): A210731.
For a guide to related arrays, see A213500.

Examples

			Northwest corner (the array is read by falling antidiagonals):
1....3....7....14...26...46
2....5....11...21...38...66
3....7....15...28...50...86
4....9....19...35...62...106
5....11...23...42...74...126
6....13...27...49...86...146

Links

Crossrefs

Cf. A000045, A213500.

Programs

  • GAP
    Flat(List([1..12], n-> List([1..n], k-> Fibonacci(k+3) + n*Fibonacci(k+2) -(n+k+2) ))); # G. C. Greubel, Jul 08 2019
  • Magma
    [[Fibonacci(k+3) + n*Fibonacci(k+2) -(n+k+2): k in [1..n]]: n in [1..12]]; // G. C. Greubel, Jul 08 2019
  • Mathematica
    (* First program *)
b[n_]:= Fibonacci[n]; c[n_]:= n;
T[n_, k_]:= Sum[b[k-i] c[n+i], {i, 0, k-1}]
TableForm[Table[T[n, k], {n, 1, 10}, {k, 1, 10}]]
Flatten[Table[T[n-k+1, k], {n, 12}, {k, n, 1, -1}]] (* A213579 *)
r[n_]:= Table[T[n, k], {k, 40}]
d = Table[T[n, n], {n, 1, 40}] (* A213580 *)
s[n_]:= Sum[T[i, n+1-i], {i, 1, n}]
s1 = Table[s[n], {n, 1, 50}] (* A053808 *)
(* Second program *)
Table[Fibonacci[n-k+4] +k*Fibonacci[n-k+3] -(n+3), {n, 12}, {k, n}]//Flatten (* G. C. Greubel, Jul 08 2019 *)
  • PARI
    t(n,k) = fibonacci(n-k+4) + k*fibonacci(n-k+3) - (n+3);
for(n=1,12, for(k=1,n, print1(t(n,k), ", "))) \\ G. C. Greubel, Jul 08 2019
  • Sage
    [[fibonacci(k+3) + n*fibonacci(k+2) -(n+k+2) for k in (1..n)] for n in (1..12)] # G. C. Greubel, Jul 08 2019

Formula

T(n,k) = 3*T(n,k-1) - 2*T(n,k-2) - T(n,k-3) + T(n,k-4).
G.f. for row n: f(x)/g(x), where f(x) = n - (n-1)*x and g(x) = (1-x-x^2) *(1-x)^2.
T(n, k) = Fibonacci(k+3) + n*Fibonacci(k+2) - (n+k+2). - G. C. Greubel, Jul 08 2019

A210730 a(n) = a(n-1) + a(n-2) + n + 2 with n>1, a(0)=a(1)=0.

Original entry on oeis.org

0, 0, 4, 9, 19, 35, 62, 106, 178, 295, 485, 793, 1292, 2100, 3408, 5525, 8951, 14495, 23466, 37982, 61470, 99475, 160969, 260469, 421464, 681960, 1103452, 1785441, 2888923, 4674395, 7563350, 12237778, 19801162, 32038975, 51840173, 83879185, 135719396
Offset: 0

Views

Author

Alex Ratushnyak, May 10 2012

Keywords

Comments

Deleting the 0's leaves row 4 of the convolution array A213579. - Clark Kimberling, Jun 20 2012

Links

Crossrefs

Cf. A033818: a(n)=a(n-1)+a(n-2)+n-5, a(0)=a(1)=0 (except first 2 terms and sign).
Cf. A002062: a(n)=a(n-1)+a(n-2)+n-4, a(0)=a(1)=0 (except the first term and sign).
Cf. A065220: a(n)=a(n-1)+a(n-2)+n-3, a(0)=a(1)=0.
Cf. A001924: a(n)=a(n-1)+a(n-2)+n-1, a(0)=a(1)=0 (except the first term).
Cf. A023548: a(n)=a(n-1)+a(n-2)+n, a(0)=a(1)=0 (except first 2 terms).
Cf. A023552: a(n)=a(n-1)+a(n-2)+n+1, a(0)=a(1)=0 (except first 2 terms).
Cf. A210731: a(n)=a(n-1)+a(n-2)+n+3, a(0)=a(1)=0.
Cf. A000045, A213579.

Programs

  • GAP
    F:=Fibonacci;; List([0..40], n-> F(n+3)+3*F(n+1)-n-5); # G. C. Greubel, Jul 08 2019
  • Magma
    I:=[0, 0, 4, 9]; [n le 4 select I[n] else 3*Self(n-1)-2*Self(n-2)-Self(n-3)+Self(n-4): n in [1..37]]; // Bruno Berselli, May 10 2012
  • Magma
    F:=Fibonacci; [F(n+3)+3*F(n+1)-n-5: n in [0..40]]; // G. C. Greubel, Jul 08 2019
  • Mathematica
    RecurrenceTable[{a[0]==a[1]==0, a[n]==a[n-1] +a[n-2] +n+2}, a, {n, 40}] (* Bruno Berselli, May 10 2012 *)
LinearRecurrence[{3,-2,-1,1},{0,0,4,9},40] (* Harvey P. Dale, Jul 24 2013 *)
With[{F=Fibonacci}, Table[F[n+3]+2*F[n+1]-n-5, {n, 40}]] (* G. C. Greubel, Jul 08 2019 *)
  • PARI
    concat(vector(2), Vec(x^2*(4-3*x)/((1-x)^2*(1-x-x^2)) + O(x^50))) \\ Colin Barker, Mar 11 2017
  • PARI
    vector(40, n, n--; f=fibonacci; f(n+3)+3*f(n+1)-n-5) \\ G. C. Greubel, Jul 08 2019
  • Sage
    f=fibonacci; [f(n+3)+3*f(n+1)-n-5 for n in (0..40)] # G. C. Greubel, Jul 08 2019

Formula

G.f.: x^2*(4-3*x)/((1-x)^2*(1-x-x^2)). - Bruno Berselli, May 10 2012
a(n) = A210677(n)-1. - Bruno Berselli, May 10 2012
a(0)=0, a(1)=0, a(2)=4, a(3)=9, a(n) = 3*a(n-1)-2*a(n-2)-a(n-3)+a(n-4). - Harvey P. Dale, Jul 24 2013
a(n) = -5 + (2^(-1-n)*((1-sqrt(5))^n*(-7+5*sqrt(5)) + (1+sqrt(5))^n*(7+5*sqrt(5)))) / sqrt(5) - n. - Colin Barker, Mar 11 2017
a(n) = Fibonacci(n+3) + 3*Fibonacci(n+1) - n - 5. - G. C. Greubel, Jul 08 2019
Showing 1-2 of 2 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.