2016年1月12日火曜日

160112

Ruby


単位分数の和が1 / 2(3)

p + q > n という条件をはずすとどうなるか調べてみた。
すなわち、
n を1 より大きい自然数とし、
互いに素となる正の整数p, q (ただし、0 < p < q ≦ n)で、
Σ(p, q) 1 / (pq)
を求めてみた。

def f(n)
  ary = []
  s = 0
  (2..n).each{|i|
    (1..i - 1).each{|j|
      s += 1r / (i * j) if i.gcd(j) == 1
    }
    ary << s
  }
  ary
end

p f(50)

出力結果
[(1/2), (1/1), (4/3), (7/4), (39/20), (23/10), (527/210), (463/168), (1469/504), (334/105), (3812/1155), (2803/792), (5723/1560), (229373/60060), (710887/180180), (66375/16016), (523689/123760), (2254127/510510), (131201039/29099070), (98005643/21162960), (1587557/335920), (15799073/3233230), (368539631/74364290), (232942821/45762640), (968698307/187210800), (70723666817/13385572200), (215198406259/40156716600), (534876161/97349616), (172178330879/31054527504), (943823276783/166363540200), (69160421412893/12033629407800), (16525794496429/2831442213600), (10032379128731/1698865328160), (1604617834569/267413986840), (622305953361/102851533400), (6691804682987/1085740999200), (583649404869259/93735639597600), (112299040878679/17809771523544), (339761515094597/53429314570632), (11511237187650241/1780977152354400), (1488661859008151/228903019581600), (16439677922891971/2489320337949900), (36411729017455/5468238801116), (4078482098102365/606736757410784), (2775287404425615383/409547311252279200), (5393065707566515283/784965679900201800), (1219770607484002541/176523382561289400), (619261627853593128913/88544128692742763040), (1454173851705467756797/206602966949733113760)]

0 件のコメント:

コメントを投稿

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