|
%I #23 Apr 23 2024 06:04:51
%S 1,1,2,1,12,20,1,30,300,488,1,56,1400,13664,22160,1,90,4200,102480,
%T 997200,1616672,1,132,9900,450912,10969200,106700352,172976960,1,182,
%U 20020,1465464,66546480,1618288672,15740903360,25518205568
%N Triangle read by rows: T(n, k) = (-8)^k*binomial(2*n, 2*k)*Euler(2*k, 1/2).
%F Triangle T(n, k), 0 <= k <=n, read by rows, given by [1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [2, 8, 18, 32, 50, 72, 98, ...] where DELTA is the operator defined in A084938. - _Philippe Deléham_, Apr 21 2024
%F T(n, k) = binomial(2*n, 2*k) * 2^k * abs(Euler(2*k)) = A086645(n, k) * A000079(k) * A000364(k). - _Philippe Deléham_, Apr 23 2024
%e Triangle starts:
%e [0] 1;
%e [1] 1, 2;
%e [2] 1, 12, 20;
%e [3] 1, 30, 300, 488;
%e [4] 1, 56, 1400, 13664, 22160;
%e [5] 1, 90, 4200, 102480, 997200, 1616672;
%e [6] 1, 132, 9900, 450912, 10969200, 106700352, 172976960;
%e [7] 1, 182, 20020, 1465464, 66546480, 1618288672, 15740903360, 25518205568;
%p T := (n, k) -> (-8)^k*binomial(2*n, 2*k)*euler(2*k, 1/2):
%p seq(print(seq(T(n, k), k = 0..n)), n = 0..7);
%t Table[(-8)^k*Binomial[2*n, 2*k]*EulerE[2*k, 1/2], {n, 0, 10}, {k, 0, n}] (* _Paolo Xausa_, Apr 17 2024 *)
%o (SageMath)
%o def DelehamDelta(R, S, dim):
%o ring = PolynomialRing(ZZ, 'x')
%o x = ring.gen()
%o A = [R(k) + x * S(k) for k in range(dim)]
%o C = [ring(0)] + [ring(1) for i in range(dim)]
%o for k in range(1, dim + 1):
%o for n in range(k - 1, 0, -1):
%o C[n] = C[n-1] + C[n+1] * A[n-1]
%o yield list(C[1])
%o def A371637_triangle(dim):
%o a = lambda n: 1 - n % 2
%o b = lambda n: 2*(n + 1)^2
%o for row in DelehamDelta(a, b, dim): print(row)
%o A371637_triangle(8) # _Peter Luschny_, Apr 21 2024
%Y Cf. A001105, A002939 (column 1), A012816 (main diagonal), A371683 (row sums), A371684 (alternating row sums).
%Y Cf. A081658, A084938, A086645, A000079, A000364.
%K nonn
%O 0,3
%A _Peter Luschny_, Apr 02 2024
|
