A306511 Number of permutations p of [n] having at least one index i with |p(i)-i| = 1.

0, 0, 1, 4, 19, 99, 603, 4248, 34115, 307875, 3085203, 33993870, 408482695, 5316309607, 74499953255, 1118421967520, 17907571955927, 304619809031127, 5486197279305911, 104289196264058030, 2086706157642260387, 43838287730208552691, 964790364323910060691

Offset

0,4

Links

Alois P. Heinz, Table of n, a(n) for n = 0..450

Wikipedia, Permutation

Formula

a(n) = n! - A078480(n).

Maple

a:= proc(n) option remember; `if`(n<5, [0$2, 1, 4, 19][n+1],
      (2*(n^3-8*n^2+20*n-14)*a(n-1)-(n-4)*(n-1)*(n^2-5*n+7)*
       a(n-2)-(n-2)*(n^2-7*n+13)*a(n-3)+(n^4-12*n^3+53*n^2
       -102*n+71)*a(n-4)+(n-4)*(n^2-5*n+7)*a(n-5))/(n^2-7*n+13))
    end:
seq(a(n), n=0..23);

Mathematica

a[n_] := n! - Sum[Sum[(-1)^k (i-k)! Binomial[2i-k, k], {k, 0, i}],
     {i, 0, n}];
a /@ Range[0, 23] (* Jean-François Alcover, May 03 2021, after Vaclav Kotesovec in A078480 *)

Crossrefs

Column k=1 of A306506.

Cf. A000142, A078480, A180191.

Sequence in context: A211855 A327115 A370024 * A177249 A083882 A007564

Adjacent sequences: A306508 A306509 A306510 * A306512 A306513 A306514

Keyword

nonn

Author

Alois P. Heinz, Feb 20 2019

Status

approved