cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

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

This page as a plain text file.
%I A338818 #13 Jun 10 2025 21:42:34
%S A338818 2,13,59,118259,182603
%N A338818 Primes p such that the sum of squares of primes < p is divisible by p.
%C A338818 Primes prime(n) such that A024450(n) is divisible by prime(n).
%C A338818 a(6) > 8 * 10^9 if it exists. - _Robert Israel_, Jun 10 2025
%e A338818 a(3) = 59 = prime(17) is in the sequence because A024450(17)=16756 is divisible by 59.
%p A338818 P:= select(isprime, [2,seq(i,i=3..10^6,2)]):
%p A338818 S:= ListTools:-PartialSums(map(`^`,P,2)):
%p A338818 P[select(t -> S[t] mod P[t]=0, [$1..nops(P)])];
%o A338818 (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
%Y A338818 Cf. A007506, A024450.
%K A338818 nonn,more
%O A338818 1,1
%A A338818 _J. M. Bergot_ and _Robert Israel_, Nov 10 2020