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

A213500 Rectangular array T(n,k): (row n) = b**c, where b(h) = h, c(h) = h + n - 1, n >= 1, h >= 1, and ** = convolution.

Original entry on oeis.org

1, 4, 2, 10, 7, 3, 20, 16, 10, 4, 35, 30, 22, 13, 5, 56, 50, 40, 28, 16, 6, 84, 77, 65, 50, 34, 19, 7, 120, 112, 98, 80, 60, 40, 22, 8, 165, 156, 140, 119, 95, 70, 46, 25, 9, 220, 210, 192, 168, 140, 110, 80, 52, 28, 10, 286, 275, 255, 228, 196, 161, 125, 90
Offset: 1

Views

Author

Clark Kimberling, Jun 14 2012

Keywords

Comments

Principal diagonal: A002412.
Antidiagonal sums: A002415.
Row 1: (1,2,3,...)**(1,2,3,...) = A000292.
Row 2: (1,2,3,...)**(2,3,4,...) = A005581.
Row 3: (1,2,3,...)**(3,4,5,...) = A006503.
Row 4: (1,2,3,...)**(4,5,6,...) = A060488.
Row 5: (1,2,3,...)**(5,6,7,...) = A096941.
Row 6: (1,2,3,...)**(6,7,8,...) = A096957.
...
In general, the convolution of two infinite sequences is defined from the convolution of two n-tuples: let X(n) = (x(1),...,x(n)) and Y(n)=(y(1),...,y(n)); then X(n)**Y(n) = x(1)*y(n)+x(2)*y(n-1)+...+x(n)*y(1); this sum is the n-th term in the convolution of infinite sequences:(x(1),...,x(n),...)**(y(1),...,y(n),...), for all n>=1.
...
In the following guide to related arrays and sequences, row n of each array T(n,k) is the convolution b**c of the sequences b(h) and c(h+n-1). The principal diagonal is given by T(n,n) and the n-th antidiagonal sum by S(n). In some cases, T(n,n) or S(n) differs in offset from the listed sequence.
b(h)........ c(h)........ T(n,k) .. T(n,n) .. S(n)
h .......... h .......... A213500 . A002412 . A002415
h .......... h^2 ........ A212891 . A213436 . A024166
h^2 ........ h .......... A213503 . A117066 . A033455
h^2 ........ h^2 ........ A213505 . A213546 . A213547
h .......... h*(h+1)/2 .. A213548 . A213549 . A051836
h*(h+1)/2 .. h .......... A213550 . A002418 . A005585
h*(h+1)/2 .. h*(h+1)/2 .. A213551 . A213552 . A051923
h .......... h^3 ........ A213553 . A213554 . A101089
h^3 ........ h .......... A213555 . A213556 . A213547
h^3 ........ h^3 ........ A213558 . A213559 . A213560
h^2 ........ h*(h+1)/2 .. A213561 . A213562 . A213563
h*(h+1)/2 .. h^2 ........ A213564 . A213565 . A101094
2^(h-1) .... h .......... A213568 . A213569 . A047520
2^(h-1) .... h^2 ........ A213573 . A213574 . A213575
h .......... Fibo(h) .... A213576 . A213577 . A213578
Fibo(h) .... h .......... A213579 . A213580 . A053808
Fibo(h) .... Fibo(h) .... A067418 . A027991 . A067988
Fibo(h+1) .. h .......... A213584 . A213585 . A213586
Fibo(n+1) .. Fibo(h+1) .. A213587 . A213588 . A213589
h^2 ........ Fibo(h) .... A213590 . A213504 . A213557
Fibo(h) .... h^2 ........ A213566 . A213567 . A213570
h .......... -1+2^h ..... A213571 . A213572 . A213581
-1+2^h ..... h .......... A213582 . A213583 . A156928
-1+2^h ..... -1+2^h ..... A213747 . A213748 . A213749
h .......... 2*h-1 ...... A213750 . A007585 . A002417
2*h-1 ...... h .......... A213751 . A051662 . A006325
2*h-1 ...... 2*h-1 ...... A213752 . A100157 . A071238
2*h-1 ...... -1+2^h ..... A213753 . A213754 . A213755
-1+2^h ..... 2*h-1 ...... A213756 . A213757 . A213758
2^(n-1) .... 2*h-1 ...... A213762 . A213763 . A213764
2*h-1 ...... Fibo(h) .... A213765 . A213766 . A213767
Fibo(h) .... 2*h-1 ...... A213768 . A213769 . A213770
Fibo(h+1) .. 2*h-1 ...... A213774 . A213775 . A213776
Fibo(h) .... Fibo(h+1) .. A213777 . A001870 . A152881
h .......... 1+[h/2] .... A213778 . A213779 . A213780
1+[h/2] .... h .......... A213781 . A213782 . A005712
1+[h/2] .... [(h+1)/2] .. A213783 . A213759 . A213760
h .......... 3*h-2 ...... A213761 . A172073 . A002419
3*h-2 ...... h .......... A213771 . A213772 . A132117
3*h-2 ...... 3*h-2 ...... A213773 . A214092 . A213818
h .......... 3*h-1 ...... A213819 . A213820 . A153978
3*h-1 ...... h .......... A213821 . A033431 . A176060
3*h-1 ...... 3*h-1 ...... A213822 . A213823 . A213824
3*h-1 ...... 3*h-2 ...... A213825 . A213826 . A213827
3*h-2 ...... 3*h-1 ...... A213828 . A213829 . A213830
2*h-1 ...... 3*h-2 ...... A213831 . A213832 . A212560
3*h-2 ...... 2*h-1 ...... A213833 . A130748 . A213834
h .......... 4*h-3 ...... A213835 . A172078 . A051797
4*h-3 ...... h .......... A213836 . A213837 . A071238
4*h-3 ...... 2*h-1 ...... A213838 . A213839 . A213840
2*h-1 ...... 4*h-3 ...... A213841 . A213842 . A213843
2*h-1 ...... 4*h-1 ...... A213844 . A213845 . A213846
4*h-1 ...... 2*h-1 ...... A213847 . A213848 . A180324
[(h+1)/2] .. [(h+1)/2] .. A213849 . A049778 . A213850
h .......... C(2*h-2,h-1) A213853
...
Suppose that u = (u(n)) and v = (v(n)) are sequences having generating functions U(x) and V(x), respectively. Then the convolution u**v has generating function U(x)*V(x). Accordingly, if u and v are homogeneous linear recurrence sequences, then every row of the convolution array T satisfies the same homogeneous linear recurrence equation, which can be easily obtained from the denominator of U(x)*V(x). Also, every column of T has the same homogeneous linear recurrence as v.

