概率基础 · 计数方法

组合数学:中奖概率是怎么算出来的

「一千三百九十八万分之一」这个数字,不是估出来的,也不是官方公布后才知道的—— 它是一条纯粹的算术恒等式,用一支笔就能在纸上算完。本文把这个过程完整走一遍。

发布 2026-05-06 更新 2026-08-27 约 3500 字 · 阅读约 14 分钟 难度:入门

一、计数的两条基本原理

整个组合数学,说到底建立在两条几乎不言自明的原理上。

加法原理。如果一件事有 m 种做法或 n 种做法(两类不重叠),总共 m + n 种。

乘法原理。如果一件事分两步完成,第一步 m 种选择,每种选择下第二步都有 n 种,总共 m × n 种。

下面所有公式都只是这两条的反复应用。

二、排列:顺序算数

n 个不同的球里依次取出 k 个排成一列,有多少种排法?

第一个位置有 n 种选择;选完后剩 n−1 个,第二个位置有 n−1 种;以此类推,第 k 个位置有 n−k+1 种。由乘法原理:

P(n, k) = n × (n−1) × (n−2) × ⋯ × (n−k+1) = n! / (n−k)! k 个因子,从 n 递减。

代入 n = 49, k = 6

P(49, 6) = 49 × 48 × 47 × 46 × 45 × 44 = 10,068,347,520 一百亿种有序结果。但摇号的结果并不看顺序。

三、组合:顺序不算数

摇号装置先摇出 3 再摇出 17,和先摇出 17 再摇出 3,是同一组号码。所以上面那一百亿把同一组号码重复数了很多次——具体重复了多少次?

6 个确定的号码,能排成多少种顺序?由排列公式取 n = k = 6

6! = 6 × 5 × 4 × 3 × 2 × 1 = 720

也就是说,每一组号码在那一百亿里都被数了整整 720 遍。除掉即可:

C(n, k) = P(n, k) / k! = n! / [ k! · (n−k)! ] 读作「n 选 k」,也写作 ⁿCₖ(nk),称为二项式系数。

「先多数、再除掉重复」

这是组合计数最核心的手法。直接数无序的东西很难,但数有序的东西容易;数完之后除以每个对象被重复计数的次数就行。后面超几何分布的推导用的还是这一招。

四、把 C(49,6) 算到底

不用计算器,用约分的方式手算。

C(49, 6) = (49 × 48 × 47 × 46 × 45 × 44) / (6 × 5 × 4 × 3 × 2 × 1)

分母 720 = 6 × 5 × 4 × 3 × 2。逐项约分:

C(49,6) 的逐步约分过程
步骤 操作 结果
148 ÷ 68
245 ÷ 59
344 ÷ 411
446 ÷ 223
547 ÷ 147
6剩下 49 ÷ 3 除不尽,改用 8 ÷ 3 也不行 → 把 3 留到最后

约完后分子是 49 × 8 × 47 × 23 × 9 × 11,分母剩 3。把 9 ÷ 3 = 3

C(49, 6) = 49 × 8 × 47 × 23 × 3 × 11

逐步相乘:

49 × 8 = 392
392 × 47 = 18,424
18,424 × 23 = 423,752
423,752 × 3 = 1,271,256
1,271,256 × 11 = 13,983,816

验算:10,068,347,520 ÷ 720 = 13,983,816。一致。

这就是「1398万分之一」的来源

在 49 个号码中随机抽 6 个,共有 13,983,816 种等可能的号码组合。买一注,命中全部 6 个的概率就是 1 / 13,983,816 ≈ 0.00000007153,约 七千一百五十三万分之五,即 7.15 × 10⁻⁸

注意这个数字里没有任何经验成分,也不依赖历史数据。它只依赖两个参数:总号码数 49 和抽取个数 6。换成别的 m 选 n 规则,把数代进去重算即可——本站的组合概率计算器就是干这个的。

4.1 对称性

从公式可以直接看出 C(n, k) = C(n, n−k):选出 6 个等价于剩下 43 个,两种描述一一对应。所以 C(49, 6) = C(49, 43)

五、帕斯卡恒等式与三角形

二项式系数满足一条递推关系:

