卡特兰数

最新推荐文章于 2023-12-24 19:11:10 发布

L_heart

最新推荐文章于 2023-12-24 19:11:10 发布

阅读量256

点赞数

分类专栏：数学文章标签： ACM OJ题目

本文链接：https://blog.csdn.net/qq1195638709/article/details/76269305

版权

数学专栏收录该内容

2 篇文章 0 订阅

订阅专栏

卡特兰数(Catalan number)

【简介】

卡特兰数又称卡塔兰数，英文名Catalan number，是组合数学中一个常出现在各种计数问题中出现的数列。以比利时的数学家欧仁·查理·卡塔兰 (1814–1894)的名字来命名，其前几项为 : 1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900, 2674440, 9694845, 35357670, 129644790, 477638700, 1767263190, 6564120420, 24466267020, 91482563640, 343059613650, 1289904147324, 4861946401452, ...

【原理】

令h(0)=1,h(1)=1，catalan数满足递推式： h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)*h(0) (n>=2)

例如：h(2)=h(0)*h(1)+h(1)*h(0)=1*1+1*1=2

另类递推式：h(n)=h(n-1)*(4*n-2)/(n+1)(C为组合数)

递推关系的另类解为：h(n)=c(2n,n)-c(2n,n-1)(n=0,1,2,...)

【应用】

实质上为递推式的应用

1. 由n个+1和n个-1构成2n项 a_1,a_2,...,a_n 其部分和满足 a_1 + a_2 + ... + a_k /ge 0 /quad , /quad 0 /le k /le 2n 的序列个数等于第n个Catalan数 C_n 。

假设不满足条件的序列个数为 U_n ，那么就有 $C_n + U_n = {2n /choose n}$ 。而对于不满足的序列，必然存在某一个奇数位使 a_1 + a_2 + ... + a_k < 0 ，而且+1的个数恰好比-1的个数少一个，此时我们将前k项中的+1变为-1，将-1变为+1,那么就得到一个有(n+1)个+1和(n-1)个-1的序列，即一个不满足条件的序列对应于一个由n+1个+1和n-1个-1组成的排列。所以 $U_n = {2n/choose n + 1}$ 。那么 $C_n = {2n/choose n} - U_n = {2n/choose n} - {2n/choose n + 1}}}$ 。

2. 我们可以将应用1变换形式：将-1看成右括号，+1看成左括号，就变成了左括号和右括号各有n个时，合法括号表达式的个数。比如2个左括号和2个右括号组成的合法表达式有 C_2 = 2 种，是()()和(())。

3. 再次转换形式，对于n+1个数连乘，乘法顺序有 C_n 种，比如三个数连乘a*b*c，n=2，相当于在式子上加括号，有2种乘法顺序，分别是(ab)c和a(bc)。四个数连乘a*b*c*d，n=3，有5种乘法顺序，分别是a(b(cd))，(ab)(cd)，((ab)c)d，(a(bc))d和a((bc)d)。同时对于n为3得到合法的括号序列有5个，分别是：((()))，()(())，()()()，(())()和(()())，也即与n+1个数连乘的顺序数一一对应关系。

4. n个节点构造二叉树的所有可能形态数为 C_n 。考虑随便取一个节点作为根，那么他左边和右边的儿子节点个数就确定了，假定根节点标号为x，那么左子树的标号就是从1到x-1,共x-1个，右子树的标号就是从x+1到n，共n-x个，那么将x从1取到n，就获得了所有的情况数 $/begin{displaymath}C_n = /sum_{i = 0}^{n - 1}C_i/,C_{n - i - 1}/end{displaymath}$ 。

5. n个非叶节点的满二叉树的形态数（对称后得到的二叉树除非自己本身对称，否则算是不同）

对于满二叉树，实际上就是将应用4中的每个子节点的空儿子上都加上叶子，就形成了对应的满二叉树，那么n个非叶节点形成的满二叉树的形态数即为Cn

6. 对于一个n*n的正方形网格，每次只能向右或者向上移动一格，那么从左下角到右上角所有在副对角线右下方的路径总数为 C_n 。

可以将一条水平边记为+1,垂直边记为-1,那么就组成了一个n个+1和n个-1的序列，并且保证前k步中水平边数不小于垂直边数，换句话说前k个元素的和非负。

7. 对凸n+2边形进行不同的三角形分割（只连接顶点对形成n个三角形）数为Cn

8. n个数入栈后的出栈的排列总数是 C_n 。例如1,2,3入栈的出栈排序有123，132，213，231和321五种。

【常规分析】

首先，我们设f（n）=序列个数为n的出栈序列种数。（我们假定，最后出栈的元素为k，显然，k取不同值时的情况是相互独立的，也就是求出每种k最后出栈的情况数后可用加法原则，由于k最后出栈，因此，在k入栈之前，比k小的值均出栈，此处情况有f(k-1)种，而之后比k大的值入栈，且都在k之前出栈，因此有f(n-k)种方式，由于比k小和比k大的值入栈出栈情况是相互独立的，此处可用乘法原则，f(n-k)*f(k-1)种，求和便是Catalan递归式。

首次出空之前第一个出栈的序数k将1~n的序列分成两个序列，其中一个是1~k-1，序列个数为k-1，另外一个是k+1~n，序列个数是n-k。

此时，我们若把k视为确定一个序数，那么根据乘法原理，f（n）的问题就等价于——序列个数为k-1的出栈序列种数乘以序列个数为n - k的出栈序列种数，即选择k这个序数的f（n）=f（k-1）×f（n-k）。而k可以选1到n，所以再根据加法原理，将k取不同值的序列种数相加，得到的总序列种数为：f（n）=f（0）f（n-1）+f（1）f（n-2）+···+f（n-1）f（0）。f（n）=h（n）= C（2n,n）/（n+1）= c（2n,n）-c（2n,n-1）（n=0，1，2，……）

