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

**１３２人目の素数さん** · 2024/08/09(金) 06:22:27.38

【質問者必読！！】
まず>>1-5をよく読んでね

このスレは医者・東大卒の人物専用スレです。
その他の人が書き込むことは許されません。

数学@5ch掲示板用掲示板での数学記号の書き方例と一般的な記号の使用例
http://mathmathmath.dotera.net/

・まずは教科書、参考書、web検索などで調べるようにしましょう。（特に基本的な公式など）
・問題の写し間違いには気をつけましょう。
・長い分母分子を含む分数はきちんと括弧でくくりましょう。
　　（× x+1/x+2　；　 ○((x+1)/(x+2)) ）
・丸文字、顔文字、その他は環境やブラウザによりうまく表示できない場合があります。
　どうしても画像を貼る場合はPCから直接見られるところに見やすい画像を貼ってください。
　ピクトはPCから見られないことがあるので避けてください。
・質問者は名前を騙られたくない場合、トリップを付けましょう。
　（トリップの付け方は　名前(N)に　俺！#oretrip　←適当なトリ)
・質問者は回答者がわかるように問題を書くようにしましょう。
　でないと放置されることがあります。
　（変に省略するより全文書いた方がいい、また説明なく習慣的でない記号を使わないように）
・質問者は何が分からないのか、どこまで考えたのかを明記しましょう。
　それがない場合、放置されることがあります。
　（特に、自分でやってみたのに合わないので教えてほしい、みたいなときは必ず書くように）
・回答者も節度ある回答を心がけてください。
・出題スレではありません。出題は該当スレにお願いします。
・970くらいになったら次スレを立ててください。
※前スレ
高校数学の質問スレ Part437
https://rio2016.5ch.net/test/read.cgi/math/1721071007/

**１３２人目の素数さん** · 2024/12/22(日) 12:35:08.21

>>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)

**１３２人目の素数さん** · 2024/12/22(日) 12:40:56.83

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

**１３２人目の素数さん** · 2024/12/22(日) 12:42:11.71

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

**１３２人目の素数さん** · 2024/12/22(日) 13:49:41.40

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

**１３２人目の素数さん** · 2024/12/22(日) 14:57:09.36

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

**１３２人目の素数さん** · 2024/12/22(日) 15:12:10.82

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

**１３２人目の素数さん** · 2024/12/22(日) 16:37:45.31

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

**１３２人目の素数さん** · 2024/12/22(日) 17:42:21.06

やはり、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]

**１３２人目の素数さん** · 2024/12/22(日) 17:55:17.31

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

**１３２人目の素数さん** · 2024/12/22(日) 17:55:22.37

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

**１３２人目の素数さん** · 2024/12/22(日) 18:41:54.83

>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)

**１３２人目の素数さん** · 2024/12/22(日) 19:03:14.31

>>463
死ね

**１３２人目の素数さん** · 2024/12/22(日) 20:14:10.65

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

**１３２人目の素数さん** · 2024/12/22(日) 21:03:18.86

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
}

**１３２人目の素数さん** · 2024/12/22(日) 22:19:23.50

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

**１３２人目の素数さん** · 2024/12/22(日) 22:32:49.57

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

**１３２人目の素数さん** · 2024/12/23(月) 08:08:07.50

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

**１３２人目の素数さん** · 2024/12/28(土) 09:03:58.56

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

**１３２人目の素数さん** · 2024/12/28(土) 10:20:03.06

>>470

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

**１３２人目の素数さん** · 2024/12/28(土) 12:26:13.98

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

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

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

**１３２人目の素数さん** · 2024/12/28(土) 12:27:13.74

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()))

**１３２人目の素数さん** · 2024/12/28(土) 12:43:10.85

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()))

**１３２人目の素数さん** · 2024/12/28(土) 14:17:39.00

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

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

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

**１３２人目の素数さん** · 2024/12/28(土) 14:21:31.42

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

**１３２人目の素数さん** · 2024/12/28(土) 18:11:22.23

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

**１３２人目の素数さん** · 2024/12/28(土) 19:00:53.56

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

**１３２人目の素数さん** · 2024/12/29(日) 07:54:27.11

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)]

**１３２人目の素数さん** · 2024/12/29(日) 11:24:48.82

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

