高校数学の質問スレ（医者・東大卒専用） Part438

>>411
生まれてから死ぬまで裕福な私立医大医者が羨ましくて仕方のない様子

　 　 　 　 　 　 定 　 公 　 こ 　 o 　 無 　 な
　 　 　 　 　 理 　 式 　 れ 　 n 　 料 　 い
　 　 　 　 は 　 は 　 も 　 　 　 で 　 ね
　 　 　 ツ 　 プ 　 ツ 　 R 　 使 　 。
　 　 ー 　 ロ 　 ー 　 　 　 え 　
　 ル 　 グ 　 ル 　 W 　 る 　
、 　 ラ 　 ） 　 o 　 時 　
　 ム 　 P 　 l 　 代 　
（ 　 y 　 f 　 。 　
　 t 　 r 　 使 　
h 　 a 　 わ 　
　 m 　 な 　
も 　 い 　
　 手 　
は 　
　

"
注射液8本のうち2本が致死的な毒薬である場合
ネズミ1: 注射液1, 2, 3, 4
ネズミ2: 注射液1, 2, 5, 6
ネズミ3: 注射液1, 3, 5, 7
に記載された投与法で実験したら全部のネズミが死亡してしまった。
可能な２つの毒薬の組み合わせをすべて列挙せよ。
"
rm(list=ls())

mat=matrix(
c(1, 2, 3, 4,
1, 2, 5, 6,
1, 3, 5, 7),nrow=3,by=TRUE)
mat
cm=t(combn(8,2)) ; cm

lost_rats=\(y) (apply(mat,1, \(x) y[1] %in% x) | apply(mat,1, \(x) y[2] %in% x)) |> which() |> sum()
lost=apply(cm,1,lost_rats)
re=cbind(cm,lost)
tmp1=re[re[,3]==6,1:2] ; tmp1
tmp2=cbind(1:7,8)
rbind(tmp1,tmp2) |> unique()

"
注射液8本のうち2本が致死的な毒薬である場合
ネズミ1: 注射液1, 2, 3, 4
ネズミ2: 注射液1, 2, 5, 6
ネズミ3: 注射液1, 3, 5, 7
に記載された投与法で実験したら全部のネズミが死亡してしまった。
可能な２つの毒薬の組み合わせをすべて列挙せよ。

この投与法で２個の毒薬が特定できる場合をすべて求めて
そのときの毒薬の組み合わせと合わせて列挙せよ。
"
rm(list=ls())

mat=matrix(
c(1, 2, 3, 4,
1, 2, 5, 6,
1, 3, 5, 7),nrow=3,by=TRUE)
cm=t(combn(8,2)) ; cm

lost=\(y){
tmp1=(apply(mat,1, \(x) y[1] %in% x) | apply(mat,1, \(x) y[2] %in% x))
tmp2=as.integer(tmp1)
tmp2 %*% c(1,2,4)|> as.vector() -> tmp3
tmp3
}

i=apply(cm,1,lost) ; i
cbind(cm,i)

lapply(1:7,\(x) cm[i==x,])

>>413
>>413
朝からレス乞食ご苦労様
今日も高校生に相手にされてなくて哀れだね

縦横変換して喜んでるとか小学生かよwww

"
[1年もので500万円お預け入れの場合] 500万円×金利年0.80%（税引前）×1年＝40,000円
40,000円×15.315%（国税）＝6,126円　40,000円×5%（地方税）＝2,000円　40,000円-6,126円-2,000円＝31,874円

[5年もので500万円お預け入れの場合] 500万円×金利年1.20%（税引前）×5年＝300,000円
300,000円×15.315%（国税）＝45,945円　300,000円×5%（地方税）＝15,000円　300,000円-45,945円-15,000円＝239,055円
"
rm(list=ls())
par(bty='l')

(1-20.315/100)
f=\(x,y=5e6) y * ( (1+6*(x/12)/100)^(5*12/6) - 1) * (1-20.315/100)
f(1.2)
a=239055
curve(f(x),0,1.5)
abline(h=a,lty=3)

uniroot(\(x) f(x)-a,c(0,1.5))

>>417
大学でていたら斜め変換プログラムくらい書けるだろう。
ただし、Fランシリツは除く

"
[1年もので500万円お預け入れの場合] 500万円×金利年0.80%（税引前）×1年＝40,000円
40,000円×15.315%（国税）＝6,126円　40,000円×5%（地方税）＝2,000円　40,000円-6,126円-2,000円＝31,874円

