prog.sf 825 B

123456789101112131415161718192021222324252627282930313233343536
  1. #!/usr/bin/ruby
  2. # Least Lucas number with n Lucas divisors
  3. # https://oeis.org/A356123
  4. # Known terms:
  5. # 1, 2, 4, 18, 1364, 1860498, 2537720636, 6440026026380244498, 8784200221406821330636, 77162173529763648886126034136172445632164498
  6. include("../../../factordb/auto.sf")
  7. var table = Hash()
  8. for k in (1..1e6) {
  9. var t = divisors(k.lucas).count{.is_lucas}
  10. if (!table.has(t)) {
  11. table{t} = k
  12. say("a(#{t}) = #{k.lucas}")
  13. }
  14. }
  15. __END__
  16. a(1) = 1
  17. a(2) = 3
  18. a(3) = 4
  19. a(4) = 18
  20. a(5) = 1364
  21. a(6) = 1860498
  22. a(7) = 2537720636
  23. a(8) = 6440026026380244498
  24. a(9) = 8784200221406821330636
  25. a(10) = 77162173529763648886126034136172445632164498
  26. a(13) = 677807981804381412733419097483256033462649109976189806684020291364
  27. a(11) = 4365101043708483494615466932242949707161871659736799144058331102381689400753867700636