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 A168552 #7 Apr 01 2022 09:15:15 %S A168552 1,1,1,1,3,1,1,11,11,1,1,43,140,43,1,1,159,1244,1244,159,1,1,551,8779, %T A168552 19954,8779,551,1,1,1819,54249,236347,236347,54249,1819,1,1,5811, %U A168552 309742,2353021,4440834,2353021,309742,5811,1,1,18167,1684634,21025310,67447952,67447952,21025310,1684634,18167,1 %N A168552 Expansion of g.f. (1/2)*( a*(1+x)^n + b*(1-x)^(n+2)*LerchPhi(x, -n-1, 1) + c*2^(n+1)*(1-x)^(n+1)*LerchPhi(x, -n, 1/2) ), where a = 3, b = -3, and c = 1, read by rows. %H A168552 G. C. Greubel, <a href="/A168552/b168552.txt">Rows n = 0..50 of the triangle, flattened</a> %F A168552 G.f.: (1/2)*( a*(1+x)^n + b*(1-x)^(n+2)*LerchPhi(x, -n-1, 1) + c*2^(n+1)*(1 - x)^(n+1)*LerchPhi(x, -n, 1/2) ), where a = 3, b = -3, and c = 1. %F A168552 From _G. C. Greubel_, Mar 31 2022: (Start) %F A168552 T(n, k) = (1/2)*( a*binomial(n,k) + sum( (-1)^(k-j)*(b*binomial(n+2, k-j)*(j+1)^(n+1) + 2*c*binomial(n+1,k-j)*(2*j+1)^n) for j in (0..k)) ), with a = 3, b = -3, and c = 1. %F A168552 T(n, n-k) = T(n, k). (End) %e A168552 Triangle begins as: %e A168552 1; %e A168552 1, 1; %e A168552 1, 3, 1; %e A168552 1, 11, 11, 1; %e A168552 1, 43, 140, 43, 1; %e A168552 1, 159, 1244, 1244, 159, 1; %e A168552 1, 551, 8779, 19954, 8779, 551, 1; %e A168552 1, 1819, 54249, 236347, 236347, 54249, 1819, 1; %e A168552 1, 5811, 309742, 2353021, 4440834, 2353021, 309742, 5811, 1; %e A168552 1, 18167, 1684634, 21025310, 67447952, 67447952, 21025310, 1684634, 18167, 1; %t A168552 p[x_, n_, a_, b_, c_]= (1/2)*(a*(1+x)^n + b*(1-x)^(n+2)*LerchPhi[x,-n-1,1] + c*2^(n+1)*(1-x)^(n+1)*LerchPhi[x,-n,1/2]); %t A168552 Table[CoefficientList[p[x,n,3,-3,1], x], {n,0,10}]//Flatten (* modified by _G. C. Greubel_, Mar 31 2022 *) %o A168552 (Sage) %o A168552 def A168552(n,k,a,b,c): return (1/2)*( a*binomial(n,k) + sum( (-1)^(k-j)*(b*binomial(n+2, k-j)*(j+1)^(n+1) + 2*c*binomial(n+1,k-j)*(2*j+1)^n) for j in (0..k)) ) %o A168552 flatten([[A168552(n,k,3,-3,1) for k in (0..n)] for n in (0..12)]) # _G. C. Greubel_, Mar 31 2022 %Y A168552 Cf. A001263, A168517, A168518, A168549, A168551. %K A168552 nonn,tabl %O A168552 0,5 %A A168552 _Roger L. Bagula_, Nov 29 2009 %E A168552 Edited by _G. C. Greubel_, Mar 31 2022