卡特兰数、dp-p1044

题目背景

栈是计算机中经典的数据结构，简单的说，栈就是限制在一端进行插入删除操作的线性表。

栈有两种最重要的操作，即 pop（从栈顶弹出一个元素）和 push（将一个元素进栈）。

栈的重要性不言自明，任何一门数据结构的课程都会介绍栈。宁宁同学在复习栈的基本概念时，想到了一个书上没有讲过的问题，而他自己无法给出答案，所以需要你的帮忙。

题目描述

宁宁考虑的是这样一个问题：一个操作数序列，1,2,\ldots ,n1,2,…,n（图示为 1 到 3 的情况），栈 A 的深度大于 nn。

现在可以进行两种操作，

1. 将一个数，从操作数序列的头端移到栈的头端（对应数据结构栈的 push 操作）
2. 将一个数，从栈的头端移到输出序列的尾端（对应数据结构栈的 pop 操作）

使用这两种操作，由一个操作数序列就可以得到一系列的输出序列，下图所示为由 1 2 3 生成序列 2 3 1 的过程。

（原始状态如上图所示）

你的程序将对给定的 nn，计算并输出由操作数序列 1,2,\ldots,n1,2,…,n 经过操作可能得到的输出序列的总数。

输入格式

输入文件只含一个整数 nn（1 \leq n \leq 181≤n≤18）。

输出格式

输出文件只有一行，即可能输出序列的总数目。

输入输出样例

输入 

3

输出 

5

方法一、卡特兰数

出入栈种类问题：基本卡特兰数模型

h(0)=h(1)=1

直接利用卡特兰数的通项 h(n)=h(n-1)*(4*n-2)/(n+1)即h(n)=C(2n,n)/(n+1) (n=0,1,2,...)

或者h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)*h(0) (n>=2)

然后写递推代码就行

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <ctype.h>
#include <stdlib.h>
#define  max(x,y) (((x)>(y)) ? (x) : (y))
#define  min(x,y) (((x)<(y)) ? (x) : (y))
#define PI acos(-1)
#define ll long long

ll c[20], n;
int main(void) {
	scanf("%lld", &n);
	c[1] = 1;
	for (int i = 2; i <= n; i++)
		c[i] = c[i - 1] * (4 * i - 2) / (i + 1);
	printf("%lld", c[n]);
	return 0;
}

或者

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <ctype.h>
#include <stdlib.h>
#define  max(x,y) (((x)>(y)) ? (x) : (y))
#define  min(x,y) (((x)<(y)) ? (x) : (y))
#define PI acos(-1)
#define ll long long
int n;
int c[20];

int main(void) {
	c[0] = c[1] = 1;
	scanf("%d", &n);
	for (int i = 2; i <= n; i++)
		for (int j = 0; j < i; j++)
			c[i] += c[j] * c[i - j - 1];
	printf("%d", c[n]);
	return 0;
}

 方法二、递归+记忆化搜索

基本思路:用二维数组存储栈内和未入栈元素，搭配记忆化，减少递归次数

状态表示：dp[i][j]表示栈内有i个元素，还有j个元素未入栈，值为情况数（i表示出栈数也可，因为元素总数是确定的）。

结束条件：j=0时，返回1；

特殊情况：i要始终大于0才能有出栈情况，当i=0时只能入栈。

状态转移：i大于0： dp[i][j]=dp[i-1][j]+dp[i+1][j-1]

                  i等于0： dp[i][j]=dp[i+1][j-1];

代码如下

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <ctype.h>
#include <stdlib.h>
#define  max(x,y) (((x)>(y)) ? (x) : (y))
#define  min(x,y) (((x)<(y)) ? (x) : (y))
#define PI acos(-1)
#define ll long long
int n;

int m[20][20];//记忆化
ll dp(int i, int j) {
	if (j == 0)
		return 1;//结束条件
	if (i == 0)
		return m[i][j] ? m[i][j] : (m[i][j] = dp(1, j - 1));
	else
		return m[i][j] ? m[i][j] : (m[i][j] = dp(i - 1, j) + dp(i + 1, j - 1));
}

int main(void) {
	scanf("%d", &n);
	printf("%lld", dp(0, n));
	return 0;
}