Penta-power prime seeds: Difference between revisions
Content added Content deleted
(added FreeBASIC) |
(→{{header|FreeBASIC}}: correct program) |
||
Line 283: | Line 283: | ||
#Include "gmp.bi" |
#Include "gmp.bi" |
||
#Define sieve_max |
#Define sieve_max 21000000 |
||
Dim As Mpz_ptr n2 = Allocate (Len(__mpz_struct)) |
Dim As Mpz_ptr n2 = Allocate (Len(__mpz_struct)) |
||
Line 307: | Line 307: | ||
Next |
Next |
||
Dim As |
Dim As LongInt n = -1, count, k |
||
Dim As LongInt si = 15 |
Dim As LongInt si = 15 |
||
Print "The first |
Print "The first thirty penta-power prime seeds are:" |
||
While count < |
While count < 30 |
||
n += |
n += 2 |
||
k = n +1 |
k = n +1 |
||
' n ^ 0 = 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_ui_pow_ui(n4, n , 4) |
||
Mpz_add_ui(n4, n4, k) |
Mpz_add_ui(n4, n4, k) |
||
Line 324: | Line 325: | ||
Mpz_add_ui(n2, n2, k) |
Mpz_add_ui(n2, n2, k) |
||
If Mpz_probab_prime_p(n2, si) >= 1 Then ' if prime then print n |
If Mpz_probab_prime_p(n2, si) >= 1 Then ' if prime then print n |
||
Print Using "########"; n; |
Print Using "##########"; n; |
||
count += 1 |
count += 1 |
||
If (count Mod 10) = 0 Then Print |
If (count Mod 10) = 0 Then Print |
||
Line 332: | Line 333: | ||
Dim As ULongInt m = 1, million = 1000000 |
Dim As ULongInt m = 1, million = 1000000 |
||
n = -1 : count = 0 |
|||
Print !"\n\nFirst |
Print !"\n\nFirst penta-power prime seed greater than:" |
||
While m < 11 |
While m < 11 |
||
n += |
n += 2 |
||
k = n +1 |
k = n +1 |
||
If sieve(n + k) Then ' skip if n + k is not prime |
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_ui_pow_ui(n4, n , 4) |
||
Mpz_add_ui(n4, n4, k) |
Mpz_add_ui(n4, n4, k) |
||
Line 357: | Line 358: | ||
Wend |
Wend |
||
Mpz_clear(n4) : Mpz_clear(n3) : Mpz_clear(n2) |
Mpz_clear(n4) : Mpz_clear(n3) : Mpz_clear(n2) |
||
Line 366: | Line 367: | ||
End</syntaxhighlight> |
End</syntaxhighlight> |
||
{{out}} |
{{out}} |
||
<pre>The first |
<pre>The first thirty penta-power prime seeds are: |
||
1 |
1 5 69 1665 2129 25739 29631 62321 77685 80535 |
||
82655 126489 207285 211091 234359 256719 366675 407945 414099 628859 |
|||
644399 770531 781109 782781 923405 1121189 1158975 1483691 1490475 1512321 |
|||
77685 79646 80535 82655 85251 86996 91014 96566 97739 105939 |
|||
108240 108681 119754 122436 123164 126489 140636 150480 153179 163070 |
|||
First |
First penta-power prime seed greater than: |
||
1 million is 1, |
1 million is 1,121,189 at index 26 |
||
2 million is 2, |
2 million is 2,066,079 at index 39 |
||
3 million is 3, |
3 million is 3,127,011 at index 47 |
||
4 million is 4, |
4 million is 4,059,525 at index 51 |
||
5 million is 5, |
5 million is 5,279,175 at index 59 |
||
6 million is 6, |
6 million is 6,320,601 at index 63 |
||
7 million is 7, |
7 million is 7,291,361 at index 68 |
||
8 million is 8, |
8 million is 8,334,915 at index 69 |
||
9 million is 9, |
9 million is 9,100,671 at index 71 |
||
10 million is 10, |
10 million is 10,347,035 at index 72</pre> |
||
=={{header|Go}}== |
=={{header|Go}}== |