[组合计数.pptVIP

组合计数问题 何亮 roba269@ 基础知识 加法原理 乘法原理 排列 组合 组合数递推式 在很多复杂问题中可以预处理出 在要求答案mod某数时常用 可重排列组合 可重排列 比如求{a,a,b,b,c}的不同全排列数： 5! / (2! * 2!) 可重组合 给定n种物品，要求取出m个(可以重复取) 容斥原理 Inclusion-Exclusion Principle 容斥原理 容斥原理的应用1——欧拉函数 定义phi(p)为比p小的与p互素的数 比如, phi(3)=2, phi(6)=2 设n的素因子有p1, p2, p3, … pk 欧拉函数 包含p1, p2…的个数为n/p1, n/p2… 包含p1*p2, p2*p3…的个数为n/(p1*p2)… 欧拉函数 在实际代码中可以用类似素数筛法求出 for (i = 1 ; i MAXN ; i++) phi[i] = i; for (i = 2 ; i MAXN ; i++) if (phi[i] == i) { for (j = i ; j MAXN ; j += i) { phi[j] /= i; phi[j] *= i - 1; } } 错位排列 容斥原理应用2——错位排列 若一个1..n的排列{a1,a2,a3…,an}满足a1≠1, a2≠2,…，则称其为一个(全)错位排列 求所有错位排列的个数 错位排列 所有排列数n! 有大于1个正确的排列数n*(n-1)! 有大于2个正确的排列数C(n,2)*(n-2)! 故全错的排列数为 错位排列 递推公式 f(n) = (n-1)*(f(n-1)+f(n-2)) 解释： 设把元素a放到位置B，则有两种情况 若元素b放到位置A，则变为f(n-2) 若b放在A,B之外，则变为f(n-1) 元素a可错放的位置有n-1种，故得上式 Catalan数 闭形式 递推式 Catalan数 1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 有多种解释方式 不穿过对角线的路径数 完全加括号方案数 满二叉树数目 多边形三角剖分方案数 Catalan数 – 形式1 不穿过对角线的路径数 一条不合法的路径 对应 一条(0,0)到(n-1,n+1)的路径 Catalan数 – 形式1 全部路径有C(2n,n)条 (0,0)到(n-1,n+1)的路径有C(2n,n-1)条 故合法路径有C(2n,n) - C(2n,n-1)条 化简可得 Catalan数 – 形式1 类似问题有 合法括号序列方案数 排队找零钱问题 Catalan数 – 形式2 Catalan数 – 形式2 任选一边(不妨选边0-1)，设其所在的三角形的顶点为点k，则多边形被分成一个k边形和一个n+1-k边形 定义k=2也为一种合法情况 则有f(n)=sigma(f(k)*f(n+1-k)) (2=k=n-1) f(n)=Catalan(n-2) Catalan数 – 形式2 N个中间结点的满二叉树数为Catalan(n) N个元素完全加括号的方案数为Catalan(n-1) 例题 POJ 3597 Polygon Division 求把一个N边形划分成三角形和四边形的方案数。(3=N=5000) 例题 确定一条起始边，分情况考虑 连出三角形——同Catalan数 连出四边形 分使用1/2/3条对角线情况 2条对角线又有相邻与不相邻两种情况 对于使用3条对角线的情况，直接写的复杂度是O(n^3)，可以通过预处理降到O(n^2) Super Catalan 任意加括号的合法方案数 不必完全加括号 UVA 10312 当四个元素时，有11种方案 xxxx, (xx)xx, x(xx)x, xx(xx), (xxx)x, x(xxx), ((xx)x)x, (x(xx))x, (xx)(xx), x((xx)x), x(x(xx)) Super Catalan 容易发现 S(N) = sigma( S(P1)*S(P2)*S(P3)*...*S(Pn) ) for all P1+P2+..+Pn = N 利用生成函数可以得出(过程略): S(n) = (3(2n-3)S(n-1) – (n-3)S(n-2)) / n S(1) = S(2) = 1 详见Brualdi《组合数学》8th chapter 第二类Stirling数 n个可区分的元素划分成k个集合的方案数 如n=4, k=2，则有7种方案