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