%I #15 Sep 08 2022 08:45:24
%S 0,3614400,1306371456000,2432274637386240000,
%T 15509750490368582860800000,265241692266421512138485760000000,
%U 10332925158674345473855915900600320000000,815905363532798455769292988741440720076800000000,119621339682330952236606797649198078512534126592000000000
%N If X_1,...,X_n is a partition of a 5n-set X into 5-blocks then a(n) is equal to the number of permutations f of X such that f(X_i) <> X_i, (i=1,...,n).
%H G. C. Greubel, <a href="/A116221/b116221.txt">Table of n, a(n) for n = 1..89</a>
%H Milan Janjic, <a href="http://www.pmfbl.org/janjic/">Enumerative Formulas for Some Functions on Finite Sets</a>
%F a(n) = Sum_{j=0..n} (-120)^j*binomial(n,j)*(5*n-5*k)!.
%p a:=n->sum((-120)^i*binomial(n,i)*(5*n-5*i)!,i=0..n).
%t Table[Sum[(-5!)^j*Binomial[n, j]*(5*n-5*j)!, {j,0,n}], {n,1,20}] (* _G. C. Greubel_, May 11 2019 *)
%o (PARI) {a(n) = sum(j=0,n, (-120)^j*binomial(n,j)*(5*(n-j))!)}; \\ _G. C. Greubel_, May 11 2019
%o (Magma) [(&+[(-120)^j*Binomial(n, j)*Factorial(5*n-5*j): j in [0..n]]): n in [1..20]]; // _G. C. Greubel_, May 11 2019
%o (Sage) [sum((-120)^j*binomial(n,j)*factorial(5*n-5*j) for j in (0..n)) for n in (1..20)] # _G. C. Greubel_, May 11 2019
%o (GAP) List([1..20], n-> Sum([0..n], j-> (-120)^j*Binomial(n,j)* Factorial(5*n-5*j))) # _G. C. Greubel_, May 11 2019
%Y Cf. A116218, A116219, A116220.
%K nonn
%O 1,2
%A _Milan Janjic_, Apr 09 2007
|