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 A156645 #9 Sep 08 2022 08:45:41 %S A156645 1,1,1,1,36,1,1,1225,1225,1,1,41616,1416100,41616,1,1,1413721, %T A156645 1634261476,1634261476,1413721,1,1,48024900,1885939157025, %U A156645 64069586905104,1885939157025,48024900,1,1,1631432881,2176372249076025,2511659716192658889,2511659716192658889,2176372249076025,1631432881,1 %N A156645 Triangle T(n, k, m) = b(n,m)/(b(k,m)*b(n-k,m)), where b(n, k) = Product_{j=1..n} (1 - ChebyshevT(j, k+1)^2), b(n, 0) = n!, and m = 2, read by rows. %H A156645 G. C. Greubel, <a href="/A156645/b156645.txt">Rows n = 0..25 of the triangle, flattened</a> %F A156645 T(n, k, m) = b(n,m)/(b(k,m)*b(n-k,m)), where b(n, k) = Product_{j=1..n} (1 - ChebyshevT(j, k+1)^2), b(n, 0) = n!, and m = 2. %F A156645 From _G. C. Greubel_, Jul 03 2021: (Start) %F A156645 T(n, k, m) = b(n,m)/(b(k,m)*b(n-k,m)), where b(n, k) = (1/2^n)*Product_{j=1..n} (1 - ChebyshevT(2*j, k+1)), b(n, 0) = n!, and m = 2. %F A156645 T(n, k, m) = Product_{j=1..n-k} ( (1 - ChebyshevT(2*j+2*k, m+1))/(1 - ChebyshevT(2*j, m+1)) ) with m = 2. (End) %e A156645 Triangle begins as: %e A156645 1; %e A156645 1, 1; %e A156645 1, 36, 1; %e A156645 1, 1225, 1225, 1; %e A156645 1, 41616, 1416100, 41616, 1; %e A156645 1, 1413721, 1634261476, 1634261476, 1413721, 1; %e A156645 1, 48024900, 1885939157025, 64069586905104, 1885939157025, 48024900, 1; %t A156645 (* First program *) %t A156645 b[n_, k_]:= b[n,k]= If[k==0, n!, Product[1 -ChebyshevT[j, k+1]^2, {j,n}]]; %t A156645 T[n_, k_, m_]= b[n,m]/(b[k,m]*b[n-k,m]); %t A156645 Table[T[n, k, 2], {n,0,12}, {k,0,n}]//Flatten (* modified by _G. C. Greubel_, Jul 03 2021 *) %t A156645 (* Second program *) %t A156645 T[n_, k_, m_]:= T[n, k, m]= Product[(1 - ChebyshevT[2*j+2*k, m+1])/(1 - ChebyshevT[2*j, m+1]), {j, n-k}]; %t A156645 Table[T[n,k,2], {n,0,12}, {k,0,n}]//Flatten (* _G. C. Greubel_, Jul 03 2021 *) %o A156645 (Magma) %o A156645 b:= func< n, k | n eq 0 select 1 else k eq 0 select Factorial(n) else (&*[1 - Evaluate(ChebyshevT(j), k+1)^2 : j in [1..n]]) >; %o A156645 T:= func< n,k,m | b(n,m)/(b(k,m)*b(n-k,m)) >; %o A156645 [T(n,k,2): k in [0..n], n in [0..12]]; // _G. C. Greubel_, Jul 03 2021 %o A156645 (Sage) %o A156645 def b(n, k): return factorial(n) if (k==0) else product( 1 - chebyshev_T(j, k+1)^2 for j in (1..n) ) %o A156645 def T(n, k, m): return b(n,m)/(b(k,m)*b(n-k,m)) %o A156645 flatten([[T(n, k, 2) for k in (0..n)] for n in (0..12)]) # _G. C. Greubel_, Jul 03 2021 %Y A156645 Cf. A007318 (m=0), A173585 (m=1), this sequence (m=2), A156646 (m=10). %Y A156645 Cf. A123583, A156647. %K A156645 nonn,tabl %O A156645 0,5 %A A156645 _Roger L. Bagula_, Feb 12 2009 %E A156645 Edited by _G. C. Greubel_, Jul 03 2021