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”).

A292490
p-INVERT of the odd positive integers, where p(S) = 1 - S - 7 S^2.
1
1, 11, 68, 365, 2019, 11328, 63321, 353483, 1974124, 11026373, 61584323, 343956104, 1921047729, 10729356747, 59925127764, 334691142941, 1869302113507, 10440343236752, 58310941508105, 325675681470731, 1818949357172988, 10159115194159989, 56740239146359107
OFFSET
0,2
COMMENTS
Suppose s = (c(0), c(1), c(2), ...) is a sequence and p(S) is a polynomial. Let S(x) = c(0)*x + c(1)*x^2 + c(2)*x^3 + ... and T(x) = (-p(0) + 1/p(S(x)))/x. The p-INVERT of s is the sequence t(s) of coefficients in the Maclaurin series for T(x). Taking p(S) = 1 - S gives the "INVERT" transform of s, so that p-INVERT is a generalization of the "INVERT" transform (e.g., A033453).
See A292480 for a guide to related sequences.
FORMULA
G.f.: -(((1 + x) (1 + 5 x + 8 x^2))/(-1 + 5 x + 17 x^3 + 7 x^4)).
a(n) = 5*a(n-1) + 17*a(n-3) + 7*a(n-4) for n >= 5.
MATHEMATICA
z = 60; s = x (x + 1)/(1 - x)^2; p = 1 - s - 7 s^2;
Drop[CoefficientList[Series[s, {x, 0, z}], x], 1] (* A005408 *)
Drop[CoefficientList[Series[1/p, {x, 0, z}], x], 1] (* A292490 *)
PROG
(PARI) x='x+O('x^99); Vec(((1+x)*(1+5*x+8*x^2))/(1-5*x-17*x^3-7*x^4)) \\ Altug Alkan, Oct 03 2017
CROSSREFS
Sequence in context: A111931 A066433 A038741 * A169731 A212057 A293550
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Oct 03 2017
STATUS
approved