prog.sf 357 B

1234567891011121314151617181920
  1. #!/usr/bin/ruby
  2. # Numbers k such that (23^k - 2^k)/21 is prime.
  3. # https://oeis.org/A375161
  4. # Known terms:
  5. # 5, 11, 197, 4159
  6. # Next term, if it exists, is greater than 14956.
  7. for k in (14956 .. 1e6) {
  8. var v = (23**k - (1 << k))/21
  9. say "Testing: #{k} (#{v.len} digits)"
  10. if (v.is_prob_prime) {
  11. die "Found new term: #{k}\n"
  12. }
  13. }