OFFSET
1,1
COMMENTS
Corresponding primes A122209(a(n)) = A024450(a(n)^2) are listed in A122210(n) = {239087, 29194283, 13459558559, 2330212120559, ...}. All a(n) are of the form 6*m, where m = {1, 2, 5, 11, 26, 30, 38, 56, 61, 93, 125, 140, 149, 163, 233, 235, 251, 260, 263, 277, 299, 300, ...}. Because A122209(2*m-1) is an even number and A122209(3*m-1) == A122209(3*m+1) == 0 (mod 3) for m >= 1. [Edited by Jinyuan Wang, Mar 23 2020]
MATHEMATICA
s=0; Do[p=Prime[n]; k=Sqrt[n]; s=s+p*p; If[PrimeQ[s]&&IntegerQ[k], Print[{k, n, s}]], {n, 1, 10^7}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Alexander Adamchuk, Aug 26 2006
EXTENSIONS
More terms from Jinyuan Wang, Mar 23 2020
STATUS
approved