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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A173779 Triangle T(n, k, q) = round( c(n,q)/(c(k,q)*c(n-k,q)) ), where c(n,q) = Product_{j=0..n} v(j, q)((1)), v(n, q) = M*v(n-1, q), v(0, q) = {1, 1, 1}, M = {{0, 1, 0}, {0, 0, 1}, {q^3, q^3, 0}}, and q = 4, read by rows. 4
 1, 1, 1, 1, 1, 1, 1, 128, 128, 1, 1, 128, 16384, 128, 1, 1, 8256, 1056768, 1056768, 8256, 1, 1, 16384, 135266304, 135266304, 135266304, 16384, 1, 1, 536576, 8791261184, 567036346368, 567036346368, 8791261184, 536576, 1, 1, 1576960, 846158888960, 108308337786880, 6985887787253760, 108308337786880, 846158888960, 1576960, 1 (list; table; graph; refs; listen; history; text; internal format)
 OFFSET 0,8 LINKS G. C. Greubel, Rows n = 0..50 of the triangle, flattened FORMULA T(n, k, q) = round( c(n,q)/(c(k,q)*c(n-k,q)) ), where c(n,q) = Product_{j=0..n} v(j, q)((1)), v(n, q) = M*v(n-1, q), v(0, q) = {1, 1, 1}, M = {{0, 1, 0}, {0, 0, 1}, {q^3, q^3, 0}}, and q = 4. T(n, k, q) = round( (1/f(k,q))*Product_{j=0..n-k} f(j+k,q)/f(j,q) ), where f(n, q) = 6*q*f(n-2, q) + 8*q^3*f(n-3, q), and f(0,q) = f(1,q) = f(2,q) = 1, and q = 4. - G. C. Greubel, Jul 06 2021 EXAMPLE Triangle begins as: 1; 1, 1; 1, 1, 1; 1, 128, 128, 1; 1, 128, 16384, 128, 1; 1, 8256, 1056768, 1056768, 8256, 1; 1, 16384, 135266304, 135266304, 135266304, 16384, 1; 1, 536576, 8791261184, 567036346368, 567036346368, 8791261184, 536576, 1; MATHEMATICA f[n_, q_]:= f[n, q] = If[n<3, 1, q^3*f[n-2, q] + q^3*f[n-3, q]]; c[n_, q_]:= Product[f[j, q], {j, 0, n}]; T[n_, k_, q_]:= Round[c[n, q]/(c[k, q]*c[n-k, q])]; Table[T[n, k, 4], {n, 0, 12}, {k, 0, n}]//Flatten (* G. C. Greubel, Jul 06 2021 *) PROG (Magma) function f(n, k) if n lt 3 then return 1; else return k^3*f(n-2, k) + k^3*f(n-3, k); end if; return f; end function; T:= func< n, k, q | Round( (&*[f(j+k, q)/f(j, q): j in [0..n-k]])/f(k, q) ) >; [T(n, k, 4): k in [0..n], n in [0..12]]; // G. C. Greubel, Jul 06 2021 (Sage) @CachedFunction def f(n, q): return 1 if (n<3) else q^3*f(n-2, q) + q^3*f(n-3, q) def T(n, k, q): return round(product( f(j+k, q)/f(j, q) for j in (0..n-k))/f(k, q)) flatten([[T(n, k, 4) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Jul 06 2021 CROSSREFS Cf. A173778 (q=2), this sequence (q=4). Cf. A173747, A173749, A173779. Sequence in context: A188829 A336776 A172532 * A187353 A248773 A217847 Adjacent sequences: A173776 A173777 A173778 * A173780 A173781 A173782 KEYWORD nonn,tabl,less AUTHOR Roger L. Bagula, Feb 24 2010 EXTENSIONS Definition corrected and edited by G. C. Greubel, Jul 06 2021 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.

Last modified April 16 08:27 EDT 2024. Contains 371698 sequences. (Running on oeis4.)