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 20050000
#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 ULongInt n, count, k
Dim As LongInt n = -1, count, k
Dim As LongInt si = 15
Dim As LongInt si = 15


Print "The first fifty quad-power prime seeds are:"
Print "The first thirty penta-power prime seeds are:"
While count < 50
While count < 30
n += 1
n += 2
k = n +1
k = n +1
If sieve(n + k) Then ' skip if n + k is not prime
' 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 quad-power prime seed greater than:"
Print !"\n\nFirst penta-power prime seed greater than:"
While m < 11
While m < 11
n += 1
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 fifty quad-power prime seeds are:
<pre>The first thirty penta-power prime seeds are:
1 2 5 6 69 131 426 1665 2129 2696
1 5 69 1665 2129 25739 29631 62321 77685 80535
5250 7929 9689 13545 14154 14286 16434 19760 25739 27809
82655 126489 207285 211091 234359 256719 366675 407945 414099 628859
29631 36821 41819 46619 48321 59030 60500 61955 62321 73610
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 quad-power prime seed greater than:
First penta-power prime seed greater than:
1 million is 1,009,286 at index 141
1 million is 1,121,189 at index 26
2 million is 2,015,496 at index 234
2 million is 2,066,079 at index 39
3 million is 3,005,316 at index 319
3 million is 3,127,011 at index 47
4 million is 4,004,726 at index 383
4 million is 4,059,525 at index 51
5 million is 5,023,880 at index 452
5 million is 5,279,175 at index 59
6 million is 6,000,554 at index 514
6 million is 6,320,601 at index 63
7 million is 7,047,129 at index 567
7 million is 7,291,361 at index 68
8 million is 8,005,710 at index 601
8 million is 8,334,915 at index 69
9 million is 9,055,151 at index 645
9 million is 9,100,671 at index 71
10 million is 10,023,600 at index 701</pre>
10 million is 10,347,035 at index 72</pre>


=={{header|Go}}==
=={{header|Go}}==
Retrieved from "https://rosettacode.org/wiki/Penta-power_prime_seeds"