【派遣s1】用R语言计算先锋存活的概率
本帖最后由 咸鱼鱼 于 2025-5-3 16:44 编辑先知：远征期间将有先知每日掷出生死骰子预测先锋的生死，累计结果。
若掷出数量 生(1)＞死(-1)，则先锋必将生还凯旋。
若掷出数量 生(1)＜死(-1)，则先锋必将英勇牺牲。
若掷出数量 生(1)= 死(-1)，则次日再掷1次定生死。
起止：先知预测有效期: 4月25日~5月7日
规则如上，可以简单翻译为
一个先知每天占卜先锋的生死，有三种情况，生、死、不明，概率相等。持续X日，可以暂时认为7日或者14日（毕竟7天后才选先锋，前面7天一直出死是不是大家就未战先怯了？）
当生的次数大于死的次数，则先锋生存。当死的次数大于生时，则先锋死亡。如果次数相同，则额外占卜若干次，直到生死次数不等。求先锋生存的概率
如果占卜仅有概率相等的两种情况（生、死），其他规则不变，求先锋生存的概率
此外在最开始的设计中，是d-1会出现-1、0、1三种情况的，我们也不知道最终结果是会保留三种情况，还是只保留生死，因此都必须计算
三种情况的代码
compute_prob_case1 <- function(X) {
# 生成所有可能的s和d组合
dat <- expand.grid(s = 0:X, d = 0:X)
dat <- dat
# 计算每个组合的概率
dat$prob <- apply(dat, 1, function(row) {
s <- row[['s']]
d <- row[['d']]
counts <- c(s, d, X - s - d)
dmultinom(counts, prob = rep(1/3, 3))
})
P_gt <- sum(dat$prob)
P_eq <- sum(dat$prob)
return(P_gt + 0.5 * P_eq)
}
X <- 7
result_case1 <- compute_prob_case1(X)
print(result_case1)# 输出：0.5
两种情况的代码
compute_prob_case2 <- function(X) {
# 由于对称性，生存概率始终为0.5
return(0.5)
}
X <- 7
result_case2 <- compute_prob_case2(X)
print(result_case2)# 输出：0.5
最终的结论
问题1分析：在每次占卜中，生、死、不明的概率均为1/3。当生死次数相等时，进入额外占卜阶段，此时每次占卜的生/死概率相等，故生存概率为0.5。由于生与死的初始概率对称，总体生存概率为0.5。
问题2分析：
每次占卜生或死的概率各为1/2。当生死次数相等时，额外占卜的生存概率为0.5。同样由于对称性，无论占卜天数X是奇数还是偶数，最终生存概率均为0.5。
两个问题的生存概率均为0.5，这可以通过对称性和概率的公平性解释，R代码验证了这一点。
匪夷所思的反直觉事件
为什么无论怎么样都是0.5呢，我觉得可以简单理解成抛硬币，你抛多少枚硬币，是不是无论每个面都是0.5的概率，那么在1次定胜负的概率是0.5的情况下，两次的情况是不是，正正，正反，反正，正正。还是0.5？。当你抛第三次的时候，在第二次的基础上
正正和反反，第三次硬币都不会改变第二次的结果。正反，反正则会基于第三次的结果发生改变，依旧是0.5
那么三种情况下，你就可以理解成不明的时候，是硬币立起来了，无需计算其中。还是只要看正反情况
所以，无论持续占卜多少天，先锋生存的概率都是0.5
死去的高数还在追我……
生存的概率都是0.5那就只能看运气了惹{:6_169:} 一半的存活概率咩(´×ω×`)感觉要遭，而且还没那个复活的勋章的……说不定压根不会被抽到欸～那没事咯，喝茶ing 一半的概率）好残酷的LV7级用户随从远征.jpg大佬们加油啊） 如果这个概率模型复杂一点就可以计算马尔可夫链了，不过看起来版主不打算折磨大家的大脑{:6_169:}
好像还可以通过转经筒复活来着，要是一死就能复活那应该能提升不少概率 毕竟最后生死是统计生死数量，而不是全部都要roll生才能活。只要这个概率不涉及条件和状态的改变，那在这个赛博空间里，这个概率是不会随着次数递增受到影响的。【真实世界扔骰子就不一样了【当然，赛博伪随机数什么的就不讨论了
话说比起apply，这种不复杂的function我更喜欢直接向量计算。感觉这样简单一点
compute_prob_case1 <- function(X) {
s_vec <- rep(0:X, each = X + 1)
d_vec <- rep(0:X, times = X + 1)
valid_idx <- s_vec + d_vec <= X
s <- s_vec
d <- d_vec
r <- X - s - d
probs <- dmultinom(x = cbind(s, d, r), prob = rep(1/3, 3))
P_gt <- sum(probs)
P_eq <- sum(probs)
return(P_gt + 0.5 * P_eq)
}
感觉还好，主要是损失在可接受范围之内，而存活奖励足够吸引人~ 只有一个先锋，感觉对于个人来说，选上的概率也不到1％ 其实不用算的，因为无论怎么算，只要生死状态是对称的，就一定是50% 很有感觉，先知一票一票的投入预知先锋生死的结果池，这个过程想想就很刺激了{:6_179:} 感觉主要问题是还不一定被选为先锋，到现在也还没具体公布到底是怎么选先锋的，怀疑我到s100都不会被选中呜呜 幸好被选上成先锋的概率也很低惹 哈哈哈，确定不是概率论吗，不过。。。。 虽然大数据是一半存活概率~
但到底还是得看运气~希望运气偏向好那边一些惹吧;P 生死的概率一样，所以不管多少次累计都是一样 反正一切都是機率，看看誰的運氣比較好了 还好结论一半一半，符合我的直觉，我还以为会出什么反直觉的结论 直接无脑选择，然后看命，嘻嘻 :funk:没仔细看这次活动喵，只知道要发帖omg，还要死随从吗oh no:'( 生或死呢, 不過有人把需要靈魂來升級的勛章放上去就精彩了:loveliness:
页:
[1]
2