OFFSET
1,1
LINKS
Peter Luschny, Table of n, a(n) for n = 1..1200
EXAMPLE
1157 is in this sequence because 1157 = (5 + 8)*(5^2 + 8^2) = 13*89.
PROG
(Julia) # Returns the terms less than or equal to b^3.
using Nemo
function A349202List(b)
b3 = b^3; R = Int[]
for n in 1:b
for k in 0:n
a = (n + k) * (n^2 + k^2)
a > b3 && break
isprime(n+k) && isprime(n^2 + k^2) && push!(R, a)
end end
sort(R) end
A349202List(34) |> println
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Nov 11 2021
STATUS
approved