Catalan Number 卡特兰数

1. 其前几项为 :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, ...
2. 令h(0)=1,h(1)=1，catalan数满足递推式 [1] ：
   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(3)=h(0)*h(2)+h(1)*h(1)+h(2)*h(0)=1*2+1*1+2*1=5
   另类递推式 [2] ：
   h(n)=h(n-1)*(4*n-2)/(n+1);
   递推关系的解为：
   h(n)=C(2n,n)/(n+1) (n=0,1,2,...)
   递推关系的另类解为：
   h(n)=c(2n,n)-c(2n,n+1)(n=0,1,2,...)
3. 推导：对于每一个数来说，必须进栈一次、出栈一次。我们把进栈设为状态‘1’，出栈设为状态‘0’。n个数的所有状态对应n个1和n个0组成的2n位二进制数。由于等待入栈的操作数按照1‥n的顺序排列、入栈的操作数b大于等于出栈的操作数a(a≤b)，因此输出序列的总数目=由左而右扫描由n个1和n个0组成的2n位二进制数，1的累计数不小于0的累计数的方案种数。显然，不符合要求的方案数为c(2n,n+1)。由此得出输出序列的总数目=c(2n,n)-c(2n,n+1)=c(2n,n)/(n+1)=h(n+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位数必对应一个不符合要求的数。

 

实际题目中，对 初始值的理解是关键：

 

1. 应用：出栈次序，3个元素，f(3)=h(3)=5,f(n)=h(n)
2. 应用：凸三角形划分:在一个凸多边形中，通过若干条互不相交的对角线，把这个多边形划分成了若干个三角形。任务是键盘上输入凸多边形的边数n，求不同划分的方案数f（n）。比如当n=6时，f(6)=h(4)=14，f(n)=h(n-2).因为凸多边形的任意一条边必定属于某一个三角形，所以我们以某一条边为基准，以这条边的两个顶点为起点P1和终点Pn（P即Point），将该凸多边形的顶点依序标记为P1、P2、……、Pn，再在该凸多边形中找任意一个不属于这两个点的顶点Pk（2<=k<=n-1），来构成一个三角形，用这个三角形把一个凸多边形划分成两个凸多边形，其中一个凸多边形，是由P1，P2，……，Pk构成的凸k边形（顶点数即是边数），另一个凸多边形，是由Pk，Pk+1，……，Pn构成的凸n-k+1边形。
3. 应用：
   买票找零f(n) = h(n)
   有2n个人排成一行进入剧场。入场费5元。其中只有n个人有一张5元钞票，另外n人只有10元钞票，剧院无其它钞票，问有多少中方法使得只要有10元的人买票，售票处就有5元的钞票找零？(将持5元者到达视作将5元入栈，持10元者到达视作使栈中某5元出 栈)

   小兔的棋盘f(n) = h(n)

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

   最后上一个求Catalan数的codes

    

   from functools import lru_cache
   class Solution:
       def catalan(self, n):
           @lru_cache(None)
           def c(n,m):
               if (m == n):
                   return 1
               elif (m == 1):
                   return n
               return c(n-1,m-1)+c(n-1,m)
           return c(2*n,n)-c(2*n,n-1)
   
       @lru_cache(None)
       def catalan2(self,n):
           if (n == 0 or n == 1):
               return 1
           res = 0
           for i in range(0,n):
               res += self.catalan2(i)*self.catalan2(n-i-1)
           return res
   s = Solution()
   print(s.catalan2(3))