login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A203003 Symmetric matrix based on A007598(n+1), by antidiagonals. 3
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 (list; table; graph; refs; listen; history; text; internal format)
OFFSET
1,2
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.
LINKS
EXAMPLE
Northwest corner:
1....4.....9....25....64
4....17....40...109...281
9....40....98...265...685
25...109...265..1865
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) *)
CROSSREFS
Sequence in context: A061886 A059815 A202670 * A339427 A319646 A214826
KEYWORD
nonn,tabl
AUTHOR
Clark Kimberling, Dec 27 2011
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 00:26 EDT 2024. Contains 371264 sequences. (Running on oeis4.)