login
This site is supported by donations to The OEIS Foundation.

 

Logo

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A160828 a(n) = 4*n^4 + 24*n^3 + 84*n^2 + 144*n + 98. 1
98, 354, 978, 2258, 4578, 8418, 14354, 23058, 35298, 51938, 73938, 102354, 138338, 183138, 238098, 304658, 384354, 478818, 589778, 719058, 868578, 1040354, 1236498, 1459218, 1710818, 1993698, 2310354, 2663378, 3055458, 3489378, 3968018, 4494354 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,1

COMMENTS

Sums of 4 consecutive fourth powers.

Subsequence of A217844. - Michel Marcus, Jun 30 2013

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 0..10000

Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5,1).

FORMULA

a(n) = Sum_{i=0..3} A000583(n+i) = Sum_{j=n..n+3} j^4 = A160827(n) + (n+3)^4.

G.f.: (18*x^4 - 72*x^3 + 188*x^2 - 136*x + 98)/(1-x)^5. - Maksym Voznyy (voznyy(AT)mail.ru), Aug 14 2009

a(0)=98, a(1)=354, a(2)=978, a(3)=2258, a(4)=4578, a(n)=5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Harvey P. Dale, Mar 25 2012

E.g.f.: 2*(49 + 128*x + 92*x^2 + 24*x^3 + 2*x^4)*exp(x). - G. C. Greubel, Apr 30 2018

MAPLE

A000583 := proc(n) n^4 ; end: A160828 := proc(n) add(A000583(i), i=n..n+3) ; end: seq(A160828(n), n=0..40) ; # R. J. Mathar, May 29 2009

MATHEMATICA

Table[4n^4+24n^3+84n^2+144n+98, {n, 0, 40}] (* or *) LinearRecurrence[ {5, -10, 10, -5, 1}, {98, 354, 978, 2258, 4578}, 40] (* Harvey P. Dale, Mar 25 2012 *)

CoefficientList[Series[(18*x^4 -72*x^3 +188*x^2 -136*x +98)/(1-x)^5, {x, 0, 50}], x] (* G. C. Greubel, Apr 30 2018 *)

PROG

(MAGMA) [4*n^4+24*n^3+84*n^2+144*n+98: n in [0..40]]; // Vincenzo Librandi, Aug 27 2011

(Python)

A160828_list, m = [], [96, 0, 80, 80, 98]

for _ in range(10**2):

    A160828_list.append(m[-1])

    for i in range(4):

        m[i+1] += m[i] # Chai Wah Wu, Jan 23 2016

(PARI) x='x+O('x^50); Vec((18*x^4 -72*x^3 +188*x^2 -136*x +98)/(1-x)^5) \\ G. C. Greubel, Apr 30 2018

CROSSREFS

Sequence in context: A202371 A195751 A072607 * A158129 A071319 A263040

Adjacent sequences:  A160825 A160826 A160827 * A160829 A160830 A160831

KEYWORD

nonn,easy

AUTHOR

Al Hakanson (hawkuu(AT)gmail.com), May 27 2009

EXTENSIONS

Edited and corrected by R. J. Mathar, May 29 2009

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 January 18 23:05 EST 2019. Contains 319282 sequences. (Running on oeis4.)