|
%I #25 Jul 16 2021 01:33:30
%S 1,0,1,0,1,0,1,1,-1,0,1,4,-5,1,0,1,11,-14,1,2,0,1,26,-24,-29,36,-9,0,
%T 1,57,1,-244,281,-104,9,0,1,120,225,-1259,1401,-454,-83,50,0,1,247,
%U 1268,-5081,4621,911,-3422,1723,-267,0,1,502,5278,-16981,5335,30871
%N Triangle of coefficients of (1 - x)^n*B_n(x/(1 - x)), where B_n(x) is the n-th Bell polynomial.
%C The n-th row consists of the coefficients in the expansion of Sum_{j=0..n} A048993(n,j)*x^j*(1 - x)^(n - j), where A048993 is the triangle of Stirling numbers of second kind.
%D M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, pp. 824-825.
%H G. C. Greubel, <a href="/A122753/b122753.txt">Rows n = 0..50 of the irregular triangle, flattened</a>
%H M. Abramowitz and I. A. Stegun, eds., <a href="http://www.convertit.com/Go/ConvertIt/Reference/AMS55.ASP">Handbook of Mathematical Functions</a>, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
%H Peter Steinbach, <a href="http://www.jstor.org/stable/2691048">Golden fields: a case for the heptagon</a>, Math. Mag. 70 (1997), no. 1, 22-31.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/BellPolynomial.html">Bell Polynomial</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/StirlingNumberoftheSecondKind.html">Stirling Number of the Second Kind</a>
%F From _Franck Maminirina Ramaharo_, Oct 10 2018: (Start)
%F E.g.f.: exp((x/(1 - x))*(exp((1 - x)*y) - 1).
%F T(n,1) = A000295(n-1). (End)
%e Triangle begins:
%e 1;
%e 0, 1;
%e 0, 1;
%e 0, 1, 1, -1;
%e 0, 1, 4, -5, 1;
%e 0, 1, 11, -14, 1, 2;
%e 0, 1, 26, -24, -29, 36, -9;
%e 0, 1, 57, 1, -244, 281, -104, 9;
%e 0, 1, 120, 225, -1259, 1401, -454, -83, 50;
%e 0, 1, 247, 1268, -5081, 4621, 911, -3422, 1723, -267;
%e ... reformatted and extended. - _Franck Maminirina Ramaharo_, Oct 10 2018
%t Table[CoefficientList[Sum[StirlingS2[m, n]*x^n*(1-x)^(m-n), {n,0,m}], x], {m,0,10}]//Flatten
%o (Maxima)
%o P(x, n) := expand(sum(stirling2(n, j)*x^j*(1 - x)^(n - j), j, 0, n))$
%o T(n, k) := ratcoef(P(x, n), x, k)$
%o tabf(nn) := for n:0 thru nn do print(makelist(T(n, k), k, 0, hipow(P(x, n), x)))$ /* _Franck Maminirina Ramaharo_, Oct 10 2018 */
%o (Sage)
%o def p(n,x): return sum( stirling_number2(n, j)*x^j*(1-x)^(n-j) for j in (0..n) )
%o def T(n): return ( p(n,x) ).full_simplify().coefficients(sparse=False)
%o flatten([T(n) for n in (0..12)]) # _G. C. Greubel_, Jul 15 2021
%Y Cf. A000110, A048993, A088996, A122610.
%Y Cf. A123018, A123019, A123021, A123027, A123199, A123202, A123217, A123221.
%K sign,tabf
%O 0,12
%A _Roger L. Bagula_ and _Gary W. Adamson_, Sep 21 2006
%E Edited, new name, and offset corrected by _Franck Maminirina Ramaharo_, Oct 10 2018
|
