Cubic special primes: Difference between revisions
→{{header|Julia}}
(Realize in F#) |
|||
Line 99:
<pre>
Same as Wren example.
</pre>
=={{header|Julia}}==
<lang julia>using Primes
function cubicspecialprimes(N = 15000)
pmask = primesmask(1, N)
cprimes, maxidx = [2], isqrt(N)
while (n = cprimes[end]) < N
for i in 1:maxidx
q = n + i * i * i
if q > N
return cprimes
elseif pmask[q] # got next qprime
push!(cprimes, q)
break
end
end
end
end
println("Cubic special primes < 16000:")
foreach(p -> print(rpad(p[2], 6), p[1] % 10 == 0 ? "\n" : ""), enumerate(cubicspecialprimes()))
</lang>{{out}}
<pre>
Cubic special primes < 16000:
2 3 11 19 83 1811 2027 2243 2251 2467
2531 2539 3539 3547 4547 5059 10891 12619 13619 13627
13691 13907 14419
</pre>
| |||