150728(2)

Ruby


ラマヌジャン予想(1)

先ほどのコードを使って、|τ(p)| / p^5.5 の最大値を記録してみる。

require 'prime'

# m次以下を取り出す

def mul(f_ary, b_ary, m)

  s1, s2 = f_ary.size, b_ary.size

  ary = Array.new(s1 + s2 - 1, 0)

  (0..s1 - 1).each{|i|

    (0..s2 - 1).each{|j|

      ary[i + j] += f_ary[i] * b_ary[j]

    }

  }

  ary[0..m]

end

# m次以下を取り出す

def power(ary, n, m)

  return [1] if n == 0

  mul(ary, power(ary, n - 1, m), m)

end

def A000594(n)

  ary = [1]

  # 無限積のうち必要なところだけ取り出す

  (1..n - 1).each{|i|

    b_ary = Array.new(i + 1, 0)

    b_ary[0], b_ary[-1] = 1, -1

    ary = mul(ary, b_ary, n)

  }

  # 24乗してx倍

  power(ary, 24, n).unshift(0)[1..n]

end

# |t(p)| / p^5.5 の最大値

def ramanujan_conjecture(ary)

  max = 0

  r_max_ary = []

  p_ary = Prime.each(ary.size).to_a

  p_ary.each{|p|

    t = ary[p - 1]

    u = t.abs / p ** 5.5

    if max < u

      r_max_ary.push([p, t, u])

      max = u

    end

  }

  r_max_ary

end

ary = A000594(1000)

p ary

p ramanujan_conjecture(ary)

実行結果

