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 A173740 #14 Feb 14 2021 01:37:04 %S A173740 1,1,1,1,4,1,1,5,5,1,1,6,8,6,1,1,7,12,12,7,1,1,8,17,22,17,8,1,1,9,23, %T A173740 37,37,23,9,1,1,10,30,58,72,58,30,10,1,1,11,38,86,128,128,86,38,11,1, %U A173740 1,12,47,122,212,254,212,122,47,12,1,1,13,57,167,332,464,464,332,167,57,13,1 %N A173740 Triangle T(n,k) = binomial(n,k) + 2 for 1 <= k <= n - 1, n >= 2, and T(n,0) = T(n,n) = 1 for n >= 0, read by rows. %C A173740 For n >= 1, row n sums to A131520(n). %H A173740 G. C. Greubel, <a href="/A173740/b173740.txt">Rows n = 0..100 of the triangle, flattened</a> %F A173740 From _Franck Maminirina Ramaharo_, Dec 08 2018:(Start) %F A173740 T(n,k) = A007318(n,k) + 2*(1 - A103451(n,k)). %F A173740 T(n,k) = 3*A007318(n,k) - 2*A132044(n,k). %F A173740 n-th row polynomial is 1 - (-1)^(2^n) + (1 + x)^n + 2*(x - x^n)/(1 - x). %F A173740 G.f.: (1 - (1 + x)*y + 3*x*y^2 - 2*(x + x^2)*y^3)/((1 - y)*(1 - x*y)*(1 - y - x*y)). %F A173740 E.g.f.: (2 - 2*x + 2*x*exp(y) - 2*exp(x*y) + (1 - x)*exp((1 + x)*y))/(1 - x). (End) %F A173740 Sum_{k=0..n} T(n, k) = 2^n + 2*(n - 1 + [n=0]) = 2*A100314(n). - _G. C. Greubel_, Feb 13 2021 %e A173740 Triangle begins: %e A173740 1; %e A173740 1, 1; %e A173740 1, 4, 1; %e A173740 1, 5, 5, 1; %e A173740 1, 6, 8, 6, 1; %e A173740 1, 7, 12, 12, 7, 1; %e A173740 1, 8, 17, 22, 17, 8, 1; %e A173740 1, 9, 23, 37, 37, 23, 9, 1; %e A173740 1, 10, 30, 58, 72, 58, 30, 10, 1; %e A173740 1, 11, 38, 86, 128, 128, 86, 38, 11, 1; %e A173740 1, 12, 47, 122, 212, 254, 212, 122, 47, 12, 1; %e A173740 ... %t A173740 T[n_, m_] = Binomial[n, m] + 2*If[m*(n - m) > 0, 1, 0]; %t A173740 Flatten[Table[T[n, m], {n, 0, 10}, {m, 0, n}]] %o A173740 (Maxima) T(n,k) := if k = 0 or k = n then 1 else binomial(n, k) + 2$ %o A173740 create_list(T(n, k), n, 0, 12, k, 0, n); /* _Franck Maminirina Ramaharo_, Dec 08 2018 */ %o A173740 (Sage) %o A173740 def T(n, k): return 1 if (k==0 or k==n) else binomial(n, k) + 2 %o A173740 flatten([[T(n,k) for k in (0..n)] for n in (0..12)]) # _G. C. Greubel_, Feb 13 2021 %o A173740 (Magma) %o A173740 T:= func< n,k | k eq 0 or k eq n select 1 else Binomial(n,k) + 2 >; %o A173740 [T(n,k): k in [0..n], n in [0..12]]; // _G. C. Greubel_, Feb 13 2021 %Y A173740 Cf. A100314, A103451, A131520, A156050. %Y A173740 Sequences of the form binomial(n, k) + q: A132823 (q=-2), A132044 (q=-1), A007318 (q=0), A132735 (q=1), this sequence (q=2), A173741 (q=4), A173742 (q=6). %K A173740 nonn,tabl,easy %O A173740 0,5 %A A173740 _Roger L. Bagula_, Feb 23 2010 %E A173740 Edited and name clarified by _Franck Maminirina Ramaharo_, Dec 08 2018