|
|
A061890
|
|
Squares that are the sum of initial primes.
|
|
9
|
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
|
|
EXAMPLE
|
100 = 10^2 = 2 + 3 + 5 + 7 + 11 + 13 + 17 + 19 + 23, so 100 is in the sequence.
|
|
MAPLE
|
N:= 10^7: # to use primes up to N
P:= select(isprime, [2, seq(2*i+1, i=1..floor(N/2))]):
S:= ListTools:-PartialSums(P):
|
|
MATHEMATICA
|
s[n_] := Sum[Prime[i], {i, 1, n}]; t := Table[s[n], {n, 20000}]; Select[t, IntegerQ[Sqrt[#]] &] (* Carlos Eduardo Olivieri, Feb 24 2015 *)
|
|
PROG
|
(PARI) lista() = {s = 0; forprime(p=2, , s += p; if (issquare(s), print1(s, ", ")); ); } \\ Michel Marcus, Mar 10 2015
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|