Examples

			Northwest corner (the array is read by southwest falling antidiagonals):
  1,  4, 10, 20,  35,  56,  84, ...
  2,  7, 16, 30,  50,  77, 112, ...
  3, 10, 22, 40,  65,  98, 140, ...
  4, 13, 28, 50,  80, 119, 168, ...
  5, 16, 34, 60,  95, 140, 196, ...
  6, 19, 40, 70, 110, 161, 224, ...
T(6,1) = (1)**(6) = 6;
T(6,2) = (1,2)**(6,7) = 1*7+2*6 = 19;
T(6,3) = (1,2,3)**(6,7,8) = 1*8+2*7+3*6 = 40.
		

Crossrefs

Cf. A000027.

Programs

  • Mathematica
    b[n_] := 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}]]
    r[n_] := Table[t[n, k], {k, 1, 60}]  (* A213500 *)
  • PARI
    t(n,k) = sum(i=0, k - 1, (k - i) * (n + i));
    tabl(nn) = {for(n=1, nn, for(k=1, n, print1(t(k,n - k + 1),", ");); print(););};
    tabl(12) \\ Indranil Ghosh, Mar 26 2017
    
  • Python
    def t(n, k): return sum((k - i) * (n + i) for i in range(k))
    for n in range(1, 13):
        print([t(k, n - k + 1) for k in range(1, n + 1)]) # Indranil Ghosh, Mar 26 2017

