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!)
A078123 Square of infinite lower triangular matrix A078122. 4
1, 2, 1, 5, 6, 1, 23, 51, 18, 1, 239, 861, 477, 54, 1, 5828, 32856, 25263, 4347, 162, 1, 342383, 3013980, 3016107, 699813, 39285, 486, 1, 50110484, 690729981, 865184724, 253656252, 19053063, 354051, 1458, 1, 18757984046, 406279238154 (list; table; graph; refs; listen; history; text; internal format)
OFFSET
0,2
LINKS
FORMULA
M(1, j) = A078125(j), M(j+1, j)=2*3^j.
EXAMPLE
Square of A078122 = A078123 as can be seen by 4 X 4 submatrix:
[1,_0,_0,0]^2=[_1,_0,_0,_0]
[1,_1,_0,0]___[_2,_1,_0,_0]
[1,_3,_1,0]___[_5,_6,_1,_0]
[1,12,_9,1]___[23,51,18,_1]
MAPLE
S:= proc(i, j) option remember;
add(M(i, k)*M(k, j), k=0..i)
end:
M:= proc(i, j) option remember; `if`(j=0 or i=j, 1,
add(S(i-1, k)*M(k, j-1), k=0..i-1))
end:
seq(seq(S(n, k), k=0..n), n=0..10); # Alois P. Heinz, Feb 27 2015
MATHEMATICA
S[i_, j_] := S[i, j] = Sum[M[i, k]*M[k, j], {k, 0, i}]; M[i_, j_] := M[i, j] = If[j == 0 || i == j, 1, Sum[S[i-1, k]*M[k, j-1], {k, 0, i-1}]]; Table[Table[S[n, k], {k, 0, n}], {n, 0, 10}] // Flatten (* Jean-François Alcover, Mar 06 2015, after Alois P. Heinz *)
CROSSREFS
Sequence in context: A178121 A302595 A113345 * A342968 A323312 A231774
KEYWORD
nonn,tabl
AUTHOR
Paul D. Hanna, Nov 18 2002
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 28 10:31 EDT 2024. Contains 371240 sequences. (Running on oeis4.)