A089041 - OEIS

%I #36 Apr 12 2019 10:20:13

%S 1,0,3,26,453,11844,439975,22056222,1436236809,117923229512,

%T 11921584264011,1455483251191650,211163237294447053,

%U 35913642489947449356,7077505637217289437423,1599980633296779087784934,411293643476907595937924625,119299057697083019137937718672

%N Inverse binomial transform of squares of factorial numbers.

%C a(n) enumerates (ordered) lists of n two-tuples such that all numbers from 1 to n appear as the first as well as the second tuple entry and the j-th list member is not the tuple (j,j), for every j=1,..,n. Called coincidence-free 2-tuple lists of length n. See the Charalambides reference for this combinatorial interpretation.

%D Ch. A. Charalambides, Enumerative Combinatorics, Chapman & Hall/CRC, Boca Raton, Florida, 2002, p. 187, Exercise 13.(a), for r=2.

%H Alois P. Heinz, <a href="/A089041/b089041.txt">Table of n, a(n) for n = 0..200</a>

%H Roland Bacher, <a href="http://www.combinatorics.org/ojs/index.php/eljc/article/view/v19i3p7">Counting Packings of Generic Subsets in Finite Groups</a>, Electr. J. Combinatorics, 19 (2012), #P7. - From _N. J. A. Sloane_, Feb 06 2013

%F G.f.: hypergeom([1, 1, 1], [], x/(1+x))/(1+x).

%F E.g.f.: exp(-x)* hypergeom([1, 1], [], x).

%F a(n) = n^2*a(n-1) + n*(n-1)*a(n-2) + (-1)^n. - _Vladeta Jovovic_, Jul 15 2004

%F a(n) = Sum_{j=0..n} ((-1)^(n-j))*binomial(n,j)*(j!)^2. See the Charalambides reference a(n)=B_{n,2}.

%F a(n) = (n-1)*(n+1)*a(n-1) + (n-1)*(2*n-1)*a(n-2) + (n-2)*(n-1)*a(n-3). - _Vaclav Kotesovec_, Aug 13 2013

%F a(n) ~ 2*Pi*exp(-2*n)*n^(2*n+1). - _Vaclav Kotesovec_, Aug 13 2013

%F G.f.: Sum_{k>=0} (k!)^2*x^k/(1 + x)^(k+1). - _Ilya Gutkovskiy_, Apr 12 2019

%e 2-tuple combinatorics: a(1)=0 because the only list of 2-tuples with numbers 1 is [(1,1)] and this is a coincidence for j=1.

%e 2-tuple combinatorics: the a(2)=3 coincidence free 2-tuple lists of length n=2 are [(1,2),(2,1)], [(2,1),(1,2)] and [(2,2),(1,1)]. The list [(1,1),(2,2)] has two coincidences (j=1 and j=2).

%p a:= proc(n) option remember;

%p        `if`(n<2, 1-n, n^2*a(n-1)+n*(n-1)*a(n-2)+(-1)^n)

%p     end:

%p seq(a(n), n=0..20);  # _Alois P. Heinz_, Jun 21 2013

%t Table[n!Sum[(-1)^k(n-k)!/k!,{k,0,n}],{n,0,15}] (* _Geoffrey Critzer_, Jun 17 2013 *)

%Y Cf. A001044, A046662 (binomial transform of squares of factorial numbers).

%Y (-1)^n times the polynomials in A099599 evaluated at -1.

%K easy,nonn

%O 0,3

%A _Vladeta Jovovic_, Dec 03 2003

%E Charalambides reference and comments with combinatorial examples from _Wolfdieter Lang_, Jan 21 2008