【卡特兰数】HDU2067小兔的棋盘

  
  

   
   

    
    Problem Description
   
   
   
   

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

    
     
   
   
   
   

   
   

    
    Input
   
   
   
   

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

    
     
   
   
   
   

   
   

    
    Output
   
   
   
   

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

    
     
   
   
   
   

   
   

    
    Sample Input
   
   
   
   

    
    
     
     

      
      1
3
12
-1
     
     
   
   
   
   

    
     
   
   
   
   

   
   

    
    Sample Output
   
   
   
   

    
    
     
     

      
      1 1 2
2 3 10
3 12 416024
     
     
     
     


     
     
   
   
  
  

带限制条件的路径总数

首先我们来看一个问题：

在一个平面直角坐标系中，只能往右或往上走一个单位长度，问有多少种不同的路径可以从左下角 (1,1) (1,1) 走到右上角  (n,n) (n,n)，并且要求路径不能经过直线  y=x y=x 上方的点，下图中的路径都是合法的（图片来源 Wikipedia）

如果没有限制条件，那么从左下角走到右上角一共有  2n 2n 步，有  n n 步是往右，另外  n n 步是往上，那么路径方案数就是  2n 2n 步中选择  n n 步往右，一共有  (2nn) (2nn) （即  Cn2n C2nn）种方案

那么我们考虑一下这里面有多少种方案是不合法的

首先对于每一种不合法的方案，它的路径一定与  y=x+1 y=x+1 有交。我们找到它与  y=x+1 y=x+1 的第一个交点，然后将这个点后面部分的路径关于  y=x+1 y=x+1 做一个对称。由于原来路径到达  (n,n) (n,n)，新的对称之后的路径就会到达  (n−1,n+1) (n−1,n+1)。这样我们把每一种不合法方案都对应到了一条从  (1,1) (1,1) 到  (n−1,n+1) (n−1,n+1) 的路径，现在再来看是否每一条这样的路径都能对应到一种不合法方案，如果是，那么这就建立了一个一一映射的关系，也就是它们的方案总数相同。这是肯定的，因为每一条这样的路径必定与  y=x+1 y=x+1 有交，那么对称回去，就得到一条不合法方案

由于从  (1,1) (1,1) 到  (n−1,n+1) (n−1,n+1) 的路径有  (n−1+n+1n−1) (n−1+n+1n−1) 条，那么合法的方案就是

。

我们把这个方案数记为  Cn Cn，这就是著名的 Catalan 数

我们来看看它的前几项（ n n 从  0 0 开始）

1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900, 2674440, 9694845, 35357670, 129644790, 477638700, 1767263190

代码：

根据卡特兰数公式1：h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)h(0) (n>=2)

#include<iostream>
using namespace std;
long long a[40];
//  公式1：h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)h(0) (n>=2)
void Catalan1(int i,int j)
{
    for(;i<36;i++){
        j=0;a[i]=0;
        for(;j<=i;j++){
            a[i]+=a[j]*a[i-j-1];
        }
    }
}
//  公式2：h(n)=h(n-1)*(4*n-2)/(n+1);
void Catalan2()
{
    for(int i=2;i<36;i++)
        a[i]=a[i-1]*1.0/(i+1)*(4*i-2);      //  因为浮点数的问题，从a[33]开始就会用误差了；所以这个公式过不了
}
int main()
{
    a[0]=1;a[1]=1;
    Catalan1(1,0);
    //Catalan2();
    int t=1,n;
    while(cin>>n&&n!=-1){
        cout<<t++<<' '<<n<<' '<<2*a[n]<<endl;           //  对称，所以是a[n]*2；
    }
    return 0;
}