Formula

T(n,k) = 4*T(n,k-1) - 6*T(n,k-2) + 4*T(n,k-3) - T(n,k-4).
T(n,k) = 2*T(n-1,k) - T(n-2,k).
G.f. for row n: x*(n - (n - 1)*x)/(1 - x)^4.

A067418 Triangle A067330 with rows read backwards.

Original entry on oeis.org

1, 2, 1, 5, 3, 2, 10, 7, 5, 3, 20, 15, 12, 8, 5, 38, 30, 25, 19, 13, 8, 71, 58, 50, 40, 31, 21, 13, 130, 109, 96, 80, 65, 50, 34, 21, 235, 201, 180, 154, 130, 105, 81, 55, 34, 420, 365, 331, 289, 250, 210, 170, 131, 89, 55, 744, 655, 600, 532, 469, 404, 340, 275, 212, 144, 89, 1308, 1164, 1075, 965
Offset: 0

Views

Author

Wolfdieter Lang, Feb 15 2002

Keywords

Comments

The column m (without leading 0's) gives the convolution of Fibonacci numbers F(n+1) := A000045(n+1), n>=0, with those with m-shifted index: a(n+m,m)=sum(F(k+1)*F(m+n+1-k),k=0..n), n>=0, m=0,1,...
The row polynomials p(n,x) := sum(a(n,m)*x^m,m=0..n) are generated by A(z)*(A(z)-x*A(x*z))/(1-x), with A(x) := 1/(1-x-x^2) (g.f. for Fibonacci F(n+1), n>=0).
The columns give A001629(n+2), A023610, A067331-4, A067430-1, A067977-8 for m= 0..9, respectively. Row sums give A067988.

Examples

			{1}; {2,1}; {5,3,2}; {10,7,5,3}; ...; p(2,n)=5+3*x+2*x^2.
		

Programs

  • Mathematica
    Reverse /@ Table[Sum[Fibonacci[k + 1] Fibonacci[n - k + 1], {k, 0, m}], {n, 0, 11}, {m, 0, n}] // Flatten (* Michael De Vlieger, Apr 11 2016 *)

Formula

a(n, m)=A067330(n, n-m), n>=m>=0, else 0.
a(n, m)= (((3*(n-m)+5)*F(n-m+1)+(n-m+1)*F(n-m))*F(m+1)+((n-m)*F(n-m+1)+2*(n-m+1)*F(n-m))*F(m))/5.
G.f. for column m=0, 1, ...: (x^m)*(F(m+1)+F(m)*x)/(1-x-x^2)^2, with F(m) := A000045(m) (Fibonacci).
a(n, m) = ((-1)^m*F(n-2*m+1)-m*L(n+2)+n*L(n+2)+5*F(n)+4*F(n-1))/5, with F(-n) = (-1)^(n+1)*F(n), hence a(n, m) = (2*(n-m+1)*L(n+2)-A067990(n, m))/5, n>=m>=0. - Ehren Metcalfe, Apr 11 2016

A104509 Matrix inverse of triangle A104505, which is the right-hand side of triangle A084610 of coefficients in (1 + x - x^2)^n.

Original entry on oeis.org

1, 1, -1, 3, -2, 1, 4, -6, 3, -1, 7, -12, 10, -4, 1, 11, -25, 25, -15, 5, -1, 18, -48, 60, -44, 21, -6, 1, 29, -91, 133, -119, 70, -28, 7, -1, 47, -168, 284, -296, 210, -104, 36, -8, 1, 76, -306, 585, -699, 576, -342, 147, -45, 9, -1, 123, -550, 1175, -1580, 1485, -1022, 525, -200, 55, -10, 1, 199, -979, 2310, -3454, 3641
Offset: 0

Views

Author

Paul D. Hanna, Mar 11 2005

Keywords

Comments

Riordan array ( (1 + x^2)/(1 - x - x^2), -x/(1 - x - x^2) ) belonging to the hitting time subgroup of the Riordan group (see Peart and Woan). - Peter Bala, Jun 29 2015
The sums of absolute values along steep diagonals in this triangle are: 1, 1, 3, 4 + |-1|, 7 + |-2|, 11 + |-6|, 18 + |-12| + 1, ... and these are the tribonacci numbers A000213 that begin with 1, 1, 1, 3. To see this, replace the y in the g.f. A(x,y) = (1 + x^2)/(1-x-x^2 + x*y) with y=-x^2, multiply by x, and add 1, to obtain the g.f. (1 - x^2)/(1-x-x^2-x^3) for A000213. - Noah Carey and Greg Dresden, Nov 02 2021

Examples

			Rows begin:
   1;
   1,   -1;
   3,   -2,   1;
   4,   -6,   3,   -1;
   7,  -12,  10,   -4,   1;
  11,  -25,  25,  -15,   5,   -1;
  18,  -48,  60,  -44,  21,   -6,   1;
  29,  -91, 133, -119,  70,  -28,   7,  -1;
  47, -168, 284, -296, 210, -104,  36,  -8, 1;
  76, -306, 585, -699, 576, -342, 147, -45, 9, -1; ...
		

Crossrefs

Leftmost column is A000204 (Lucas numbers). Other columns include: A045925, A067988. Row sums are: {1,0,2,0,2,0,2,...}. Absolute row sums form: A099425. Antidiagonal sums are: {1,1,2,2,2,2,2,...}. Absolute antidiagonal sums are: A084214.

Programs

  • Maple
    S:= series((1 + x^2)/(1-x-x^2 + x*y),x, 20):
    for n from 0 to 19 do R[n]:= coeff(S,x,n) od:
    seq(seq(coeff(R[n],y,j),j=0..n), n=0..19); # Robert Israel, Jun 30 2015
  • Mathematica
    nmax = 11;
    T[n_, k_] := Coefficient[(1 + x - x^2)^n, x, n + k];
    M = Table[T[n, k], {n, 0, nmax}, {k, 0, nmax}] // Inverse;
    Table[M[[n+1, k+1]], {n, 0, nmax}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 27 2019 *)
  • PARI
    { T(n,k) = my(X=x+x*O(x^n),Y=y+y*O(y^k)); polcoeff(polcoeff((1 + X^2)/(1-X-X^2 + X*Y),n,x),k,y); }
    
  • PARI
    { tabl(nn) = my(m = matrix(nn, nn, n, k, n--; k--; if((nMichel Marcus, Jun 30 2015
    
  • PARI
    { A104509(n,k) = if(n==0, k==0, (-1)^k * sum(i=0, (n-k)\2, n/(n-i) * binomial(n-k-i,i) * binomial(n-i,k) )); } \\ Max Alekseyev, Oct 11 2021

Formula

For n>=1, a(n,k) = (-1)^k * Sum_{i=0..[(n-k)/2]} n/(n-i) * binomial(n-i,i) * binomial(n-2*i,k) = (-1)^k * Sum_{i=0..[(n-k)/2]} n/(n-i) * binomial(n-k-i,i) * binomial(n-i,k). - Max Alekseyev, Oct 11 2021
G.f.: A(x, y) = (1 + x^2)/(1-x-x^2 + x*y).
G.f. for column k: g_k(x) = -(x^2+1)*x^k/(x^2+x-1)^(k+1). - Robert Israel, Jun 30 2015
G.f. for row n>=1 is the Lucas polynomial L_n(1-x). - Max Alekseyev, Oct 11 2021
Showing 1-3 of 3 results.