A157117 - OEIS

%I #7 Sep 08 2022 08:45:41

%S 2,1,1,1,8,1,1,131,131,1,1,8204,29216,8204,1,1,2097187,44136233,

%T 44136233,2097187,1,1,2147483736,846839476071,503464582368,

%U 846839476071,2147483736,1,1,8796093022411,150092195453359483,288159195861579519,288159195861579519,150092195453359483,8796093022411,1

%N Triangle T(n, k) = f(n, k) + f(n, n-k), where f(n, k) = A008292(n*k + 1, n-k) if k <= n otherwise A008292(n*(n-k), k), read by rows.

%H G. C. Greubel, <a href="/A157117/b157117.txt">Rows n = 0..25 of the triangle, flattened</a>

%F T(n, k) = f(n, k) + f(n, n-k), where f(n, k) = A008292(n*k + 1, n-k) if k <= n otherwise A008292(n*(n-k), k).

%F T(n, n-k) = T(n, k).

%e   2;

%e   1,          1;

%e   1,          8,            1;

%e   1,        131,          131,            1;

%e   1,       8204,        29216,         8204,            1;

%e   1,    2097187,     44136233,     44136233,      2097187,          1;

%e   1, 2147483736, 846839476071, 503464582368, 846839476071, 2147483736, 1;

%t f[n_, k_]:= If[k<=n, Eulerian[n*k+1,n-k], Eulerian[n*(n-k)+1,k]];

%t T[n_, k_]:= f[n,k] + f[n,n-k];

%t Table[T[n, k], {n,0,10}, {k,0,n}]//Flatten (* modified by _G. C. Greubel_, Jan 11 2022 *)

%o (Magma)

%o Eulerian:= func< n,k | (&+[(-1)^j*Binomial(n+1,j)*(k-j+1)^n: j in [0..k+1]]) >;

%o f:= func< n,k | k le n select Eulerian(n*k+1,n-k) else Eulerian(n*(n-k)+1, k) >;

%o A157117:= func< n,k | f(n,k) + f(n,n-k) >;

%o [A157117(n,k): k in [0..n], n in [0..12]]; // _G. C. Greubel_, Jan 11 2022

%o (Sage)

%o def Eulerian(n,k): return sum((-1)^j*binomial(n+1,j)*(k-j+1)^n for j in (0..k+1))

%o def f(n,k): return Eulerian(n*k+1,n-k) if (k<n+1) else Eulerian(n*(n-k)+1, k)

%o def A157117(n,k): return f(n,k) + f(n,n-k)

%o flatten([[A157117(n,k) for k in (0..n)] for n in (0..12)]) # _G. C. Greubel_, Jan 11 2022

%Y Cf. A008292, A157114, A157118.

%K nonn,tabl

%O 0,1

%A _Roger L. Bagula_, Feb 23 2009

%E Edited by _G. C. Greubel_, Jan 11 2022