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

A203002 Array: row n shows the coefficients of the characteristic polynomial of the n-th principal submatrix of the symmetric matrix A203001; by antidiagonals.

Original entry on oeis.org

1, -1, 1, -3, 1, 1, -14, 21, -1, 1, -29, 162, -120, 1, 1, -48, 540, -1736, 844, -1, 1, -71, 1267, -8091, 17022, -5664, 1, 1, -98, 2475, -24908, 105503, -158690, 39045, -1, 1, -129, 4312, -60994, 408508, -1250056, 1416673
Offset: 1

Views

Author

Clark Kimberling, Dec 27 2011

Keywords

Comments

Let p(n)=p(n,x) be the characteristic polynomial of the n-th principal submatrix. The zeros of p(n) are positive, and they interlace the zeros of p(n+1).

Examples

			Top of the array:
1...-1
1...-3....1
1...-14...21....-1
1...-29...162...-120...1
		

Crossrefs

Programs

  • Mathematica
    f[k_] := Fibonacci[k]^2;
    U[n_] := NestList[Most[Prepend[#, 0]] &, #, Length[#] - 1] &[Table[f[k], {k, 1, n}]];
    L[n_] := Transpose[U[n]];
    F[n_] := CharacteristicPolynomial[L[n].U[n], x];
    c[n_] := CoefficientList[F[n], x]
    TableForm[Flatten[Table[F[n], {n, 1, 10}]]]
    Table[c[n], {n, 1, 12}]
    Flatten[%]
    TableForm[Table[c[n], {n, 1, 10}]]

A202605 Array: row n shows the coefficients of the characteristic polynomial of the n-th principal submatrix of the Fibonacci self-fusion matrix (A202453).

Original entry on oeis.org

1, -1, 1, -3, 1, 1, -6, 9, -1, 1, -9, 26, -24, 1, 1, -12, 52, -96, 64, -1, 1, -15, 87, -243, 326, -168, 1, 1, -18, 131, -492, 1003, -1050, 441, -1, 1, -21, 184, -870, 2392, -3816, 3265, -1155, 1, 1, -24, 246, -1404, 4871, -10500, 13710
Offset: 1

Views

Author

Clark Kimberling, Dec 21 2011

Keywords

Comments

Let p(n)=p(n,x) be the characteristic polynomial of the n-th principal submatrix. The zeros of p(n) are positive and interlace the zeros of p(n+1). (See the references and examples.)
Following is a guide to sequences (f(n)) for symmetric matrices (self-fusion matrices) and characteristic polynomials. Notation: F(k)=A000045(k) (Fibonacci numbers); floor(n*tau)=A000201(n) (lower Wythoff sequence); "periodic x,y" represents the sequence (x,y,x,y,x,y,...).
f(n)........ symmetric matrix.. char. polynomial
1............... A087062....... A202672
n............... A115262....... A202673
n^2............. A202670....... A202671
2n-1............ A202674....... A202675
3n-2............ A202676....... A202677
n(n+1)/2........ A185957....... A202678
2^n-1........... A202873....... A202767
2^(n-1)......... A115216....... A202868
floor(n*tau).... A202869....... A202870
F(n)............ A202453....... A202605
F(n+1).......... A202874....... A202875
Lucas(n)........ A202871....... A202872
F(n+2)-1........ A202876....... A202877
F(n+3)-2........ A202970....... A202971
(F(n))^2........ A203001....... A203002
(F(n+1))^2...... A203003....... A203004
C(2n,n)......... A115255....... A203005
(-1)^(n+1)...... A003983....... A076757
periodic 1,0.... A203905....... A203906
periodic 1,0,0.. A203945....... A203946
periodic 1,0,1.. A203947....... A203948
periodic 1,1,0.. A203949....... A203950
periodic 1,0,0,0 A203951....... A203952
periodic 1,2.... A203953....... A203954
periodic 1,2,3.. A203955....... A203956
...
In the cases listed above, the zeros of the characteristic polynomials are positive. If more general symmetric matrices are used, the zeros are all real but not necessarily positive - but they do have the interlace property. For a guide to such matrices and polynomials, see A202605.

Examples

			The 1st principal submatrix (ps) of A202453 is {{1}} (using Mathematica matrix notation), with p(1) = 1-x and zero-set {1}.
...
The 2nd ps is {{1,1},{1,2}}, with p(2) = 1-3x+x^2 and zero-set {0.382..., 2.618...}.
...
The 3rd ps is {{1,1,2},{1,2,3},{2,3,6}}, with p(3) = 1-6x+9x^2-x^3 and zero-set {0.283..., 0.426..., 8.290...}.
  ...
Top of the array A202605:
  1,   -1;
  1,   -3,    1;
  1,   -6,    9,   -1;
  1,   -9,   26,  -24,    1;
  1,  -12,   52,  -96,   64,   -1;
  1,  -15,   87, -243,  326, -168,    1;
		

Crossrefs

Programs

  • Mathematica
    f[k_] := Fibonacci[k];
    U[n_] := NestList[Most[Prepend[#, 0]] &, #, Length[#] - 1] &[Table[f[k], {k, 1, n}]];
    L[n_] := Transpose[U[n]];
    F[n_] := CharacteristicPolynomial[L[n].U[n], x];
    c[n_] := CoefficientList[F[n], x]
    TableForm[Flatten[Table[F[n], {n, 1, 10}]]]
    Table[c[n], {n, 1, 12}]
    Flatten[%]
    TableForm[Table[c[n], {n, 1, 10}]]

A115255 "Correlation triangle" of central binomial coefficients A000984.

Original entry on oeis.org

1, 2, 2, 6, 5, 6, 20, 14, 14, 20, 70, 46, 41, 46, 70, 252, 160, 134, 134, 160, 252, 924, 574, 466, 441, 466, 574, 924, 3432, 2100, 1672, 1534, 1534, 1672, 2100, 3432, 12870, 7788, 6118, 5506, 5341, 5506, 6118, 7788, 12870, 48620, 29172, 22692, 20152, 19174
Offset: 0

Views

Author

Paul Barry, Jan 18 2006

Keywords

Comments

Row sums are A033114. Diagonal sums are A115256. T(2n,n) is A115257. Corresponds to the triangle of antidiagonals of the correlation matrix of the sequence array for C(2n,n).
Let s=(1,2,6,20,...), (central binomial coefficients), and let T be the infinite square matrix whose n-th row is formed by putting n-1 zeros before the terms of s. Let T' be the transpose of T. Then A115255 represents the matrix product M=T'*T. M is the self-fusion matrix of s, as defined at A193722. See A203005 for characteristic polynomials of principal submatrices of M, with interlacing zeros. - Clark Kimberling, Dec 27 2011

Examples

			Triangle begins:
  1;
  2, 2;
  6, 5, 6;
  20, 14, 14, 20;
  70, 46, 41, 46, 70;
  252, 160, 134, 134, 160, 252;
Northwest corner (square format):
  1    2    6    20    70
  2    5    14   46    160
  6    14   41   134   466
  20   46   134  441   1534
		

Crossrefs

Programs

  • Mathematica
    s[k_] := Binomial[2 k - 2, k - 1];
    U = NestList[Most[Prepend[#, 0]] &, #, Length[#] - 1] &[Table[s[k], {k, 1, 15}]];
    L = Transpose[U]; M = L.U; TableForm[M]
    m[i_, j_] := M[[i]][[j]]; (* A115255 in square format *)
    Flatten[Table[m[i, n + 1 - i], {n, 1, 12}, {i, 1, n}]]
    f[n_] := Sum[m[i, n], {i, 1, n}] + Sum[m[n, j], {j, 1, n - 1}]; Table[f[n], {n, 1, 12}]
    Table[Sqrt[f[n]], {n, 1, 12}]  (* A006134 *)
    Table[m[1, j], {j, 1, 12}]     (* A000984 *)
    Table[m[j, j], {j, 1, 12}]     (* A115257 *)
    Table[m[j, j + 1], {j, 1, 12}] (* 2*A082578 *)
    (* Clark Kimberling, Dec 27 2011 *)

Formula

G.f.: 1/(sqrt(1-4*x)*sqrt(1-4*x*y)*(1-x^2*y)) (format due to Christian G. Bower).
T(n, k) = Sum_{j=0..n} [j<=k]*C(2*k-2*j, k-j)*[j<=n-k]*C(2*n-2*k-2*j, n-k-j).

A203003 Symmetric matrix based on A007598(n+1), by antidiagonals.

Original entry on oeis.org

1, 4, 4, 9, 17, 9, 25, 40, 40, 25, 64, 109, 98, 109, 64, 169, 281, 265, 265, 281, 169, 441, 740, 685, 723, 685, 740, 441, 1156, 1933, 1802, 1865, 1865, 1802, 1933, 1156, 3025, 5065, 4709, 4910, 4819, 4910, 4709, 5065, 3025, 7921, 13256, 12337, 12827
Offset: 1

Views

Author

Clark Kimberling, Dec 27 2011

Keywords

Comments

Let s=A007598(n+1) (squared Fibonacci numbers, beginning with F(2)), and let T be the infinite square matrix whose n-th row is formed by putting n-1 zeros before the terms of s. Let T' be the transpose of T. Then A203003 represents the matrix product M=T'*T. M is the self-fusion matrix of s, as defined at A193722. See A203004 for characteristic polynomials of principal submatrices of M, with interlacing zeros.

Examples

			Northwest corner:
1....4.....9....25....64
4....17....40...109...281
9....40....98...265...685
25...109...265..1865
		

Crossrefs

Programs

  • Mathematica
    s[k_] := Fibonacci[k + 1]^2;
    U = NestList[Most[Prepend[#, 0]] &, #, Length[#] - 1] &[Table[s[k], {k, 1, 15}]];
    L = Transpose[U]; M = L.U; TableForm[M]
    m[i_, j_] := M[[i]][[j]];  (* A203003 *)
    Flatten[Table[m[i, n + 1 - i], {n, 1, 12}, {i, 1, n}]]
    f[n_] := Sum[m[i, n], {i, 1, n}] + Sum[m[n, j], {j, 1, n - 1}]; Table[f[n], {n, 1, 12}]
    Table[Sqrt[f[n]], {n, 1, 12}]    (* A119996 *)
    Table[m[1, j], {j, 1, 12}]       (* A007598(n+1) *)
Showing 1-4 of 4 results.