OFFSET
0,6
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (16,-101,298,-278,-656,1854,-676,-2581,2744,935,-2422,460,696,-288).
FORMULA
a(n) = A219836(n,n-4).
G.f.: 8 * x^5 * (1-3*x+5*x^2-29*x^3+63*x^4-46*x^5+33*x^6-36*x^7)/((1-x^2)^4 * (1-2*x)^3 * (1-3*x)^2 * (1-4*x)).
a(n) = 16*a(n-1) - 101*a(n-2) + 298*a(n-3) - 278*a(n-4) - 656*a(n-5) + 1854*a(n-6) - 676*a(n-7) - 2581*a(n-8) + 2744*a(n-9) + 935*a(n-10) - 2422*a(n-11) + 460*a(n-12) + 696*a(n-13) - 288*a(n-14) for n > 13.
EXAMPLE
For n=5, the a(5) = 8 derangements are: [2, 3, 4, 5, 1] (2<3, 3<4, 4<5), [2, 3, 5, 1, 4] (2<3, 3<5, 1<4), [2, 4, 5, 1, 3] (2<4, 4<5, 1<3), [2, 5, 1, 3, 4] (2<5, 1<3, 3<4), [3, 4, 5, 1, 2] (3<4, 4<5, 1<2), [3, 5, 1, 2, 4] (3<5, 1<2, 2<4), [4, 5, 1, 2, 3] (4<5, 1<2, 2<3), and [5, 1, 2, 3, 4] (1<2, 2<3, 3<4).
PROG
(PARI) my(N=30, x='x+O('x^N)); concat([0, 0, 0, 0, 0], Vec(8*x^5*(1-3*x+5*x^2-29*x^3+63*x^4-46*x^5+33*x^6-36*x^7)/((1-x^2)^4*(1-2*x)^3*(1-3*x)^2*(1-4*x))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Apr 25 2026
STATUS
approved
