login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo

Please make a donation to keep the OEIS running. We are now in our 56th year. In the past year we added 10000 new sequences and reached almost 9000 citations (which often say "discovered thanks to the OEIS").
Other ways to donate

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A290929 p-INVERT of the positive integers, where p(S) = (1 - S)(1 - S^2). 2
1, 4, 13, 39, 114, 330, 948, 2703, 7655, 21554, 60389, 168468, 468199, 1296826, 3581185, 9862749, 27096216, 74277342, 203200986, 554869701, 1512575195, 4116813032, 11188568267, 30367047720, 82316338381, 222875101936, 602784607477, 1628612506131 (list; graph; refs; listen; history; text; internal format)
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 A290890 for a guide to related sequences.

LINKS

Clark Kimberling, Table of n, a(n) for n = 0..1000

Index entries for linear recurrences with constant coefficients, signature (7, -18, 23, -18, 7, -1)

FORMULA

a(n) = 7*a(n-1) - 18*a(n-2) + 23*a(n-3) - 18*a(n-4) + 7*a(n-5) - a(n-6).

G.f.: (1 - 3*x + 3*x^2 - 3*x^3 + x^4) / ((1 - 3*x + x^2)^2*(1 - x + x^2)). - Colin Barker, Aug 19 2017

MATHEMATICA

z = 60; s = x/(1 - x)^2; p = (1 - s)(1 - s^2);

Drop[CoefficientList[Series[s, {x, 0, z}], x], 1] (* A000027 *)

Drop[CoefficientList[Series[1/p, {x, 0, z}], x], 1]  (* A290929 *)

LinearRecurrence[{7, -18, 23, -18, 7, -1}, {1, 4, 13, 39, 114, 330}, 40] (* Vincenzo Librandi, Aug 20 2017 *)

PROG

(PARI) Vec((1 - 3*x + 3*x^2 - 3*x^3 + x^4) / ((1 - 3*x + x^2)^2*(1 - x + x^2)) + O(x^30)) \\ Colin Barker, Aug 19 2017

(MAGMA) I:=[1, 4, 13, 39, 114, 330]; [n le 6 select I[n] else 7*Self(n-1)-18*Self(n-2)+23*Self(n-3)-18*Self(n-4)+7*Self(n-5)-Self(n-6): n in [1..40]]; // Vincenzo Librandi, Aug 20 2017

CROSSREFS

Cf. A000027, A033453, A290890.

Sequence in context: A266429 A105693 A215404 * A121192 A133409 A000746

Adjacent sequences:  A290926 A290927 A290928 * A290930 A290931 A290932

KEYWORD

nonn,easy

AUTHOR

Clark Kimberling, Aug 19 2017

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified December 3 20:35 EST 2020. Contains 338920 sequences. (Running on oeis4.)