G.f.: (1 -x +x^2 -x^3 +x^4)/((1+x^2)*(1+x^4)*(1+x)^2*(1-x)^3). - R. J. Mathar, Dec 18 2014
a(n) = floor((2*n+7+5*(-1)^n)^2/128). - Hoang Xuan Thanh, Jun 20 2025
MAPLE
seq(coeff(series((1+x^5)/((1-x^2)^2*(1-x^8)), x, n+1), x, n), n = 0 .. 70); # modified by G. C. Greubel, Jul 30 2019
MATHEMATICA
CoefficientList[Series[(1+x^5)/((1-x^2)^2*(1-x^8)), {x, 0, 70}], x] (* G. C. Greubel, Jul 30 2019 *)
PROG
(PARI) my(x='x+O('x^70)); Vec((1+x^5)/((1-x^2)^2*(1-x^8))) \\ G. C. Greubel, Jul 30 2019
(Magma) R<x>:=PowerSeriesRing(Integers(), 70); Coefficients(R!( (1+x^5)/((1-x^2)^2*(1-x^8)) )); // G. C. Greubel, Jul 30 2019
(SageMath) ((1+x^5)/((1-x^2)^2*(1-x^8))).series(x, 70).coefficients(x, sparse=False) # G. C. Greubel, Jul 30 2019
(GAP) a:=[1, 0, 2, 0, 3, 1, 4, 2, 6, 3, 8];; for n in [12..70] do a[n]:=a[n-1]+a[n-2]-a[n-3] +a[n-8]-a[n-9]-a[n-10]+a[n-11]; od; a; # G. C. Greubel, Jul 30 2019
