OFFSET
1,1
PROG
(Julia)
function isA299505(n)
n % 5 >= 2 && return false
n == 5 && return true
K = Int(floor(5.383*log(n)^1.161))
M = Int(floor(2*sqrt(n/3)))
for k in 3:K
for y in 2:M, x in 1:y
n == x^4+y*x^3+y^2*x^2+y^3*x+y^4 && return true
end end
return false
end
A299505list(upto) = [n for n in 1:upto if isA299505(n)]
println(A299505list(18000))
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Mar 02 2018
STATUS
approved