Carmichael 3 strong pseudoprimes: Difference between revisions

(→‎{{header|Julia}}: A new entry for Julia)
Line 1,146:
{{works with|Ruby|1.9}}
<lang ruby># Generate Charmichael Numbers
 
#
# Nigel_Galloway
# November 30th., 2012.
#
require 'prime'
 
IntegerPrime.each_primeeach(61) {do |p|
(2...p).each {do |h3|
g = h3 + p
(1...g).each {do |d|
next if (g*(p-1)) % d != 0 or (-1*p*p) % h3 != d % h3
q = 1 + ((p - 1) * g / d)
next if notunless q.prime?
r = 1 + (p * q / h3)
next if notunless r.prime? or notand (q * r) % (p - 1) == 1
puts "#{p} Xx #{q} Xx #{r}"
}end
}end
puts ""
}end</lang>
 
{{out}}
<pre style="height:30ex;overflow:scroll">
3 Xx 11 Xx 17
 
5 Xx 29 Xx 73
5 Xx 17 Xx 29
5 Xx 13 Xx 17
 
7 Xx 19 Xx 67
7 Xx 31 Xx 73
7 Xx 13 Xx 31
7 Xx 23 Xx 41
7 Xx 73 Xx 103
7 Xx 13 Xx 19
 
 
13 Xx 61 Xx 397
13 Xx 37 Xx 241
13 Xx 97 Xx 421
13 Xx 37 Xx 97
13 Xx 37 Xx 61
 
17 Xx 41 Xx 233
17 Xx 353 Xx 1201
 
19 Xx 43 Xx 409
19 Xx 199 Xx 271
 
23 Xx 199 Xx 353
 
29 Xx 113 Xx 1093
29 Xx 197 Xx 953
 
31 Xx 991 Xx 15361
31 Xx 61 Xx 631
31 Xx 151 Xx 1171
31 Xx 61 Xx 271
31 Xx 61 Xx 211
31 Xx 271 Xx 601
31 Xx 181 Xx 331
 
37 Xx 109 Xx 2017
37 Xx 73 Xx 541
37 Xx 613 Xx 1621
37 Xx 73 Xx 181
37 Xx 73 Xx 109
 
41 Xx 1721 Xx 35281
41 Xx 881 Xx 12041
41 Xx 101 Xx 461
41 Xx 241 Xx 761
41 Xx 241 Xx 521
41 Xx 73 Xx 137
41 Xx 61 Xx 101
 
43 Xx 631 Xx 13567
43 Xx 271 Xx 5827
43 Xx 127 Xx 2731
43 Xx 127 Xx 1093
43 Xx 211 Xx 757
43 Xx 631 Xx 1597
43 Xx 127 Xx 211
43 Xx 211 Xx 337
43 Xx 433 Xx 643
43 Xx 547 Xx 673
43 Xx 3361 Xx 3907
 
47 Xx 3359 Xx 6073
47 Xx 1151 Xx 1933
47 Xx 3727 Xx 5153
 
53 Xx 157 Xx 2081
53 Xx 79 Xx 599
53 Xx 157 Xx 521
 
59 Xx 1451 Xx 2089
 
61 Xx 421 Xx 12841
61 Xx 181 Xx 5521
61 Xx 1301 Xx 19841
61 Xx 277 Xx 2113
61 Xx 181 Xx 1381
61 Xx 541 Xx 3001
61 Xx 661 Xx 2521
61 Xx 271 Xx 571
61 Xx 241 Xx 421
61 Xx 3361 Xx 4021
</pre>
 
Anonymous user