|
|
A027616
|
|
Number of permutations of n elements containing a 2-cycle.
|
|
4
|
|
|
0, 0, 1, 3, 9, 45, 285, 1995, 15855, 142695, 1427895, 15706845, 188471745, 2450132685, 34301992725, 514529890875, 8232476226975, 139952095858575, 2519137759913775, 47863617438361725, 957272348112505425, 20102719310362613925, 442259824841726816925, 10171975971359716789275
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,4
|
|
LINKS
|
|
|
FORMULA
|
E.g.f.: (1 - exp(-x^2/2)) / (1-x).
a(n) = n! * ( 1 - Sum_{k=0..floor(n/2)} (-1)^k / (2^k * k!) ).
a(n) + A000266(n) = n!. - Yuval Dekel (dekelyuval(AT)hotmail.com), Nov 09 2003
|
|
MAPLE
|
S:= series((1-exp(-x^2/2))/(1-x), x, 101):
|
|
MATHEMATICA
|
nn=30; Table[n!, {n, 0, nn}]-Range[0, nn]!CoefficientList[Series[Exp[-x^2/2]/(1-x), {x, 0, nn}], x] (* Geoffrey Critzer, Oct 20 2012 *)
|
|
PROG
|
(PARI)
a(n) = n! * (1 - sum(k=0, floor(n/2), (-1)^k / (2^k * k!) ) );
(PARI)
N=33; x='x+O('x^N);
v=Vec( 'a0 + serlaplace( (1-exp(-x^2/2))/(1-x) ) );
v[1]-='a0; v
(Magma)
A027616:= func< n | Factorial(n)*(1- (&+[(-1/2)^j/Factorial(j): j in [0..Floor(n/2)]]) ) >;
(SageMath)
def A027616(n): return factorial(n)*(1-sum((-1/2)^k/factorial(k) for k in (0..(n//2))))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Joe Keane (jgk(AT)jgk.org)
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|