OFFSET
2,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 2..450
Wikipedia, Permutation
EXAMPLE
a(2) = 1: (1)(2).
a(3) = 6: (1)(23), (23)(1), (2)(13), (13)(2), (3)(12), (12)(3).
a(4) = 19: (123)(4), (4)(123), (132)(4), (4)(132), (124)(3), (3)(124), (142)(3), (3)(142), (134)(2), (2)(134), (143)(2), (2)(143), (1)(234), (234)(1), (1)(243), (243)(1), (12)(34), (13)(24), (14)(23).
MAPLE
a:= n-> 2*add(binomial(n, k)*(k-1)!*(n-k-1)!, k=1..n/2)-
`if`(n::even, 3/2*binomial(n, n/2)*(n/2-1)!^2, 0):
seq(a(n), n=2..25);
# second Maple program:
a:= proc(n) option remember; `if`(n<5, [0, 1, 6, 19][n],
((2*n-1)*(n-1)*a(n-1)+(n-2)*(2*n^2-5*n-1)*a(n-2)
-(n-3)^2*((2*n^2-5*n+4)*a(n-3)+(n-4)^2*a(n-4)))/(2*n))
end:
seq(a(n), n=2..25);
MATHEMATICA
Table[(n-1)!*(2*HarmonicNumber[n] - (3 + (-1)^n)/n), {n, 2, 25}] (* Vaclav Kotesovec, Apr 29 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Apr 27 2017
STATUS
approved