[過去ログ]
臨床統計もおもしろいですよ、その2 (1002レス)
上
下
前
次
1-
新
このスレッドは過去ログ倉庫に格納されています。
次スレ検索
歴削→次スレ
栞削→次スレ
過去ログメニュー
902
: 2020/02/01(土)22:05
ID:hIisy8jC(2/2)
AA×
[240|
320
|
480
|
600
|
100%
|
JPG
|
べ
|
レス栞
|
レス消
]
902: [sage] 2020/02/01(土) 22:05:31 ID:hIisy8jC H(n) = Σ[k=1,2,...,n] 1/k とする。H(n)を既約分数で表したときの分子の整数をf(n)と表す。 f <- function(n){ # Σ 1/kを既約分数表示する if(n==1){ cat(n, ':' ,1,'\n') invisible(1) }else{ GCD <- function(a,b){ # ユークリッドの互除法 r = a%%b # a=bq+r ⇒ a%%b=b%%rで最大公約数表示 while(r!=0){a = b ; b = r ; r = a%%b} b } nn = 1:n # nn : 1 2 3 ... n nume=list() # 分子の容器となるlist length(nume)=n # を設定 Nume=0 library(gmp) for(i in nn) Nume = Nume + prod.bigz(nn[-i]) # nnからi番目の要素を除いて乗算して総和を分子に Deno=factorialZ(n) # 分母 gcd = GCD(Nume,Deno) # 約分するため最大公約数を計算 NUME=Nume/gcd DENO=Deno/gcd ratio=as.numeric(Nume/Deno) RE=list(NUME,DENO,ratio) cat(n, ':' ,as.character(NUME),'/',as.character(DENO),'\n') invisible(RE) }} http://egg.5ch.net/test/read.cgi/hosp/1540905566/902
とするを既約分数で表したときの分子の整数をと表す を既約分数表示する ユークリッドの互除法 で最大公約数表示 分子の容器となる を設定 から番目の要素を除いて乗算して総和を分子に 分母 約分するため最大公約数を計算
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 100 レスあります
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
AAサムネイル
ぬこの手
ぬこTOP
0.027s