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

高校数学の質問スレ（医者・東大卒専用） Part438 (991ﾚｽ)
上下前次1-新
抽出解除ﾚｽ栞

65: １３２人目の素数さん [sage] 2024/08/16(金)17:29:04.80 ID:NAPgnpHZ(1/2)
問題　油分け算のソルバーを作れ。（言語は問わない）

R言語でのソルバー
https://i.imgur.com/bqkmHbB.png

これをWolframに移植する気にならんなぁ。
Wolframでないと分数での厳密解が出せないというわけでもないから。

黒瀬のスパイスをつかってチキンジャーキーを作る方が楽しい。
焼き肉のタレでつくったら嫁から評判が悪かったので作り直し。

301: １３２人目の素数さん [sage] 2024/11/04(月)07:56:10.80 ID:pGz/mTDl(1)
>>299
f(x)=0とa≦x≦<b
から存在する(∃;Exists)xは示せるが、
f(x)=0とa<x<b
から任意の(∀;for All)xは示せないから

338: １３２人目の素数さん [sage] 2024/11/23(土)16:01:29.80 ID:qyJPztJt(3/4)
年賀状の３等は１００枚に３枚が当たる確率である。
(1) １００枚の年賀状を無作為に購入したときの当たりの枚数の期待値を求めよ。
(2) ３等が当たっている年賀状の枚数が奇数の確率をＰとする。Pは0.5より大きいか小さいか、直感で即答せよ。
(3) (2)の確率を分数で求めよ。
(4) (3)をシミュレーションで検証してみよ。

p=3/100;
Table[Binomial[100,2n-1] p^(2n-1) (1-p)^(100-(2n-1)),{n,1,50}] // Total
% // N
Table[Boole@OddQ@RandomVariate[BinomialDistribution[100,3/100]],10^6] // Mean

575(1): １３２人目の素数さん [sage] 01/16(木)04:27:48.80 ID:jkCW5V3d(1/2)
やっぱりFランなんだろうな。
中央値の区間推定すらできないとは。

688: １３２人目の素数さん [sage] 02/06(木)17:18:59.80 ID:Kqo8KOCs(6/6)
しまった、隔離スレのつもりでレスしてたけど違った
皆様申し訳ない
もうスレ違いの書き込みはしません

783: １３２人目の素数さん [sage] 04/10(木)08:52:41.80 ID:LedCjmFj(1)
やっぱりこのポンコツには反応したらあかんな
ゴミみたいな問題をいつまでもいつまでもダラダラダラダラ

800(1): １３２人目の素数さん [sage] 04/29(火)18:09:40.80 ID:pY4WJf3b(3/4)
library(RcppAlgos)
library(fmsb)
library(matrixStats)

N <- 50

# Generate all combinations
cm <- comboGeneral(0:N, 3, repetition = FALSE)

# Pre-allocate
n_vec <- rep(N, 3)
success <- cm
failure <- matrix(n_vec, nrow = nrow(cm), ncol = 3, byrow = TRUE) - cm

# Modified fast Fisher function - correct implementation
fast_fisher <- function(success, failure) {
# Initialize p-value matrix
pvals <- matrix(NA_real_, nrow = nrow(success), ncol = 3)

# Perform pairwise comparisons
for (i in 1:nrow(success)) {
# 1 vs 2
pvals[i,1] <- fisher.test(matrix(c(success[i,1], failure[i,1],
success[i,2], failure[i,2]), nrow = 2))$p.value

# 1 vs 3
pvals[i,2] <- fisher.test(matrix(c(success[i,1], failure[i,1],
success[i,3], failure[i,3]), nrow = 2))$p.value

# 2 vs 3
pvals[i,3] <- fisher.test(matrix(c(success[i,2], failure[i,2],
success[i,3], failure[i,3]), nrow = 2))$p.value
}

# Bonferroni adjustment
pmin(pvals * 3, 1) # Cap at 1 after adjustment
}

# Run with timing
system.time({
# Overall Fisher tests
overall_p <- apply(cbind(success, failure), 1, function(x) {
fisher.test(matrix(x, nrow = 2))$p.value
})

# Pairwise Fisher tests
pairwise_p <- fast_fisher(success, failure)
min_pairwise_p <- rowMins(pairwise_p, na.rm = TRUE)

# Filter condition
keep <- overall_p > 0.05 & min_pairwise_p < 0.05
result <- cm[keep, ]
})

# Print first few results
head(result)
nrow(result) # Number of qualifying combinations

865: １３２人目の素数さん [sage] 06/04(水)07:49:17.80 ID:n75lIIio(1)
set.seed(123)
solve = function(x,k=1e5){
f = function(x) runif(1,x-0.5,x+0.5)
y=replicate(k,mean(sapply(x,f)))
quantile(y,p=c(0.025,0.975))
}
solve(c(9,10,11,11,12))

set.seed(123) # 再現性のため
solve2 <- function(x, k = 1e5) {
# 各x_iに対して一様乱数を生成し、平均を計算（ブートストラップ）
bootstrap_means <- replicate(k, {
x_true <- runif(length(x), x - 0.5, x + 0.5)
mean(x_true)
})
# 95%信頼区間を計算
quantile(bootstrap_means, probs = c(0.025, 0.975))
}

# 実行例
x_rounded <- c(9, 10, 11, 11, 12)
solve2(x_rounded)

924: １３２人目の素数さん [] 06/22(日)13:28:47.80 ID:DRKE+jFQ(1/5)
医者東大卒専用が聞いて呆れる
日本語読めないチンパンが発狂してるだけのスレ

934(1): １３２人目の素数さん [sage] 06/22(日)17:04:39.80 ID:AY7cZjkg(9/10)
傑作質問が頭の中に次々と浮かんできます
明日から1日1問、質問してもいいですか？

上下前次1-新書関写板覧索設栞歴

ｽﾚ情報赤ﾚｽ抽出画像ﾚｽ抽出歴の未読ｽﾚ AAｻﾑﾈｲﾙ

ぬこの手ぬこTOP 0.039s