login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A078630 Numerators of coefficients of asymptotic expansion of probability p(n) (see A002816) in powers of 1/n. 4

%I

%S 1,-4,0,20,58,796,7858,40324,140194,2444744,40680494,-7117319032,

%T -149539443124,-223750776484,-4960419494993024,-46146161037854692,

%U -689434674121075448,-132496988938839119444,-9686633414582239854958,-442788087926096759821484

%N Numerators of coefficients of asymptotic expansion of probability p(n) (see A002816) in powers of 1/n.

%H Vaclav Kotesovec, <a href="/A078630/b078630.txt">Table of n, a(n) for n = 0..50</a>

%H B. Aspvall and F. M. Liang, <a href="http://www-db.stanford.edu/TR/cstr8x.html">The dinner table problem</a>, Technical Report CS-TR-80-829, Computer Science Department, Stanford, California, 1980.

%e p(n) = exp(-2)*(1 - 4/n + 20/(3n^3) + 58/(3n^4) + ...).

%t t = 15;

%t y[n_]:=(1+Sum[Subscript[p,k]/n^k,{k,1,t}]);

%t mul=1;start=9; If[t>9,mul=n^(t-9);start=t];

%t w=Apart[Expand[mul*Simplify[

%t y[n]*n*(n-1)*(n-2)*(n-3)*(n-4)*(n-5)*(n-6)*(n-7)*(n-8)*(n-9)*(n-10)

%t -((3*n-30)*y[n-11]

%t +(6*n-45)*y[n-10]*(n-10)

%t +(5*n+18)*y[n-9]*(n-9)*(n-10)

%t -(8*n-139)*y[n-8]*(n-8)*(n-9)*(n-10)

%t -(26*n-204)*y[n-7]*(n-7)*(n-8)*(n-9)*(n-10)

%t -(4*n-30)*y[n-6]*(n-6)*(n-7)*(n-8)*(n-9)*(n-10)

%t +(26*n-148)*y[n-5]*(n-5)*(n-6)*(n-7)*(n-8)*(n-9)*(n-10)

%t +(8*n-74)*y[n-4]*(n-4)*(n-5)*(n-6)*(n-7)*(n-8)*(n-9)*(n-10)

%t -(9*n-18)*y[n-3]*(n-3)*(n-4)*(n-5)*(n-6)*(n-7)*(n-8)*(n-9)*(n-10)

%t -(2*n-15)*y[n-2]*(n-2)*(n-3)*(n-4)*(n-5)*(n-6)*(n-7)*(n-8)*(n-9)*(n-10)

%t +(n+2)*y[n-1]*(n-1)*(n-2)*(n-3)*(n-4)*(n-5)*(n-6)*(n-7)*(n-8)*(n-9)*(n-10))],n],n];

%t sol=Solve[Table[Coefficient[w,n,j]==0,{j,start,start-t+1,-1}]];

%t asympt=y[n]/.sol[[1]];

%t Table[Numerator[Coefficient[asympt,n,-j]],{j,0,t}] (* _Vaclav Kotesovec_, Apr 06 2012 *)

%Y Cf. A078631, A089222.

%K sign

%O 0,2

%A _N. J. A. Sloane_, Dec 13 2002

%E Terms a(5)-a(19) from Vaclav Kotesovec, Apr 06 2012 (terms a(5)-a(7) were wrong, see A089222 for more information)

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 13 22:55 EDT 2020. Contains 336473 sequences. (Running on oeis4.)