OFFSET
0,4
COMMENTS
In other words, there are no two objects which the permutation leaves fixed and no two objects that the permutation swaps.
The limit as n -> infinity of a(n)/n! = 2/exp(3/2) or approximately 0.4462603203. - Les Reid, Jun 04 2012
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..450
Hannah Jackson, Kathryn Nyman and Les Reid, Properties of generalized derangement graphs, Involve, Vol. 6 (2013), No. 1, 25-33; DOI: 10.2140/involve.2013.6.25.
FORMULA
E.g.f.: (1+x)*exp(-x)*exp(-x^2/2)/(1-x).
a(n) = (n-1)*a(n-1) - a(n-2) + (n-2)*n*a(n-3) for n > 2, a(n) = (n+1)*(2-n)/2 for n < 3. - Alois P. Heinz, Feb 19 2019
EXAMPLE
a(3)=2 because we have 312 and 231.
MAPLE
g:=(1+x)*exp(-x)*exp(-(1/2)*x^2)/(1-x): gser:=series(g, x=0, 25): seq(factorial(n)*coeff(gser, x, n), n=0..20);
# second Maple program:
a:= proc(n) option remember; `if`(n<3, (n+1)*(2-n)/2,
(n-1)*a(n-1)-a(n-2)+(n-2)*n*a(n-3))
end:
seq(a(n), n=0..23); # Alois P. Heinz, Feb 19 2019
MATHEMATICA
With[{nn=20}, CoefficientList[Series[((1+x)Exp[-x]Exp[-x^2/2])/(1-x), {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Nov 17 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Jono Henshaw (jjono(AT)hotmail.com), Apr 22 2008, corrected Apr 30 2008
STATUS
approved