高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
449: １３２人目の素数さん [sage] 2024/12/22(日) 11:08:27.94 ID:iBxkrbTc(2/3) AAS
>>445
煽って乞食行為とかテンプレいらないぞ
少しはオリジナリティのある返しをしてくれないとwww
450: １３２人目の素数さん [sage] 2024/12/22(日) 11:09:55.89 ID:iBxkrbTc(3/3) AAS
その論理なら、クズコード書き込んでるやつはwolfram使えないの確定だぞwww
自白かよwww
さすが中卒www
451: １３２人目の素数さん [sage] 2024/12/22(日) 12:28:31.05 ID:PT9UhnTA(1) AAS
>>447
なに自分で自演してんの？馬鹿なの？
それとも自分がFラン卒だって自白したかっただけ？
452(3): １３２人目の素数さん [sage] 2024/12/22(日) 12:33:00.32 ID:rWRciXK9(3/8) AAS
>>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
最小最大のときの形状が判明したので
あとは一般解を求めればいい。
453(1): １３２人目の素数さん [sage] 2024/12/22(日) 12:35:08.21 ID:rWRciXK9(4/8) AAS
>>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)
454: １３２人目の素数さん [sage] 2024/12/22(日) 12:40:56.83 ID:rWRciXK9(5/8) AAS
>>453
算出のために作った関数の動作確認の図
https://i.imgur.com/0mvHN8w.png
455: １３２人目の素数さん [sage] 2024/12/22(日) 12:42:11.71 ID:rWRciXK9(6/8) AAS
毒薬が２本のときのWolframのコード まだぁ？
456: １３２人目の素数さん [sage] 2024/12/22(日) 13:49:41.40 ID:pkRHj4OZ(1) AAS
>>452
rが1/2に近いとき、明らかに長方形の時最小なるような…
457: １３２人目の素数さん [sage] 2024/12/22(日) 14:57:09.36 ID:S8dVAroW(1) AAS
別スレを潰したのに飽き足らずこっちのスレも潰す気かよ
458: １３２人目の素数さん [sage] 2024/12/22(日) 15:12:10.82 ID:rKRhCRuS(1/2) AAS
>>452
横軸にｒ、縦軸にPQRSの面積の最大値（青）最小値（赤）をグラフ化。
https://i.imgur.com/iio0LJI.png
Rのコードが読めれば、単純作業。（但し、Fラン卒は除く）。
459: １３２人目の素数さん [sage] 2024/12/22(日) 16:37:45.31 ID:+EIC8tFQ(1) AAS
なんか発狂して書き込んでるみたいwww
クズ問題読んで貰えると思ってるのウケるwww
自意識過剰すぎだろwwww
460(2): １３２人目の素数さん [sage] 2024/12/22(日) 17:42:21.06 ID:rKRhCRuS(2/2) AAS
やはり、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]
461: １３２人目の素数さん [sage] 2024/12/22(日) 17:55:17.31 ID:yf8QUsfH(1) AAS
>>460
クズはいろんなスレに関係ないことしか書き込まないから総スカンにされてバカにされてるアンタのことだよ
462: １３２人目の素数さん [sage] 2024/12/22(日) 17:55:22.37 ID:gaGlxwVG(2/3) AAS
>>460
やはりスレタイ読めない
保育園卒確定！
463(2): １３２人目の素数さん [sage] 2024/12/22(日) 18:41:54.83 ID:rWRciXK9(7/8) AAS
>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)
464: １３２人目の素数さん [sage] 2024/12/22(日) 19:03:14.31 ID:bA66atg/(1) AAS
>>463
死ね
465: １３２人目の素数さん [sage] 2024/12/22(日) 20:14:10.65 ID:gaGlxwVG(3/3) AAS
>>463
スレタイに則った書き込みすればいいのに
それもできないのでホイ卒確定
466(1): １３２人目の素数さん [sage] 2024/12/22(日) 21:03:18.86 ID:rWRciXK9(8/8) AAS
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
}
467: １３２人目の素数さん [sage] 2024/12/22(日) 22:19:23.50 ID:ODwmJOyZ(1) AAS
https://i.imgur.com/Q34Vduy.jpeg
>>466
468: １３２人目の素数さん [sage] 2024/12/22(日) 22:32:49.57 ID:DH2naynY(1) AAS
類は何とやら
キチガイはキチガイをよぶ
469: １３２人目の素数さん [sage] 2024/12/23(月) 08:08:07.50 ID:5pyB+nba(1) AAS
キチガイが複数人いるより、キチガイの自演の方が自然
470(1): １３２人目の素数さん [sage] 2024/12/28(土) 09:03:58.56 ID:pTaNTUpL(1) AAS
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
471: １３２人目の素数さん [sage] 2024/12/28(土) 10:20:03.06 ID:byAnFSyo(1/2) AAS
>>470
137：卵の名無しさん：[sage]：2024/12/27(金) 16:30:06.83 ID:UPA/6UxQ
レス乞食じゃなきゃ然るべきスレでやるもんね
472: １３２人目の素数さん [sage] 2024/12/28(土) 12:26:13.98 ID:ABHDoLY/(1/3) AAS
横綱、大関、関脇で巴戦を行う。
２力士が対戦して，勝った方が残りの 1 力士と対戦をする。
これを繰り返して，２連勝した力士が優勝をする。
勝利確率は通算成績に従うとする（巴戦での成績も組み入れて計算する）
巴戦前の通算成績は
横綱 対 大関は 横綱の１０勝９負
大関 対 関脇は 大関の ８勝７負
横綱 対 関脇は 横綱の ６勝５負
とする。
大関 対 関脇で巴戦を始めるときにそれぞれの力士の優勝する確率を算出しなさい。
473: １３２人目の素数さん [sage] 2024/12/28(土) 12:27:13.74 ID:ABHDoLY/(2/3) AAS
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()))
474: １３２人目の素数さん [sage] 2024/12/28(土) 12:43:10.85 ID:ABHDoLY/(3/3) AAS
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()))
475: １３２人目の素数さん [sage] 2024/12/28(土) 14:17:39.00 ID:1sBW2m9k(1) AAS
"
横綱、大関、関脇で巴戦を行う。
２力士が対戦して，勝った方が残りの 1 力士と対戦をする。
これを繰り返して，２連勝した力士が優勝をする。
勝利確率は通算成績に従うとする（巴戦での成績も組み入れて計算する）
巴戦前の通算成績は
横綱 対 大関は 横綱の１０勝９負
大関 対 関脇は 大関の ８勝７負
横綱 対 関脇は 横綱の ６勝５負
とする。
(1) 大関 対 関脇で巴戦を始めるときにそれぞれの力士の優勝する確率を算出しなさい。
(2) 巴戦の試合数の期待値を求めなさい。
476(1): １３２人目の素数さん [sage] 2024/12/28(土) 14:21:31.42 ID:2IDswPeZ(1) AAS
尿瓶ジジイのレス乞食はいつになったら願いが届くのかな？w
477: １３２人目の素数さん [sage] 2024/12/28(土) 18:11:22.23 ID:wnhYoY4/(1) AAS
>>476
あなたが返信してる時点で目的を達してる
完全スルーができないならスレを見ないで欲しい
478: １３２人目の素数さん [sage] 2024/12/28(土) 19:00:53.56 ID:byAnFSyo(2/2) AAS
尿瓶ジジイはそんなにバカにされたいのかねぇww
479(1): １３２人目の素数さん [sage] 2024/12/29(日) 07:54:27.11 ID:/ygSVCno(1) AAS
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)]
480: １３２人目の素数さん [sage] 2024/12/29(日) 11:24:48.82 ID:8wG+CYHS(1) AAS
バク頭脳のFimoseくんには解決できない問題
>479のバグを指摘せよ。
481: １３２人目の素数さん [sage] 2024/12/29(日) 11:25:58.82 ID:bEHNzlJY(1/2) AAS
バグ頭脳って自己紹介ですか？w
482(3): １３２人目の素数さん [sage] 2024/12/29(日) 12:22:06.12 ID:90INKYGd(1) AAS
Rの単純なコードも読めないのが確定。
やっぱり、Fラン卒みたいだね。
483: １３２人目の素数さん [sage] 2024/12/29(日) 12:47:08.60 ID:lQcX1F1j(1) AAS
>>482
日本語も読めないのが確定
やっぱりチンパンだね
484: １３２人目の素数さん [sage] 2024/12/29(日) 14:05:39.63 ID:bEHNzlJY(2/2) AAS
尿瓶ジジイ医者でも東大卒でもないのにこのスレばっか伸びてるのほんと図々しいなw
485: １３２人目の素数さん [sage] 2024/12/29(日) 23:58:05.11 ID:/vJ/3Q0O(1) AAS
尿瓶ジジイまたダンマリかよ？
486: １３２人目の素数さん [sage] 2024/12/30(月) 03:14:24.14 ID:DAgtZ3Zg(1/5) AAS
# 遷移行列
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))
487: １３２人目の素数さん [sage] 2024/12/30(月) 03:17:30.79 ID:DAgtZ3Zg(2/5) AAS
(* 遷移行列の定義 *)
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}
488: １３２人目の素数さん [sage] 2024/12/30(月) 03:31:36.28 ID:DAgtZ3Zg(3/5) AAS
# 遷移行列
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))
489: １３２人目の素数さん [sage] 2024/12/30(月) 03:34:00.29 ID:DAgtZ3Zg(4/5) AAS
(* 遷移行列の定義 *)
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}
490(1): １３２人目の素数さん [sage] 2024/12/30(月) 10:41:19.52 ID:dmM3LvbP(1/2) AAS
>>482
自分の問題すら碌に解けないチンパンジーだったみたい
491(1): １３２人目の素数さん [sage] 2024/12/30(月) 21:13:17.75 ID:DAgtZ3Zg(5/5) AAS
>>490
類似問題のコードを弄るだけなのに
クレクレ レス乞食かよ。
492(1): １３２人目の素数さん [sage] 2024/12/30(月) 21:39:27.87 ID:dmM3LvbP(2/2) AAS
>>491
クレクレレス乞食って自分で出した問題すら解けずに逃走するアンタのことだろ
493: １３２人目の素数さん [sage] 2024/12/31(火) 07:15:10.51 ID:41bOMap+(1/4) AAS
(*
応用問題
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]
494(1): １３２人目の素数さん [sage] 2024/12/31(火) 07:15:43.90 ID:41bOMap+(2/4) AAS
>>492
レス乞食のために答を算出するRのコードを底辺シリツ医スレにアップしておいたから
実行してみ。
495(1): １３２人目の素数さん [sage] 2024/12/31(火) 10:02:09.35 ID:6B6ig4bh(1) AAS
>>494
これ答えじゃないよな？
レス乞食さん
496: １３２人目の素数さん [sage] 2024/12/31(火) 14:55:01.89 ID:41bOMap+(3/4) AAS
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]
497(2): １３２人目の素数さん [sage] 2024/12/31(火) 14:56:07.64 ID:41bOMap+(4/4) AAS
>>495
実行すれば答がでるだろ。
証言をすべてみたす組み合わせのすべてで嘘つきになるのを列挙すればいい。
Rコード読めないの？
498: １３２人目の素数さん [sage] 2024/12/31(火) 18:47:36.81 ID:B87FjDFd(1) AAS
>>497
結局グダグダ言い訳してるが自分の問題すら答えられないアホってことね
499: １３２人目の素数さん [sage] 2024/12/31(火) 23:12:32.74 ID:PjhLbT7v(1) AAS
>>497
さすが答えの意味分かってないwww
500(1): １３２人目の素数さん [sage] 2025/01/01(水) 08:56:34.48 ID:trf8HhUX(1) AAS
ハイポサイクロイド
https://i.imgur.com/Y699ttV.gif
501: １３２人目の素数さん [sage] 2025/01/01(水) 11:12:05.21 ID:vYaE/f7n(1) AAS
>>500
年明けても日本語読めないみたいだな
502: １３２人目の素数さん [sage] 2025/01/03(金) 19:11:24.08 ID:KUn71WA0(1/2) AAS
https://i.imgur.com/PiWG2FN.gif
503: １３２人目の素数さん [sage] 2025/01/03(金) 20:33:58.18 ID:gmyxpwxe(1/2) AAS
高校生にすら相手にされてないみたいだね
実に哀れ
504: １３２人目の素数さん [sage] 2025/01/03(金) 22:47:45.54 ID:KUn71WA0(2/2) AAS
https://i.imgur.com/XfJfU12.gif
505: １３２人目の素数さん [sage] 2025/01/03(金) 23:25:35.89 ID:rOKue8hC(1) AAS
>>482
荒らすの止めてくれないか
506: １３２人目の素数さん [sage] 2025/01/03(金) 23:25:35.96 ID:gmyxpwxe(2/2) AAS
さっさと精神科行け
507(1): １３２人目の素数さん [sage] 2025/01/05(日) 12:23:25.71 ID:tuC7O64g(1) AAS
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}]
508: １３２人目の素数さん [sage] 2025/01/05(日) 12:47:54.51 ID:/2bPdu/v(1) AAS
>>507
まだバカにされたりないのか尿瓶ジジイは
しかし正規のスレじゃもうバカにすらしてもらえないみたいww
509(1): １３２人目の素数さん [sage] 2025/01/07(火) 19:53:08.02 ID:rJfoFeSd(1) AAS
【問題】
固定円の半径を２、回転ローターの円の半径を３とするときに
ルーローの三角の軌跡の内部でローターが移動できるには三角の頂点とロータ円の中心の距離が
回転ローターの円の半径の何倍以上であることが必要ですか？
> 2/3*(2+sqrt(3))
[1] 2.488034
510: １３２人目の素数さん [sage] 2025/01/07(火) 22:35:53.49 ID:zsZUtuIb(1) AAS
>>509
高校数学で回転ローターとやらが定義されてるのかよwww
中卒以下の低脳でまともに数学用語使えないゴミなのはしょうがないけど、もう少し努力しろよwww
511: １３２人目の素数さん [sage] 2025/01/08(水) 02:59:50.29 ID:0uIf2Djd(1) AAS
ルーローって魯肉のことだな
512: １３２人目の素数さん [sage] 2025/01/08(水) 04:11:37.72 ID:X/URcmTf(1/3) AAS
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)
513(2): １３２人目の素数さん [sage] 2025/01/08(水) 19:26:21.90 ID:SKdwPI7I(1/4) AAS
コイントスの終了条件を
Aくんは、 表が2回連続して出たら終了
Bくんは、 表裏の順に出たら終了（裏表の順では終了しない）
することになった。
(1) Aくん、Bくんについて、コイントスの回数の期待値を求めよ。
(2) Aくんはの方が少ない回数で終了する確率を求めよ。小数解でよい。
514(1): １３２人目の素数さん [sage] 2025/01/08(水) 19:29:44.12 ID:X/URcmTf(2/3) AAS
>>513
（２）回数の差の分布
https://i.imgur.com/jDZubdl.png
515: １３２人目の素数さん [sage] 2025/01/08(水) 20:16:25.91 ID:X/URcmTf(3/3) AAS
>>514
https://i.imgur.com/cmdTCLY.png
516(1): １３２人目の素数さん [sage] 2025/01/08(水) 20:27:25.85 ID:bjQ79ehj(1/2) AAS
日本語すら読めないチンパンジーが数学やった気になってるのか
517(1): １３２人目の素数さん [sage] 2025/01/08(水) 20:27:57.47 ID:SKdwPI7I(2/4) AAS
>>513
1000万回のシミュレーション
https://i.imgur.com/JKpHsKb.png
https://i.imgur.com/A1NiFF0.png
518(1): １３２人目の素数さん [sage] 2025/01/08(水) 20:28:49.92 ID:SKdwPI7I(3/4) AAS
>>516
罵倒しか投稿できないのがFimoseくん。FランにちなんでPhimoseくんから昇格ｗ
519: １３２人目の素数さん [sage] 2025/01/08(水) 20:30:20.34 ID:SKdwPI7I(4/4) AAS
>>517
期待値の差よりも
回数の差の分布の方が直感に近いな。
520: １３２人目の素数さん [sage] 2025/01/08(水) 20:33:00.33 ID:bjQ79ehj(2/2) AAS
チンパンジーが図星で発狂中w
521: １３２人目の素数さん [sage] 2025/01/09(木) 00:40:37.23 ID:bJYGsz4u(1) AAS
>>518
他人がどうこうなんてお前が無能な事とは全く関係ないんだから、まずは国語の勉強してスレタイくらい読めるようになれよwwww
汚い言葉勉強してる場合じゃないぞwww
522(1): １３２人目の素数さん [sage] 2025/01/09(木) 06:57:53.73 ID:0PPPIAFh(1/2) AAS
やはり、下手な考え休むに似たりだな。
オッペンハイマーでの会話
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.
523(1): １３２人目の素数さん [sage] 2025/01/09(木) 07:00:06.18 ID:0PPPIAFh(2/2) AAS
Fimoseくんはロータリーエンジンの動くシェーマの動画を作れないの？
完成した動画をみてもハウジングが何かわからんようだし。
やっぱり、弄れるのはforeskinだけみたいだね。
524(1): １３２人目の素数さん [sage] 2025/01/09(木) 08:14:29.03 ID:jj9pV73T(1) AAS
>>522
アンタのことじゃん
525: １３２人目の素数さん [sage] 2025/01/09(木) 08:14:47.34 ID:o9+qLk7l(1/2) AAS
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
526: １３２人目の素数さん [sage] 2025/01/09(木) 08:16:10.89 ID:o9+qLk7l(2/2) AAS
nを10000未満の正整数とする。
1/n+1/x=1/y を満たす{x,y}の正整数解の組み合わせの数が最大なnはいくつか？
527: １３２人目の素数さん [sage] 2025/01/09(木) 12:34:42.59 ID:1fujxzuR(1/2) AAS
こういう報道に接するとZ3が必須だな。
https://www.373news.com/_news/storyid/155006/
528(1): １３２人目の素数さん [sage] 2025/01/09(木) 12:35:57.90 ID:1fujxzuR(2/2) AAS
>>524
映画みてないの？
英語がわからないFラン？
529: １３２人目の素数さん [sage] 2025/01/09(木) 17:28:07.32 ID:6xzpzcRA(1) AAS
>>528
下手な考え休むに似たりって完全にアンタのことなんだけどまさか自覚ないんか？ww
530: １３２人目の素数さん [sage] 2025/01/10(金) 02:04:06.88 ID:sfHDQvoL(1) AAS
>>523
他人がどうこうとか自分の無能を誤魔化すばっかりwww
自分で定義分かってるとか思い込んでるから何にも理解出来ないんだよwww
思い込みじゃなくて自分できちんと日本語使って数学的な定義述べてみろよwww
出来ないんだろwww
531(1): １３２人目の素数さん [sage] 2025/01/10(金) 02:06:10.84 ID:5hjYx106(1/2) AAS
Fimoseくんはロータリーエンジンの動くシェーマの動画を作れないの？
完成した動画をみてもハウジングが何かわからんようだし。
やっぱり、弄れるのはforeskinだけ。
532: １３２人目の素数さん [sage] 2025/01/10(金) 02:14:15.07 ID:H6m+2x8m(1/2) AAS
結局できなくて発狂かよ
自分で出したクソの処理すらできない模様
533: １３２人目の素数さん [sage] 2025/01/10(金) 02:18:01.93 ID:MdFFAMIP(1) AAS
>>531
汚い言葉で誤魔化すばっかりかよwww
他人がどうこうなんて関係なくお前が無能な事を指摘して貰ってるんだから誤魔化さずに自分で努力してみろよwww
534(1): １３２人目の素数さん [sage] 2025/01/10(金) 02:31:21.32 ID:H6m+2x8m(2/2) AAS
毎回毎回同じことしか言えないし知性のかけらもないのに自称医者東大卒とは笑わせるねww
535(2): １３２人目の素数さん [sage] 2025/01/10(金) 05:36:36.58 ID:5hjYx106(2/2) AAS
>>534
俺は医科歯科卒
理１を蹴って入学。
あんたは動画すら作れないFランなんだろ。
このスレは出禁だぞ。
536: １３２人目の素数さん [sage] 2025/01/10(金) 07:02:22.75 ID:nzxEBzqG(1) AAS
>>535
医科歯科卒wwww
日本語すらまともに使えない中卒以下の低脳が何言ってんの？
Fラン受験出来るのが羨ましいのかよwww
537: １３２人目の素数さん [sage] 2025/01/10(金) 08:29:22.85 ID://fsJZeZ(1/2) AAS
医科歯科ってスレタイ読めなくても卒業できるの？
ルールに則らないことを推奨してるの？
医師国家試験でルールに則らないことを是とする回答をしたの？
538: １３２人目の素数さん [sage] 2025/01/10(金) 08:30:21.85 ID://fsJZeZ(2/2) AAS
質問スレで出題ばかりするし質問には回答しないしで、
世の中のいなくなったほうが世界か良くなる人物の中の一人
539: １３２人目の素数さん [sage] 2025/01/10(金) 09:16:14.22 ID:zQblG/AR(1/2) AAS
>>535
まあ匿名掲示板ならいくらでも経歴盛れるもんなそんな主張意味ないよ
本当に医科歯科卒ならそれ相応の能力を示せよ
540(1): １３２人目の素数さん [sage] 2025/01/10(金) 11:01:30.43 ID:QR184n+C(1/2) AAS
今日は内視鏡が少なくて10:30に終了。予約の検査終了したら帰っていい契約なので( ・∀・)ｲｲ!!
541: １３２人目の素数さん [sage] 2025/01/10(金) 11:16:40.37 ID:QR184n+C(2/2) AAS
熟練職人不足で新年から内視鏡週3回になった。
スタッフが可愛いので続けられる。
542: １３２人目の素数さん [sage] 2025/01/10(金) 11:31:20.92 ID:zQblG/AR(2/2) AAS
本物の医者がこんな数学板の糞スレにいるわけないじゃん忙しいはずなのに妄想乙
543: １３２人目の素数さん [sage] 2025/01/10(金) 11:56:51.07 ID:02D377cs(1) AAS
>>540
これが何かの証拠になるとか思ってんの？
説得力皆無で今まで誰一人信じて貰ってないのが現実だろwwww
何の意味もない事繰り返し書き込んじゃうって知能低すぎん？
544: １３２人目の素数さん [sage] 2025/01/10(金) 18:50:35.59 ID:nUmCuFV/(1) AAS
尿瓶ジジイ急にダンマリかよ
545(1): １３２人目の素数さん [sage] 2025/01/11(土) 07:06:37.50 ID:5YoQBcPi(1/2) AAS
A君とB君がそれぞれ同じ規格立方体のサイコロを投げます。
終了条件は
A君は、 同じ数の目が2回連続して出たら終了
B君は、 前の目の数より１多い数がでたら修了（但し、６→１でも終了するとする）
【問題】
(1) A君の投げた回数の期待値を求めよ
(2) B君の投げた回数の期待値を求めよ
(3) A君の投げた回数の方が多い確率をもとめよ
(4) 二人の投げた回数が同数で終了する確率をもとめよ
546: １３２人目の素数さん [sage] 2025/01/11(土) 08:13:55.14 ID:5YoQBcPi(2/2) AAS
(*
コイントスの終了条件を
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}]
547: １３２人目の素数さん [sage] 2025/01/11(土) 09:15:51.78 ID:EiIJz6LX(1/2) AAS
>>545
結局誰にも信じられずに哀れだね
548(1): １３２人目の素数さん [sage] 2025/01/11(土) 19:01:33.32 ID:WI+T/Kwt(1/3) AAS
(*
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
上下前次1-新書関写板覧索設栞歴
あと 443 ﾚｽあります
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.030s