[5年もので500万円お預け入れの場合] 500万円×金利年1.20%（税引前）×5年＝300,000円
300,000円×15.315%（国税）＝45,945円　300,000円×5%（地方税）＝15,000円　300,000円-45,945円-15,000円＝239,055円
"
rm(list=ls())
par(bty='l')

(1-20.315/100)
f=\(x,y=5e6) y * ( (1+6*(x/12)/100)^(5*12/6) - 1) * (1-20.315/100)
f(1.2)
a=239055
curve(f(x),0,1.5)
abline(h=a,lty=3)

uniroot(\(x) f(x)-a,c(0,1.5))

ネズミと試薬の組み合わせ行列から単射になるかを判定する関数はできたが
総当たりするには数が多すぎる。

"
行：ネズミ　列：混合試薬　から　単射　か否かを判定する。
試薬数が不揃いなときはNAを入れる
[,1] [,2] [,3] [,4]
[1,] 1 2 3 8
[2,] 2 4 5 6
[3,] 1 2 4 5
[4,] 1 3 4 6
[5,] 4 6 7 NA
[6,] 1 2 5 NA
"

rm(list=ls())

rats <- function(num, N=2, digit = 6){
r=num%%N
q=num%/%N
while(q > 0 | digit > 1){
r=append(q%%N,r)
q=q%/%N
digit=digit-1
}
which(r==1)
}

cm=t(combn(8,2))
calc82=\(mat,verbose=FALSE){
lost=\(y){
m=nrow(mat)
tmp1=(apply(mat,1, \(x) y[1] %in% x) | apply(mat,1, \(x) y[2] %in% x))
tmp2=as.integer(tmp1)
tmp2 %*% c(2^(0:(m-1)))|> as.vector() -> tmp3
tmp3
}
i=apply(cm,1,lost)
re=sapply(1:max(i),\(x) cm[i==x,])
if(verbose) print(re)
all(unlist(lapply(re,length))<=2)
}


mat=t(replicate(6,sort(sample(8,4)))) ; mat
calc82(mat,TRUE)

id:qk9aTmTZ

正規の高校数学スレじゃもう高校生にバカにすらしてもらえなくて哀れだね
そうやってリアルでも誰にも相手にしてもらえず見捨てられた老害なんだね

[,1] [,2] [,3]
[1,] 1 2 3
[2,] 1 4 5
[3,] 2 4 6
[4,] 3 5 6
[5,] 1 6 7
[6,] 2 5 8

> mat
[,1] [,2] [,3] [,4]
[1,] 6 1 8 NA
[2,] NA 5 8 3
[3,] 5 2 NA 1
[4,] NA 6 7 5
[5,] 3 NA 4 6
[6,] 3 1 7 NA

尿瓶ジジイダンマリきたw

答を出すのに小中学生レベル以上の基礎知識が必要かを検討。
「
注射液が８本あり、そのうち２本が致死的な毒薬であることがわかっている。
注射液を混合して６匹のネズミに注射して８本のうちどの２本が毒薬かを特定したい。
複数の注射液を混合して注射をしても毒性は変わらないものとする。
各々のネズミに注射できるのは１回だけである。
毒薬を特定する混合法を例示せよ。
」
という問題に下記の例を使う。

ネズミ６匹への混合注射液の組み合わせは下記であったとする。

[1,] 1 2 7
[2,] 1 3 6
[3,] 1 4 5
[4,] 2 3 5
[5,] 2 4 6
[6,] 3 4 7

(1)　ネズミ1,2,3のみが死亡したときの毒薬の組み合わせを求めよ。
(2)　ネズミ4,5,6のみが死亡したときの毒薬の組み合わせを求めよ。

>>425
今度はこっちを荒らすのかよ。死ねクソ野郎

東大卒なら算出できるんじゃないの?
Wolfram,Pyth9n,C,Rのどれかは理工系卒なら使えて当然。
Fランは違うのか?

>>427
東大卒なのにスレタイも読めないのかよ？
ほんとは小学校すら通ってないんだろ？

