Carmichael number/Code

From Citizendium, the Citizens' Compendium
Jump to: navigation, search
This article is developing and not approved.
Main Article
Talk
Related Articles  [?]
Bibliography  [?]
External Links  [?]
Citable Version  [?]
Code [?]
 
A collection of code samples relating to the topic of Carmichael number.

Following Source Code generates Carmichael numbers, using Korselts Theorem:

/* Ein Rexx-Programm */
p.1  =   3
p.2  =   5
p.3  =   7
p.4  =  11
p.5  =  13
p.6  =  17
p.7  =  19
p.8  =  23
p.9  =  29
p.10 =  31
p.11 =  37
p.12 =  41
p.13 =  43
p.14 =  47
p.15 =  53
p.16 =  59
p.17 =  61
p.18 =  67
p.19 =  71
p.20 =  73
p.21 =  79
p.22 =  83
p.23 =  89
p.24 =  97
p.25 = 101
p.26 = 103
p.27 = 107
p.28 = 109
p.29 = 113
p.30 = 127
p.31 = 131
p.32 = 137
p.33 = 139
p.34 = 149
p.35 = 151
p.36 = 157
p.37 = 163
p.38 = 167
p.39 = 173
p.40 = 179
p.41 = 181
p.42 = 191
p.43 = 193
p.44 = 197
p.45 = 199
p.46 = 211
p.47 = 223
p.48 = 227
p.49 = 229
p.50 = 233
p.51 = 239
p.52 = 241
p.53 = 251
p.54 = 257
p.55 = 263
p.56 = 269
p.57 = 271
p.58 = 277
p.59 = 281
p.60 = 283
p.61 = 293
p.62 = 307
p.63 = 311
p.64 = 313
p.65 = 317
p.66 = 331
p.67 = 337
p.68 = 347
p.69 = 349
p.70 = 353
p.71 = 359
p.72 = 367
p.73 = 373
p.74 = 379
p.75 = 383
p.76 = 389
p.77 = 397
p.78 = 401
p.79 = 409
p.80 = 419
p.81 = 421
p.82 = 431
p.83 = 433
p.84 = 439
p.85 = 443
p.86 = 449
p.87 = 457
p.88 = 461
p.89 = 463
p.90 = 467
p.91 = 479
p.92 = 487
p.93 = 491
p.94 = 499
p.95 = 503
p.96 = 509
p.97 = 521
p.98 = 523
p.99 = 541
p.100 = 547
p.101 = 557
p.102 = 563
p.103 = 569
p.104 = 571
p.105 = 577
p.106 = 587
p.107 = 593
p.108 = 599
p.109 = 601
p.110 = 607
p.111 = 613
p.112 = 617
p.113 = 619
p.114 = 631
p.115 = 641
p.116 = 643
p.117 = 647
p.118 = 653
p.119 = 659
p.120 = 661
p.121 = 673
p.122 = 677
p.123 = 683
p.124 = 691
p.125 = 701
p.126 = 709
p.127 = 719
p.128 = 727
p.129 = 733
p.130 = 739
p.131 = 743
p.132 = 751
p.133 = 757
p.134 = 761
p.135 = 769
p.136 = 773
p.137 = 787
p.138 = 797
p.139 = 809
p.140 = 811
p.141 = 821
p.142 = 823
p.143 = 827
p.144 = 829
p.145 = 839
p.146 = 853
p.147 = 857
p.148 = 859
p.149 = 863
p.150 = 877
p.151 = 881
p.152 = 883
p.153 = 887
p.154 = 907
p.155 = 911
p.156 = 919
p.157 = 929
p.158 = 937
p.159 = 941
p.160 = 947
p.161 = 953
p.162 = 967
p.163 = 971
p.164 = 977
p.165 = 983
p.166 = 991
p.167 = 997
p.168 = 1009
p.169 = 1013
p.170 = 1019
p.171 = 1021
p.172 = 1031
p.173 = 1033
p.174 = 1039
p.175 = 1049
p.176 = 1051
p.177 = 1061
p.178 = 1063
p.179 = 1069
p.180 = 1087
p.181 = 1091
p.182 = 1093
p.183 = 1097
p.184 = 1103
p.185 = 1109
p.186 = 1117
p.187 = 1123
p.188 = 1129
p.189 = 1151
p.190 = 1153
p.191 = 1163
p.192 = 1171
p.193 = 1181
p.194 = 1187
p.195 = 1193
p.196 = 1201
p.197 = 1213
p.198 = 1217
p.199 = 1223
p.200 = 1229
p.201 = 1231
p.202 = 1237
p.203 = 1249
p.204 = 1259
p.205 = 1277
p.206 = 1279
p.206 = 1283
p.207 = 1289
p.208 = 1291
p.209 = 1297
p.210 = 1301
p.211 = 1303
p.212 = 1307
p.213 = 1319
p.214 = 1321
p.215 = 1327
p.216 = 1361
p.217 = 1367
p.218 = 1373
p.219 = 1381
p.220 = 1399
p.221 = 1409
p.222 = 1423
p.223 = 1427
p.224 = 1429
p.225 = 1433
p.226 = 1439
p.227 = 1447
p.228 = 1451
p.229 = 1453
p.230 = 1459
p.231 = 1471
p.232 = 1481
p.233 = 1483
p.234 = 1487
p.235 = 1489
p.236 = 1493
p.237 = 1499
p.238 = 1511
p.239 = 1523
p.240 = 1531
p.241 = 1543
p.242 = 1549
p.243 = 1553
p.244 = 1559
p.245 = 1567
p.246 = 1571
p.247 = 1579
p.248 = 1583
p.249 = 1597
p.250 = 1601
p.251 = 1607
p.252 = 1609
p.253 = 1613
p.254 = 1619
p.255 = 1621
p.256 = 1627
p.257 = 1637
p.258 = 1657
p.259 = 1663
p.260 = 1667
p.261 = 1669
p.262 = 1693
p.263 = 1697
p.264 = 1699
p.265 = 1709
p.266 = 1721
p.267 = 1723
p.268 = 1733
p.269 = 1741
p.270 = 1747
p.271 = 1753
p.272 = 1759
p.273 = 1777
p.274 = 1783
p.275 = 1787
p.276 = 1789
p.277 = 1801
p.278 = 1811
p.279 = 1823
p.280 = 1831
p.281 = 1847
p.282 = 1861
p.283 = 1867
p.284 = 1871
p.285 = 1873
p.286 = 1877
p.287 = 1879
p.288 = 1889
p.289 = 1901
p.290 = 1907
p.291 = 1913
p.292 = 1931
p.293 = 1933
p.294 = 1949
p.295 = 1951
p.296 = 1973
p.297 = 1979
p.298 = 1987
p.299 = 1993
p.300 = 1997
p.301 = 1999
p.302 = 2003
p.303 = 2011
p.304 = 2017
p.305 = 2027
p.306 = 2029
p.307 = 2039
p.308 = 2053
p.309 = 2063
p.310 = 2069
p.311 = 2081
p.312 = 2083
p.313 = 2087
p.314 = 2089
p.315 = 2099
p.316 = 2111
p.317 = 2113
p.318 = 2129
p.319 = 2131
p.320 = 2137
p.321 = 2141
p.322 = 2143
p.323 = 2153
p.324 = 2161
p.325 = 2179
p.326 = 2203
p.327 = 2207
p.328 = 2213
p.329 = 2221
p.330 = 2237
p.331 = 2239
p.332 = 2243
p.333 = 2251
p.334 = 2267
p.335 = 2269
p.336 = 2273
p.337 = 2281
p.338 = 2287
p.339 = 2293
p.340 = 2297
p.341 = 2309
p.342 = 2311
p.343 = 2333
p.344 = 2339
p.345 = 2341
p.346 = 2347
p.347 = 2351
p.348 = 2357
p.349 = 2371
p.350 = 2377
p.351 = 2381
p.352 = 2383
p.353 = 2389
p.354 = 2393
p.355 = 2399
p.356 = 2411
p.357 = 2417
p.358 = 2423
p.359 = 2437
p.360 = 2441
p.361 = 2447
p.362 = 2459
p.363 = 2467
p.364 = 2473
p.365 = 2477
p.366 = 2503
p.367 = 2521
p.368 = 2531
p.369 = 2539
p.370 = 2543
p.371 = 2549
p.372 = 2551
p.373 = 2557
p.374 = 2579
p.375 = 2591
p.376 = 2593
p.377 = 2609
p.378 = 2617
p.379 = 2621
p.380 = 2633
p.381 = 2647
p.382 = 2657
p.383 = 2659
p.384 = 2663
p.385 = 2671
p.386 = 2677
p.387 = 2683
p.388 = 2687
p.389 = 2689
p.390 = 2693
p.391 = 2699
p.392 = 2707
p.393 = 2711
p.394 = 2713
p.395 = 2719
p.396 = 2729
p.397 = 2731
p.398 = 2741
p.399 = 2749
p.400 = 2753
p.401 = 2767
p.402 = 2777
p.403 = 2789
p.404 = 2791
p.405 = 2797
p.406 = 2801
p.407 = 2803
p.408 = 2819
p.409 = 2833
p.410 = 2837
p.411 = 2843
p.412 = 2851
p.413 = 2857
p.414 = 2861
p.415 = 2879
p.416 = 2887
p.417 = 2897
p.418 = 2903
p.419 = 2909
p.420 = 2917
p.421 = 2927
p.422 = 2939
p.423 = 2953
p.424 = 2957
p.425 = 2963
p.426 = 2969
p.427 = 2971
p.428 = 2999
p.429 = 3001
p.430 = 3011
p.431 = 3019
p.432 = 3023
p.433 = 3037
p.434 = 3041
p.435 = 3049
p.436 = 3061
p.437 = 3067
p.438 = 3079
p.439 = 3083
p.440 = 3089
p.441 = 3109
p.442 = 3119
p.443 = 3121
p.444 = 3137
p.445 = 3163
p.446 = 3167
p.447 = 3169
p.448 = 3181
p.449 = 3187
p.450 = 3191
p.451 = 3203
p.452 = 3209
p.453 = 3217
p.454 = 3221
p.455 = 3229
p.456 = 3251
p.457 = 3253
p.458 = 3257
p.459 = 3259
p.460 = 3271
p.461 = 3299
p.462 = 3301
p.463 = 3307
p.464 = 3313
p.465 = 3319
p.466 = 3323
p.467 = 3329
p.468 = 3331
p.469 = 3343
p.470 = 3347
p.471 = 3359
p.472 = 3361
p.473 = 3371
p.474 = 3373
p.475 = 3389
p.476 = 3391
p.477 = 3407
p.478 = 3413
p.479 = 3433
p.480 = 3449
p.481 = 3457
p.482 = 3461
p.483 = 3463
p.484 = 3467
p.485 = 3469
p.486 = 3491
p.487 = 3499
p.488 = 3511
p.489 = 3517
p.490 = 3527
p.491 = 3529
p.492 = 3533
p.493 = 3539
p.494 = 3541
p.495 = 3547
p.496 = 3557
p.497 = 3559
p.498 = 3571
p.499 = 3581
p.500 = 3583
p.501 = 3593
p.502 = 3607
p.503 = 3613
p.504 = 3617
p.505 = 3623
p.506 = 3631
p.507 = 3637
p.508 = 3643
p.509 = 3659
p.510 = 3671
p.511 = 3673
p.512 = 3677
p.513 = 3691
p.514 = 3697
p.515 = 3701
p.516 = 3709
p.517 = 3719
p.518 = 3727
p.519 = 3733
p.520 = 3739
p.521 = 3761
p.522 = 3767
p.523 = 3769
p.524 = 3779
p.525 = 3793
p.526 = 3797
p.527 = 3803
p.528 = 3821
p.529 = 3823
p.530 = 3833
p.531 = 3847
p.532 = 3851
p.533 = 3853
p.534 = 3863
p.535 = 3877
p.536 = 3881
p.537 = 3889
p.538 = 3907
p.539 = 3911
p.540 = 3917
p.541 = 3919
p.542 = 3923
p.543 = 3929
p.544 = 3931
p.545 = 3943
p.546 = 3947
p.547 = 3967
p.548 = 3989
p.549 = 4001
p.550 = 4003
p.551 = 4007
p.552 = 4013
p.553 = 4019
p.554 = 4021
p.555 = 4027
p.556 = 4049
p.557 = 4051
p.558 = 4057
i=558
do a=1 to (i-2)
  do b=a+1 to (i-1)
    do c=b+1 to i
      t = 0
      z=p.a*p.b*p.c
      ax=p.a-1
      bx=p.b-1
      cx=p.c-1
      zax=(z/p.a)-1
      zbx=(z/p.b)-1
      zcx=(z/p.c)-1
      pz=(zax // ax) + (zbx // bx) + (zcx // cx)
      if (pz = 0) then t = 1
      if (t = 1) then do
        say z||"="||p.a||"*"||p.b||"*"||p.c
        lineout("rcrmn___.txt",z||" = "||p.a||"*"||p.b||"*"||p.c)
      end
      t=0
    end
  end
end

A Source Code, generating Carmichael Numbers with three prime factors

do a=1 to (i-3)
  do b=a+1 to (i-2)
    do c=b+1 to (i-1)
      do d=c+1 to i
        t = 0
        z=p.a*p.b*p.c*p.d
        ax=p.a-1
        bx=p.b-1
        cx=p.c-1
        dx=p.d-1
        zax=(z/p.a)-1
        zbx=(z/p.b)-1
        zcx=(z/p.c)-1
        zdx=(z/p.d)-1
        pz=(zax // ax) + (zbx // bx) + (zcx // cx) + (zdx // dx)
        if (pz = 0) then t = 1
        if (t = 1) then do
          say z||"="||p.a||"*"||p.b||"*"||p.c||"*"||p.d
          lineout("rcrmn4_.txt",z||" = "||p.a||"*"||p.b||"*"||p.c||"*"||p.d)
        end
        t=0
      end
    end
  end
end

Main routine expanded to Calculate Carmichael numbers with four Prime factors