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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A060891 a(n) = n^6 - n^3 + 1. 5
1, 1, 57, 703, 4033, 15501, 46441, 117307, 261633, 530713, 999001, 1770231, 2984257, 4824613, 7526793, 11387251, 16773121, 24132657, 34006393, 47039023, 63992001, 85756861, 113369257, 148023723, 191089153, 244125001, 308898201 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

COMMENTS

a(n) = Phi_18(n) where Phi_k(x) is the k-th cyclotomic polynomial.

LINKS

Harry J. Smith, Table of n, a(n) for n = 0..1000

Index to values of cyclotomic polynomials of integer argument

Index entries for linear recurrences with constant coefficients, signature (7, -21, 35, -35, 21, -7, 1).

FORMULA

G.f.: (1-6*x+71*x^2+290*x^3+309*x^4+52*x^5+3*x^6)/(1-x)^7. [Colin Barker, Apr 22 2012]

MAPLE

with (combinat):seq(fibonacci(3, n^3)-n^3, n=0..30); # Zerinvary Lajos, May 25 2008

# Alternative

A060891 := proc(n)

        numtheory[cyclotomic](18, n) ;

end proc:

seq(A060891(n), n=0..20) ; # R. J. Mathar, Feb 11 2014

PROG

(PARI) { for (n=0, 1000, write("b060891.txt", n, " ", n^6 - n^3 + 1); ) } \\ Harry J. Smith, Jul 14 2009

(PARI) a(n) = polcyclo(18, n); \\ Michel Marcus, Dec 16 2017

CROSSREFS

Sequence in context: A195606 A232875 A084220 * A241967 A098995 A210156

Adjacent sequences:  A060888 A060889 A060890 * A060892 A060893 A060894

KEYWORD

nonn,easy

AUTHOR

N. J. A. Sloane, May 05 2001

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 26 23:05 EST 2020. Contains 331289 sequences. (Running on oeis4.)