**１３２人目の素数さん** · 2024/12/29(日) 11:25:58.82

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

**１３２人目の素数さん** · 2024/12/29(日) 12:22:06.12

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

**１３２人目の素数さん** · 2024/12/29(日) 12:47:08.60

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

**１３２人目の素数さん** · 2024/12/29(日) 14:05:39.63

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

**１３２人目の素数さん** · 2024/12/29(日) 23:58:05.11

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

**１３２人目の素数さん** · 2024/12/30(月) 03:14:24.14

# 遷移行列
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))

**１３２人目の素数さん** · 2024/12/30(月) 03:17:30.79

(* 遷移行列の定義 *)
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}

**１３２人目の素数さん** · 2024/12/30(月) 03:31:36.28

# 遷移行列
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))

**１３２人目の素数さん** · 2024/12/30(月) 03:34:00.29

(* 遷移行列の定義 *)
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}

**１３２人目の素数さん** · 2024/12/30(月) 10:41:19.52

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

**１３２人目の素数さん** · 2024/12/30(月) 21:13:17.75

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

**１３２人目の素数さん** · 2024/12/30(月) 21:39:27.87

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

**１３２人目の素数さん** · 2024/12/31(火) 07:15:10.51

(*
応用問題
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]

**１３２人目の素数さん** · 2024/12/31(火) 07:15:43.90

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

**１３２人目の素数さん** · 2024/12/31(火) 10:02:09.35

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

**１３２人目の素数さん** · 2024/12/31(火) 14:55:01.89

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]

**１３２人目の素数さん** · 2024/12/31(火) 14:56:07.64

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

**１３２人目の素数さん** · 2024/12/31(火) 18:47:36.81

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

**１３２人目の素数さん** · 2024/12/31(火) 23:12:32.74

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

**１３２人目の素数さん** · 2025/01/01(水) 08:56:34.48

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

**１３２人目の素数さん** · 2025/01/01(水) 11:12:05.21

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

**１３２人目の素数さん** · 2025/01/03(金) 19:11:24.08

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

**１３２人目の素数さん** · 2025/01/03(金) 20:33:58.18

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

**１３２人目の素数さん** · 2025/01/03(金) 22:47:45.54

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

**１３２人目の素数さん** · 2025/01/03(金) 23:25:35.89

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

**１３２人目の素数さん** · 2025/01/03(金) 23:25:35.96

さっさと精神科行け

**１３２人目の素数さん** · 2025/01/05(日) 12:23:25.71

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}]

**１３２人目の素数さん** · 2025/01/05(日) 12:47:54.51

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

**１３２人目の素数さん** · 2025/01/07(火) 19:53:08.02

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

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

**１３２人目の素数さん** · 2025/01/07(火) 22:35:53.49

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

**１３２人目の素数さん** · 2025/01/08(水) 02:59:50.29

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

**１３２人目の素数さん** · 2025/01/08(水) 04:11:37.72

rm(list=ls())

JugRiddle <- function(
big=5,
small=3,
end=4){

# starting from the bigger jug
movebig <- function(xy){ # start from c(big,0)
x=xy[1] ; y=xy[2]
# x==big
if(x==big) re=c(big-(small-y),small)
# x==0
if(x==0) re=c(big,y)
# y==small
if(y==small) re=c(x,0)
# y==0
if(y==0 & x!=big){
if(x>=small) re=c(x-small,small)
else re=c(0,x)
}
return(re)
}

STATUS=status=c(big,0)
i=1
while(!identical(status,c(0,0))){#
i=i+1
status=movebig(status)
STATUS=rbind(STATUS,status)
}
rownames(STATUS)=1:nrow(STATUS)
colnames(STATUS)=c(paste0(big,'L'),paste0(small,'L'))
(Bigger=STATUS)

# starting from the smaller jug
movesmall <- function(xy){ # start from c(0,small)
x=xy[1] ; y=xy[2]
if(y==small){
if(x<=(big-small)) re=c(x+small,0)
else re=c(big, small-(big-x))
}
if(y==0) re=c(x,small)
if(x==big) re=c(0,y)
if(x==0) re=c(y,0)
return(re)
}

STATUS=status=c(0,small)
i=1
while(!identical(status,c(0,0))){ #
i=i+1
status=movesmall(status)
STATUS=rbind(STATUS,status)
}
rownames(STATUS)=1:nrow(STATUS)
colnames(STATUS)=c(paste0(big,'L'),paste0(small,'L'))
(Smaller=STATUS)
if(all(end != c(Bigger,Smaller))) return(NA)

min_Bigger=min(which(apply(Bigger,1,function(x) end %in% x)))
min_Smaller=min(which(apply(Smaller,1,function(x) end %in% x)))

list(Bigger=as.matrix(Bigger)[1:min_Bigger,],
Smaller=as.matrix(Smaller)[1:min_Smaller,],
min_Bigger=min_Bigger,min_Smaller=min_Smaller)
}