【非常规分析】

对于每一个数来说，必须进栈一次、出栈一次。我们把进栈设为状态‘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)=c(2n,n)/(n+1)=h(n)。

9. 对于集合 /{1,2,...,2n/} ，将集合元素两两分为n个子集，若任意两个子集都不交叉，那么我们称此划分为一个不交叉划分。此时不交叉的划分数就是 C_n 。考虑将每个子集中较小的数用左括号代替，较大的用右括号代替，那么带入原来的1至2n的序列中就形成了合法括号问题。例如集合{1,2,3,4,5,6}的不交叉划分有五个：{{1,2},{3,4},{5,6}}，{{1,2},{3,6},{4,5}}，{{1,4},{2,3},{5,6}}，{{1,6},{2,3},{4,5}}和{{1,6},{2,5},{3,4}}。

10. n层的阶梯切割为n个矩形的切法数也是 C_n 。如下图所示：

考虑先绘制如下图片，即n为5的时候的阶梯：

注意到每个切割出来的矩形都必需包括一块标示为*的小正方形，那么此时枚举每个*与#标示的两角作为矩形，剩下的两个小阶梯就是我们的两个更小的子问题了，于是我们的 C_5 = C_0 * C_4 + C_1 * C_3 + C_2 * C_2 + C_1 * C_3 + C_0 * C_4 。

11. 在一个2*n的格子中填入1到2n这些数值使得每个格子内的数值都比其右边和上边的所有数值都小的情况数也是 C_n 。

12. 平面上连接可以形成凸包的2n个点分成2个一组连成n条线段，两两线段之间不相交的情况总数是 C_n ，这里实际上和应用7本质上是一样的。

13. n+m个人排队买票，并且满足 n /ge m ，票价为50元，其中n个人各手持一张50元钞票，m个人各手持一张100元钞票，初始时候售票窗口没有钱，问有多少种排队的情况数能够让大家都买到票。

如果m=n的话那么这就是初始的Catalan数问题，也就是将手持50元的人看成是+1，手持100元的人看成是-1，任前k个数值的和都非负的序列数。

对于n>m的情况，假设大家都可以买到票的情况数是 $D_{n+m}$ ，无法让每个人都买到的情况数是 $U_{n + m}$ ，那么就有 $D_{n + m} + U_{n +m} = {n + m /choose n}$ ，假设最早买不到票的人编号是k，他手持的是100元并且售票处没有钱，那么将前k个人的钱从50元变成100元，从100元变成50元，这时候就有n+1个人手持50元，m-1个手持100元的，所以就得到 $U_{n + m} = {n + m /choose n + 1}$ ，同时 $D_{n + m} = {n + m /choose n} - {n + m /choose n + 1}$ 。

参考： http://blog.csdn.net/zhangmh93425/article/details/44677891

https://baike.baidu.com/item/%E5%8D%A1%E7%89%B9%E5%85%B0%E6%95%B0/6125746?fr=aladdin

【实现】

#define max 36
__int64 a[max];
void catalan()
{
	int i,j;
	memset(a,0,sizeof(a));
	a[0]=a[1]=1;
	for(i=2;i<max;i++)
		for(j=0;j<i;j++)
		{  
			a[i]+=a[j]*a[i-j-1];
		}  
}

小兔的棋盘

Problem Description

小兔的叔叔从外面旅游回来给她带来了一个礼物，小兔高兴地跑回自己的房间，拆开一看是一个棋盘，小兔有所失望。不过没过几天发现了棋盘的好玩之处。从起点(0，0)走到终点(n,n)的最短路径数是C(2n,n),现在小兔又想如果不穿越对角线(但可接触对角线上的格点)，这样的路径数有多少?小兔想了很长时间都没想出来，现在想请你帮助小兔解决这个问题，对于你来说应该不难吧!

Input

每次输入一个数n(1<=n<=35)，当n等于－1时结束输入。

Output

对于每个输入数据输出路径数，具体格式看Sample。

Sample Input

Sample Output

   
   
    
    1 1 2
2 3 10
3 12 416024

代码1

#include<iostream>
using namespace std;
#define max 71
__int64 c[71][71];
__int64 a[36];
void p()				//打杨辉三角组合数
{
	int i,j;
	c[0][0]=1;
	for(i=1;i<max;i++)
	{
		c[i][0]=c[i][i]=1;
		for(j=1;j<i;j++)
		{
			c[i][j]=c[i-1][j]+c[i-1][j-1];
		}
	}
}
void f()				//打卡特兰数
{
	int i;
	p();
	for(i=1;i<36;i++)
	{
		a[i-1]=c[2*i][i]-c[2*i][i-1];
	}
}
int main()
{
	int n,count=0;
	f();
	while(cin>>n&&n+1)
	{
		printf("%d %d %I64d\n",++count,n,2*a[n-1]);
	}
	return 0;
}

代码2

#include<iostream>
using namespace std;
#define max 36
__int64 a[max];
void catalan()
{
	int i,j;
	memset(a,0,sizeof(a));
	a[0]=a[1]=1;
	for(i=2;i<max;i++)
		for(j=0;j<i;j++)
		{  
			a[i]+=a[j]*a[i-j-1];
		}  
}
int main()
{
	int n;
	int count=0;
	catalan();
	while(cin>>n&&n+1)
	{
		printf("%d %d %I64d\n",++count,n,2*a[n]);
	}
	return 0;
}