【卡特兰数+高精度】树屋阶梯

暑假期间，小龙报名了一个模拟野外生存作战训练班来锻炼体魄，训练的第 一个晚上，教官就给他们出了个难题。由于地上露营湿气重，必须选择在高处的 树屋露营。小龙分配的树屋建立在一颗高度为 N+1 尺（ N 为正整数）的大树上， 正当他发愁怎么爬上去的时候，发现旁边堆满了一些空心四方钢材（如图 1.1 ）， 经过观察和测量，这些钢材截面的宽和高大小不一，但都是 1 尺的整数倍，教官 命令队员们每人选取 N 个空心钢材来搭建一个总高度为 N 尺的阶梯来进入树屋， 该阶梯每一步台阶的高度为 1 尺，宽度也为 1 尺。如果这些钢材有各种尺寸，且 每种尺寸数量充足，那么小龙可以有多少种搭建方法？（注：为了避免夜里踏空， 钢材空心的一面绝对不可以向上。）

 

   以树屋高度为4尺、阶梯高度N=3尺为例，小龙一共有如图1.2所示的5种

   搭 建方法：

  

 

Input

一个正整数 N(1≤N≤500)，表示阶梯的高度

Output

一个正整数，表示搭建方法的个数。（注：搭建方法个数可能很大。）

Sample Input

3

Sample Output

5

Hint

 

1  ≤N≤500

对于每一个给出的阶梯的高度，我们可以观察上图，然后将包含左上角的矩形删    掉, 以后，原图最多被分成了两部分的新的阶梯（也可以是一个）,且高度是i和n-i-1那么这两个部分相对独立，运用乘法原理f[n]=∑f[i]*f[n-i-1]（与公式2相同）所以得出是卡特兰数，h(n)=h(n-1)*(4*n-2)/(n+1) 或h(n)=C(2n,n)/(n+1) (n=0,1,2,...)

#include <math.h>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#define maxn 10001
int n,len=1,ans[10000]={1,1};
void mult(int x)           //高精度乘法
{
    int i;
    for( i=1;i<=len;i++)
        ans[i]=ans[i]*x;
    for( i=1;i<=len;i++)
        if(ans[i]>9)
        {
            ans[i+1]=ans[i+1]+ans[i]/10;
            ans[i]=ans[i]%10;
        }
    while(ans[len+1])
    {
        len++;
        ans[len+1]=ans[len+1]+ans[len]/10;
        ans[len]=ans[len]%10;
    }
}
void divide(int x)          //高精度除法
{
    int a[maxn],b[maxn],tmp=0,i;
    for( i=1;i<=len;i++)
        a[i]=ans[len-i+1];         //先转置
    for( i=1;i<=len;i++)
    {
        b[i]=(tmp*10+a[i])/x;          //注意前几位可能为0
        tmp=(tmp*10+a[i])%x;
    }
    //for( i=1;i<=len;i++)
     //   printf("%d",b[i]);        //若n=10，b数组为016796
    //printf("\n");
    for( i=1;i<=len;i++)
        ans[i]=b[len-i+1];        //转置，再从len开始输出
    while(!ans[len])
        len--;
}
void Catalan(int k)       //h(n)=h(n-1)*(4*n-2)/(n+1)
{                         //h(n)=C(2n,n)/(n+1) (n=0,1,2,...)
    mult(4*k-2);
    divide(k+1);
}
int main()
{
    int i;

    ans[1]=1;
    scanf("%d",&n);
    for(i=2;i<=n;i++)
        Catalan(i);        //卡特兰数：1, 1, 2, 5, 14, 42, 132, 429, 1430...
    for(i=len;i>=1;i--)    //h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)*h(0) (n>=2)
        printf("%d",ans[i]);
    return 0;
}