- 688 名前：卵の名無しさん mailto:sage [2018/08/24(金) 07:53:38.25 ID:y2eRJo9I.net]
- >>649
結局はプログラムで出すのだけれど
もっとエレガントな解法が数学板にC言語で投稿されたのでこれをRに移植してみた。
a=355
b=113
f <- function(k){
dk=abs(floor(k*pi)/k-pi)
dk1=abs(floor(k*pi+1)/k-pi)
min(dk,dk1)
}
d=abs(a/b-pi)
k=1
while(f(k)>=d){
if(f(k)<d) break
k=k+1
}
k
> k
[1] 16604
瞬時、答が返ってきた。
数学板は切れ者が多いなぁ。
このスレのド底辺シリツ医大卒は