JugRiddle(5,3,4)
JugRiddle(10,5,3)

**１３２人目の素数さん** · 2025/01/08(水) 19:26:21.90

コイントスの終了条件を
Aくんは、　表が2回連続して出たら終了
Bくんは、　表裏の順に出たら終了（裏表の順では終了しない）
することになった。
(1) Aくん、Bくんについて、コイントスの回数の期待値を求めよ。
(2) Aくんはの方が少ない回数で終了する確率を求めよ。小数解でよい。

**１３２人目の素数さん** · 2025/01/08(水) 19:29:44.12

>>513
（２）回数の差の分布
https://i.imgur.com/jDZubdl.png

**１３２人目の素数さん** · 2025/01/08(水) 20:16:25.91

>>514
https://i.imgur.com/cmdTCLY.png

**１３２人目の素数さん** · 2025/01/08(水) 20:27:25.85

日本語すら読めないチンパンジーが数学やった気になってるのか

**１３２人目の素数さん** · 2025/01/08(水) 20:27:57.47

>>513
1000万回のシミュレーション
https://i.imgur.com/JKpHsKb.png
https://i.imgur.com/A1NiFF0.png

**１３２人目の素数さん** · 2025/01/08(水) 20:28:49.92

>>516
罵倒しか投稿できないのがFimoseくん。FランにちなんでPhimoseくんから昇格ｗ

**１３２人目の素数さん** · 2025/01/08(水) 20:30:20.34

>>517
期待値の差よりも
回数の差の分布の方が直感に近いな。

**１３２人目の素数さん** · 2025/01/08(水) 20:33:00.33

チンパンジーが図星で発狂中w

**１３２人目の素数さん** · 2025/01/09(木) 00:40:37.23

>>518
他人がどうこうなんてお前が無能な事とは全く関係ないんだから、まずは国語の勉強してスレタイくらい読めるようになれよwwww
汚い言葉勉強してる場合じゃないぞwww

**１３２人目の素数さん** · 2025/01/09(木) 06:57:53.73

やはり、下手な考え休むに似たりだな。

オッペンハイマーでの会話

Hahn and Strassmann in Germany.

They split the uranium nucleus.

How?

Bombard it with neutrons.

It's a nuclear fission.

They did it,
they split the atom.

It's not possible.

I'm gonna try to reproduce it.

See? Can't be done.

Very elegant. Quite clear.

**１３２人目の素数さん** · 2025/01/09(木) 07:00:06.18

Fimoseくんはロータリーエンジンの動くシェーマの動画を作れないの？
完成した動画をみてもハウジングが何かわからんようだし。
やっぱり、弄れるのはforeskinだけみたいだね。

**１３２人目の素数さん** · 2025/01/09(木) 08:14:29.03

>>522
アンタのことじゃん

**１３２人目の素数さん** · 2025/01/09(木) 08:14:47.34

limit = 10000

max_count = 0
best_n = 0

for (n in 1:(limit - 1)) {
count = 0
for (y in 1:(n - 1)) {
if ((n*y) %% (n-y) ==0 & n!=y){ # x が整数になる条件
count = count + 1
}
}
if (count > max_count) {
max_count = count
best_n = n
}
}

best_n
max_count

**１３２人目の素数さん** · 2025/01/09(木) 08:16:10.89

nを10000未満の正整数とする。
1/n+1/x=1/y を満たす{x,y}の正整数解の組み合わせの数が最大なnはいくつか？

**１３２人目の素数さん** · 2025/01/09(木) 12:34:42.59

