160413

Ruby


Zagier's problems(10)

m = 5, n = 196 のときの計算結果を残しておく。

def f(n)
  return 1 if n <= 1
  (1..n).inject(:*)
end

def u(m, n)
  i = (0..n).inject(0){|s, k| s + f(m * n + k).to_r / (f(m * n - m * k) * f((m + 1) * k + 1))}
  j = i.denominator
  return i.to_i if j == 1
  i
end

(190..200).each{|i| p u(5, i)}

出力結果
(14096922886878560361335788817154728199798428416254259239690126227412717494279242744885446032320006077757571420585333094476837266931412679357223627408073913710663264928402597278964770272621778099851913481302243430074505739640018230274303208289030498671883402934119131578330068074908960836562780185014295700732017844/5705)
(124986505056204611409009555132533466163901413328921120429053202544167993963952358227335809779686077744167837827204340115494700960030565447866788610852515859477162518308957099770037857047229094294860960786374202220231297888906516948302554246934333006937003041061039838632339276873281118924934440104513661400483694808/1147)
(27703965204891092190861560363395388039986658437725672403625445307986004215898496854365631405464858723553599285347435970756381438683067982505528628800625855347124348948156445916183446125518040658290849780845314245998005214264657113911688456155021139123548857215083084441870785851658653894097623547062189581304563768132/5765)
(1228148091233815050052958419210126612524844315005562951783519373592877523355790891395756160550693916540899220202237801122612252644741224118051699654375408165019080020843100055139056546466390902769491377489958122318550666429362907760791683976744901000112970574765397994203294380055335653394969971676475314692575636516479/5795)
(54445193056153811412056862596370167876223873867082957601678911375022525503157093226353797320000307546461886455892884704571408839383761689187751753837083021804115367916578903837399928129550224273788323022818256933895270089423673457676885432114780252007469796438030281822903616814653961900360170376807076804109794516249413/5825)
(2413616947402411616676741451517267513528638802008443525312198879373539437285699820243147938542481339297423963802129347301094141782063408826084261649771864134686014205444511491443301872478590618968853867861245510959115122584957961052517803983012396993907227859643363768391025750830073959352608287561703251067127503684705229/5855)
18181540761540712867149487897053111408341145153247021925742297078013378875918440895079511284629222125948897573866467613788397445534812254595922136833868755211587284257245954094057967606737047389560734422923576883191584992444803051004333753987804237556089739248567296230384275913141809451316512176176377053312906111739000
(4743358565932963848095491665985122216650207829948956324154504038186089524765496114993915291969756350781023801019272206655357377755866659520010906834744364101420092177060177807314020268372200539994761360652698124737441342421527858536862221380708754404696710971290364464165568190402315407515194885444978408637736425531193116524/5915)
(42055689326087115762783402823125733699826295093344726808790573319464007329213101151235893565233316965204569718371598113197970569412479945738957000670791549732822027499394896866257172662185531178825309459166679779923773314959820184133293532950727048075258277284458499371731392536785755174046448700085118409447913062966683916919/1189)
(9321881216165652840273527917467630283950547748371610933737455857725642086131629624711182672950147237203302293120429774901540426519471243432702750308438234254139399789921208901254874816723162779833911372678185485081404386352193067105351817762135529682005488932497168730911200884508027219942175578078870222502568940488317547673957/5975)
(413249530804383173485275814472324437620635643432036662479799372337212013144150979354876931502961993685566541310681669294760175877220158686677328394523077300538238183248829631436826842120471524918631177941266315869588930609048606478273028565442188335846010022617201120108944030409468018235323923914203958076324506866736290311266247/6005)