Quad-power prime seeds: Difference between revisions
Content added Content deleted
m (Reformatting.) |
(julia example) |
||
Line 553: | Line 553: | ||
2 million is 2015496 at index 234 |
2 million is 2015496 at index 234 |
||
3 million is 3005316 at index 319 |
3 million is 3005316 at index 319 |
||
</pre> |
|||
=={{header|Julia}}== |
|||
{{trans|Python}} |
|||
<syntaxhighlight lang="julia">using Primes, Formatting |
|||
isquadpowerprime(x) = all(isprime, [2x + 1, x^2 + x + 1, x^3 + x + 1, x^4 + x + 1]) |
|||
const qpprimes = filter(isquadpowerprime, Int128(1):10_100_000) |
|||
foreach(n -> print(lpad(qpprimes[n], 9), n % 10 == 0 ? "\n" : ""), 1:50) |
|||
for j in 1_000_000:1_000_000:10_000_000 |
|||
for p in qpprimes |
|||
if p > j |
|||
println("The first quad-power prime seed over ", format(j, commas = true), |
|||
" is ", format(p, commas = true)) |
|||
break |
|||
end |
|||
end |
|||
end |
|||
</syntaxhighlight>{{out}} |
|||
<pre> |
|||
1 2 5 6 69 131 426 1665 2129 2696 |
|||
5250 7929 9689 13545 14154 14286 16434 19760 25739 27809 |
|||
29631 36821 41819 46619 48321 59030 60500 61955 62321 73610 |
|||
77685 79646 80535 82655 85251 86996 91014 96566 97739 105939 |
|||
108240 108681 119754 122436 123164 126489 140636 150480 153179 163070 |
|||
The first quad-power prime seed over 1,000,000 is 1,009,286 |
|||
The first quad-power prime seed over 2,000,000 is 2,015,496 |
|||
The first quad-power prime seed over 3,000,000 is 3,005,316 |
|||
The first quad-power prime seed over 4,000,000 is 4,004,726 |
|||
The first quad-power prime seed over 5,000,000 is 5,023,880 |
|||
The first quad-power prime seed over 6,000,000 is 6,000,554 |
|||
The first quad-power prime seed over 7,000,000 is 7,047,129 |
|||
The first quad-power prime seed over 8,000,000 is 8,005,710 |
|||
The first quad-power prime seed over 9,000,000 is 9,055,151 |
|||
The first quad-power prime seed over 10,000,000 is 10,023,600 |
|||
</pre> |
</pre> |
||