Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A090582
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A090582 Numerator Q(m,n) of probability P(m,n)=Q(m,n)/n^m to see each card at least once if m>=n cards are drawn with replacement from a deck of n cards, written in a two-dimensional array read by antidiagonals with Q(m,m) as first row and Q(m,1) as first column. +0
6
1, 2, 1, 6, 6, 1, 24, 36, 14, 1, 120, 240, 150, 30, 1, 720, 1800, 1560, 540, 62, 1, 5040, 15120, 16800, 8400, 1806, 126, 1, 40320, 141120, 191520, 126000, 40824, 5796, 254, 1, 362880, 1451520, 2328480, 1905120, 834120, 186480, 18150, 510, 1, 3628800 (list; table; graph; listen)
OFFSET

1,2

COMMENT

The sequence is given as a matrix with the first row containing the cases #draws=size_of_deck. The second row contains #draws=1+size_of_deck. If "mn" indicates m cards drawn from a deck with n cards then the locations in the matrix are:

11 22 33 44 55 66 77 ...

21 32 43 54 65 76 87 ...

31 42 53 64 75 86 97 ...

41 52 63 74 85 .. .. ...

read by antidiagonals ->:

11, 22, 21, 33, 32, 31, 44, 43, 42, 41, 55, 54, 53, 52, ....

The probabilities are given by Q(m,n)/n^m:

.(m,n):.....11..22..21..33..32..31..44..43..42..41...55...54..53..52..51

.....Q:......1...2...1...6...6...1..24..36..14...1..120..240.150..30...1

...n^m:......1...4...1..27...8...1.256..81..16...1.3125.1024.243..32...1

%.Success:.100..50.100..22..75.100...9..44..88.100....4...23..62..94.100

P(n,n) = n!/n^n which can be approximated by sqrt(pi*(2n+1/3))/e^n (Gosper's approximation to n!).

Triangle T(n,k), 1<=k<=n, read by rows given by [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] DELTA [[0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, ...] where DELTA is the operator defined in A084938 . - Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Nov 10 2006

REFERENCES

See A073593.

LINKS

Kevin Buhr, Michael Press and DMB, Collecting a deck of cards on the street. Thread in NG sci.math.

Hugo Pfoertner, "Hit all" probabilities: Program and results.

FORMULA

Q(m, n) = sum_(k=0..n-1) (-1)^k * C(n, k) * (n-k)^m where C(n, k) is the binomial coefficient "n choose k".

Generated by Stirling numbers of the second kind multiplied by a factorial. - Victor Adamchik, Oct 05 2005

Formulae from Tom Copeland (tcjpn(AT)msn.com), Oct 07 2008: (Start)

G(x,t) = 1/ {1 + [1-exp(x t)]/t} = 1 + 1 x + (2 + t) x^2/2! + (6 + 6t + t^2) x^3/3! + ...

gives row polynomials of A090582-- reverse f-polynomials for the permutohedra (see A019538).

G(x,t-1) = 1 + 1 x + (1 + t) x^2 / 2! + (1 + 4t + t^2) x^3 / 3! + ...

gives row polynomials for A008292, the h-polynomials for permutohedra.

G[(t+1)x,-1/(t+1)] = 1 + (1+ t) x + (1 + 3t + 2 t^2) x^2 / 2! + ...

gives row polynomials of A028246. (End)

MATHEMATICA

In[1]:= Table[Table[k! StirlingS2[n, k], {k, n, 1, -1}], {n, 1, 6}] - Victor Adamchik, Oct 05 2005

PROGRAM

FORTRAN program given at Pfoertner link.

CROSSREFS

Cf. A073593 first m>=n giving at least 50% probability, A085813 ditto for 95%, A055775 n^n/n!, A090583 Gosper's approximation to n!.

Reflected version of A019538.

Sequence in context: A105278 A008297 A048999 this_sequence A079641 A075181 A052121

Adjacent sequences: A090579 A090580 A090581 this_sequence A090583 A090584 A090585

KEYWORD

frac,nonn,tabl

AUTHOR

Hugo Pfoertner (hugo(AT)pfoertner.org), Jan 11 2004

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 December 10 12:37 EST 2009. Contains 170569 sequences.


AT&T Labs Research