Uva 10007 / HDU 1131 - Count the Trees （卡特兰数）

 Count the Trees 

Another common social inability is known as ACM (Abnormally Compulsive Meditation). This psychological disorder is somewhat common among programmers. It can be described as the temporary (although frequent) loss of the faculty of speech when the whole power of the brain is applied to something extremely interesting or challenging.

Juan is a very gifted programmer, and has a severe case of ACM (he even participated in an ACM world championship a few months ago). Lately, his loved ones are worried about him, because he has found a new exciting problem to exercise his intellectual powers, and he has been speechless for several weeks now. The problem is the determination of the number of different labeled binary trees that can be built using exactlyn different elements.

For example, given one element A, just one binary tree can be formed (using A as the root of the tree). With two elements, A and B, four different binary trees can be created, as shown in the figure.

 

If you are able to provide a solution for this problem, Juan will be able to talk again, and his friends and family will be forever grateful.

Input

The input will consist of several input cases, one per line. Each input case will be specified by the numbern ( 1 ≤ n ≤ 300 ) of different elements that must be used to form the trees. A number 0 will mark the end of input and is not to be processed.

Output

For each input case print the number of binary trees that can be built using the n elements, followed by a newline character.

Sample Input

1
2
10
25
0

Sample Output

1
4
60949324800
75414671852339208296275849248768000000

---

Miguel Revilla 
2000-08-21

 

解题思路：result(n) = Canlatan(n)*n!

这题完全是水过去的，给你的测试数据前两个还是可以算出来，但第三个数据以后就很大了，因为是Cantalan的一个系列，所以我试出了10的sampleOutput刚好是:10的阶乘*Cantalan数，然后计算了25的输出，发现也是如此，在原有的代码上稍加改进在杭电过了之后交到了Uva，Uva的数据比较大 n<=300

 #include<cstdio>
 #include<cstring>
 #define SIZE 200
 #define MAXN 302
 #define BASE 10000

 using namespace std;

 int Canlan[MAXN][SIZE];
 int temp[SIZE];

 void A(int n)
 {
     memcpy(temp, Canlan[n], sizeof(temp));
     int  e = ;
     for(int fac = n; fac>; --fac)
     {
         for(int i = SIZE-, remain = ; i>=; --i)
         {
             e = temp[i] * fac + remain;
             temp[i] = e % BASE;
             remain = e / BASE;
         }
     }
     memcpy(Canlan[n], temp, sizeof(int)*SIZE);
 }

 void multiply(int elem)
 {
     for(int i=SIZE-, remain = ; i >= ; --i)
     {
         remain = temp[i] * elem + remain;
         temp[i] = remain % BASE;
         remain = remain / BASE;
     }
     return;
 }

 void divide(int elem)
 {
     int i;
     for(i=; i<SIZE && temp[i] == ; i++);
     for(int remain = ; i < SIZE; ++i)
     {
         remain += temp[i];
         temp[i] = remain / elem;
         remain = (remain % elem) * BASE;
     }
     return;
 }

 void Cantalan()
 {
     memset(Canlan[], , sizeof(int)*SIZE);
     Canlan[][SIZE-] = ;
     for(int i=; i<MAXN; ++i)
     {
         memcpy(temp, Canlan[i-], sizeof(int)*SIZE);
         multiply(*i-);
         divide(i+);
         memcpy(Canlan[i], temp, sizeof(int)*SIZE);
     }
     return;
 }

 int main()
 {
     int n;
     Cantalan();
     for(int i=; i<MAXN; ++i) A(i);
     while(scanf("%d", &n) != EOF && n)
     {
         int i;
         for(i=; i<SIZE && Canlan[n][i] == ; i++);
         printf("%d", Canlan[n][i]);
         if(i+ == SIZE)
         {
             printf("\n");
             continue;
         }
         for(i=i+; i<SIZE; ++i)
             printf("%04d", Canlan[n][i]);
         printf("\n");
     }
     return ;
 }