[1, -24, 252, -1472, 4830, -6048, -16744, 84480, -113643, -115920, 534612, -3709

44, -577738, 401856, 1217160, 987136, -6905934, 2727432, 10661420, -7109760, -42

19488, -12830688, 18643272, 21288960, -25499225, 13865712, -73279080, 24647168,

128406630, -29211840, -52843168, -196706304, 134722224, 165742416, -80873520, 16

7282496, -182213314, -255874080, -145589976, 408038400, 308120442, 101267712, -1

7125708, -786948864, -548895690, -447438528, 2687348496, 248758272, -1696965207,

 611981400, -1740295368, 850430336, -1596055698, 1758697920, 2582175960, -141453

3120, 2686677840, -3081759120, -5189203740, -1791659520, 6956478662, 1268236032,

 1902838392, 2699296768, -2790474540, -3233333376, -15481826884, 10165534848, 46

98104544, 1940964480, 9791485272, -9600560640, 1463791322, 4373119536, -64258047

00, -15693610240, -8951543328, 3494159424, 38116845680, 4767866880, 1665188361,

-7394890608, -29335099668, 6211086336, -33355661220, 411016992, 32358470760, 451

64021760, -24992917110, 13173496560, 9673645072, -27442896384, -13316478336, -64

496363904, 51494658600, -49569988608, 75013568546, 40727164968, -60754911516, 37

534859200, 81742959102, 41767088832, -225755128648, -48807306240, -20380127040,

38305336752, 90241258356, 107866805760, 73482676310, -61972223040, -45917755128,

 -16528605184, -85146862638, -64480268160, 90047003760, -189014559360, 656558795

34, 124540889760, 115632958896, 102825676800, 498319933, -166955487888, 77646351

384, 77785143296, -359001100500, -45668121408, -262717201024, 338071388160, -431

5678416, 66971388960, 631528759932, -198311113728, -178514816480, 371563845216,

-353937956400, -583413304320, -297198746214, -112754509056, 596793577940, 119045

821440, 677211820992, -234995646528, -308865667656, -112181096448, 620204022900,

 -35130991728, -427635232164, 268217998208, -1115433620850, 154219312800, -82444

7297848, 900676761600, 784811057562, 214837039872, -255232501440, 214308444672,

1315116754406, -914804296320, -402206035896, -950091448320, -312162946368, -3996

4520664, -357832759588, -453553290624, 650708341920, 704042392032, 2754833892216

, -356462346240, -1458379197393, 800535869280, -1211595753060, 25209042176, -950

387449578, -776603298240, 426959023400, 527734751232, -1307679342480, 5998300106

40, 1681384224780, 807974455680, -996774496018, -232167481728, 1753032622824, 15

74983618560, -880090306620, 319595480064, -3691995187608, -3955776986112, 122698

4915520, -1235871806400, 2762403350592, 680222785536, 5442387685442, -1800325645

104, -703199584080, 2497932784704, -2876091504354, 1458117876384, 728391402200,

-2154174528000, -3901420374768, -1961831018448, -2150040612720, 2561714781696, 1

488221734860, 5418123087552, -2118677359896, -570305978368, 5699723069040, 48912

3048960, -6793168439188, 2349393987456, 2467454288544, -2165790200544, -82717169

640, -6190616678400, 884806004992, -1763584231440, 368875413144, -3800963013120,

 3989820497292, 1102026123072, 7334863021472, 3293650354176, 2897808426675, 2043

524703312, -1359839565924, -3954789780480, -11824411223170, -2161128090240, -225

5788918656, 10847792102400, -17563353448518, -1575741108816, 12979893235680, 763

8507905280, 9605445111360, -2775191013504, -7139577462960, 1201502453760, -23130

6909358, -11959678392, 13400796651732, -10239936590464, -8196341949810, -1863512

433216, -6159507467960, -4464190832640, -7392445116336, 8616026412000, 129830535

45252, -2800978113024, 9966916930464, 6305212824576, -8405626627440, -1364187309

6704, 23961192565506, 103576281984, 3050979729616, 4107578522880, -1459251465309

0, -15156690238368, -24273728464488, 11381333483520, -7708949021340, 42843555955

20, -6298215111720, 22789249173248, 25837706543670, 8494510953600, -376793236052

8, -6817096065024, 2437758558144, 7132769909136, -13632191675700, -6915609888768

, -16418932005874, -14323045870560, 6005256141024, -6832194969600, 2103572290708

2, -16253083703808, 16713176326532, -14413066320384, 12976653967200, 74127760237

44, -5159168680848, 22354294505472, 13420028104723, -14884896549600, 18903419273

592, -2154700825984, -23926858987458, 10263245571936, -25063854064200, -15393380

766720, -39175875516960, 26770406900400, -10770926678736, 9458784518400, 2867528

54752, 19786735148352, 20599225693704, 10524271493120, 33599791937460, -18835465

381488, 15311092828556, 13176671778816, -56890292419296, 6125580034560, 49875160

575912, -12299441172480, -99480832756438, -31562802105744, 9190709433360, -56107

996840960, 83369248359366, 9652944861504, 68647725277560, 13037603389440, 227407

97105712, 7491910712832, -73627062866280, -2451157267392, 14731871253050, 858798

6230112, 18517634430120, 26030014940160, -44996963217024, -15617000206080, -6358

4021925868, 43181266711296, 20707267642902, -66116013413184, -74777223849720, -4

165208506368, 121001428335986, 35001100737432, -21457009384776, 49099533315840,

-28250591730816, 29078298073440, 61522344410800, -1446779811840, 22691844947520,

 22809298789872, -155661561078204, -47631668958720, -25643022194650, -1024701656

1600, 42336109121040, -105161550594048, 24909815245602, 31384304219520, 47292873

863760, 36789573985920, 29139505641792, -40353221394720, 157584150853560, -46370

707891200, -2824382481819, 23922587904432, 125576623116, -14239605545984, 707011

2085260, -42072782947776, -177901220129584, 18403444948992, -35015731390206, 211

22167358880, 26724356607312, 19601856110592, -55161734023378, 88607884502592, -9

0468277326000, 227027200942080, -74185389602940, -29447637972480, 14646311632298

0, -75800137459200, -66204734658048, -66297680414208, 231449571733632, 851939898

16320, -43235954274240, -130617304450608, 1946216834244, -110419972899712, -1498

71571611810, 16876790017920, -128749205976048, -143359620687360, 159145247502864

, 69026196104496, 184104364634400, 89431229751552, 208110680273846, -17481393652

800, -44985733752960, -25171202969600, -133407937691598, 93634088994432, 3052950

6193984, -120325635798144, 8042859783630, 51600974705280, -97413424224168, -1470

20152688640, -206167580638390, -35717321636640, -74894084045928, 332311549369856

, 86888027422560, 50848256637504, -141688531396440, 113644706660352, 15039198164

0880, -136793353656960, 73403515193820, 29999547002880, 171111932338622, 1630360

42540512, -305398345130928, -134834785367040, 176095964901150, -59218902925056,

-116479278716528, -132835132300032, -77834148249312, 1985212071360, -71775829446

768, -72336417914880, 99881248225682, -21235344119808, 156291413770800, -1081664

99528320, 198763752966240, -8853009915456, -29031220908760, 218142225100800, 192

848217019101, -95755691935008, 328369848718692, 67590935548416, -120715789641300

, -176036712515328, -281089272454200, -45197025083392, -612368143631550, -695474

02240200, 164724885738504, 125336181803136, -207760719057696, 32636149582176, 46

723705697760, 226970543923200, 303483032911706, 283785869356080, 506060490060720

, -132549189534720, -729307946668938, 54138934047744, 122188164073712, 126754807

111680, -64318590362880, 421520482764432, -617380683662484, -96645454674048, 259

227709345696, -311517437656320, 331409422110312, -438383931955200, -915560900529

6, -230530682672640, -271857947399500, -170211715494912, 181380557687814, 171349

859111040, 1050837984850080, -239423044976640, 105271555603732, 5551365824592, -

78665062484736, -733526941376, 362315536077180, -321619119641568, -2199099717618

64, 587683317365760, -90173855416176, 196712206795440, -483863128068108, -114295

429237248, -886767711942420, 147828179231040, -293446222622280, -52163393486848,

 -163948629394368, 177418682792064, -108877719272500, 528449619936000, 694218140

838432, -311593285086048, 506588355787752, 160751787356160, 394818492462660, -23

9206006331136, -367511557743036, 386719719907328, 85753393288710, 20173503905856

0, -24509721895568, -364965248630784, -781259049093600, -575068621572144, -10903

97271369840, 6352678628352, 1436688754143552, -73223513510784, -239497637293656,

 -235739289139200, 927574652509722, 350220351674160, -21818651341228, -929610334

619904, 107593673896800, 582569483147712, 364931430558912, 132989157310464, -605

238167047943, 185014776512160, 589716680624820, 262773809858560, -17801288792019

6, 151157162681280, 435865277859480, -1307904735160320, 423708824644560, -620104

957048080, -907217963244684, 520996671820800, -1695266465052058, 90430376652672,

 -251187172996536, 1358440752807936, 354921326577300, -58506205395456, 752144751

323996, 437476554427008, -790555104585666, 327172600216800, 1368997013214600, 39

6895871877120, -638228464065920, 394054368140976, -221782757268240, -87848014672

7680, 187489032934806, -144126147384576, 9894172288504, -79833163038720, -930382

787277216, -504857349769968, 244971285916812, -996855800500224, -411259346541540

, -401116231836768, -27881913916584, 827184675778560, 1352425875836970, -3114396

95212800, 1432228477996772, 454650262789632, 696125644349184, 123820048340352, -

475388987464200, -306756182605824, -877659324192574, -322080674513352, 137148169

6731384, -912940321708800, 491186908840992, -453682062566208, -853270528819176,

123661090882560, 317117898149220, 574244615698992, -2434250501782764, 6294790617

45408, -563383208178560, 601532497540800, -724775059097208, -179869321928704, -3

03318212278158, 940221012407040, 558507191467680, 1641918289891200, 183554633354

400, 258502240289664, -1701984261329400, -542851981056000, 2339216023326602, -68

82068514048, 1759401252578412, 1213586422432256, 2406885276390, -494381416648896

, -2496072486898144, -2097168523591680, -541810234405440, -806395006499040, -155

2583345382048, -1155241876731264, 2473009012516862, -367466227885344, 3750318771

84720, -756226380349440, 24336777156666, 1365367018063104, 4225450805148020, 375

702242119680, -1366161820349760, -1197003853821888, 418481404089840, -1437171065

48736, -488895969711875, 2387539986154512, 1436330213398080, -1935851862485632,

1258353120405276, -220577026400640, -4263261111420568, 3220111123046400, -171187

8446675376, -2000861960624784, -1268924080945920, 592047284838912, 9804012847617

66, -1647545406661440, -1112733760765896, 1632884804812800, 1008295928583042, -5

45779130537088, 303982416160892, 459503857053696, -20844726749280, 1767049508790

720, 3435828424131096, 140675112737280, -2774210589848880, -353564910073200, 222

971113257984, 526729822113536, -1185387156769098, -444423226322880, 305028391047

1560, 304156780634112, -166349637206046, 1079927117208576, -2265102374509140, -9

57842679306240, -5330121520762738, 1526016526220832, 1005434765317584, -24782292

19952640, -862226563598400, -496974423429648, 2393919729693360, -405511548934195

2, 1848385481410944, 1794653372393280, 3719016970449144, 829999889252352, 474119

8635421922, -2904034280063664, 1868559748713000, 2146734178562496, -141307163864

1474, 514968225234624, -1256027191734224, -2817886259865600, -342679570612848, 6

78014201539584, -3031157627077068, 1783468948504320, -1435469944213620, -1476536

265859200, -2979751628238840, -16905402892288, 922102026851124, -544604278740480

, -2747313442193908, 1398970325778816, 1017280238423904, 3735877465876896, 28825

12981450200, 2733643609804800, -2127859436502828, 615432532671600, -442596506902

6536, -628483682444800, 5727465442629702, -1016066618904960, -1942652670145880,

1443076443734016, 3270933095391360, -597835565894448, -1368704107203888, 1924903

992130560, 698326457818910, -1135028972730240, -4331712693612240, -2111401637452

800, -985169554365696, -699348135403008, -1491821174778480, -2474997578876160, -

1799173520665920, -3782019620485440, 9709787093595120, -541834695843840, 3780043

874082112, 67785179563656, -58289341158216, 1467252058138496, -3274269549861750,

 -3013838954784, 2464686886685576, 817229535682560, 3082017633650397, -169682690

046240, 118269009151272, -2580464020796928, 7912852398754982, 4269629283110016,

-2065478171352120, -3667249129586688, -8276770434109008, 840377553364944, -84069

36179151460, 1295492931344640, -1552195881925920, -641384558575488, 136286979447

1992, -1124976089825280, -5387544388705500, 1323881616561072, 3333728731570524,

5434616916158976, -1510999629912864, 2171238655824000, 6817220120892752, 2652778

444947456, 3271729493403504, 1780449350470560, -3982080448605840, -1806121795645

440, -667048975014994, -3515114791751520, 2511663066476928, 4350268758528000, -7

744079729695638, 1588913631793152, -1230393932134640, -4066257732071424, 3790637

408024460, -5554789721607168, 2998000190340120, -3437752020369408, 2524111226919

170, 1037662902581760, 6038220526507512, -8011194672970624, -11145268112466978,

-46709204021856, 1347459830544800, 6337146270766080, 768846891863232, 3596917718

683440, 3285001442747640, 1035109787765760, 5234645524234464, 3089980943425152,

-9409519712300400, -1675135446577152, 6352013923780980, -3819485940068736, 13227

104255899436, 4233606694409088, -6116979573050976, -4418504751225600, 1425699068

010672, -5132574924871680, -4019022069226556, -4994656326572304, -19426551533776

80, -1072192144038400, 2302477758246246, 1079657610071040, -18558651348375264, 5

015858304614400, 2840270079131730, 3201790504598352, 782560406237064, 5742890791

658496, -1507747030957440, -732708148655616, 6511102049004840, 6905645184936960,

 5604717670694010, -193028634807120, -5085164185024588, 3164859781923840, -94951

8954853056, 2337922181380032, -1728332228810040, -1717908208386048, -18258436578

5360, 4948021935321360, -1099342046917296, -2190662393713920, 279110773799022, 1

797458017102272, -13526517861603928, -19071793268183040, -3435312302276400, -208

5312658141440, 272543929352676, 3118693073766912, 18045917610367430, 34005247535

14560, -4137570865480248, -1559486316150784, 11719129719838338, -360940755938112

0, 13305847699403280, -8389992357626880, 3872298735325440, -1761684364651680, -7

961834707368360, -1721713132339200, 4287752862251071, -4106686376126928, 5301002

172584664, 9999543942484736, -7043971523408190, 7329560283142272, -8343868958152

, -1575524037500928, 4211720434286064, -4226303157627600, -3397052374923408, -36

32092712736768, -14982644001390898, 2795502689196672, -5852007487279800, 7623581

505914880, -22256130175851894, 1868019557983488, 5442366024533060, 1217596737100

80, -1300110507573696, 1722619906722432, 10810987426382112, 14414456987320320, -

4590371381461740, -2397149957416368, 3381847082390196, -1302434439348224, 203777

23102676160, -3750993930499200, 8944439700308392, 6207816494668800, -85247669702

73078, -4770330071189760, 6011114426772000, -542984608147968, -28102437759471274

, 696749301810240, -6029568464839416, 2548958848450560, 4222089494633682, -46283

57208458424, 516091995437132, -5873015772013824, -6316091224178400, -78808763692

48608, 5719060466578536, -3879131953213440, 4398936813945856, 2897178951391200,

890229680050932, -10796918367606784, 28650951002224320, 6746142538900800, 812108

5805687400, -5660667323351040, -2714273523041472, 14696835447157200, -6785413121

403840, -4265574004065600, 11022255310711932, -3953397257724096, 72261719397504,

 -7193206955658240, -4814420815766940, 4986257257384704, -84377784300844, 200168

3841040128, -9289515101228586, -1121368936746240, -11009058285505488, 2652116416

266240, -15682896303708816, -7283592789880944, 8467147568239920, 174055333205062

40, -10574317556301408, -12145451761457280, -4863513801242680, 7607170877644800,

 3858395392796112, 17503390720054512, -5656913118074736, 3320521288261632, 46462

98291681650, -2932515937769088, 25655490084944664, -25258393596395520, 357534397

6258530, 1543646168709120, -18092058797213940, 25853256276218496, 12568540465129

824, 14817136407899616, -17832336756146640, 5546608703032320, 38637331724206586,

 -6221465024296704, -25069169854622376, -19106402842920960, -34899727878407658,

-7953826130647488, 5744373208966224, -5122449823088640, 5926337141961600, 219734

616127104, -28512291118415604, -14139215203921920, -845687870789636, 65245907375

88000, 21009050586560232, 9768672367534080, 40033434235820202, -4353133384507536

, 13342408183359360, 10509458025477120, 17299226769945120, -25220111636401920, 4

976295806607216, 3285476054138880, -22616076492128607, -2526517334489568, -10255

287323350908, 340483770574976, 26286732520684860, 1887961499633664, 184953043083

97016, 42098067939840, -18554019842302560, -8695572865852320, -21487665696030828

, -19725972671349504, -9992711669027360, 5277839322284736, 3712431555768600, 686

6990520492032, -8738799717564174, 2164172529988224, -13361513402011320, 12065015

350120320, -8350791783897330, 11612715073634592, 11892426157508232, 655956376492

0320, -13891521966029820, 21282425086618080, -11339234730690048, 906679499283712

0, -319279232436576, 7042709342934720, 23440857788319392, 10394584268931072, -16

023173525318736, 3934767105464832, 3518130472626000, 10881679211246592, -2140041

5987399554, 2613065262540000, 13352424013671120, -30328412970240000]

[[2, -24, 0.5303300858899106], [3, 252, 0.5987336124929452], [5, 4830, 0.6912133

33204735], [11, 534612, 1.0008729094970374], [17, -6905934, 1.1796504994162533],

 [47, 2687348496, 1.7091720053051764], [103, -225755128648, 1.9188140482390927]]