SharedMadeleine.sagewsOpen in CoCalc
Author: Paul Zeitz
Views : 9
# define the primes P = Primes()
P
Set of all prime numbers: 2, 3, 5, 7, ...
P.next(99)
101
P.unrank(8)
23
for k in range(20): print k, P.unrank(k)
0 2 1 3 2 5 3 7 4 11 5 13 6 17 7 19 8 23 9 29 10 31 11 37 12 41 13 43 14 47 15 53 16 59 17 61 18 67 19 71
#random prime out of first n n = 1000 d= randint(2,20) for _ in range(1000): p=P.unrank(randint(1,n)) a = p+d b = p + 2*d if a.is_prime()*b.is_prime()==1: print p, a, b
1531 1549 1567 3371 3389 3407 2861 2879 2897 163 181 199 1583 1601 1619 1913 1931 1949 1213 1231 1249 4621 4639 4657 5 23 41 2693 2711 2729 233 251 269 383 401 419 4603 4621 4639 3911 3929 3947 2251 2269 2287 2423 2441 2459 163 181 199 1601 1619 1637 4603 4621 4639 1213 1231 1249 2671 2689 2707 2423 2441 2459 53 71 89 7603 7621 7639 5413 5431 5449 3691 3709 3727 911 929 947 673 691 709 443 461 479 4373 4391 4409 1531 1549 1567 4373 4391 4409 7481 7499 7517 53 71 89 5843 5861 5879 4621 4639 4657 2671 2689 2707 5003 5021 5039 233 251 269 421 439 457 71 89 107 673 691 709 4373 4391 4409 4951 4969 4987 2731 2749 2767 733 751 769 3691 3709 3727 61 79 97 2063 2081 2099 5413 5431 5449 2203 2221 2239 1013 1031 1049 2203 2221 2239 7481 7499 7517 2423 2441 2459 1583 1601 1619 5003 5021 5039 1871 1889 1907 3001 3019 3037 421 439 457 43 61 79 7541 7559 7577 641 659 677 4373 4391 4409 5413 5431 5449 5261 5279 5297 691 709 727 4093 4111 4129 673 691 709 1453 1471 1489 4603 4621 4639 7901 7919 7937 7193 7211 7229 5381 5399 5417 3001 3019 3037 61 79 97 2251 2269 2287 23 41 59 6673 6691 6709 1283 1301 1319 6983 7001 7019 2063 2081 2099 751 769 787 61 79 97 7541 7559 7577 2801 2819 2837 11 29 47 1013 1031 1049 7541 7559 7577 2521 2539 2557 4621 4639 4657 53 71 89