generate_psp_from_primorials_2.sf 379 B

12345678910111213141516171819
  1. #!/usr/bin/ruby
  2. # Most of these pseudoprimes are also super-pseudoprimes.
  3. var j = 2 # even multiple
  4. var min = 2**64
  5. for n in (1..1000) {
  6. #var t = n.pn_primorial
  7. #var t = n.factorial
  8. var t = n.prime.consecutive_lcm
  9. for k in (n .. 1000) {
  10. var u = ((k*t + 1) * (j*k*t + 1))
  11. if (u > min && u.is_psp) {
  12. say u
  13. }
  14. }
  15. }