卡特兰数

卡特兰数

转:http://www.mathoe.com/dispbbs.asp?boardid=89&ID=34522

卡特兰式:

Cn=1/(n+1)*C(2n,n)//C(2n,n)为组合数,图片显示不出来...囧

其递推关系试有:

Cn+1=(4n+2)/(n+2)*Cn或

Cn=(4*n-2)/n+1*Cn-1

另外:

令C(1)＝1,C(0)=1

catalan数满足递归式：

h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)h(0) (其中n>=2)

另类递归式：

h(n)=((4*n-2)/(n+1))*h(n-1);

该递推关系的解为：

h(n)=C(2n,n)/(n+1)= 1/(n+1)*c(2n,n)。 (n=1,2,3,...)

卡特兰数的应用

a. 括号化问题

矩阵链乘： P=a1×a2×a3×……×an，依据乘法结合律，不改变其顺序，只用括号表示成对的乘积，试问有几种括号化的方案？(h(n)种)

b. 出栈次序问题

　 一个栈(无穷大)的进栈序列为1，2，3，…，n，有多少个不同的出栈序列?

　　分析：对于每一个数来说，必须进栈一次、出栈一次。我们把进栈设为状态‘1’，出栈设为状态‘0’。n个数的所有状态对应n个1和n个0组成的2n位二进制数。由于等待入栈的操作数按照1‥n的顺序排列、入栈的操作数b大于等于出栈的操作数a(a≤b)，因此输出序列的总数目=由左而右扫描由n个1和n个0组成的2n位二进制数，1的累计数不小于0的累计数的方案种数。

　　在2n位二进制数中填入n个1的方案数为c(2n,n),不填1的其余n位自动填0。从中减去不符合要求（由左而右扫描，0的累计数大于1的累计数）的方案数即为所求。

　　不符合要求的数的特征是由左而右扫描时，必然在某一奇数位2m+1位上首先出现m+1个0的累计数和m个1的累计数，此后的2(n-m)-1位上有n-m个 1和n-m-1个0。如若把后面这2(n-m)-1位上的0和1互换，使之成为n-m个0和n-m-1个1，结果得1个由n+1个0和n-1个1组成的2n位数，即一个不合要求的数对应于一个由n+1个0和n-1个1组成的排列。

　　反过来，任何一个由n+1个0和n-1个1组成的2n位二进制数，由于0的个数多2个，2n为偶数，故必在某一个奇数位上出现0的累计数超过1的累计数。同样在后面部分0和1互换，使之成为由n个0和n个1组成的2n位数，即n+1个0和n-1个1组成的2n位数必对应一个不符合要求的数。

　　因而不合要求的2n位数与n＋1个0，n－1个1组成的排列一一对应。

　　显然，不符合要求的方案数为c(2n,n+1)。由此得出 输出序列的总数目=c(2n,n)-c(2n,n+1)=1/(n+1)*c(2n,n)。（这个公式的下标是从h(0)=1开始的）

类似：有2n个人排成一行进入剧场。入场费5元。其中只有n个人有一张5元钞票，另外n人只有10元钞票，剧院无其它钞票，问有多少中方法使得只要有10元的人买票，售票处就有5元的钞票找零？(将持5元者到达视作将5元入栈，持10元者到达视作使栈中某5元出栈)

c.凸多边形的三角剖分问题

　　求将一个凸多边形区域分成三角形区域的方法数。

　　类似：一位大城市的律师在她住所以北n个街区和以东n个街区处工作。每天她走2n个街区去上班。如果她从不穿越（但可以碰到）从家到办公室的对角线，那么有多少条可能的道路？

　　类似：在圆上选择2n个点,将这些点成对连接起来使得所得到的n条线段不相交的方法数?

前几个卡特兰数：

规定C0＝1，而

C₁＝1，C₂＝2，C₃＝5，C₄＝14，C₅＝42，

C₆＝132，C₇＝429，C₈＝1430，C₉＝4862，C₁₀＝16796，

C₁₁＝58786，C₁₂＝208012，C₁₃＝742900，C₁₄＝2674440，C₁₅＝9694845

相关题目:

http://acm.fzu.edu.cn/problem.php?pid=1581

http://acm.fzu.edu.cn/problem.php?pid=1902

1902_java代码:

import java.util.*; import java.math.BigInteger; public class Main { static BigInteger a,b,c,d,res,t1,t2,e=BigInteger.valueOf(1); public static void main(String[] args) { Scanner in = new Scanner(System.in); BigInteger ans,m=BigInteger.valueOf(20100501); int t,n,i,tmp,id=1; t=in.nextInt(); while(t>0) { t--; n=in.nextInt(); for(i=0;i<n;i++) tmp=in.nextInt(); ans=BigInteger.valueOf(0); for(i=1;i<=n/2;i++) ans=ans.add(catalan(n,n-i)).mod(m); System.out.println("Case #" + id + ": " + ans); id++; } } public static BigInteger catalan(int n,int m)//计算卡特兰数 { a=b=c=d=e; int j; for(j=1;j<=m;j++) { a=a.multiply(BigInteger.valueOf(n+1-j)); b=b.multiply(BigInteger.valueOf(j)); } t1=a.divide(b); for(j=1;j<=m+1;j++) { c=c.multiply(BigInteger.valueOf(n+1-j)); d=d.multiply(BigInteger.valueOf(j)); } t2=c.divide(d); res=t1.subtract(t2); return res; } }