%I #47 Apr 26 2024 05:54:23
%S 1,0,1,346,748521,3993445276,45131501617225,964363228180815366,
%T 35780355973270898382001,2158610844939711892526650456,
%U 201028342764877992289387752167601,27708893753238763155350683269145066450,5459844285803153226360263675364357481841881
%N Number of derangements of a multiset comprising 4 repeats of an n-element set.
%C Previous name was: Card-matching numbers (Dinner-Diner matching numbers).
%C A deck has n kinds of cards, 4 of each kind. The deck is shuffled and dealt in to n hands with 4 cards each. A match occurs for every card in the j-th hand of kind j. A(n) is the number of ways of achieving no matches. The probability of no matches is A(n)/((4n)!/4!^n).
%C Number of fixed-point-free permutations of n distinct letters (ABCD...), each of which appears 4 times: 1111, 11112222, 111122223333, 1111222233334444, etc. If there is only one letter of each type we get A000166 - _Zerinvary Lajos_, Nov 05 2006
%C a(n) is the maximal number of totally mixed Nash equilibria in games of n players, each with 5 pure options. [_Raimundas Vidunas_, Jan 22 2014]
%D F. N. David and D. E. Barton, Combinatorial Chance, Hafner, NY, 1962, Ch. 7 and Ch. 12.
%D R.D. McKelvey and A. McLennan, The maximal number of regular totally mixed Nash equilibria, J. Economic Theory, 72:411-425, 1997.
%D S. G. Penrice, Derangements, permanents and Christmas presents, The American Mathematical Monthly 98(1991), 617-620.
%D J. Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, pp. 174-178.
%D R. P. Stanley, Enumerative Combinatorics Volume I, Cambridge University Press, 1997, p. 71.
%H Michael De Vlieger, <a href="/A059074/b059074.txt">Table of n, a(n) for n = 0..100</a>
%H Shalosh B. Ekhad, Christoph Koutschan, and Doron Zeilberger, <a href="https://arxiv.org/abs/2101.10147">There are EXACTLY 1493804444499093354916284290188948031229880469556 Ways to Derange a Standard Deck of Cards (ignoring suits) [and many other such useful facts]</a>, arXiv:2101.10147 [math.CO], 2021.
%H Shalosh B. Ekhad, <a href="https://sites.math.rutgers.edu/~zeilberg/tokhniot/oMultiDer1.txt">Terms and recurrences for multiset derangements</a>.
%H S. Even and J. Gillis, <a href="https://doi.org/10.1017/S0305004100052154">Derangements and Laguerre polynomials</a>, Mathematical Proceedings of the Cambridge Philosophical Society, Volume 79, Issue 1, January 1976, pp. 135-143.
%H F. F. Knudsen and I. Skau, <a href="http://www.jstor.org/stable/2691467">On the Asymptotic Solution of a Card-Matching Problem</a>, Mathematics Magazine 69 (1996), 190-197.
%H Barbara H. Margolius, <a href="http://www.jstor.org/stable/3219303">The Dinner-Diner Matching Problem</a>, Mathematics Magazine, 76 (2003), 107-118.
%H Barbara H. Margolius, <a href="http://academic.csuohio.edu/bmargolius/homepage/dinner/dinner/cardentry.htm">Dinner-Diner Matching Probabilities</a>
%H Raimundas Vidunas, <a href="http://arxiv.org/abs/1401.5400">MacMahon's master theorem and totally mixed Nash equilibria</a>, arXiv preprint arXiv:1401.5400 [math.CO], 2014.
%H Raimundas Vidunas, <a href="https://doi.org/10.1007/s00026-017-0344-2">Counting derangements and Nash equilibria</a> Ann. Comb. 21, No. 1, 131-152 (2017).
%H <a href="/index/Ca#cardmatch">Index entries for sequences related to card matching</a>
%F G.f.: sum(coeff(R(x, n, k), x, j)*(t-1)^j*(n*k-j)!, j=0..n*k) where n is the number of kinds of cards, k is the number of cards of each kind (3 in this case) and R(x, n, k) is the rook polynomial given by R(x, n, k)=(k!^2*sum(x^j/((k-j)!^2*j!))^n (see Stanley or Riordan). coeff(R(x, n, k), x, j) indicates the coefficient for x^j of the rook polynomial.
%F From _Jeremy Tan_, Apr 25 2024: (Start)
%F a(n) = Integral_{x=0..oo} exp(-x)*L_4(x)^n dx, where L_n(x) is the Laguerre polynomial of degree n (Even and Gillis).
%F D-finite with recurrence 3*(128*n^3 - 560*n^2 + 840*n - 537)*a(n) - n*(4096*n^6 - 24064*n^5 + 62720*n^4 - 92992*n^3 + 75248*n^2 - 38670*n + 4179)*a(n-1) - 2*n*(n-1)*(18432*n^5 - 99072*n^4 + 197120*n^3 - 191776*n^2 + 144568*n - 92531)*a(n-2) + 48*n*(n-1)*(n-2)*(768*n^4 - 2976*n^3 + 3104*n^2 - 2438*n + 1583)*a(n-3) + 288*n*(n-1)*(n-2)*(n-3)*(128*n^3 - 176*n^2 + 104*n - 129)*a(n-4) = 8192*n^6 - 28672*n^5 + 23680*n^4 - 7904*n^3 + 1416*n^2 + 14382*n - 1611 (Ekhad).
%F a(n) ~ A014608(n)/exp(4) ~ n^(4*n)*(32/3)^n*sqrt(8*Pi*n)/exp(4*n+4). (End)
%e There are 346 ways of achieving zero matches when there are 4 cards of each kind and 3 kinds of card so A(3)=346.
%p p := (x,k)->k!^2*sum(x^j/((k-j)!^2*j!),j=0..k); R := (x,n,k)->p(x,k)^n; f := (t,n,k)->sum(coeff(R(x,n,k),x,j)*(t-1)^j*(n*k-j)!,j=0..n*k); seq(f(0,n,4)/4!^n,n=0..18);
%t p[x_, k_] := k!^2*Sum[x^j/((k - j)!^2*j!), {j, 0, k}]; r[x_, n_, k_] := p[x, k]^n; f[t_, n_, k_] := Sum[Coefficient[r[x, n, k], x, j]*(t - 1)^j*(n*k - j)!, {j, 0, n*k}]; Table[f[0, n, 4]/4!^n, {n, 0, 18}] // Flatten (* _Jean-François Alcover_, Oct 21 2013, after Maple *)
%t Table[Integrate[Exp[-x] LaguerreL[4, x]^n, {x, 0, Infinity}], {n, 0, 16}] (* _Jeremy Tan_, Apr 25 2024 *)
%t rec = 3*(128*n^3 - 560*n^2 + 840*n - 537)*a[n] - n*(4096*n^6 - 24064*n^5 + 62720*n^4 - 92992*n^3 + 75248*n^2 - 38670*n + 4179)*a[n-1] - 2*n*(n-1)*(18432*n^5 - 99072*n^4 + 197120*n^3 - 191776*n^2 + 144568*n - 92531)*a[n-2] + 48*n*(n-1)*(n-2)*(768*n^4 - 2976*n^3 + 3104*n^2 - 2438*n + 1583)*a[n-3] + 288*n*(n-1)*(n-2)*(n-3)*(128*n^3 - 176*n^2 + 104*n - 129)*a[n-4] == 8192*n^6 - 28672*n^5 + 23680*n^4 - 7904*n^3 + 1416*n^2 + 14382*n - 1611;
%t RecurrenceTable[{rec, a[0] == 1, a[1] == 0, a[2] == 1, a[3] == 346}, a, {n, 0, 16}] (* _Jeremy Tan_, Apr 25 2024 *)
%o (Python)
%o def A059074(n):
%o l = [1, 0, 1, 346]
%o for k in range(4, n+1):
%o num = (((((8192*k-28672)*k+23680)*k-7904)*k+1416)*k+14382)*k-1611 \
%o + k*((((((4096*k-24064)*k+62720)*k-92992)*k+75248)*k-38670)*k+4179)*l[-1] \
%o + 2*k*(k-1)*(((((18432*k-99072)*k+197120)*k-191776)*k+144568)*k-92531)*l[-2] \
%o - 48*k*(k-1)*(k-2)*((((768*k-2976)*k+3104)*k-2438)*k+1583)*l[-3] \
%o - 288*k*(k-1)*(k-2)*(k-3)*(((128*k-176)*k+104)*k-129)*l[-4]
%o r = num // (3*(((128*k-560)*k+840)*k-537))
%o l.append(r)
%o return l[n] # _Jeremy Tan_, Apr 25 2024
%Y Cf. A000166, A008290, A059056-A059071, A371252.
%K nonn,nice
%O 0,4
%A Barbara Haas Margolius (margolius(AT)math.csuohio.edu)
%E Name changed by _Jeremy Tan_, Apr 25 2024