|
|
A061789
|
|
a(n) = Sum_{k=1..n} prime(k)^prime(k).
|
|
6
|
|
|
4, 31, 3156, 826699, 285312497310, 303160419089563, 827240565046755853740, 1979246896225360344977719, 20880469979094808259715377888286, 2567686153182091604540923022990731504371755
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
a(3) = 2^2 + 3^3 + 5^5 = 3156.
|
|
MATHEMATICA
|
|
|
PROG
|
(PARI) a=n=0; forprime (p=2, 383, write("b061789.txt", n++, " ", a+=p^p)) \\ Harry J. Smith, Jul 28 2009
(Python)
from itertools import accumulate, count, islice
from sympy import prime
def A061789_gen(): # generator of terms
yield from accumulate(((p:=prime(k))**p for k in count(1)))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|