login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A147605
Expansion of g.f.: 1/((1 - x^2 - x^3 - x^4 - x^5 - x^6 - x^7)*(1 + x + x^2 + x^3 + x^4 + x^5 - x^7)).
4
1, -1, 1, 0, 1, 1, 3, 3, 3, 10, 11, 21, 32, 52, 77, 128, 206, 320, 517, 817, 1297, 2060, 3290, 5220, 8298, 13205, 20980, 33360, 53056, 84366, 134114, 213263, 339086, 539123, 857240, 1363034, 2167197, 3445840, 5478951, 8711511, 13851359
OFFSET
0,7
LINKS
FORMULA
G.f.: 1/(1 + x - x^3 - 2*x^4 - 3*x^5 - 5*x^6 - 7*x^7 - 5*x^8 - 3*x^9 - 2*x^10 - x^11 + x^13 + x^14).
G.f.: -1/(x^7*f(x)*f(1/x)), where f(x) = -1 + x^2 + x^3 + x^4 + x^5 + x^6 + x^7. - G. C. Greubel, Oct 24 2022
MATHEMATICA
f[x_]= -1 +x^2 +x^3 +x^4 +x^5 +x^6 +x^7;
CoefficientList[Series[-1/(x^7*f[x]*f[1/x]), {x, 0, 50}], x] (* G. C. Greubel, Oct 24 2022 *)
PROG
(PARI) Vec(1/(1 +x -x^3 -2*x^4 -3*x^5 -5*x^6 -7*x^7 -5*x^8 -3*x^9 -2*x^10 -x^11 + x^13 +x^14) + O(x^40)) \\ Jinyuan Wang, Mar 10 2020
(Magma) R<x>:=PowerSeriesRing(Integers(), 50); Coefficients(R!( 1/((1-x^2-x^3-x^4-x^5-x^6-x^7)*(1+x+x^2+x^3+x^4+x^5-x^7)) )); // G. C. Greubel, Oct 24 2022
(SageMath)
def A147605_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( 1/((1-x^2-x^3-x^4-x^5-x^6-x^7)*(1+x+x^2+x^3+x^4+x^5-x^7)) ).list()
A147605_list(50) # G. C. Greubel, Oct 24 2022
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Roger L. Bagula, Nov 08 2008
EXTENSIONS
Definition corrected by N. J. A. Sloane, Nov 09 2008
STATUS
approved