OFFSET
1,1
COMMENTS
This is a subset of A335410.
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000 (first 102 terms from Chai Wah Wu)
EXAMPLE
Prime(23)=83. Prime(22)=79. 2*(83^2 - 79^2) = 6^4.
Prime(378)=2593. Prime(377)=2591. 2*(2593^2 - 2591^2) = 12^4.
MATHEMATICA
Select[Prime@Range[2, 500000], IntegerQ@Sqrt[Sqrt[2(#^2 - NextPrime[#, -1]^2)]]&] (* a modification of Giovanni Resta's program for A335410 *)
PROG
(PARI) isok(p) = isprime(p) && ispower(2*(p^2-precprime(p-1)^2), 4); \\ Michel Marcus, Jun 08 2020
(PARI) lista(nn) = {my(pp=2); forprime(p=3, nn, if (ispower(2*(p^2 - pp^2), 4), print1(p, ", ")); pp = p; ); } \\ Michel Marcus, Jun 08 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Jeff Brown, Jun 08 2020
EXTENSIONS
More terms from Amiram Eldar, Jun 08 2020
More terms from Giovanni Resta, Jun 08 2020
STATUS
approved