ACMSTEP 2.3.8 小兔的棋盘 //高精度 卡特兰数

原题链接


卡特兰数有一个重要的应用：
　　求将一个凸多边形区域分成三角形区域的方法数。

　　类似：一位大城市的律师在她住所以北n个街区和以东n个街区处工作。每天她走2n个街区去上班。如果她从不穿越（但可以碰到）从家到办公室的对角线，那么有多少条可能的道路？

　　类似：在圆上选择2n个点,将这些点成对连接起来使得所得到的n条线段不相交的方法数

此题恰好为这个应用不过并不是原本的卡特兰数 恰好为所有卡特兰数的2倍。和之前How Many Trees那道题一样（那道题的应用时卡特兰数的二叉树计数）只改了几个地方，一次AC

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <vector>
#include <queue>
#include <algorithm>

using namespace std;

int a[105][80];    //catalan number
int b[105];         //the length of catalan

void catalan()  
{
    int i, j, len, carry, temp;
    a[1][0] = b[1] = 2;
    len = 1;
    for(i = 2; i <= 100; i++)
    {
        for(j = 0; j < len; j++)    //multi
            a[i][j] = a[i-1][j]*(4*(i-1)+2);
        carry = 0;
        for(j = 0; j < len; j++)    
        {
            temp = a[i][j] + carry;
            a[i][j] = temp % 10;
            carry = temp / 10;
        }
        while(carry)    //carried processing
        {
            a[i][len++] = carry % 10;
            carry /= 10;
        }
        carry = 0;
        for(j = len-1; j >= 0; j--) //div
        {
            temp = carry*10 + a[i][j];
            a[i][j] = temp/(i+1);
            carry = temp%(i+1);
        }
        while(!a[i][len-1]) 	//drop the zeros
            len --;
        b[i] = len;
    }
}

int main()
{
    int i, n, cases=1;
    catalan();
    while(cin>>n)
    {
    	if (n==-1) break;
    	if (n==1)
		{
			cout<<"1 1 2"<<endl;
			cases++;
			continue;
		}
		cout<<cases<<" "<<n<<" ";
        for(i = b[n]-1; i>=0; i--)
        {
            cout<<a[n][i];
        }
       cout<<endl;
       cases++;
    }

    return 0;
}