1 | module math |
2 | |
3 | const ( |
4 | log_sqrt_2pi = 9.18938533204672741780329736e-1 |
5 | bernoulli = [ |
6 | /* |
7 | Bernoulli numbers B(2),B(4),B(6),...,B(20). Only B(2),...,B(10) currently |
8 | * used. |
9 | */ |
10 | 1.0 / (6.0 * 2.0 * 1.0), |
11 | -1.0 / (30.0 * 4.0 * 3.0), |
12 | 1.0 / (42.0 * 6.0 * 5.0), |
13 | -1.0 / (30.0 * 8.0 * 7.0), |
14 | 5.0 / (66.0 * 10.0 * 9.0), |
15 | -691.0 / (2730.0 * 12.0 * 11.0), |
16 | 7.0 / (6.0 * 14.0 * 13.0), |
17 | -3617.0 / (510.0 * 16.0 * 15.0), |
18 | 43867.0 / (796.0 * 18.0 * 17.0), |
19 | -174611.0 / (330.0 * 20.0 * 19.0), |
20 | ] |
21 | factorials_table = [ |
22 | // 0! |
23 | 1.000000000000000000000e+0, |
24 | // 1! |
25 | 1.000000000000000000000e+0, |
26 | // 2! |
27 | 2.000000000000000000000e+0, |
28 | // 3! |
29 | 6.000000000000000000000e+0, |
30 | // 4! |
31 | 2.400000000000000000000e+1, |
32 | // 5! |
33 | 1.200000000000000000000e+2, |
34 | // 6! |
35 | 7.200000000000000000000e+2, |
36 | // 7! |
37 | 5.040000000000000000000e+3, |
38 | // 8! |
39 | 4.032000000000000000000e+4, |
40 | // 9! |
41 | 3.628800000000000000000e+5, |
42 | // 10! |
43 | 3.628800000000000000000e+6, |
44 | // 11! |
45 | 3.991680000000000000000e+7, |
46 | // 12! |
47 | 4.790016000000000000000e+8, |
48 | // 13! |
49 | 6.227020800000000000000e+9, |
50 | // 14! |
51 | 8.717829120000000000000e+10, |
52 | // 15! |
53 | 1.307674368000000000000e+12, |
54 | // 16! |
55 | 2.092278988800000000000e+13, |
56 | // 17! |
57 | 3.556874280960000000000e+14, |
58 | // 18! |
59 | 6.402373705728000000000e+15, |
60 | // 19! |
61 | 1.216451004088320000000e+17, |
62 | // 20! |
63 | 2.432902008176640000000e+18, |
64 | // 21! |
65 | 5.109094217170944000000e+19, |
66 | // 22! |
67 | 1.124000727777607680000e+21, |
68 | // 23! |
69 | 2.585201673888497664000e+22, |
70 | // 24! |
71 | 6.204484017332394393600e+23, |
72 | // 25! |
73 | 1.551121004333098598400e+25, |
74 | // 26! |
75 | 4.032914611266056355840e+26, |
76 | // 27! |
77 | 1.088886945041835216077e+28, |
78 | // 28! |
79 | 3.048883446117138605015e+29, |
80 | // 29! |
81 | 8.841761993739701954544e+30, |
82 | // 30! |
83 | 2.652528598121910586363e+32, |
84 | // 31! |
85 | 8.222838654177922817726e+33, |
86 | // 32! |
87 | 2.631308369336935301672e+35, |
88 | // 33! |
89 | 8.683317618811886495518e+36, |
90 | // 34! |
91 | 2.952327990396041408476e+38, |
92 | // 35! |
93 | 1.033314796638614492967e+40, |
94 | // 36! |
95 | 3.719933267899012174680e+41, |
96 | // 37! |
97 | 1.376375309122634504632e+43, |
98 | // 38! |
99 | 5.230226174666011117600e+44, |
100 | // 39! |
101 | 2.039788208119744335864e+46, |
102 | // 40! |
103 | 8.159152832478977343456e+47, |
104 | // 41! |
105 | 3.345252661316380710817e+49, |
106 | // 42! |
107 | 1.405006117752879898543e+51, |
108 | // 43! |
109 | 6.041526306337383563736e+52, |
110 | // 44! |
111 | 2.658271574788448768044e+54, |
112 | // 45! |
113 | 1.196222208654801945620e+56, |
114 | // 46! |
115 | 5.502622159812088949850e+57, |
116 | // 47! |
117 | 2.586232415111681806430e+59, |
118 | // 48! |
119 | 1.241391559253607267086e+61, |
120 | // 49! |
121 | 6.082818640342675608723e+62, |
122 | // 50! |
123 | 3.041409320171337804361e+64, |
124 | // 51! |
125 | 1.551118753287382280224e+66, |
126 | // 52! |
127 | 8.065817517094387857166e+67, |
128 | // 53! |
129 | 4.274883284060025564298e+69, |
130 | // 54! |
131 | 2.308436973392413804721e+71, |
132 | // 55! |
133 | 1.269640335365827592597e+73, |
134 | // 56! |
135 | 7.109985878048634518540e+74, |
136 | // 57! |
137 | 4.052691950487721675568e+76, |
138 | // 58! |
139 | 2.350561331282878571829e+78, |
140 | // 59! |
141 | 1.386831185456898357379e+80, |
142 | // 60! |
143 | 8.320987112741390144276e+81, |
144 | // 61! |
145 | 5.075802138772247988009e+83, |
146 | // 62! |
147 | 3.146997326038793752565e+85, |
148 | // 63! |
149 | 1.982608315404440064116e+87, |
150 | // 64! |
151 | 1.268869321858841641034e+89, |
152 | // 65! |
153 | 8.247650592082470666723e+90, |
154 | // 66! |
155 | 5.443449390774430640037e+92, |
156 | // 67! |
157 | 3.647111091818868528825e+94, |
158 | // 68! |
159 | 2.480035542436830599601e+96, |
160 | // 69! |
161 | 1.711224524281413113725e+98, |
162 | // 70! |
163 | 1.197857166996989179607e+100, |
164 | // 71! |
165 | 8.504785885678623175212e+101, |
166 | // 72! |
167 | 6.123445837688608686152e+103, |
168 | // 73! |
169 | 4.470115461512684340891e+105, |
170 | // 74! |
171 | 3.307885441519386412260e+107, |
172 | // 75! |
173 | 2.480914081139539809195e+109, |
174 | // 76! |
175 | 1.885494701666050254988e+111, |
176 | // 77! |
177 | 1.451830920282858696341e+113, |
178 | // 78! |
179 | 1.132428117820629783146e+115, |
180 | // 79! |
181 | 8.946182130782975286851e+116, |
182 | // 80! |
183 | 7.156945704626380229481e+118, |
184 | // 81! |
185 | 5.797126020747367985880e+120, |
186 | // 82! |
187 | 4.753643337012841748421e+122, |
188 | // 83! |
189 | 3.945523969720658651190e+124, |
190 | // 84! |
191 | 3.314240134565353266999e+126, |
192 | // 85! |
193 | 2.817104114380550276949e+128, |
194 | // 86! |
195 | 2.422709538367273238177e+130, |
196 | // 87! |
197 | 2.107757298379527717214e+132, |
198 | // 88! |
199 | 1.854826422573984391148e+134, |
200 | // 89! |
201 | 1.650795516090846108122e+136, |
202 | // 90! |
203 | 1.485715964481761497310e+138, |
204 | // 91! |
205 | 1.352001527678402962552e+140, |
206 | // 92! |
207 | 1.243841405464130725548e+142, |
208 | // 93! |
209 | 1.156772507081641574759e+144, |
210 | // 94! |
211 | 1.087366156656743080274e+146, |
212 | // 95! |
213 | 1.032997848823905926260e+148, |
214 | // 96! |
215 | 9.916779348709496892096e+149, |
216 | // 97! |
217 | 9.619275968248211985333e+151, |
218 | // 98! |
219 | 9.426890448883247745626e+153, |
220 | // 99! |
221 | 9.332621544394415268170e+155, |
222 | // 100! |
223 | 9.332621544394415268170e+157, |
224 | // 101! |
225 | 9.425947759838359420852e+159, |
226 | // 102! |
227 | 9.614466715035126609269e+161, |
228 | // 103! |
229 | 9.902900716486180407547e+163, |
230 | // 104! |
231 | 1.029901674514562762385e+166, |
232 | // 105! |
233 | 1.081396758240290900504e+168, |
234 | // 106! |
235 | 1.146280563734708354534e+170, |
236 | // 107! |
237 | 1.226520203196137939352e+172, |
238 | // 108! |
239 | 1.324641819451828974500e+174, |
240 | // 109! |
241 | 1.443859583202493582205e+176, |
242 | // 110! |
243 | 1.588245541522742940425e+178, |
244 | // 111! |
245 | 1.762952551090244663872e+180, |
246 | // 112! |
247 | 1.974506857221074023537e+182, |
248 | // 113! |
249 | 2.231192748659813646597e+184, |
250 | // 114! |
251 | 2.543559733472187557120e+186, |
252 | // 115! |
253 | 2.925093693493015690688e+188, |
254 | // 116! |
255 | 3.393108684451898201198e+190, |
256 | // 117! |
257 | 3.969937160808720895402e+192, |
258 | // 118! |
259 | 4.684525849754290656574e+194, |
260 | // 119! |
261 | 5.574585761207605881323e+196, |
262 | // 120! |
263 | 6.689502913449127057588e+198, |
264 | // 121! |
265 | 8.094298525273443739682e+200, |
266 | // 122! |
267 | 9.875044200833601362412e+202, |
268 | // 123! |
269 | 1.214630436702532967577e+205, |
270 | // 124! |
271 | 1.506141741511140879795e+207, |
272 | // 125! |
273 | 1.882677176888926099744e+209, |
274 | // 126! |
275 | 2.372173242880046885677e+211, |
276 | // 127! |
277 | 3.012660018457659544810e+213, |
278 | // 128! |
279 | 3.856204823625804217357e+215, |
280 | // 129! |
281 | 4.974504222477287440390e+217, |
282 | // 130! |
283 | 6.466855489220473672507e+219, |
284 | // 131! |
285 | 8.471580690878820510985e+221, |
286 | // 132! |
287 | 1.118248651196004307450e+224, |
288 | // 133! |
289 | 1.487270706090685728908e+226, |
290 | // 134! |
291 | 1.992942746161518876737e+228, |
292 | // 135! |
293 | 2.690472707318050483595e+230, |
294 | // 136! |
295 | 3.659042881952548657690e+232, |
296 | // 137! |
297 | 5.012888748274991661035e+234, |
298 | // 138! |
299 | 6.917786472619488492228e+236, |
300 | // 139! |
301 | 9.615723196941089004197e+238, |
302 | // 140! |
303 | 1.346201247571752460588e+241, |
304 | // 141! |
305 | 1.898143759076170969429e+243, |
306 | // 142! |
307 | 2.695364137888162776589e+245, |
308 | // 143! |
309 | 3.854370717180072770522e+247, |
310 | // 144! |
311 | 5.550293832739304789551e+249, |
312 | // 145! |
313 | 8.047926057471991944849e+251, |
314 | // 146! |
315 | 1.174997204390910823948e+254, |
316 | // 147! |
317 | 1.727245890454638911203e+256, |
318 | // 148! |
319 | 2.556323917872865588581e+258, |
320 | // 149! |
321 | 3.808922637630569726986e+260, |
322 | // 150! |
323 | 5.713383956445854590479e+262, |
324 | // 151! |
325 | 8.627209774233240431623e+264, |
326 | // 152! |
327 | 1.311335885683452545607e+267, |
328 | // 153! |
329 | 2.006343905095682394778e+269, |
330 | // 154! |
331 | 3.089769613847350887959e+271, |
332 | // 155! |
333 | 4.789142901463393876336e+273, |
334 | // 156! |
335 | 7.471062926282894447084e+275, |
336 | // 157! |
337 | 1.172956879426414428192e+278, |
338 | // 158! |
339 | 1.853271869493734796544e+280, |
340 | // 159! |
341 | 2.946702272495038326504e+282, |
342 | // 160! |
343 | 4.714723635992061322407e+284, |
344 | // 161! |
345 | 7.590705053947218729075e+286, |
346 | // 162! |
347 | 1.229694218739449434110e+289, |
348 | // 163! |
349 | 2.004401576545302577600e+291, |
350 | // 164! |
351 | 3.287218585534296227263e+293, |
352 | // 165! |
353 | 5.423910666131588774984e+295, |
354 | // 166! |
355 | 9.003691705778437366474e+297, |
356 | // 167! |
357 | 1.503616514864999040201e+300, |
358 | // 168! |
359 | 2.526075744973198387538e+302, |
360 | // 169! |
361 | 4.269068009004705274939e+304, |
362 | // 170! |
363 | 7.257415615307998967397e+306, |
364 | ] |
365 | log_factorials_table = [ |
366 | // 0! |
367 | 0.000000000000000000000e+0, |
368 | // 1! |
369 | 0.000000000000000000000e+0, |
370 | // 2! |
371 | 6.931471805599453094172e-1, |
372 | // 3! |
373 | 1.791759469228055000812e+0, |
374 | // 4! |
375 | 3.178053830347945619647e+0, |
376 | // 5! |
377 | 4.787491742782045994248e+0, |
378 | // 6! |
379 | 6.579251212010100995060e+0, |
380 | // 7! |
381 | 8.525161361065414300166e+0, |
382 | // 8! |
383 | 1.060460290274525022842e+1, |
384 | // 9! |
385 | 1.280182748008146961121e+1, |
386 | // 10! |
387 | 1.510441257307551529523e+1, |
388 | // 11! |
389 | 1.750230784587388583929e+1, |
390 | // 12! |
391 | 1.998721449566188614952e+1, |
392 | // 13! |
393 | 2.255216385312342288557e+1, |
394 | // 14! |
395 | 2.519122118273868150009e+1, |
396 | // 15! |
397 | 2.789927138384089156609e+1, |
398 | // 16! |
399 | 3.067186010608067280376e+1, |
400 | // 17! |
401 | 3.350507345013688888401e+1, |
402 | // 18! |
403 | 3.639544520803305357622e+1, |
404 | // 19! |
405 | 3.933988418719949403622e+1, |
406 | // 20! |
407 | 4.233561646075348502966e+1, |
408 | // 21! |
409 | 4.538013889847690802616e+1, |
410 | // 22! |
411 | 4.847118135183522387964e+1, |
412 | // 23! |
413 | 5.160667556776437357045e+1, |
414 | // 24! |
415 | 5.478472939811231919009e+1, |
416 | // 25! |
417 | 5.800360522298051993929e+1, |
418 | // 26! |
419 | 6.126170176100200198477e+1, |
420 | // 27! |
421 | 6.455753862700633105895e+1, |
422 | // 28! |
423 | 6.788974313718153498289e+1, |
424 | // 29! |
425 | 7.125703896716800901007e+1, |
426 | // 30! |
427 | 7.465823634883016438549e+1, |
428 | // 31! |
429 | 7.809222355331531063142e+1, |
430 | // 32! |
431 | 8.155795945611503717850e+1, |
432 | // 33! |
433 | 8.505446701758151741396e+1, |
434 | // 34! |
435 | 8.858082754219767880363e+1, |
436 | // 35! |
437 | 9.213617560368709248333e+1, |
438 | // 36! |
439 | 9.571969454214320248496e+1, |
440 | // 37! |
441 | 9.933061245478742692933e+1, |
442 | // 38! |
443 | 1.029681986145138126988e+2, |
444 | // 39! |
445 | 1.066317602606434591262e+2, |
446 | // 40! |
447 | 1.103206397147573954291e+2, |
448 | // 41! |
449 | 1.140342117814617032329e+2, |
450 | // 42! |
451 | 1.177718813997450715388e+2, |
452 | // 43! |
453 | 1.215330815154386339623e+2, |
454 | // 44! |
455 | 1.253172711493568951252e+2, |
456 | // 45! |
457 | 1.291239336391272148826e+2, |
458 | // 46! |
459 | 1.329525750356163098828e+2, |
460 | // 47! |
461 | 1.368027226373263684696e+2, |
462 | // 48! |
463 | 1.406739236482342593987e+2, |
464 | // 49! |
465 | 1.445657439463448860089e+2, |
466 | // 50! |
467 | 1.484777669517730320675e+2, |
468 | // 51! |
469 | 1.524095925844973578392e+2, |
470 | // 52! |
471 | 1.563608363030787851941e+2, |
472 | // 53! |
473 | 1.603311282166309070282e+2, |
474 | // 54! |
475 | 1.643201122631951814118e+2, |
476 | // 55! |
477 | 1.683274454484276523305e+2, |
478 | // 56! |
479 | 1.723527971391628015638e+2, |
480 | // 57! |
481 | 1.763958484069973517152e+2, |
482 | // 58! |
483 | 1.804562914175437710518e+2, |
484 | // 59! |
485 | 1.845338288614494905025e+2, |
486 | // 60! |
487 | 1.886281734236715911873e+2, |
488 | // 61! |
489 | 1.927390472878449024360e+2, |
490 | // 62! |
491 | 1.968661816728899939914e+2, |
492 | // 63! |
493 | 2.010093163992815266793e+2, |
494 | // 64! |
495 | 2.051681994826411985358e+2, |
496 | // 65! |
497 | 2.093425867525368356464e+2, |
498 | // 66! |
499 | 2.135322414945632611913e+2, |
500 | // 67! |
501 | 2.177369341139542272510e+2, |
502 | // 68! |
503 | 2.219564418191303339501e+2, |
504 | // 69! |
505 | 2.261905483237275933323e+2, |
506 | // 70! |
507 | 2.304390435657769523214e+2, |
508 | // 71! |
509 | 2.347017234428182677427e+2, |
510 | // 72! |
511 | 2.389783895618343230538e+2, |
512 | // 73! |
513 | 2.432688490029827141829e+2, |
514 | // 74! |
515 | 2.475729140961868839366e+2, |
516 | // 75! |
517 | 2.518904022097231943772e+2, |
518 | // 76! |
519 | 2.562211355500095254561e+2, |
520 | // 77! |
521 | 2.605649409718632093053e+2, |
522 | // 78! |
523 | 2.649216497985528010421e+2, |
524 | // 79! |
525 | 2.692910976510198225363e+2, |
526 | // 80! |
527 | 2.736731242856937041486e+2, |
528 | // 81! |
529 | 2.780675734403661429141e+2, |
530 | // 82! |
531 | 2.824742926876303960274e+2, |
532 | // 83! |
533 | 2.868931332954269939509e+2, |
534 | // 84! |
535 | 2.913239500942703075662e+2, |
536 | // 85! |
537 | 2.957666013507606240211e+2, |
538 | // 86! |
539 | 3.002209486470141317540e+2, |
540 | // 87! |
541 | 3.046868567656687154726e+2, |
542 | // 88! |
543 | 3.091641935801469219449e+2, |
544 | // 89! |
545 | 3.136528299498790617832e+2, |
546 | // 90! |
547 | 3.181526396202093268500e+2, |
548 | // 91! |
549 | 3.226634991267261768912e+2, |
550 | // 92! |
551 | 3.271852877037752172008e+2, |
552 | // 93! |
553 | 3.317178871969284731381e+2, |
554 | // 94! |
555 | 3.362611819791984770344e+2, |
556 | // 95! |
557 | 3.408150588707990178690e+2, |
558 | // 96! |
559 | 3.453794070622668541074e+2, |
560 | // 97! |
561 | 3.499541180407702369296e+2, |
562 | // 98! |
563 | 3.545390855194408088492e+2, |
564 | // 99! |
565 | 3.591342053695753987760e+2, |
566 | // 100! |
567 | 3.637393755555634901441e+2, |
568 | // 101! |
569 | 3.683544960724047495950e+2, |
570 | // 102! |
571 | 3.729794688856890206760e+2, |
572 | // 103! |
573 | 3.776141978739186564468e+2, |
574 | // 104! |
575 | 3.822585887730600291111e+2, |
576 | // 105! |
577 | 3.869125491232175524822e+2, |
578 | // 106! |
579 | 3.915759882173296196258e+2, |
580 | // 107! |
581 | 3.962488170517915257991e+2, |
582 | // 108! |
583 | 4.009309482789157454921e+2, |
584 | // 109! |
585 | 4.056222961611448891925e+2, |
586 | // 110! |
587 | 4.103227765269373054205e+2, |
588 | // 111! |
589 | 4.150323067282496395563e+2, |
590 | // 112! |
591 | 4.197508055995447340991e+2, |
592 | // 113! |
593 | 4.244781934182570746677e+2, |
594 | // 114! |
595 | 4.292143918666515701285e+2, |
596 | // 115! |
597 | 4.339593239950148201939e+2, |
598 | // 116! |
599 | 4.387129141861211848399e+2, |
600 | // 117! |
601 | 4.434750881209189409588e+2, |
602 | // 118! |
603 | 4.482457727453846057188e+2, |
604 | // 119! |
605 | 4.530248962384961351041e+2, |
606 | // 120! |
607 | 4.578123879812781810984e+2, |
608 | // 121! |
609 | 4.626081785268749221865e+2, |
610 | // 122! |
611 | 4.674121995716081787447e+2, |
612 | // 123! |
613 | 4.722243839269805962399e+2, |
614 | // 124! |
615 | 4.770446654925856331047e+2, |
616 | // 125! |
617 | 4.818729792298879342285e+2, |
618 | // 126! |
619 | 4.867092611368394122258e+2, |
620 | // 127! |
621 | 4.915534482232980034989e+2, |
622 | // 128! |
623 | 4.964054784872176206648e+2, |
624 | // 129! |
625 | 5.012652908915792927797e+2, |
626 | // 130! |
627 | 5.061328253420348751997e+2, |
628 | // 131! |
629 | 5.110080226652360267439e+2, |
630 | // 132! |
631 | 5.158908245878223975982e+2, |
632 | // 133! |
633 | 5.207811737160441513633e+2, |
634 | // 134! |
635 | 5.256790135159950627324e+2, |
636 | // 135! |
637 | 5.305842882944334921812e+2, |
638 | // 136! |
639 | 5.354969431801695441897e+2, |
640 | // 137! |
641 | 5.404169241059976691050e+2, |
642 | // 138! |
643 | 5.453441777911548737966e+2, |
644 | // 139! |
645 | 5.502786517242855655538e+2, |
646 | // 140! |
647 | 5.552202941468948698523e+2, |
648 | // 141! |
649 | 5.601690540372730381305e+2, |
650 | // 142! |
651 | 5.651248810948742988613e+2, |
652 | // 143! |
653 | 5.700877257251342061414e+2, |
654 | // 144! |
655 | 5.750575390247102067619e+2, |
656 | // 145! |
657 | 5.800342727671307811636e+2, |
658 | // 146! |
659 | 5.850178793888391176022e+2, |
660 | // 147! |
661 | 5.900083119756178539038e+2, |
662 | // 148! |
663 | 5.950055242493819689670e+2, |
664 | // 149! |
665 | 6.000094705553274281080e+2, |
666 | // 150! |
667 | 6.050201058494236838580e+2, |
668 | // 151! |
669 | 6.100373856862386081868e+2, |
670 | // 152! |
671 | 6.150612662070848845750e+2, |
672 | // 153! |
673 | 6.200917041284773200381e+2, |
674 | // 154! |
675 | 6.251286567308909491967e+2, |
676 | // 155! |
677 | 6.301720818478101958172e+2, |
678 | // 156! |
679 | 6.352219378550597328635e+2, |
680 | // 157! |
681 | 6.402781836604080409209e+2, |
682 | // 158! |
683 | 6.453407786934350077245e+2, |
684 | // 159! |
685 | 6.504096828956552392500e+2, |
686 | // 160! |
687 | 6.554848567108890661717e+2, |
688 | // 161! |
689 | 6.605662610758735291676e+2, |
690 | // 162! |
691 | 6.656538574111059132426e+2, |
692 | // 163! |
693 | 6.707476076119126755767e+2, |
694 | // 164! |
695 | 6.758474740397368739994e+2, |
696 | // 165! |
697 | 6.809534195136374546094e+2, |
698 | // 166! |
699 | 6.860654073019939978423e+2, |
700 | // 167! |
701 | 6.911834011144107529496e+2, |
702 | // 168! |
703 | 6.963073650938140118743e+2, |
704 | // 169! |
705 | 7.014372638087370853465e+2, |
706 | // 170! |
707 | 7.065730622457873471107e+2, |
708 | // 171! |
709 | 7.117147258022900069535e+2, |
710 | ] |
711 | ) |