2016年1月11日月曜日

160111

Ruby


逆数をとり、単位分数を加えていく数列

以下のような操作を繰り返してみた。
0 + 1 / 2 = 1 / 2
2 / 1 + 1 / 3 = 7 / 3
3 / 7 + 1 / 4 = 19 / 28
28 / 19 + 1 / 5 = 159 / 95

def f(k, n)
  s = 0
  ary = []
  (k..k + n - 1).each{|i|
    s += 1r / i
    ary << s
    s = 1r / s
  }
  ary
end

p f(1, 50)
p f(2, 50)
p f(3, 50)
p f(10, 50)
p f(100, 50)

出力結果
[(1/1), (3/2), (1/1), (5/4), (1/1), (7/6), (1/1), (9/8), (1/1), (11/10), (1/1), (13/12), (1/1), (15/14), (1/1), (17/16), (1/1), (19/18), (1/1), (21/20), (1/1), (23/22), (1/1), (25/24), (1/1), (27/26), (1/1), (29/28), (1/1), (31/30), (1/1), (33/32), (1/1), (35/34), (1/1), (37/36), (1/1), (39/38), (1/1), (41/40), (1/1), (43/42), (1/1), (45/44), (1/1), (47/46), (1/1), (49/48), (1/1), (51/50)]
[(1/2), (7/3), (19/28), (159/95), (81/106), (823/567), (5359/6584), (64615/48231), (21877/25846), (306183/240647), (354883/408244), (5662055/4613479), (1433689/1617730), (25699639/21505335), (369784999/411194224), (7360086807/6286344983), (1487830821/1635574846), (32563752895/28268785599), (23917578595/26051002316), (570988627231/502269150495), (96040578001/103816114042), (2483811200967/2208933294023), (6166467806391/6623496535912), (171753881204191/154161695159775), (24733479025789/26423674031414), (738172677873967/667803933696303), (396666996354499/421812958785124), (12629242801123095/11503342894280471), (1589909020575721/1683899040149746), (53790779265217847/49287179637847351), (1630980527676333079/1721304936486971104), (58434043431746379511/53822357413318991607), (6534235970534920741/6874593344911338766), (247145003042431777551/228698258968722225935), (104694843529832492731/109842223574414123356), (4168857115783155056903/3873709210603802231047), (419306944927223379049/438827064819279479674), (17533562472879123086335/16352970852161711782911), (26866095862373903776111/28053699956606596938136), (1177067794083244378239687/1101509930357330054820551), (107574795621521783856469/112101694674594702689494), (4927947666629093999504711/4625716211725436705828167), (1722805462665688504594579/1791980969683306908910804), (82361949098414499405580759/77526245819955982706756055), (6897106345588638381694441/7161908617253434730920066), (343506811356500070734937543/324163998242666003939638727), (1767042080778274251093066271/1832036327234667043919666896), (91536822115276959403156744175/86585061958135438303560247279), (7073263872035278199329870573/7322945769222156752252539534), (380543498102365272564209386807/360736457473799188165823399223)]
[(1/3), (13/4), (33/65), (47/22), (201/329), (2833/1608), (17305/25497), (10891/6922), (87033/119801), (169405/116044), (1677977/2202265), (663463/479422), (7854793/9951945), (167085913/125676688), (2303589609/2840460521), (659652827/511908802), (10385920065/12533403713), (10442159773/8308736052), (184925616865/219285355233), (41398375471/33622839430), (814723682361/952162635833), (2629625215817/2172596486296), (56944537373217/65740630395425), (10450893063043/8760698057418), (246989740613329/282174112702161), (166282957066813/141136994636188), (4259255801516265/4822205754937577), (661890793109527/567900773535502), (18266814772710089/20518614586395337), (674862481537360873/584538072726722848), (19964618881519214857/22270461890732908809), (2689136066319855067/2348778691943437042), (84896390284340151537/94119762321194927345), (42879109059843920197/37731729015262289572), (1438953082624548634361/1586527035214225047289), (170988865431482272663/151468745539426172038), (6078269941469102982145/6668565751827808633857), (10912836000583257933457/9725231906350564771432), (409647344160956413562169/447426276023913575271737), (43540931830306862981803/39014032777233944148778), (1721144341251366461379257/1872260068703195108217529), (695046176563569844817773/625870669545951440501548), (28859226306131384667387433/31277077945360643016799785), (2774300211328395015954967/2509497939663598666729342), (120720703375517532352234041/130392109932434565749883449), (708837997792486298705182177/643843751336093505878581552), (32257181813261068086755678225/34733061891831828636553926673), (2830256442247763999863123219/2580574545060885446940454258), (134439558240352921793826290377/144343078554635963993019284169), (45208755284505461830951651285/41366017920108591321177320116)]
[(1/10), (111/11), (27/148), (1951/351), (6865/27314), (16663/4119), (82567/266608), (4614903/1403639), (3320045/9229806), (178686359/63080855), (1440303459/3573727180), (173443479/68585879), (13903577/31535178), (739212671/319782271), (8413987175/17741104104), (18077663591/8413987175), (1401427989/2781179014), (8499251263/4204283967), (126219202339/237979035364), (7027611227895/3660356867831), (519281410057/937014830386), (29566741152023/16097723711767), (544693899928567/946135716864736), (31767172556464855/17974898697642711), (2224614976734661/3737314418407630), (2714910604510239/1589010697667615), (739744391611659/1206626935337884), (45384940999113367/27370542489631383), (3006829793919961/4777362210432986), (189323956000806415/117266361962878479), (195199137380637823/302918329601290264), (12614850651033538647/8003164632606150743), (38749751691165763317/58869303038156513686), (2570129782331895851815/1666239322720127822631), (627145950264607603699/934592648120689400660), (42683815115695630633399/28221567761907342166455), (2534737111840138696201/3711636097017011359426), (176981633671639672589223/119132644256486518721447), (655038728664776952357631/943902046248744920475856), (957270183568434246034175/655038728664776952357631), (1348368264672291274556629/1914540367136868492068350), (98989926988652584370042479/68766781498286855002388079), (21744748904731177778072323/30458439073431564421551532), (1636042019796604092120303519/1152471691950752422237833119), (788512510927620183863744345/1090694679864402728080202346), (2431068796138790809130994935/1734727524040764404500237559), (2032159390661665254309067311/2778364338444332353292565640), (17822102964665401043553923199/12870342807523879943957426303), (908801409989358427815796253/1229110549287269037486477462), (73426323817938231639517966511/53619283189372147241131978927)]
[(1/100), (10101/101), (2267/114478), (11793501/233501), (36077605/1226524104), (5152844341/151525941), (21214594087/546201500146), (58464775109709/2269961567309), (33735624931009/701577301316508), (76505661468430381/3677183117479981), (43726891308293481/765056614684303810), (9440556791251779599/539298326135619599), (1425346312621248463/21578415522861210512), (2439786300395938036319/161064133326201076319), (20801097499582860736685/278135638245136936140366), (1280255979827613220675151/95685048498081159388751), (12379721605605027709770267/148509693660003133598317516), (1270217976464750650793551/105809586372692544527951), (13755749168442470905091769/149885721222840576793639018), (61765247663274986537536799/5664132010535135078567199), (741461088927491195965600679/7411829719592998384504415880), (81597532469061845065545447469/8156071978202403155621607469), (1076638313809755030051381558687/9954898961225545097996544591218), (136169912282727978009291818475389/14714056955399985410702214635389), (1960712974752326168936366433263625/16885069123058269273152185490948236), (3380150965371257576500735284482869/392142594950465233787273286652725), (1077349345492446470075452436790331/8691816768097519482430462160098806), (1104938078893877420738744146769338693/136823366877540701699582459472372037), (18618329039219087238285298959232959429/141432074098416309854559250786475352704), (2029261765303880339830714294490505939805/266862716228806917082089285082339085149), (1468856595001951182420092854207783480367/10552161179580177767119714331350630886986), (1383801971120005238675102670261140429675533/192420213945255604897032163901219635928077), (221349340594163182521349985993571341918857/1509602150312732987645566549375789559646036), (556782369341240029194963160811561171121445/81549757061007488297339468523947336496421), (11484449815516243461038451943020504261641859/74608837491726163912125063548749196930273630), (124489845817266029278986691741039038146278789/19140749692527072435064086571700840436069765), (2727631804000947880447702465492353337451766829/16930619031148179981942190076781309187893915304), (2322222439071301605406527742984531712078918163477/373685557148129859621335237772452407230892055573), (11319224811071880326223648509889301388351611391/67310795335400046533522543274913962668954149666), (9367519776431678348485857163722930112372978414965/1573372248738991365345087142874612892980873983349), (187460926203992219997386169115239783779343131497/1070573688735048954112669390139763441414054675996), (151138351037845894749883770178821885023161052446933/26431990594762903019631449845248809512887381541077), (3904481015494178123537549648204152835853169231279867/21461645847374117054483495365392707673288869447464486), (18182957616508833827897499330765451775835275149223085/3303791628495073796839465086941975476491143195698349), (54880994679977717841420052427823324487839988369976149/290927321864141341246359989292247228413364402387569360), (50225139910797232102929391795248123001695217995918589/9462240462065123765762078004797124911696549718961405), (1431712247372305301904192780495628360109391476964283719/7332870426976395887027691202106225958247501827404113994), (119929296112544722299441644387790093802499128900596560093/23384633373747653264435148748095263215120060790416634077), (2615672049252883422524356215563103761607208273106105481/12965329309464294302642339933815145275945851773037465956), (1934449739159432734516233006354019749877539122455688532925/389735135338679629956129076118902460479474032692809716669)]

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。