login
Coefficients of certain polynomials (rising powers).
10

%I #34 Jul 18 2024 21:02:11

%S 1,2,1,6,6,2,24,36,22,6,120,240,210,100,24,720,1800,2040,1350,548,120,

%T 5040,15120,21000,17640,9744,3528,720,40320,141120,231840,235200,

%U 162456,78792,26136,5040,362880,1451520,2751840,3265920,2693880,1614816

%N Coefficients of certain polynomials (rising powers).

%C This is the unsigned triangle A048594 with rows read backwards.

%C The row polynomials p(n,y) := Sum_{m=0..n-1}a(n,m)*y^m, n>=1, are obtained from (log(x)*(-x*log(x))^n)*(d^n/dx^n)(1/log(x)), n>=1, after replacement of log(x) by y.

%C The gcd of row n is A075182(n). Row sums give A007840(n), n>=1.

%C The columns give A000142 (factorials), A001286 (Lah), 2* A075183, 2*A075184, 4*A075185, 4!*A075186, 4!*A075187 for m=0..6.

%C Coefficients T(n,k) of the differential operator expansion

%C [x^(1+y)D]^n = x^(n*y)[T(n,1)* (xD)^n / n! + y * T(n,2)* (xD)^(n-1) / (n-1)! + ... + y^(n-1) * T(n,n) * (xD)], where D = d/dx. Note that (xD)^n = Bell(n,:xD:), where (:xD:)^n = x^n * D^n and Bell(n,x) are the Bell / Touchard polynomials. See A094638. - _Tom Copeland_, Aug 22 2015

%H Vincenzo Librandi, <a href="/A075181/b075181.txt">Rows n = 1..100, flattened</a>

%H Y.-Z. Huang, J. Lepowsky and L. Zhang, <a href="https://arxiv.org/abs/math/0311235">A logarithmic generalization of tensor product theory for modules for a vertex operator algebra</a>, arXiv:math/0311235 [math.QA], 2003; Internat. J. Math. 17 (2006), no. 8, 975-1012. See page 984 eq. (3.9) MR2261644.

%H D. Lubell, <a href="http://www.jstor.org/stable/3647779">Problem 10992, problems and solutions</a>, Amer. Math. Monthly 110 (2003) p. 155. <a href="http://www.jstor.org/stable/4145207">Equal Sums of Reciprocal Products: 10992</a> (2004) pp. 827-829.

%F a(n, m) = (n-m)!*|S1(n, n-m)|, n>=m+1>=1, else 0, with S1(n, m) := A008275(n, m) (Stirling1).

%F a(n, m) = (n-m)*a(n-1, m)+(n-1)*a(n-1, m-1), if n>=m+1>=1, a(n, -1) := 0 and a(1, 0)=1, else 0.

%e Triangle starts:

%e 1;

%e 2,1;

%e 6,6,2;

%e 24,36,22,6;

%e ...

%e n=2: (x^2*log(x)^3)*(d^2/d^x^2)(1/log(x)) = 2 + log(x).

%p seq(seq(k!*abs(Stirling1(n,k)),k=n..1,-1),n=1..10); # _Robert Israel_, Jul 12 2015

%t Table[ Table[ k!*StirlingS1[n, k] // Abs, {k, 1, n}] // Reverse, {n, 1, 9}] // Flatten (* _Jean-François Alcover_, Jun 21 2013 *)

%o (PARI) {T(n, k)= if(k<0 || k>=n, 0, (-1)^k* stirling(n, n-k)* (n-k)!)} /* _Michael Somos_ Apr 11 2007 */

%Y Cf. A048594, A075178, A007840, A075182.

%Y Cf. A094638.

%K nonn,easy,tabl

%O 1,2

%A _Wolfdieter Lang_, Sep 19 2002