一、概念
以比利时的数学家欧仁·查理·卡塔兰 (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, ...
(一)括号化
矩阵连乘: P=a1×a2×a3×……×an,依据乘法结合律,不改变其顺序,只用括号表示成对的乘积,试问有几种括号化的方案?(h(n)种)
(二)出栈次序
首先,我们设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递归式。
(四)给定结点组成二叉搜索树
给定N个节点,能构成多少种不同的二叉搜索树?(能构成h(N)个)
(这个公式的下标是从h(0)=1开始的)
(五)n对括号正确匹配的数目
考虑n对括号时的任意一种配对方案,最后一个右括号有唯一的与之匹配的左括号,于是有唯一的表示A(B),其中A和B也是合法的括号匹配序列
假设S(n)为n对括号的正确配对数目,那么有递推关系S(n)=S(0)S(n-1)+S(1)S(n-2) +...+S(n-1)S(0),显然S(n)是卡特兰数。
(感谢百度百科)
int a[110][110]; //大数卡特兰数
int b[110]; //卡特兰数的长度
void Catalan(){ //求卡特兰数
int i,j,len,carry,tmp;
a[1][0]=b[1]=1;
len=1;
for(i=2;i<=100;i++){
for(j=0;j<len;j++)
a[i][j]=a[i-1][j]*(4*i-2);
carry=0;
for(j=0;j<len;j++){
tmp=carry+a[i][j];
a[i][j]=tmp%10;
carry=tmp/10;
}
while(carry){
a[i][len++]=carry%10;
carry/=10;
}
//carry=0;
for(j=len-1;j>=0;j--){
tmp=carry*10+a[i][j];
a[i][j]=tmp/(i+1);
carry=tmp%(i+1);
}
while(!a[i][len-1]) //高位零处理
len--;
b[i]=len;
}
}
(二)、卡塔兰数与逆元求解
const int mod = 1000000007;
ll Catelan[N];
//用逆元求解
ll exgcd(ll a, ll b, ll &x, ll &y)
{
ll d, t;
if (b == 0){x = 1; y = 0; return a;}
d = exgcd(b, a%b, x, y);
t = x;
x = y;
y = t - a/b*y;
return d;
}
void pre()
{
int i;
ll x, y;
Catelan[0] = 1, Catelan[1] = 1;
for(i = 2; i < N-5; i++)
{
Catelan[i] = Catelan[i-1]*(4*i-2) % mod;
exgcd(i+1, mod, x, y);
Catelan[i] = (Catelan[i]*((x+mod)%mod)) % mod;
}
}