Pierpont primes: Difference between revisions

m
Undo revision 287012 by Wherrera (talk)
(julia example)
m (Undo revision 287012 by Wherrera (talk))
Line 164:
 
250th Pierpont prime of the second kind: 4111131172000956525894875083702271
</pre>
 
=={{header|Julia}}==
<lang julia>using Primes
 
function pierponts(N, firstkind = true)
count, ret, incdec = 0, typeof(N)[], firstkind ? 1 : -1
for k2 in 0:1000
for k3 in 0:k2
n = typeof(N)(2)^k2 * typeof(N)(3)^k3 + incdec
if isprime(n)
push!(ret, n)
if length(ret) >= N
return ret
end
end
end
end
throw("Failed to find $N primes")
end
 
println("The first 50 Pierpont primes (first kind) are: ", pierponts(50))
 
println("\nThe first 50 Pierpont primes (second kind) are: ", pierponts(50, false))
 
println("\nThe 250th Pierpont prime (first kind) is: ", pierponts(BigInt(250))[250])
 
println("\nThe 250th Pierpont prime (second kind) is: ", pierponts(BigInt(250), false)[250])
</lang>{{out}}
<pre>
The first 50 Pierpont primes (first kind) are: [2, 3, 7, 5, 13, 37, 73, 17, 433, 1297, 97, 2593, 193, 577, 1153, 3457, 10369, 257, 769, 746497, 18433, 1492993, 120932353, 12289, 331777, 995329, 725594113, 1990657, 483729409, 147457, 2654209, 71663617, 65537, 1769473, 5308417, 11609505793, 1179649, 69657034753, 786433, 63700993, 1253826625537, 33853318889473, 14155777, 3439853569, 22568879259649, 28311553, 6879707137, 169869313, 123834728449, 90275517038593]
 
The first 50 Pierpont primes (second kind) are: [5, 3, 11, 7, 23, 71, 47, 431, 31, 863, 2591, 191, 15551, 127, 383, 1151, 6911, 62207, 27647, 6143, 165887, 497663, 995327, 8191, 73727, 17915903, 4353564671, 442367, 294911, 214990847, 131071, 10616831, 95551487, 5642219814911, 786431, 63700991, 417942208511, 524287, 10319560703, 185752092671, 18874367, 169869311, 9172942847, 82556485631, 60183678025727, 6115295231, 495338913791, 361102068154367, 4076863487, 10871635967]
 
The 250th Pierpont prime (first kind) is: 14847887199655905933701142860524349856292012033
 
The 250th Pierpont prime (second kind) is: 325879554450609732426911228763223042113404927
</pre>
 
4,108

edits