Penta-power prime seeds: Difference between revisions
→{{header|FreeBASIC}}: correct program
(added FreeBASIC) |
(→{{header|FreeBASIC}}: correct program) |
||
Line 283:
#Include "gmp.bi"
#Define sieve_max
Dim As Mpz_ptr n2 = Allocate (Len(__mpz_struct))
Line 307:
Next
Dim As
Dim As LongInt si = 15
Print "The first
While count <
n +=
k = n +1
If sieve(1 + k) And sieve(n + k) Then ' skip if 1 + k or n + k is not prime
Mpz_ui_pow_ui(n4, n , 4)
Mpz_add_ui(n4, n4, k)
Line 324 ⟶ 325:
Mpz_add_ui(n2, n2, k)
If Mpz_probab_prime_p(n2, si) >= 1 Then ' if prime then print n
Print Using "##########"; n;
count += 1
If (count Mod 10) = 0 Then Print
Line 332 ⟶ 333:
Dim As ULongInt m = 1, million = 1000000
n = -1 : count = 0
Print !"\n\nFirst
While m < 11
n +=
k = n +1
If sieve(1 + k) And sieve(n + k) Then ' skip if 1 + k or n + k is not prime
Mpz_ui_pow_ui(n4, n , 4)
Mpz_add_ui(n4, n4, k)
Line 357 ⟶ 358:
Wend
Mpz_clear(n4) : Mpz_clear(n3) : Mpz_clear(n2)
Line 366 ⟶ 367:
End</syntaxhighlight>
{{out}}
<pre>The first
1
First
1 million is 1,
2 million is 2,
3 million is 3,
4 million is 4,
5 million is 5,
6 million is 6,
7 million is 7,
8 million is 8,
9 million is 9,
10 million is 10,
=={{header|Go}}==
| |||