こういう報道に接するとZ3が必須だな。
https://www.373news.com/_news/storyid/155006/

**１３２人目の素数さん** · 2025/01/09(木) 12:35:57.90

>>524
映画みてないの？
英語がわからないFラン？

**１３２人目の素数さん** · 2025/01/09(木) 17:28:07.32

>>528
下手な考え休むに似たりって完全にアンタのことなんだけどまさか自覚ないんか？ww

**１３２人目の素数さん** · 2025/01/10(金) 02:04:06.88

>>523
他人がどうこうとか自分の無能を誤魔化すばっかりwww
自分で定義分かってるとか思い込んでるから何にも理解出来ないんだよwww

思い込みじゃなくて自分できちんと日本語使って数学的な定義述べてみろよwww
出来ないんだろwww

**１３２人目の素数さん** · 2025/01/10(金) 02:06:10.84

Fimoseくんはロータリーエンジンの動くシェーマの動画を作れないの？
完成した動画をみてもハウジングが何かわからんようだし。
やっぱり、弄れるのはforeskinだけ。

**１３２人目の素数さん** · 2025/01/10(金) 02:14:15.07

結局できなくて発狂かよ
自分で出したクソの処理すらできない模様

**１３２人目の素数さん** · 2025/01/10(金) 02:18:01.93

>>531
汚い言葉で誤魔化すばっかりかよwww
他人がどうこうなんて関係なくお前が無能な事を指摘して貰ってるんだから誤魔化さずに自分で努力してみろよwww

**１３２人目の素数さん** · 2025/01/10(金) 02:31:21.32

毎回毎回同じことしか言えないし知性のかけらもないのに自称医者東大卒とは笑わせるねww

**１３２人目の素数さん** · 2025/01/10(金) 05:36:36.58

>>534
俺は医科歯科卒
理１を蹴って入学。
あんたは動画すら作れないFランなんだろ。
このスレは出禁だぞ。

**１３２人目の素数さん** · 2025/01/10(金) 07:02:22.75

>>535
医科歯科卒wwww
日本語すらまともに使えない中卒以下の低脳が何言ってんの？
Fラン受験出来るのが羨ましいのかよwww

**１３２人目の素数さん** · 2025/01/10(金) 08:29:22.85

医科歯科ってスレタイ読めなくても卒業できるの？
ルールに則らないことを推奨してるの？
医師国家試験でルールに則らないことを是とする回答をしたの？

**１３２人目の素数さん** · 2025/01/10(金) 08:30:21.85

質問スレで出題ばかりするし質問には回答しないしで、
世の中のいなくなったほうが世界か良くなる人物の中の一人

**１３２人目の素数さん** · 2025/01/10(金) 09:16:14.22

>>535
まあ匿名掲示板ならいくらでも経歴盛れるもんなそんな主張意味ないよ
本当に医科歯科卒ならそれ相応の能力を示せよ

**１３２人目の素数さん** · 2025/01/10(金) 11:01:30.43

今日は内視鏡が少なくて10:30に終了。予約の検査終了したら帰っていい契約なので( ・∀・)ｲｲ!!

**１３２人目の素数さん** · 2025/01/10(金) 11:16:40.37

熟練職人不足で新年から内視鏡週3回になった。
スタッフが可愛いので続けられる。

**１３２人目の素数さん** · 2025/01/10(金) 11:31:20.92

本物の医者がこんな数学板の糞スレにいるわけないじゃん忙しいはずなのに妄想乙

**１３２人目の素数さん** · 2025/01/10(金) 11:56:51.07

>>540
これが何かの証拠になるとか思ってんの？
説得力皆無で今まで誰一人信じて貰ってないのが現実だろwwww

何の意味もない事繰り返し書き込んじゃうって知能低すぎん？

**１３２人目の素数さん** · 2025/01/10(金) 18:50:35.59

尿瓶ジジイ急にダンマリかよ

**１３２人目の素数さん** · 2025/01/11(土) 07:06:37.50

A君とB君がそれぞれ同じ規格立方体のサイコロを投げます。
終了条件は
　A君は、　同じ数の目が2回連続して出たら終了
　B君は、　前の目の数より１多い数がでたら修了（但し、６→１でも終了するとする）

