OFFSET
1,1
COMMENTS
Unlike the results of applying the partial sum operator once (A007504), twice (A014148), or thrice (A014150) to primes, this sequence begins with 4 primes. The next prime in the sequence is a(26) = 367673.
Row 4 in A254858. - Reinhard Zumkeller, Feb 08 2015
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 1000 terms from Harvey P. Dale)
EXAMPLE
a(15) = 2 + 9 + 26 + 60 + 122 + 225 + 386 + 624 + 962 + 1429 + 2056 + 2880 + 3942 + 5285 + 6956 = 24964 = 2^2 x 79^2.
MAPLE
A007504 := proc(n) option remember; add( ithprime(i), i=1..n) ; end proc:
seq(A178138(n), n=1..80) ; # R. J. Mathar, Oct 19 2010
# Alternative:
b:= proc(n, k) option remember; `if`(n=0, 0,
`if`(k=0, ithprime(n), b(n-1, k)+b(n, k-1)))
end:
a:= n-> b(n, 4):
seq(a(n), n=1..39); # Alois P. Heinz, Mar 20 2026
MATHEMATICA
Nest[Accumulate[#]&, Prime[Range[40]], 4] (* Harvey P. Dale, Sep 25 2014 *)
PROG
(Haskell)
a178138 n = a178138_list !! (n-1)
a178138_list = (iterate (scanl1 (+)) a000040_list) !! 4
-- Reinhard Zumkeller, Feb 08 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jonathan Vos Post, May 20 2010
STATUS
approved