R语言生成和为1、步长0.05的两位小数随机概率方法

你最初的代码逻辑有两个硬伤：

1. 直接对随机抽样结果做除法归一化，会产生浮点计算误差，必然出现非0.05步长的多位小数
2. 没有提前约束三个数的和匹配目标值，靠事后归一化永远没法严格满足步长+总和的双重要求

把所有数值放大20倍转为整数运算，从根源规避小数精度问题：

• 步长0.05放大20倍为整数1
• 取值范围0.05_{0.95放大20倍为整数区间1}19
• 总和要求1放大20倍为整数20
  此时问题简化为：找3个落在1~19区间的整数，三者和为20，最后将结果统一除以20，就能得到完全符合要求的两位小数概率值。

# 枚举所有满足整数约束的合法组合
valid_combinations <- list()
for (n1 in 1:19) {
  for (n2 in 1:19) {
    n3 <- 20 - n1 - n2
    if (n3 >= 1 && n3 <= 19) {
      valid_combinations[[length(valid_combinations) + 1]] <- c(n1, n2, n3)
    }
  }
}

# 随机抽取一组合法值，还原为概率
selected <- sample(valid_combinations, size = 1)[[1]]
x <- selected / 20

# 结果校验（可省略）
print(x)
print(paste0("三个数求和结果：", sum(x)))

• 输出的三个值固定为0.05的整数倍，天然是两位小数，不会出现无限小数问题
• 每个值严格落在0.05~0.95的区间内
• 三者求和严格等于1，无浮点误差
• 如果需要三个概率值互不相等，只需要在组合判断的条件里追加&& n1 != n2 && n2 != n3 && n1 != n3即可。

内容的提问来源于stack exchange，提问作者Joker312