【异界方块】R语言计算本次活动概率

先说结论
一等奖概率: 0.045208 （4.5%）
二等奖概率: 0.32307 （32.3%）
三等奖概率: 0.631722 （63.2%）
总和: 1


比起新年的活动，1-2等奖的中奖概率可以说翻倍了，我记得今年过年的中奖率从8个数字提升到15个数字，中奖率才35%
（但是我依旧什么都没有）
感兴趣的同学可以让AI跑一下发1贴和2贴的概率是多少，我就摸鱼去了

R语言线上执行网页https://rtool.cn/run/r400.html
相关代码如下

1. # 定义颜色及其概率
   
2. colors <- c("金色", "黑色", "白色", "绿色", "蓝色", "粉色",
   
3.             "橙色", "红色", "紫色", "青色")
   
4. probs  <- c(0.01, 0.04, 0.04, 0.13, 0.13, 0.13,
   
5.             0.13, 0.13, 0.13, 0.13)
   
6. 
   
7. # 生成所有有序组合（3 次抽取，每次独立）
   
8. results <- expand.grid(rep(list(colors), 3), stringsAsFactors = FALSE)
   
9. colnames(results) <- c("c1", "c2", "c3")
   
10. 
    
11. # 计算每个结果的概率（乘积）
    
12. results$prob <- apply(results, 1, function(x) {
    
13.   prod(probs[match(x, colors)])
    
14. })
    
15. 
    
16. # 判断每个结果的中奖等级（优先级：一 > 二 > 三）
    
17. award <- character(nrow(results))
    
18. 
    
19. for (i in 1:nrow(results)) {
    
20.   triple <- as.character(results[i, 1:3])
    
21.   freq <- table(triple)
    
22.   n_unique <- length(freq)
    
23.   
    
24.   # ----- 一等奖 -----
    
25.   # 条件1：含有金色（至少一次）
    
26.   # 条件2：三个颜色相同
    
27.   if ("金色" %in% triple || n_unique == 1) {
    
28.     award[i] <- "一等奖"
    
29.     next
    
30.   }
    
31.   
    
32.   # ----- 二等奖 -----
    
33.   # 条件1：两个相同 + 一个不同（且不含金色，因为金色已过滤）
    
34.   if (n_unique == 2 && max(freq) == 2) {
    
35.     award[i] <- "二等奖"
    
36.     next
    
37.   }
    
38.   # 条件2：黑、白、其他（三个不同颜色，且同时包含黑和白）
    
39.   if (n_unique == 3 && "黑色" %in% triple && "白色" %in% triple) {
    
40.     award[i] <- "二等奖"
    
41.     next
    
42.   }
    
43.   
    
44.   # ----- 三等奖 -----
    
45.   # 三个不同颜色（此时一定不含金色，也不含黑+白组合）
    
46.   if (n_unique == 3) {
    
47.     award[i] <- "三等奖"
    
48.     next
    
49.   }
    
50.   
    
51.   # 理论上不应有未覆盖的情况，但以防万一
    
52.   award[i] <- "无奖"
    
53. }
    
54. 
    
55. # 汇总概率
    
56. prob_first  <- sum(results$prob[award == "一等奖"])
    
57. prob_second <- sum(results$prob[award == "二等奖"])
    
58. prob_third  <- sum(results$prob[award == "三等奖"])
    
59. 
    
60. cat("一等奖概率:", prob_first, "\n")
    
61. cat("二等奖概率:", prob_second, "\n")
    
62. cat("三等奖概率:", prob_third, "\n")
    
63. cat("总和:", prob_first + prob_second + prob_third, "\n")

复制代码

中奖的概率个人感觉差不多吧，最近的都是30%+左右变化不会很大

哇，概率居然有提高吗，看来我太非了，什么感觉都没有，这下只能靠模拟爽爽了

三十多的概率已经很不错了，虽然大概率又是陪跑，但是霉运小精灵谁不想要呢（真的吗）

看样子这次的活动获奖概率还是很大的，而且三等奖也有个保底。

主要这次三等奖也有个保底呢，总体上来说就是中了血赚不中也是不亏的程度

有横向对比才知道新年福卡的中奖率原来这么低吗，而且这回还有保底，以后新年福卡估计也有了，这是好事呀.jpg

有保底就是很好的事情啦，总好过没有保底至于中奖概率什么的，根本没指望过qaq

本非酋体质并不抱太大希望，0.5概率的抽卡都能5连歪