(*
８本の注射液があり、赤,青,黄,緑,白,黒,紫,茶のラベルが貼られている。
いずれも十分な量があり、そのうち１本が致死的な毒薬であることがわかっている。
３匹のネズミに注射液を混合して注射してどれが毒薬かを特定したい。
複数の注射液を混合しても毒性は変わらないとする。
各々のネズミに注射できるのは１回だけである。

赤 青 黄 黒 の混合液
赤 白 黒 茶 の混合液
青 緑 白 黒 の混合液

を３匹に１つずつ注射すれば特定できる。
ネズミを区別しないとして、このような組み合わせは何通りあるか。
*)


s4=Subsets[Range[8],{4}];
grp=Subsets[s4,{3}];
Short[grp]
rat=Subsets[Range[3],3]

f[x_]:=(
z1=Length@Complement[Range[8],Flatten@x];
z28=Table[
dead = Intersection@@x[[rat[[i]]]];
alive= Flatten@x[[Complement[Range[3],rat[[i]]]]];
Length@Complement[dead,alive],{i,2,8}];
AllTrue[Flatten@{z1,z28},#==1&]
)
ans=Select[grp,f];
Length[ans]
Short[ans]

RをWolframに移植してカウント、数値は合致した。
理工系卒ならWolframくらい使えて普通じゃないの？
Fランシリツは違うのか？それとも私文卒かよ？

>>428
俺は医科歯科卒だよ。二期校時代に理Iは蹴って医科歯科入学。

>>431
妄想www
もしくは大学の恥さらしwww

いずれにせよ現在は痴呆か無能だろ

>>431
じゃあ何でスレタイ読めないんだよアホが

>>431
今まで散々アホ晒してきたのに今更誰が信じるんだよ？
そんなことも分からないくらい救いようのないアホだから一生クソみたいな書き込みしてんだろうなw

医者板の自分が立てたクソスレは散々論破されて全く関係ないことでしか発狂できなくなった模様、実に哀れ

https://itest.5ch.net/egg/test/read.cgi/hosp/1733201932

まあここでも同じかww

嫌われる奴って何処でも嫌われるんだな

医科歯科卒ってスレタイも読めないんだな
レベル低ｗ

理工系卒ならWolframのコードくらい読めるんじゃないのか？
Fランは例外なのか？

>>438
読めるよ。で？
クズ問題とクズコードに価値が無いのは変わらない。

>>438
これは酷い。
君は、何が悪くて総スカンされているのか全く理解をしていない。
もう、社会性が完全に欠落した悪意を持った生き物と言える。

>>438
社会はもちろん5chでも相手にされないレス乞食のゴミって生きてて何が楽しいの？

㋖印は自分が㋖印で有ることを理解しない

>>431
結局スレタイすら読めないチンパンがキーキー発狂してるだけかよ

お願いします。


１辺が１の正方形の各頂点を中心とする半径rの円を考え、各円の周上にそれぞれ１点ずつP,Q,R,Sをとるとき、
四角形PQRSの面積の取りうる値の範囲を求めよ。ただし、0<r<1/2とする。

>>439
では毒薬が２本のときのWolframのコードを書いてみ！

>>445
スレタイ読めるならスレタイに則った書き込みしてごらんよ

>>445
Wolframを使えないのが確定。
やはり、Fラン卒なのか？

>>445
クレクレ乞食www
しかも個人の識別も出来ないのかよwww
クズ中のクズwww

>>445
煽って乞食行為とかテンプレいらないぞ
少しはオリジナリティのある返しをしてくれないとwww

その論理なら、クズコード書き込んでるやつはwolfram使えないの確定だぞwww
自白かよwww
さすが中卒www

>>447
なに自分で自演してんの？馬鹿なの？
それとも自分がFラン卒だって自白したかっただけ？

>>444
>>444
Wolframで思い通りに作図するスキルはないので
R言語で最小値と最大値を求めてそのときの図を作図。
r=0.3のとき
> c(min,max)
[1] 0.3314719 2.0285281
最小のときの図
https://i.imgur.com/U9IMSgw.png
最大のときの図
https://i.imgur.com/olgz6rn.png

最小最大のときの形状が判明したので
あとは一般解を求めればいい。

>>452
おまけ、R言語でのクズコード
理工系卒ならWolframで書けるはず。
但し、Fラン卒は除く。

r=0.3
source("toolmini.R")
p=0.5+0.5i
q=-0.5+0.5i
r0=-0.5-0.5i
s=0.5-0.5i
fP=\(t) r*exp(-1i*t)+p
fQ=\(t) r*exp(-1i*t)+q
fR=\(t) r*exp(-1i*t)+r0
fS=\(t) r*exp(-1i*t)+s

f=\(tt,verbose=FALSE,...){
P=fP(tt[1])
Q=fQ(tt[2])
R=fR(tt[3])
S=fS(tt[4])
if(verbose){
Plot(-1,1,zero=FALSE)
Polygon(p,q,r0,s,Col=8)
Cir(p,r,col=8)
Cir(q,r,col=8)
Cir(r0,r,col=8)
Cir(s,r,col=8)
Polygon(P,Q,R,S,...)
pta(P);pta(Q);pta(R);pta(S)
}
ABC2S(P,Q,R)+ABC2S(P,S,R)
}

opt=optim(runif(4,-pi,pi),f)
opt=optim(opt$par,f)
tt=opt$par
min=f(tt,TRUE,Col=2)

opt=optim(runif(4,-pi,pi),f,control=list(fnscale=-1))
opt=optim(opt$par,f,control=list(fnscale=-1))
tt=opt$par
max=f(tt,TRUE,Col=4)
c(min,max)

>>453
算出のために作った関数の動作確認の図
https://i.imgur.com/0mvHN8w.png

毒薬が２本のときのWolframのコード 　まだぁ？

>>452
rが1/2に近いとき、明らかに長方形の時最小なるような…

別スレを潰したのに飽き足らずこっちのスレも潰す気かよ

>>452
横軸にｒ、縦軸にPQRSの面積の最大値（青）最小値（赤）をグラフ化。
https://i.imgur.com/iio0LJI.png
Rのコードが読めれば、単純作業。（但し、Fラン卒は除く）。

なんか発狂して書き込んでるみたいwww
クズ問題読んで貰えると思ってるのウケるwww
自意識過剰すぎだろwwww

やはり、RのコードもWolframのコードもないクズ書き込みのみ。
Fラン確定！

Clear[{R,f,a,b,g}]

(* 球の半径 *)
R = 20; (* 半径を任意の値に設定 *)

f[a_,b_]:=(
{lat1,lon1}={a,0} Degree;
{lat2,lon2}={b,90} Degree;
cosc=Sin[lat1]Sin[lat2]+Cos[lat1]Cos[lat2]Cos[lon1-lon2];
R ArcCos[cosc]
)
f[a,b]

(* 点の座標 *)
latLonToCartesian[R_, φ_, λ_] := R * {
Cos[φ Degree] * Cos[λ Degree],
Cos[φ Degree] * Sin[λ Degree],
Sin[φ Degree]
};

g[a_,b_]:=(
(* 点 pC, pA, pB を設定 *)
pC = latLonToCartesian[R, 90, 0] ;(* 点 P (90°, 経度0°) *)
pA = latLonToCartesian[R, a, 0] ; (* 点 A (緯度a°, 経度0°) *)
pB = latLonToCartesian[R, b , 90] ; (* 点 B (緯度0°, 経度90°) *)

(* 大円の法線ベクトル *)
n = Cross[pA, pB];

(* 点 P から大円への最短距離（球面上） *)
distance = ArcSin[Abs[Dot[n, pC]] / Norm[n]] * R
)
g[a,b]

>>460
クズはいろんなスレに関係ないことしか書き込まないから総スカンにされてバカにされてるアンタのことだよ

>>460
やはりスレタイ読めない
保育園卒確定！

>444に答えてあげればいいのに。
それもできないのでFラン確定！

fd=\(lat1,lon1,lat2,lon2,R){
R*acos(
sin(lat1)*sin(lat2)+cos(lat1)*cos(lat2)*cos(lon1-lon2)
)
}
fd(pi/5,0,pi/3,pi/2,1)

>>463
死ね

>>463
スレタイに則った書き込みすればいいのに
それもできないのでホイ卒確定

spherical_triangle_area <- function(R, a, b, c) {
# 直交座標への変換
to_cartesian <- function(lat, lon) {
c(cos(lat) * cos(lon), cos(lat) * sin(lon), sin(lat))
}

p <- to_cartesian(a[1], a[2])
q <- to_cartesian(b[1], b[2])
r <- to_cartesian(c[1], c[2])

# 辺のベクトルを計算
u_ab <- q - p
u_bc <- r - q
u_ca <- p - r

# 内角を計算
angle_between <- function(u, v) {
acos(sum(u * v) / (sqrt(sum(u^2)) * sqrt(sum(v^2))))
}

alpha <- angle_between(u_ab, -u_ca)
beta <- angle_between(u_bc, -u_ab)
gamma <- angle_between(u_ca, -u_bc)

# 超過角を計算
E <- alpha + beta + gamma - pi

# 球面三角形の面積を計算
R^2 * E
}

https://i.imgur.com/Q34Vduy.jpeg
>>466

類は何とやら
キチガイはキチガイをよぶ

キチガイが複数人いるより、キチガイの自演の方が自然

Wolfram Language 14.0.0 Engine for Microsoft Windows (64-bit)
Copyright 1988-2023 Wolfram Research, Inc.

In[1]:= 4 * Integrate[Sqrt[r^2-x^2]-r/2,{x,r/2,r*Sqrt[3]/2}]

2
r (-3 (-1 + Sqrt[3]) r + Pi Sqrt[r ])
Out[1]= -------------------------------------
3

In[2]:= r=6

Out[2]= 6

In[3]:= f[p_,o_] := Total[(p-o)^2] <= r^2

In[4]:= f4[p_] :=(
AllTrue[{f[p,{0,0}],f[p,{r,0}],f[p,{r,r}],f[p,{0,r}]},#==True&]
)

In[5]:= ratio=f4 /@ Tuples[Range[0,r,0.001],2] // Boole // Mean;

In[6]:= r^2 * ratio // N

Out[6]= 11.3415

In[7]:=
11.3415

>>470

137：卵の名無しさん：[sage]：2024/12/27(金) 16:30:06.83 ID:UPA/6UxQ
レス乞食じゃなきゃ然るべきスレでやるもんね

横綱、大関、関脇で巴戦を行う。
２力士が対戦して，勝った方が残りの 1 力士と対戦をする。
これを繰り返して，２連勝した力士が優勝をする。

勝利確率は通算成績に従うとする（巴戦での成績も組み入れて計算する）
巴戦前の通算成績は
横綱 対 大関は　横綱の１０勝９負
大関 対 関脇は　大関の　８勝７負
横綱 対 関脇は　横綱の　６勝５負
とする。

大関　対　関脇で巴戦を始めるときにそれぞれの力士の優勝する確率を算出しなさい。

sim=\(){
h12=c(10,9) # bout3
h23=c(8,7) # bout1
h13=c(6,5) # bout2

bout=1
last_winner=0
flg= FALSE
count=0
while(!flg){
count=count+1
if(bout==1){
i=2-rbinom(1,1,h23[1]/sum(h23))
h23[i]=h23[i]+1
winner=c(2,3)[i]
bout=c(2,3)[-i]
}
if(bout==2){
i=2-rbinom(1,1,h13[1]/sum(h13))
h13[i]=h13[i]+1
winner=c(1,3)[i]
bout=c(1,3)[-i]
}
if(bout==3){
i=2-rbinom(1,1,h12[1]/sum(h12))
h12[i]=h12[i]+1
winner=c(1,2)[i]
bout=c(1,2)[-i]
}
flg <- winner==last_winner
last_winner=winner
}
c(winner=winner,count=count)
}

t(replicate(100,sim()))

sim=function(){
h12=c(10,9) # bout3　横綱1vs大関2
h23=c(8,7) # bout1　大関2vs関脇3
h13=c(6,5) # bout2　横綱1vs関脇3
bout=1 # 初戦の試合番号
last_winner=0　# 前勝者
flg= FALSE # 連勝フラッグ
count=0 # 試合数
while(!flg){
count=count+1
if(bout==1){
i=2-rbinom(1,1,h23[1]/sum(h23)) # 勝者のindex
h23[i]=h23[i]+1 # 通算成績更新
winner=c(2,3)[i] # 勝者番号
bout=c(2,3)[-i] # 敗者番号=次の試合番号
}
if(bout==2){
i=2-rbinom(1,1,h13[1]/sum(h13))
h13[i]=h13[i]+1
winner=c(1,3)[i]
bout=c(1,3)[-i]
}
if(bout==3){
i=2-rbinom(1,1,h12[1]/sum(h12))
h12[i]=h12[i]+1
winner=c(1,2)[i]
bout=c(1,2)[-i]
}
flg <- winner==last_winner
last_winner=winner
}
c(winner=winner,count=count)
}

t(replicate(100,sim()))

"
横綱、大関、関脇で巴戦を行う。
２力士が対戦して，勝った方が残りの 1 力士と対戦をする。
これを繰り返して，２連勝した力士が優勝をする。

勝利確率は通算成績に従うとする（巴戦での成績も組み入れて計算する）
巴戦前の通算成績は
横綱 対 大関は　横綱の１０勝９負
大関 対 関脇は　大関の　８勝７負
横綱 対 関脇は　横綱の　６勝５負
とする。

(1) 大関　対　関脇で巴戦を始めるときにそれぞれの力士の優勝する確率を算出しなさい。
(2) 巴戦の試合数の期待値を求めなさい。

尿瓶ジジイのレス乞食はいつになったら願いが届くのかな？w

>>476
あなたが返信してる時点で目的を達してる
完全スルーができないならスレを見ないで欲しい

尿瓶ジジイはそんなにバカにされたいのかねぇww

rm(list=ls())
library(RcppAlgos)
pm12=permuteGeneral(0:1,12,repetition = TRUE)
colnames(pm12)=LETTERS[1:12]
f12=\(x){
A=x[1];B=x[2];C=x[3];D=x[4]
E=x[5];F=x[6];G=x[7];H=x[8]
I=x[9];J=x[10];K=x[11];L=x[12]
all(
(A==1 & A==1) | (A==0 & A==0), # A 「私は正直者である」
(B==1 & A==0) | (B==1 & A==1), # B 「Aは嘘つきである」
(C==1 & B==0) | (C==1 & B==1), # C 「Bは嘘つきである」
(D==1 & C==0) | (D==1 & C==1), # D 「Cは嘘つきである」

(E==1 & A==1) | E==0, # E 「Aは正直者である」
(F==1 & B==1) | F==0, # F 「Bは正直者である」
(G==1 & C==1) | G==0, # G 「Cは正直者である」
(H==1 & D==1) | H==0, # H 「Dは正直者である」

I==1 | (I==0 & E==1), # I 「Eは嘘つきである」
J==1 | (J==0 & F==1), # J 「Fは嘘つきである」
K==1 | (K==0 & G==1), # K 「Gは嘘つきである」
L==1 | (L==0 & H==1)) # L 「Hは嘘つきである」
}

re=pm12[apply(pm12,1,f12),]
colnames(re)[colSums(re)==nrow(re)]

バク頭脳のFimoseくんには解決できない問題
>479のバグを指摘せよ。

バグ頭脳って自己紹介ですか？w

Rの単純なコードも読めないのが確定。
やっぱり、Fラン卒みたいだね。

>>482
日本語も読めないのが確定
やっぱりチンパンだね

尿瓶ジジイ医者でも東大卒でもないのにこのスレばっか伸びてるのほんと図々しいなw

尿瓶ジジイまたダンマリかよ？

# 遷移行列
P <- matrix(c(
0, 0.5, 0.5, 0, 0, # S0: 初戦 A vs B
0, 0, 0, 0.5, 0.5, # S1: A vs C
0, 0, 0, 0.5, 0.5, # S2: B vs C
0.5, 0, 0, 0, 0.5, # S3: 再試合 A vs B
0.5, 0, 0, 0.5, 0 # S4: 再試合 B vs A
), nrow = 5, byrow = TRUE)

# 固有値と固有ベクトルの計算
eig <- eigen(t(P))
stationary <- Re(eig$vectors[, which(abs(Re(eig$values) - 1) < 1e-8)])
stationary <- stationary / sum(stationary)

# 各プレイヤーの勝率を計算
A_win <- stationary[2] + stationary[4]
B_win <- stationary[3] + stationary[5]
C_win <- stationary[1]

cat("定常分布:\n")
print(stationary)
cat(sprintf("\nAの勝率: %.3f\n", A_win))
cat(sprintf("Bの勝率: %.3f\n", B_win))
cat(sprintf("Cの勝率: %.3f\n", C_win))

(* 遷移行列の定義 *)
P = {
{0, 0.5, 0.5, 0, 0}, (* S0: 初戦 A vs B *)
{0, 0, 0, 0.5, 0.5}, (* S1: A vs C *)
{0, 0, 0, 0.5, 0.5}, (* S2: B vs C *)
{0.5, 0, 0, 0, 0.5}, (* S3: 再試合 A vs B *)
{0.5, 0, 0, 0.5, 0} (* S4: 再試合 B vs A *)
};

(* 遷移行列の転置 *)
PT = Transpose[P];

(* 固有値と固有ベクトルの計算 *)
{eigValues, eigVectors} = Eigensystem[PT];

(* 固有値1に対応する固有ベクトル *)
stationaryVector = Select[eigVectors,
Norm[PT.# - #] & /@ eigValues == 1
][[1]];

(* 正規化（確率の総和を1にする） *)
stationaryVector = stationaryVector/Total[stationaryVector];

(* 各プレイヤーの勝率を計算 *)
AWin = stationaryVector[[2]] + stationaryVector[[4]];
BWin = stationaryVector[[3]] + stationaryVector[[5]];
CWin = stationaryVector[[1]];

(* 結果表示 *)
{stationaryVector, AWin, BWin, CWin}

# 遷移行列
P <- matrix(c(
0, 0.5, 0.5, 0, 0, # S0: 初戦 A vs B
0, 0, 0, 0.5, 0.5, # S1: A vs C
0, 0, 0, 0.5, 0.5, # S2: B vs C
0.5, 0, 0, 0, 0.5, # S3: 再試合 A vs B
0.5, 0, 0, 0.5, 0 # S4: 再試合 B vs A
), nrow = 5, byrow = TRUE)

# 固有値と固有ベクトルの計算
eig <- eigen(t(P))
stationary <- Re(eig$vectors[, which(abs(Re(eig$values) - 1) < 1e-8)])
stationary <- stationary / sum(stationary)

# 各プレイヤーの勝率を計算
A_win <- stationary[2] + stationary[4]
B_win <- stationary[3] + stationary[5]
C_win <- stationary[1]

cat("定常分布:\n")
print(stationary)
cat(sprintf("\nAの勝率: %.3f\n", A_win))
cat(sprintf("Bの勝率: %.3f\n", B_win))
cat(sprintf("Cの勝率: %.3f\n", C_win))

(* 遷移行列の定義 *)
P = {
{0, 0.5, 0.5, 0, 0}, (* S0: 初戦 A vs B *)
{0, 0, 0, 0.5, 0.5}, (* S1: A vs C *)
{0, 0, 0, 0.5, 0.5}, (* S2: B vs C *)
{0.5, 0, 0, 0, 0.5}, (* S3: 再試合 A vs B *)
{0.5, 0, 0, 0.5, 0} (* S4: 再試合 B vs A *)
};

(* 遷移行列の転置 *)
PT = Transpose[P];

(* 固有値と固有ベクトルの計算 *)
{eigValues, eigVectors} = Eigensystem[PT];

(* 固有値1に対応する固有ベクトル *)
stationaryVector = Select[eigVectors,
Norm[PT.# - #] & /@ eigValues == 1
][[1]];

(* 正規化（確率の総和を1にする） *)
stationaryVector = stationaryVector/Total[stationaryVector];

(* 各プレイヤーの勝率を計算 *)
AWin = stationaryVector[[2]] + stationaryVector[[4]];
BWin = stationaryVector[[3]] + stationaryVector[[5]];
CWin = stationaryVector[[1]];

(* 結果表示 *)
{stationaryVector, AWin, BWin, CWin}

>>482
自分の問題すら碌に解けないチンパンジーだったみたい

>>490
類似問題のコードを弄るだけなのに
クレクレ　レス乞食かよ。

>>491
クレクレレス乞食って自分で出した問題すら解けずに逃走するアンタのことだろ

(*
応用問題
600×nの三乗根が正整数になるｎの最小値を求めよ
演習問題
1234567890×nの五乗根が正整数になるｎの最小値を求めよ
*)
solve[x_,m_:3] := Times@@(FactorInteger[x][[All,1]]^Mod[(m-Mod[FactorInteger[x][[All,2]],m]),m])
solve[600,3]
solve[1234567890,5]

>>492
レス乞食のために答を算出するRのコードを底辺シリツ医スレにアップしておいたから
実行してみ。

>>494
これ答えじゃないよな？
レス乞食さん

EpiCycloid[R_,r_] := Module[
{circleA,circleB,redPoint,redPointTrace,frames,theta},
(* theta:Bの中心の偏角 *)
circleA = Circle[{0, 0}, R];
circleB[theta_] := Circle[{(R+r) Cos[theta], (R+r) Sin[theta]}, r];
redPoint[theta_] := {(R+r) Cos[theta] + r Cos[(r+R)/r theta + Pi], (R+r) Sin[theta] + r Sin[(R+r)/r theta + Pi]};
(* 赤点の軌跡を保存するためのリスト *)
redPointTrace = Table[redPoint[theta], {theta, 0, 2 Pi, 2 Pi/100}];
(* 円Bの中心を円Aの円周に沿って移動させるためのアニメーション *)
frames = Table[
Show[
Graphics[{EdgeForm[Black], FaceForm[None], circleA}],
Graphics[{EdgeForm[None], FaceForm[Red], circleB[theta]}],
Graphics[{Red, PointSize[Large], Point[redPoint[theta]]}],
Graphics[{Red, Line[redPointTrace[[1 ;; Round[theta/ (2 Pi/100)]]]]}],
PlotRange -> {{-2(R+r), 2(R+r)}, {-2(R+r), 2(R+r)}},
Axes -> True,
AxesOrigin -> {0, 0}
],
{theta, 0, 2 Pi, 2 Pi/100}];
(* アニメーション表示 *)
ListAnimate[frames]
]
EpiCycloid[3,1]


HypoCycloid[R_,r_] := Module[
{circleA,circleB,redPoint,theta,frames},
(* theta 円Bの中心の偏角 *)
circleA = Circle[{0, 0}, R];
circleB[theta_] := Circle[{(R - r) Cos[theta], (R - r) Sin[theta]}, r];
redPoint[theta_] := {(R-r) Cos[theta] + r Cos[(1-R/r) theta], (R-r) Sin[theta] + r Sin[(1-R/r) theta]};
(* 円Bの中心を円Aの円周に沿って移動させるためのアニメーション *)
frames = Table[
Show[
Graphics[{EdgeForm[Black], FaceForm[None], circleA}],
Graphics[{EdgeForm[None], FaceForm[Red], circleB[theta]}],
Graphics[{Red, PointSize[Large], Point[redPoint[theta]]}],
Graphics[{Red, Line[Table[redPoint[theta], {theta, 0, 2 Pi, 2 Pi/100}]]}],
PlotRange -> {{-R, R}, {-R, R}},
Axes -> True,
AxesOrigin -> {0, 0}
],
{theta, 0, 2 Pi, 2 Pi/100}
];
(* アニメーションを表示 *)
ListAnimate[frames]
]
HypoCycloid[3,1]

>>495
実行すれば答がでるだろ。
証言をすべてみたす組み合わせのすべてで嘘つきになるのを列挙すればいい。
Rコード読めないの？

>>497
結局グダグダ言い訳してるが自分の問題すら答えられないアホってことね

>>497
さすが答えの意味分かってないwww

ハイポサイクロイド
https://i.imgur.com/Y699ttV.gif

>>500
年明けても日本語読めないみたいだな

https://i.imgur.com/PiWG2FN.gif

高校生にすら相手にされてないみたいだね
実に哀れ

https://i.imgur.com/XfJfU12.gif

>>482
荒らすの止めてくれないか

さっさと精神科行け

R=3;
r=2;
K=2; (*　ロータ中心と頂点距離/ロータギア半径 *)

apex[t_] := Module[
{p,c,d},
p = {(R-r) Cos[t] + R Cos[(t+Pi)/3], (R-r) Sin[t] + R Sin[(t+Pi)/3]};
c = {(R-r) Cos[t],(R-r) Sin[t]};
d = K(p-c)+c
]
theta=Flatten@{Range[0,6 Pi, Pi/1000]};
ListPlot[apex[#]& /@ theta,AspectRatio ->1]
x[t_]:=apex[t][[1]]
y[t_]:=apex[t][[2]]
length=Integrate[(D[x[t],t])^2+(D[y[t],t])^2,{t,0,6 Pi}]
area=(1/2) Integrate[x[t] D[y[t],t] - y[t] D[x[t],t],{t,0, 6 Pi}]

Clear[{R,r,K}]

>>507
まだバカにされたりないのか尿瓶ジジイは
しかし正規のスレじゃもうバカにすらしてもらえないみたいww

【問題】
固定円の半径を２、回転ローターの円の半径を３とするときに
ルーローの三角の軌跡の内部でローターが移動できるには三角の頂点とロータ円の中心の距離が
回転ローターの円の半径の何倍以上であることが必要ですか？

> 2/3*(2+sqrt(3))
[1] 2.488034

>>509
高校数学で回転ローターとやらが定義されてるのかよwww
中卒以下の低脳でまともに数学用語使えないゴミなのはしょうがないけど、もう少し努力しろよwww

ルーローって魯肉のことだな