OFFSET
1,1
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..385
FORMULA
a(n) = 1 + Sum_{i=1..n} Sum_{j=1..n} i^j.
a(n) = n + 1 + Sum_{j=2..n} j*(j^n - 1)/(j-1).
a(n) = A086787(n) + 1.
MAPLE
seq( n+1+add(j*(j^n-1)/(j-1), j=2..n), n=1..30); # G. C. Greubel, Dec 25 2019
MATHEMATICA
Table[Sum[i^j, {i, 1, n}, {j, 1, n}]+1, {n, 1, 20}]
PROG
(PARI) vector(30, n, n+1 + sum(j=2, n, j*(j^n-1)/(j-1)) ) \\ G. C. Greubel, Dec 25 2019
(Magma) [0] cat [n+1 + (&+[j*(j^n-1)/(j-1): j in [2..n]]): n in [2..30]]; // G. C. Greubel, Dec 25 2019
(Sage) [n+1 + sum(j*(j^n-1)/(j-1) for j in (2..n)) for n in (1..30)] # G. C. Greubel, Dec 25 2019
(GAP) List([1..30], n-> n+1 + Sum([2..n], j-> j*(j^n-1)/(j-1)) ); # G. C. Greubel, Dec 25 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Dec 14 2006
EXTENSIONS
Edited by Max Alekseyev, Jan 29 2012
STATUS
approved