Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A075263
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A075263 Triangle of coefficients of polynomials H(n,x) formed from the first (n+1) terms of the power series expansion of ( -x/ln(1-x) )^(n+1), multiplied by n!. +0
5
1, 1, -1, 2, -3, 1, 6, -12, 7, -1, 24, -60, 50, -15, 1, 120, -360, 390, -180, 31, -1, 720, -2520, 3360, -2100, 602, -63, 1, 5040, -20160, 31920, -25200, 10206, -1932, 127, -1, 40320, -181440, 332640, -317520, 166824, -46620, 6050, -255, 1, 362880, -1814400, 3780000, -4233600, 2739240, -1020600 (list; table; graph; listen)
OFFSET

0,4

COMMENT

Special values: H(n,1)=0, H(2n,2)=0, H(n,-x) ~= ( x/ln(1+x) )^(n+1), for x>0. H'(n,1) = -1/n!, where H'(n,x) = d/dx H(n,x).

The zeros of these polynomials are all positive reals >= 1. If we order the zeros of H(n,x), {r_k, k=0..(n-1)}, by magnitude so that r_0 = 1, r_k > r_(k-1), for 0<k<n, then r_(n-k) = r_k/(r_k - 1) when 0<k<n, n>1, where r_(n/2) = 2 for even n.

Also product_{k=0..(n-1)} r_k = n!, r_(n-1) ~ C 2^n.

I believe that these numbers are the coefficients of the Eulerian polynomials An(z) written in powers of z-1. That is, the sequence is: A0(1); A1(1), A1'(1); A2(1), A2'(1), A2''(1)/2!; A3(1), A3'(1), A3''(1)/2!, A3'''(1)/3!; A4(1), A4'(1), A4''(1)/2!, A4'''(1)/3!, A4''''(1)/4! etc. My convention: A0(z)=z, A1(z)=z, A2(z)=z+z^2, A3(z)=z+4z^2+z^3, A4(z)=z+11z^2+11z^3+z^4, etc. - Louis Zulli (zullil(AT)lafayette.edu), Jan 19 2005

H(n,2) gives 1,-1,0,2,0,-16,0,272,0,-7936,0,..., see A009006 . - Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Aug 20 2007

FORMULA

Generated by [1, 1, 2, 2, 3, 3, ...] DELTA [ -1, 0, -2, 0, -3, 0, ...], where DELTA is the operator defined in A084938.

T(n, k) = Sum_{i=0..n-k} (-1)^(n-i)*C(n-k, i)*(i+1)^n; n >= 0, 0<=k<=n. E.g.f.: A(x, y) = -log(1-(1-exp(-x*y))/y). - Paul D. Hanna (pauldhanna(AT)juno.com), Jul 21 2005

EXAMPLE

H(0,x) = 1

H(1,x) = (1 - 1x)/1!

H(2,x) = (2 - 3x + 1x^2)/2!

H(3,x) = (6 - 12x + 7x^2 - 1x^3)/3!

H(4,x) = (24 - 60x + 50x^2 - 15x^3 + 1x^4)/4!

H(5,x) = (120 - 360x + 390x^2 - 180x^3 + 31x^4 - 1x^5)/5!

H(6,x) = (720 - 2520x + 3360x^2 - 2100x^3 + 602x^4 - 63x^5 + 1x^5)/6!

PROGRAM

(PARI) T(n, k)=if(k<0|k>n, 0, n!*polcoeff((-x/log(1-x+x^2*O(x^n)))^(n+1), k))

(PARI) T(n, k)=sum(i=0, n-k, (-1)^(n-i)*binomial(n-k, i)*(i+1)^n) (PARI) /* Using E.g.f. A(x, y): */ {T(n, k)=local(X=x+x*O(x^n), Y=y+y^2*O(y^(k))); n!*polcoeff(polcoeff(-log(1-(1-exp(-X*Y))/y), n, x), k, y)}

(PARI) /* DELEHAM's DELTA: T(n, k) = [x^(n-k)*y^k] P(n, 0) */ {P(n, k)=if(n<0|k<0, 0, if(n==0, 1, P(n, k-1)+(x*(k\2+1)+y*(-(k\2+1)*((k+1)%2)))*P(n-1, k+1)))} {T(n, k)=polcoeff(polcoeff(P(n, 0), n-k, x), k, y)}

CROSSREFS

Cf. A075264, A028246, A084938.

Sequence in context: A105613 A135894 A130850 this_sequence A130405 A058372 A128264

Adjacent sequences: A075260 A075261 A075262 this_sequence A075264 A075265 A075266

KEYWORD

nice,sign,tabl

AUTHOR

Paul D. Hanna (pauldhanna(AT)juno.com), Sep 13 2002

EXTENSIONS

Error in one term corrected by Benoit Cloitre, Aug 20 2007

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified July 26 13:41 EDT 2008. Contains 142293 sequences.


AT&T Labs Research