卡特兰数

这是和栈有关的知识点，我们通过一道例题来解释说明这个神奇的理论。
题目背景
栈是计算机中经典的数据结构，简单的说，栈就是限制在一端进行插入删除操作的线性表。

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

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

题目描述

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

现在可以进行两种操作，

将一个数，从操作数序列的头端移到栈的头端（对应数据结构栈的pushpush操作）

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

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

（原始状态如上图所示）

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

输入格式
输入文件只含一个整数n(1≤n≤18)n(1≤n≤18)

输出格式
输出文件只有11行，即可能输出序列的总数目

输入输出样例
输入 #1 复制
3
输出 #1 复制
5

核心思想：
OK，题意我们都理解的话，我们就开始思考。
建立数组f。f[i]表示i个数的全部可能性。
f[0] = 1, f[1] = 1; //当然只有一个
设 x 为当前出栈序列的最后一个，则x有n种取值
由于x是最后一个出栈的，所以可以将已经出栈的数分成两部分
比x小
比x大
比x小的数有x-1个，所以这些数的全部出栈可能为f[x-1]
比x大的数有n-x个，所以这些数的全部出栈可能为f[n-x]
这两部分互相影响，所以一个x的取值能够得到的所有可能性为f[x-1] * f[n-x]
另外，由于x有n个取值，所以
ans = f[0]*f[n-1] + f[1]f[n-2] + … + f[n-1]f[0];
这，就是传说中的卡特兰数

AC代码：

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5+10;
int n , f[N];
int main()
{
	cin >> n;
	f[0] = 1;
	f[1] = 1;
	for(int i = 2 ; i <= n ; i++)
	{
		for(int j = 0 ; j < i ; j++)
		{
			f[i] = f[i]+f[j]*f[i-j-1];
		}
	}
	cout << f[n] << endl;
	return 0;
} 

注意：
上面的只在数据很小的时候才适用，那么当数据较大的时候，就存不下那么多的数，怎么办呢？我们就要模拟乘法和除法来进行计算。
并且，卡特兰数还有一种计算方式，就是
H【1】 = 1；H【n】 = H【n-1】（4n-2）/（n+1）{n>=2}；
并且卡特兰数一般的适用范围有以下几种情况：
1括号匹配的种数。
2在栈中的自然数出栈的种数。
3求多边形内三角形的个数。
4，n个数围城圆圈，找不相交线段的个数。
5给定n个数，求组成二叉树的种数……

下面是模拟乘除法的代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define jz 10000
int a[110][110];
void chengfa(int k,int num){
    int jw = 0;
    for(int i = 100;i >= 0;i--){
        jw += a[k-1][i] * num;
        a[k][i] = jw % jz;
        jw /= jz;
    }
}
void chufa(int k,int num){
    int div = 0;
    for(int i = 0;i <= 100;i++){
        div = div*jz + a[k][i];
        a[k][i] = div / num;
        div %= num;
    }
}
void init(){
    memset(a,0,sizeof(a));
    a[1][100] = 1;
    int num;
    for(int i = 2;i <= 100;i++){
        num = (4*i-2);
        chengfa(i,num);
        num = i+1;
        chufa(i,num);
    }
}
int main()
{
    init();
    int n;
    while(~scanf("%d",&n)){
        if(n==-1) break;
        int pos = -1;
        for(int i = 0;i <= 100;i++){
            if(a[n][i]){
                pos = i;
                break;
            }
        }
        printf("%d",a[n][pos]);
        for(int i = pos+1;i <= 100;i++){
            printf("%04d",a[n][i]);
        }
        puts("");
    }
    return 0;
}