OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (9,-36,84,-126,126,-84,36,-9,1).
FORMULA
a(n) = (n^2+n+1) * (n^6+n^3+1) and so is never prime. - Jonathan Vos Post, Dec 21 2012
G.f.: (x^8 + 162*x^7 + 3418*x^6 + 14212*x^5 + 16578*x^4 + 5482*x^3 + 466*x^2 + 1)/(1-x)^9. - Colin Barker, Nov 05 2012, edited by M. F. Hasler, Dec 31 2012
a(n) = (n^9-1)/(n-1) with a(1) = 9. - L. Edson Jeffery and M. F. Hasler, Dec 30 2012
MATHEMATICA
1 + Sum[Range[0, 30]^j, {j, 1, 8}] (* G. C. Greubel, Feb 13 2018 *)
PROG
(PARI) a(n)=n^8+n^7+n^6+n^5+n^4+n^3+n^2+n+1 \\ Charles R Greathouse IV, Oct 07 2015
(Magma) [(&+[n^j: j in [0..8]]): n in [0..30]]; // G. C. Greubel, Feb 13 2018
(Sage) [sum(n^j for j in (0..8)) for n in (0..30)] # G. C. Greubel, Apr 14 2019
CROSSREFS
Cf. similar sequences of the type a(n) = Sum_{j=0..m} n^j: A000027 (m=1), A002061 (m=2), A053698 (m=3), A053699 (m=4), A053700 (m=5), A053716 (m=6), A053717 (m=7), this sequence (m=8), A103623 (m=9), A060885 (m=10), A105067 (m=11), A060887 (m=12), A104376 (m=13), A104682 (m=14), A105312 (m=15), A269442 (m=16), A269446 (m=18).
KEYWORD
nonn,easy
AUTHOR
Douglas Winston (douglas.winston(AT)srupc.com), Mar 01 2005
EXTENSIONS
Offset corrected by N. J. A. Sloane, Dec 30 2012
STATUS
approved