A338818 Primes p such that the sum of squares of primes < p is divisible by p.

2, 13, 59, 118259, 182603

Offset

1,1

Comments

Primes prime(n) such that A024450(n) is divisible by prime(n).

Links

Table of n, a(n) for n=1..5.

Example

a(3) = 59 = prime(17) is in the sequence because A024450(17)=16756 is divisible by 59.

Maple

P:= select(isprime, [2, seq(i, i=3..10^6, 2)]):
S:= ListTools:-PartialSums(map(`^`, P, 2)):
map(t -> P[t], select(t -> S[t] mod P[t]=0, [$1..nops(P)]);

Prog

(PARI) lista(nn) = {my(s=0, list=List()); forprime(p=2, nn, s += p^2; if (!(s % p), listput(list, p)); ); Vec(list); } \\ Michel Marcus, Nov 11 2020

Crossrefs

Cf. A007506, A024450.

Sequence in context: A042061 A229736 A187560 * A290721 A354942 A205532

Adjacent sequences: A338815 A338816 A338817 * A338819 A338820 A338821

Keyword

nonn,more

Author

J. M. Bergot and Robert Israel, Nov 10 2020

Status

approved