数论及组合数学计科0907班.pptVIP

特殊的计数序列-Catalan数 此时，将前k项反号，即1--1,-1-1。那么整个序列中，有n+1个1，n-1个-1，一共有C(2n,n-1)种。这种排列和不满足条件的n个1和-1的排列是一一对应的。 那么满足条件的排列数：C(2n,n)-C(2n,n-1)= C(2n,n)/(n+1) Catalan数的一个应用 凸多边形的划分：将一个凸多边形通过其对角线的分割使其成为若干个三角形，问有多少种分法。 solution 确定一条基边，由该基边可引出n-1个三角形，且对于每个三角形，都将原多边形分成2部分。设h[i]为(i-1)边形划分总数。则有 h[n]=(h[1]*h[n-1])+(h[2]*h[n-2])+….+(h[n-1]*h[1]) 不加证明的给出(证明较复杂，大家注意应用就好）： h[n]=C(2n-2,n-1)/n.即h[n+1]=C(2n,n)/(n+1); 第二类sitirling数 问题的引入： 将P个不同的元素的集合划分成K个无差别的非空盒的划分个数 用s[p,k]表示划分个数 数论与组合数学 计科0907班 何剑 Email:jianhe25@ 主讲内容 数论 解模线性方程 解模线性方程组（中国剩余定理） 欧拉函数，欧拉定理 组合数学 鸽巢原理 全排列的生产 catalan数 polya计数 数论入门 解模线性方程： 求解方程 ax = b (%n) 其中a,b,n已知，求x 解法： 欧几里德方法求gcd： gcd(a,b) = gcd(b,a%b); //循环计算式,直到b==0. 扩展欧几里德： ax = b(%n) == ax + ny = b 的整数解 设 d=gcd(a,n)，则方程等价为 ax + ny = b; (a = a/d ; n = n/d ; b = b/d;需要d | b,不整除无解 ） 如何解ax+ny = b? ax+ny=b 的解是 ax+ny=1 的解的b倍 观察欧几里德循环过程： ax + ny = 1 与 gcd(a,n)=gcd(n,a % n); 方程 nx+(a % n)y=1 的解与ax+ny=1的解 是否可以通过某种变换使得x=f(x),y=g(y); 答案是显然的： y = x; x = x - [a/n]*y; 注意到a % n = a - [a/n]*n，自己演算一下就知道了 Ex_gcd(a,b,n) int ex_gcd(int a,int b,int x,int y) { int t,rec; if ( b==0 ) { x=1; y=0; return a; } rec = ex_gcd(b,a%b,x,y); t=x; x=y; y=t-a/b*y; return rec; } 模线性方程组 中国剩余定理： x = b[0] (mod m[0]) x = b[1] (mod m[1]) ... x = b[n-1] (mod m[n-1]) 当m[0],m[1],...m[n-1]两两互质时，方程有唯一解！解的范围在0..T-1;T=m[0]*..*m[n-1] 如何解满足中国剩余定理的方程： 记Mi=M/mi(0=in),因为gcd(Mi,mi)=1,故有二整数pi,qi满足Mi*pi+mi*qi=1.如果记ei=Mi*pi. 那么有： ei = 0(mod mj, j!=i) 或 1(mod mj, j=i) 于是解就是: 素数相关问题 欧拉函数phi(x): phi(x)={y|yx gcd(x,y)==1} 意义：比x小且和x互质的数。 群 Zx ; * *定义为 a*b % (x). 也就是一个(mod x) 下的包含乘法的群。 欧拉定理:a^phi(x) = 1 (%x) a属于整数集Z 由群论中的拉格朗日定理直接可得！ 拉格朗日定理： 有限群GS,*中任意一个元素的阶o(a) | #G. o(a)定义为a^o(a) % #G = a; 既a的周期 phi(x)怎么求： phi(x) = x * (1-1/fac[0]) * (1-1/fac[i]) *...(1-1/fac[n-1]) 简单来说就是每个x的质因子都删掉一部分数！ 组合数学 N个人坐在跟自己编号不同的位置上，有多少种方法？ 12个不同颜色的珠子串成一条项链，能够做成多少个不同的项链？ 什么是组合数学？ 第一个问题是经典的错排问题，考虑第N个人 若前N-1个人已完全错排则有第N个人和前N-1个任意一个人换即可. 若前N-1个人未完全错排，其中有一个仍在自己位置上，则第N个人跟它换即可. 若有两个人未错排，则无论如何都无法完全错排所以可得 F[N]=(N-1)*