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.

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

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