verify.sf 2.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263
  1. #!/usr/bin/ruby
  2. # a(n) = largest prime factor of the number with decimal expansion 20305070...0p_n where p_n = n-th prime.
  3. # https://oeis.org/A308899
  4. DATA.lines.each_kv {|n,line|
  5. say "Testing: #{n+1}"
  6. var p = line.nums[-1]
  7. var k = primes(prime(n+1)).join('0').to_i
  8. assert(k % p == 0)
  9. assert(p.is_prime)
  10. assert((k / p).gpf < p)
  11. }
  12. __END__
  13. 1 2
  14. 2 29
  15. 3 131
  16. 4 33287
  17. 5 17627
  18. 6 1754975809
  19. 7 59218567
  20. 8 318879703697
  21. 9 2030507011013017019023
  22. 10 14400758943354730631369
  23. 11 1016015647
  24. 12 32002443156997
  25. 13 2464082401591041689
  26. 14 4916481866859605372937116297910511
  27. 15 2030507011013017019023029031037041043047
  28. 16 5517847030789010049276008334588917177
  29. 17 54878567865216676189811595433433541703974407
  30. 18 676835670337672339674343010345680347682351019687
  31. 19 370133692634226265626294781970616577
  32. 20 4576955195368289951819122267757
  33. 21 9263128224909627760096007651
  34. 22 2825491798345128304972687
  35. 23 2520893999924289042092178403737739835508730989770043950773
  36. 24 577141136050035836784414812892484771166264394975674023
  37. 25 2030507011013017019023029031037041043047053059061067071073079083089097
  38. 26 33453922728575399447059064285237470339
  39. 27 2605679122539420040611276520344407664126070570864875068894443081659773
  40. 28 7078390162660648591303027564507832327314969886325011
  41. 29 10880835735448418808729106659755283530233
  42. 30 566547395294271356579972489809180130253116886673343686243582315021290875271719
  43. 31 3672340814862605587840224392789976548187081758921409955297392715911547640976399991700233
  44. 32 29855508690799180530426359664796915000531321031
  45. 33 4143891859210238814332712308238859271524598079716463410353222618549177571634755320430410469646965513
  46. 34 3122496557476233834596094561307115324227505351820364739558022224664403954799
  47. 35 3718314791261006233951166825065157
  48. 36 3291111452472429759816895931882352753219740329
  49. 37 2766933537764655468658721968146169913610369507647616102703629443059
  50. 38 1826610717820635648564307697794849319246252573523144839322588582804037373866989618610063
  51. 39 44797427613949353737281434417090568954360344661385909506008921161983
  52. 40 4287487697165721794684904064572954699216098893895481252886108172041785197926957988912118166019508690806034377701995099160457
  53. 41 287352875009393763297406312604453171476946890684342438221480138897422567250246925218132917080131509261866780015453
  54. 42 870386791374524000382532028387196131700129377669591865030078014157724334829456524490874097563575226284459650060893709646525437
  55. 43 83500354876663346987233791194771234357889225153327905842255232683548803811782081502976665146717