OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..400
EXAMPLE
a(3) = 29983 is a term because it is prime and A001414(29983^4-1) = 2*7+3+5+13*2+19+263+521+937+1021 = 2809 = 53^2 and 53 is a prime.
MAPLE
spf:= proc(n) local t; add(t[1]*t[2], t=ifactors(n)[2]) end proc:
filter:= proc(n) local s; s:= spf(n-1)+spf(n+1)+spf(n^2+1); issqr(s) and isprime(sqrt(s)) end proc:
select(filter, [seq(ithprime(i), i=1..10^5)]);
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Mar 23 2021
STATUS
approved