OFFSET
1,1
COMMENTS
There are 266 such primes < 10^6, the largest being 994871 -> 9^2+4^2+7^2 = 9^2+8^2+1^2 = 146.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
1487 is in the sequence because 1^2+8^2 = 4^2+7^4 = 65.
MAPLE
filter:= proc(n) local L, i;
L:= convert(n, base, 10);
add(L[i]^2*(-1)^i, i=1..nops(L)) = 0 and isprime(n)
end proc:
select(filter, [seq(i, i=3..2*10^5, 2)]); # Robert Israel, Mar 18 2026
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Zak Seidov, Nov 01 2002
STATUS
approved
