OFFSET
0,3
COMMENTS
Two n-permutations are randomly selected from S_n with replacement. a(n)/(n!)^2 is the probability that they will have the same number of inversions. - Geoffrey Critzer, May 15 2010
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..100
Eric Weisstein's World of Mathematics, q-Factorial.
FORMULA
a(n) ~ 6 * sqrt(Pi) * n^(2*n - 1/2) / exp(2*n). - Vaclav Kotesovec, Oct 22 2020
a(n) = Sum_{k>=0} A008302(n,k)^2. - R. J. Mathar, Jan 06 2022
EXAMPLE
Definition of q-factorial of n:
faq(n) = Product_{k=1..n} (1-q^k)/(1-q) for n>0, with faq(0)=1;
faq(4) = 1*(1 + q)*(1 + q + q^2)*(1 + q + q^2 + q^3) = 1 + 3*q + 5*q^2 + 6*q^3 + 5*q^4 + 3*q^5 + q^6;
then a(n) is the sum of squared coefficients of q:
a(4) = 1^2 + 3^2 + 5^2 + 6^2 + 5^2 + 3^2 + 1^2 = 106.
MATHEMATICA
Table[Total[ CoefficientList[Expand[Product[Sum[x^i, {i, 0, m}], {m, 1, n - 1}]], x]^2], {n, 0, 15}] (* Geoffrey Critzer, May 15 2010 *)
PROG
(PARI) {a(n)=local(faq_n=if(n==0, 1, prod(k=1, n, (1-q^k)/(1-q)))); sum(k=0, n*(n-1)/2, polcoeff(faq_n, k, q)^2)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 25 2007
STATUS
approved