lcm.sf 3.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180
  1. #!/usr/bin/ruby
  2. var a = [
  3. 561
  4. 62745
  5. 576480525985
  6. 1886616373665
  7. 3193231538989185
  8. 11947816523586945
  9. 101817952350880305
  10. 171800042106877185
  11. 101817952350880305
  12. 11947816523586945
  13. 101817952350880305
  14. 171800042106877185
  15. 11947816523586945
  16. 101817952350880305
  17. 171800042106877185
  18. 1177908521713261185
  19. 1068792168802099105
  20. 1114090435183709985
  21. 1177908521713261185
  22. 1177996857330634705
  23. 1254323866883045505
  24. 1344418370696745345
  25. 1355376045311269345
  26. 1841108721934863745
  27. 1873766365160431105
  28. 1974201289481769985
  29. 1975738329089659905
  30. 2102479398374245345
  31. 2107496096157739105
  32. 2166093710941883905
  33. 2325605354531800705
  34. 2337521512943074305
  35. 2435191635734599105
  36. 2435279804766429985
  37. 2731573851245420545
  38. 2950323884558198785
  39. 2978083893490447105
  40. 3077672792834337985
  41. 3174706634458860865
  42. 3239491439742459265
  43. 3297345708582907345
  44. 3368213950731767985
  45. 3390302891095115905
  46. 3399464645479365585
  47. 4400268692628910465
  48. 4442444347156637185
  49. 4448649946743866865
  50. 4521969833754314305
  51. 4682767770306971745
  52. 4757780389128977665
  53. 4997145281197228545
  54. 5092454223678002545
  55. 5473367446820370945
  56. 5935592741005669105
  57. 6088257084243875905
  58. 6143972295188340705
  59. 6687537029316101505
  60. 7258662586258825105
  61. 7626934787034827265
  62. 8104948400265498145
  63. 8638848567964719105
  64. 8684725437411304065
  65. 8932813088456120865
  66. 9142879289181859585
  67. 9422671481221099465
  68. 9689119695229916305
  69. 9718427839254355585
  70. 9744543540935384545
  71. 10186034230241370865
  72. 10885189367088417745
  73. 11356964096926592145
  74. 11595472526856600705
  75. 12275497073183645665
  76. 12376763082073211185
  77. 12525679299253814785
  78. 13294381492836585505
  79. 14631273312068834785
  80. 15131182782493999585
  81. 15157968873827062465
  82. 16888413473957767105
  83. 17215556534723568385
  84. 17445622572008652385
  85. 18329749222420733505
  86. 1886616373665
  87. 11985185775745
  88. 23239986511105
  89. 27607174936705
  90. 62139586530385
  91. 98572590097345
  92. 158353658932305
  93. 269040992399565
  94. 312220789782265
  95. 655510549443465
  96. 831047375197585
  97. 881715504450705
  98. 1590202557925345
  99. 1817671359979245
  100. 1936529152138945
  101. 2263411171923265
  102. 3193231538989185
  103. 3852971941960065
  104. 4267414072422865
  105. 5071419883911745
  106. 5560847142362785
  107. 6128613921672705
  108. 7969193810377345
  109. 10370544300097345
  110. 11947816523586945
  111. 12301576752408945
  112. 13462627333098945
  113. 16057190782234785
  114. 21904922730392065
  115. 26708253318968145
  116. 31454143858820145
  117. 52451051117566465
  118. 53502800037708385
  119. 60977817398996785
  120. 63469753433627905
  121. 75131642415974145
  122. 87966790826033665
  123. 101817952350880305
  124. 105377502147819265
  125. 107185855472437105
  126. 113536268221379185
  127. 113781032286917905
  128. 121719617715279585
  129. 122261077742785345
  130. 135964101464432605
  131. 171800042106877185
  132. 174946667573899585
  133. 179493438867778945
  134. 180950795673242145
  135. 184882434303977985
  136. 186028123265766145
  137. 202703674932642885
  138. 207897332231129905
  139. 218395364809636785
  140. 220603906882380385
  141. 220629545178715905
  142. 224811969542371905
  143. 242746392608701345
  144. 292561243007134465
  145. 337671763976568385
  146. 345319996314995065
  147. 371717470696746945
  148. 446716685971537345
  149. 462937246809774465
  150. 585360421327663345
  151. 606721976846086465
  152. 624654503446875265
  153. 755739161051095585
  154. 781506906921758865
  155. 897396829788960865
  156. 916541202603455265
  157. 932542749584987905
  158. 947087538769733505
  159. 978502408508476545
  160. ].uniq.sort
  161. var b = [
  162. 1886616373665, 23239986511105, 881715504450705, 3193231538989185, 6128613921672705, 10370544300097345, 11947816523586945, 12301576752408945, 31454143858820145, 75131642415974145, 101817952350880305, 107185855472437105, 171800042106877185, 180950795673242145, 202703674932642885, 220629545178715905, 224811969542371905, 462937246809774465, 781506906921758865, 897396829788960865, 916541202603455265, 947087538769733505, 1114090435183709985, 1177908521713261185, 1975738329089659905, 2325605354531800705, 2435279804766429985, 2978083893490447105, 3239491439742459265, 4448649946743866865, 4682767770306971745, 5092454223678002545, 5473367446820370945, 7258662586258825105, 7626934787034827265, 9718427839254355585, 10186034230241370865, 13294381492836585505
  163. ]
  164. #var f = [3, 5, 17, 23, 89].prod
  165. #var f = [5, 7, 13, 17, 19, 23].prod
  166. #var f = [3, 5, 17, 23, 29, 53, 89].prod
  167. var f = [3, 5, 17, 23, 29, 53, 89, 197].prod
  168. var v = b.map{.factor}.flat.uniq.grep{_ < 1e6 }.grep{gcd(_, f) == 1 }.lcm
  169. say f.sigma/f
  170. say f
  171. say v.sigma0.log10
  172. say v