OFFSET
0,2
LINKS
Indranil Ghosh, Rows 0..125, flattened
FORMULA
From G. C. Greubel, Aug 21 2025: (Start)
T(n, k) = binomial(n,k) + 3*binomial(n+4, k+2) - 18.
Sum_{k=0..n} (-1)^k*T(n, k) = 6*(1+(-1)^n)* floor((n+1)/2) + [n=0]. (End)
EXAMPLE
Triangle starts:
1;
13, 13;
28, 44, 28;
46, 90, 90, 46;
67, 154, 198, 154, 67;
91, 239, 370, 370, 239, 91;
...
MATHEMATICA
t[n_, k_]:=Binomial[n + 4, k + 2 ] + Binomial[n + 3, k + 1] + Binomial[n + 3, k + 2] + Binomial[n + 2, k] + Binomial[n + 2, k + 1] + Binomial[n + 2, k + 2] + Binomial[n + 1, k] + Binomial[n + 1, k + 1] + Binomial[n, k] ; T[n_, k_]:=t[n, k] - t[0, 0] + 1; Flatten[Table[T[n, k], {n, 0, 9}, {k, 0, n}]] (* Indranil Ghosh, Mar 13 2017 *)
PROG
(PARI) alias(C, binomial);
t(n, k) = C(n+4, k+2) + C(n+3, k+1) + C(n+3, k+2) + C(n+2, k) + C(n+2, k+1) + C(n+2, k+2) + C(n+1, k) + C(n+1, k+1) + C(n, k);
T(n, k) = t(n, k)-t(0, 0)+1;
tabl(nn) = for (n=0, nn, for (k=0, n, print1(T(n, k), ", ")); print());
\\ Michel Marcus, Mar 13 2017
(Magma)
A026907:= func< n, k | Binomial(n, k) + 3*Binomial(n+4, k+2) - 18 >;
[A026907(n, k): k in [0..n], n in [0..12]]; // G. C. Greubel, Aug 21 2025
(SageMath)
def A026907(n, k): return binomial(n, k) +3*binomial(n+4, k+2) -18
print(flatten([[A026907(n, k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Aug 21 2025
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
STATUS
approved