【問題】
(1) A君の投げた回数の期待値を求めよ
(2) B君の投げた回数の期待値を求めよ
(3) A君の投げた回数の方が多い確率をもとめよ
(4) 二人の投げた回数が同数で終了する確率をもとめよ　

**１３２人目の素数さん** · 2025/01/11(土) 08:13:55.14

(*
コイントスの終了条件を
Aくんは、　表が2回連続して出たら終了
Bくんは、　表裏の順に出たら終了（裏表の順では終了しない）
する。
*)

pA[k_]:=Fibonacci[k-1]/2^k (* Aの終了回数ごとの確率 *)
Sum[pA[k],{k,2,Infinity}] (* Σpmf==1 確認 *)
Sum[k pA[k],{k,2,Infinity}](*　期待値 *)

pB[k_]:=(k-1)/2^k (* Bの終了回数ごとの確率 *)
Sum[pB[k],{k,2,Infinity}](* Σpmf==1 確認 *)
Sum[k pB[k],{k,2,Infinity}](* 期待値 *)

pABd[d_] := Sum[pA[k+d] pB[k],{k,2,Infinity}] (* Aの方がd回多く終わる確率 *)
pABd[0] (* 同じ回数で終わる確率 *)

(* Aの方が多く終わる確率 *)
Sum[pABd[d],{d,1,Infinity}] // Simplify
pbA[k_]:=Sum[pB[j],{j,2,k-1}] (* A=ｋのときのB<Aの確率　*)
Sum[pA[k] pbA[k],{k,2,Infinity}]

**１３２人目の素数さん** · 2025/01/11(土) 09:15:51.78

>>545
結局誰にも信じられずに哀れだね

**１３２人目の素数さん** · 2025/01/11(土) 19:01:33.32

(*
1ヶ月間に馬に蹴られて死亡する兵士の数が母数λ=2のポアソン分布に従い、
1ヶ月間に補充される兵士の数はp=1/3の幾何分布に従うとする。
いずれの分布も定義域は非負整数。
1ヶ月後に兵士が増えている確率を算出し、シミュレーションとの合致を確認せよ。
*)

p=1/3;
lambda=2;
pDiff[d_] := Sum[PDF[GeometricDistribution[p],k+d] PDF[PoissonDistribution[lambda],k],{k,0,Infinity}]
Sum[pDiff[d],{d,1,Infinity}]
% // N

k=10^6;
po=RandomVariate[PoissonDistribution[lambda],k];
ge=RandomVariate[GeometricDistribution[p],k];
d=ge-po;
Boole[#>0& /@ d]//Mean // N

**１３２人目の素数さん** · 2025/01/11(土) 21:32:12.61

>>548
マルチポストしてるのに誰にも相手にされてないみたい
実に哀れ

**１３２人目の素数さん** · 2025/01/11(土) 23:12:46.87

COVID-19の潜伏期はμ=1.6 、σ=0.5の対数正規分布(オレンジ色)に
インフルエンザの潜伏期は形状母数k=2.0,尺度母数θ=1.0のガンマ分布(青色)に従うとする。
https://i.imgur.com/QmlK7yi.png
同時に感染したときにインフルエンザの方が先に発症する確率を求めなさい。
分布は独立でウイルス干渉はないものとする。

pdfz =\(z) integrate( \(x) dlnorm(x,1.6,0.5) * dgamma(x-z,shape=2,scale=1),
0,Inf,rel.tol = 1e-12)$value
pdfz=Vectorize(pdfz)
curve(pdfz(x),-10,30)
1-integrate(pdfz,-Inf,0,rel.tol = 1e-12)$value

k=1e7
Cov=rlnorm(k,1.6,0.5)
Flu=rgamma(k,shape=2,scale=1)
plotPost(Cov-Flu,compVal = 0,col=7)
curve(pdfz(x),add=TRUE,lwd=2)
mean(Cov>Flu)

https://i.imgur.com/NKroLQc.png

**１３２人目の素数さん** · 2025/01/11(土) 23:13:48.64

ここはRやWolframによる解法のコード置き場。
亀レスとの照合用。

**１３２人目の素数さん** · 2025/01/12(日) 08:08:57.37

>>551
ちがいます
>>5