OFFSET
1,1
COMMENTS
Numbers k with at least one divisor d < sqrt(k) such that (d + k/d)/2 is prime.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(5) = 13 is a term because 13 = 7^2 - 6^2 where 7 is prime.
MAPLE
filter:= proc(n) ormap(t -> (t+n/t)::even and t > n/t and isprime((t + n/t)/2), NumberTheory:-Divisors(n)) end proc:
select(filter, [$1..300]);
CROSSREFS
KEYWORD
nonn
AUTHOR
Will Gosnell and Robert Israel, Jan 31 2026
STATUS
approved
