OFFSET
1,1
COMMENTS
First terms absent in A064397: 2, 18, 28, 564, 1332, 3928, 12415, 13005, 16886.
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
FORMULA
a(n) >> n^2/log^2 n. - Charles R Greathouse IV, Mar 08 2025
EXAMPLE
2 is a term because prime(2) + prime(3) = 3 + 5 = 8 = 2^3 (perfect power);
7 is a term because prime(7) + prime(8) = 17 + 19 = 36 = 6^2 (perfect power);
39867 is a term because prime(39867) + prime(39868) = 478241 + 478243 = 956484 = 978^2 (perfect power).
MAPLE
ispow:= proc(n) local F;
F:= ifactors(n)[2];
igcd(F[.., 2]) > 1
end proc:
q:= 2: count:= 0: R:= NULL:
for i from 1 while count < 100 do
p:= q; q:= nextprime(q);
if ispow(p+q) then count:= count+1; R:= R, i;
fi
od:
R; # Robert Israel, Mar 09 2026
MATHEMATICA
Select[Range[16342], ResourceFunction["PerfectPowerQ"][Prime[#]+Prime[#+1]]&] (* James C. McMahon, Mar 08 2025 *)
PROG
(PARI) s=[]; for(n=1, 41530, a=prime(n)+prime(n+1); if(ispower(a), s=concat(s, n))); s
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Nov 17 2007
STATUS
approved