C(n, k) = C(n−1, k−1) + C(n−1, k)

组合证明。盯住第 n 个球。所有 k 元子集分为两类,且不重不漏:

  • 含第 n 个球:剩下要从前 n−1 个里选 k−1 个,共 C(n−1, k−1) 种。
  • 不含第 n 个球:全部从前 n−1 个里选 k 个,共 C(n−1, k) 种。

由加法原理,两类之和就是总数。∎

把这条递推排成三角形,就是帕斯卡三角(中国称杨辉三角,见于 1261 年杨辉《详解九章算法》引贾宪之法,早于帕斯卡约四百年)。每个数等于它上方两个数之和:

帕斯卡三角前 8 行(第 n 行第 k 列即 C(n,k),从 0 计数)
nC(n, 0) … C(n, n)
01
111
2121
31331
414641
515101051
61615201561
7172135352171

另外,第 n 行之和恒为 2ⁿ——因为 n 个元素的所有子集共 2ⁿ 个,按大小分类求和即得。

六、超几何分布:各等奖概率

现在算更细的问题:买一注 6 个号码,恰好命中 k 个的概率是多少?

把 49 个号码分成两堆:你选中的 6 个(记为「命中候选」),和你没选的 43 个。开奖时装置从 49 个里无放回抽 6 个。要恰好命中 k 个,意味着这 6 个开奖号里,有 k 个来自你的 6 个,有 6−k 个来自另外 43 个。

用乘法原理数一数满足条件的开奖组合数:

  • 从你的 6 个里挑 k 个:C(6, k)
  • 从另外 43 个里挑 6−k 个:C(43, 6−k)
  • 两步独立选择,相乘

除以总数 C(49, 6),得到超几何分布:

P(X = k) = C(6, k) · C(43, 6−k) / C(49, 6) 一般形式:P(X = k) = C(K, k) · C(N−K, n−k) / C(N, n),其中 N 为总数,K 为「特殊品」个数,n 为抽取个数。

6.1 逐项计算

k = 3 为例:

C(6, 3) = 20   C(43, 3) = (43 × 42 × 41) / 6 = 74,046 / 6 ...

43 × 42 × 41 算清楚:43 × 42 = 18061806 × 41 = 74,046。除以 3! = 612,341。所以

P(X = 3) = 20 × 12,341 / 13,983,816 = 246,820 / 13,983,816 ≈ 0.017650 约 1.765%,即 56.7 注里约有 1 注能命中 3 个。

七、完整分布表与验算

下表列出全部 7 种可能结果。所有分子都是精确整数,没有四舍五入。

49 选 6 规则下,一注命中 k 个号码的精确概率分布
命中 k 个 C(6, k) C(43, 6−k) 组合数(乘积) 概率 约合「多少注中 1 注」
016,096,4546,096,4540.4359652.3
16962,5985,775,5880.4130192.4
215123,4101,851,1500.1323787.6
32012,341246,8200.01765056.7
41590313,5450.0009691,032
56432580.0000184554,201
61110.0000000715313,983,816
合计13,983,8161.000000

合计行就是最好的验算

七个组合数相加,必须精确等于 C(49,6) = 13,983,816,一个不多一个不少。这是范德蒙德卷积恒等式的一个特例:

Σk C(K, k) · C(N−K, n−k) = C(N, n)

组合意义:从 N 个里选 n 个,按「其中有几个来自那 K 个」分类求和,必然遍历所有情况且不重复。任何超几何计算,都该用这条来自查——本站计算器每次运行都会执行这项校验。

命中个数的概率分布

49 选 6 规则,超几何分布精确值。注意纵轴是线性刻度——k ≥ 4 的柱子小到几乎看不见,这正是问题所在。

图里 k = 0k = 1 两根柱子加起来接近 85%。也就是说,随便买一注,最可能发生的事情是「一个都没中或只中一个」。而 k ≥ 4 的三根柱子在图上根本画不出来——它们合计只占 0.0988%,大约千分之一。

八、期望命中个数

超几何分布的期望有一个简洁的闭式:

E[X] = n · K / N = 6 × 6 / 49 = 36/49 ≈ 0.7347 平均每注命中约 0.73 个号码。

