OFFSET
0,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Sam Miner, Enumeration of several two-by-four classes, arXiv:1610.01908 [math.CO], 2016.
FORMULA
G.f.: (3 - 22*x + 54*x^2 - 54*x^3 + 25*x^4 - 4*x^5 - (1 - 6*x + 14*x^2 - 16*x^3 + 5*x^4)*sqrt(1 - 4*x))/(2 * sqrt(1 - 4*x) * (1 - 3*x + x^2)^2).
a(n) ~ 2^(2*n+4)/(25*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 07 2016
EXAMPLE
For n = 4, a(n)=22, since all permutations of length 4 except 4123 and 1324 are included.
MATHEMATICA
CoefficientList[Series[(3 -22*x +54*x^2 -54*x^3 +25*x^4 -4*x^5 -(1 -6*x + 14*x^2 -16*x^3 +5*x^4)*Sqrt[1 -4*x])/(2*Sqrt[1 -4*x]*(1 -3*x + x^2)^2), {x, 0, 40}], x] (* G. C. Greubel, Apr 09 2017 *)
PROG
(PARI) x='x+O('x^40); Vec((3-22*x+54*x^2-54*x^3+25*x^4-4*x^5-(1-6*x + 14*x^2 -16*x^3 +5*x^4)*sqrt(1-4*x))/(2*sqrt(1-4*x)*(1-3*x+x^2)^2)) \\ G. C. Greubel, Apr 09 2017
(Magma) m:=40; R<x>:=PowerSeriesRing(Rationals(), m); Coefficients(R!((3 -22*x +54*x^2 -54*x^3 +25*x^4 -4*x^5 -(1 -6*x + 14*x^2 -16*x^3 +5*x^4)*Sqrt(1 -4*x))/(2*Sqrt(1 -4*x)*(1 -3*x + x^2)^2))); // G. C. Greubel, Oct 22 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Sam Miner, Oct 05 2016
STATUS
approved