This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A174451 #9 Sep 08 2022 08:45:51 %S A174451 1,1,1,1,18,1,1,36,36,1,1,60,2160,60,1,1,90,5400,5400,90,1,1,126, %T A174451 11340,680400,11340,126,1,1,168,21168,1905120,1905120,21168,168,1,1, %U A174451 216,36288,4572288,411505920,4572288,36288,216,1,1,270,58320,9797760,1234517760,1234517760,9797760,58320,270,1 %N A174451 Triangle T(n, k, q) = n!*(n+1)!*q^k/((n-k)!(n-k+1)!) if floor(n/2) > k-1, otherwise n!*(n+1)!*q^(n-k)/(k!*(k+1)!) for q = 3, read by rows. %H A174451 G. C. Greubel, <a href="/A174451/b174451.txt">Rows n = 0..50 of the triangle, flattened</a> %F A174451 T(n, k, q) = n!*(n+1)!*q^k/((n-k)!(n-k+1)!) if floor(n/2) > k-1, otherwise n!*(n+1)!*q^(n-k)/(k!*(k+1)!) for q = 3. %F A174451 T(n, n-k, q) = T(n, k, q). %F A174451 From _G. C. Greubel_, Nov 29 2021: (Start) %F A174451 T(2*n, n, q) = q^n*(2*n+1)!*Catalan(n) for q = 3. %F A174451 T(n, k, q) = binomial(n, k)*binomial(n+1, k+1) * ( k!*(k+1)!*q^k/(n-k+1) if (floor(n/2) >= k), otherwise ((n-k)!)^2*q^(n-k) ), for q = 3. (End) %e A174451 Triangle begins as: %e A174451 1; %e A174451 1, 1; %e A174451 1, 18, 1; %e A174451 1, 36, 36, 1; %e A174451 1, 60, 2160, 60, 1; %e A174451 1, 90, 5400, 5400, 90, 1; %e A174451 1, 126, 11340, 680400, 11340, 126, 1; %e A174451 1, 168, 21168, 1905120, 1905120, 21168, 168, 1; %e A174451 1, 216, 36288, 4572288, 411505920, 4572288, 36288, 216, 1; %e A174451 1, 270, 58320, 9797760, 1234517760, 1234517760, 9797760, 58320, 270, 1; %t A174451 T[n_, k_, q_]:= If[Floor[n/2]>k-1, n!*(n+1)!*q^k/((n-k)!*(n-k+1)!), n!*(n+1)!*q^(n-k)/(k!*(k+1)!)]; %t A174451 Table[T[n, k, 3], {n,0,12}, {k,0,n}]//Flatten %o A174451 (Magma) %o A174451 F:= Factorial; // T = A174451 %o A174451 T:= func< n,k,q | Floor(n/2) gt k-1 select F(n)*F(n+1)*q^k/(F(n-k)*F(n-k+1)) else F(n)*F(n+1)*q^(n-k)/(F(k)*F(k+1)) >; %o A174451 [T(n,k,3): k in [0..n], n in [0..12]]; // _G. C. Greubel_, Nov 29 2021 %o A174451 (Sage) %o A174451 f=factorial %o A174451 def A174451(n,k,q): %o A174451 if ((n//2)>k-1): return f(n)*f(n+1)*q^k/(f(n-k)*f(n-k+1)) %o A174451 else: return f(n)*f(n+1)*q^(n-k)/(f(k)*f(k+1)) %o A174451 flatten([[A174451(n,k,3) for k in (0..n)] for n in (0..12)]) # _G. C. Greubel_, Nov 29 2021 %Y A174451 Cf. A174449 (q=1), A174450 (q=2), this sequence (q=3). %Y A174451 Cf. A000108. %K A174451 nonn,tabl,easy %O A174451 0,5 %A A174451 _Roger L. Bagula_, Mar 20 2010 %E A174451 Edited by _G. C. Greubel_, Nov 29 2021