为什么是这个式子?用示性变量。令 Iⱼ = 1 表示你的第 j 个号码被开出,否则为 0。由对称性,任一特定号码被开出的概率是 6/49,所以 E[Iⱼ] = 6/49。命中总数 X = I₁ + ⋯ + I₆,由期望的线性性

E[X] = Σ E[Iⱼ] = 6 × (6/49) = 36/49

线性性不需要独立

这些 Iⱼ 之间并不独立(无放回抽样,见《独立事件与条件概率》第八节),但期望的线性性对任意随机变量都成立,不需要独立性。

这一点在《期望值与返奖率》里会变成一个决定性的论证:正因为线性性不挑条件,任何下注组合方式——无论多复杂、无论怎样选号——的总期望都只是各注期望之和,谁也逃不掉。

方差的公式为(含有限总体修正因子):

Var(X) = n · (K/N) · (1 − K/N) · (N−n)/(N−1) = 6 × (6/49) × (43/49) × (43/48) ≈ 0.5773

最后那个因子 (N−n)/(N−1) = 43/48 ≈ 0.896 就是无放回带来的修正:方差比同参数的二项分布约 10%。这是因为无放回抽样中各次结果存在负相关,互相抵消了一部分波动。

九、这个数字到底有多大

13,983,816 这个数写出来只有八位,看着不吓人。换几个参照物:

1/13,983,816 的几个等价表述
表述方式具体数值
逐格搜索把 13,983,816 个格子排成方阵,边长约 3,740 格;在其中找一个特定格子
连续时间每秒看一注,不吃不睡不间断,需要约 162 天才能看完全部组合
抛硬币大致相当于连抛 23.7 次全部正面(2²³·⁷ ≈ 1.4 × 10⁷)
每周买一注期望等待时间约 26.9 万年(13,983,816 ÷ 52)
面积类比在一个标准足球场(约 7,140 m²)上任撒一粒芝麻,正好落在指定的 0.5 mm² 上

「买多注」能改善多少

m不同的号码,命中概率精确等于 m / 13,983,816——线性增长,没有任何加速。

买 100 注,概率从 7.15 × 10⁻⁸ 变成 7.15 × 10⁻⁶,仍然是十四万分之一。要让概率达到 50%,需要买约 699 万注。而每一注的期望回报都是负的(见《期望值与返奖率》),所以买得越多,期望损失越大——概率的线性增长和损失的线性增长是同一条直线的两种读法。

9.1 那些「缩水」「杀号」技巧为什么无效

有一类说法是:通过某种规则先排除一部分号码,把 C(49,6) 缩小到一个更小的集合,从而提高命中率。

这在数学上的问题很直接:被排除的组合,和被保留的组合,中奖概率完全一样,都是 1/13,983,816缩小候选集不会改变每注的概率,只改变你买了几注。如果最终仍然买 m 注,命中概率仍是 m/13,983,816;如果排除规则不巧把真正的中奖组合排除在外,概率反而变成 0。

唯一真实存在的效应是与他人分享奖金:如果很多人偏好同一类号码(生日集中在 1–31、连号、对称图案),那么选择冷门组合在中奖之后可以少分一些。但这不改变中奖概率本身,只影响条件期望。


延伸阅读与参考

  • Feller, W. An Introduction to Probability Theory and Its Applications, Vol. 1 (3rd ed.). 第 2 章组合分析;第 2.6 节超几何分布。
  • Graham, R. L., Knuth, D. E., & Patashnik, O. Concrete Mathematics (2nd ed.). 第 5 章二项式系数,含范德蒙德卷积的多种证明。
  • Stanley, R. P. Enumerative Combinatorics, Vol. 1 (2nd ed.). 计数理论的标准参考。
  • Ross, S. A First Course in Probability (10th ed.). 第 1 章组合分析;第 4.8.3 节超几何随机变量。
  • 杨辉《详解九章算法》(1261),其中记载贾宪「开方作法本源图」,即帕斯卡三角。
  • Vandermonde, A.-T. (1772). Mémoire sur des irrationnelles de différents ordres. 范德蒙德卷积的原始出处。