ZJU 2604 Little Brackets - 动态规划 Catalan数变形

题目描述：

定义一个括号序列的Depth为最深的嵌套层数。问n对括号，depth为k的括号序列共有多少种。

分析：

很类似Catalan数的题目。回忆经典的Catalan数：f[n]表示n对括号能够组成多少种不同的序列。那么可以得出f[n]的递推式：

f[n]=∑f[i]*f[n-i-1](0<=i<n)

对于此题可以用类似的思路得到递推式。

设f[i][j]表示i对括号能组成多少个depth为j的的序列。那么初始条件为f[0][0]=f[1][1]=1。考虑f[i][j]，设与第一个前括号配对的后括号之间有k对括号（如下图），那么可以把f[i][j]分为两个子问题：括号里的k对括号，和括号后面的i-k-1对括号。

 ( ..[k对括号].. ) ..[i-k-1对括号]..

那么可以写出以下递推式：

f[i][j] = ∑ ( f[k][j-1]*g[i-k-1][j-1] + g[k][j-1]*f[i-k-1][j] ) ,(0<=k<i)，其中g[i][j]=∑f[i][k](0<=k<=j)。

显然最后数字很大，需要大数。

1. /*
2. ZJU2604 Little Brackets
3. */
5. #include <stdio.h>
6. #include <string.h>
8. #define clr(a) memset(a,0,sizeof(a))
9. #define cpy(a,b) memcpy(a,b,sizeof(b))
10. #define N 51
12. /****************************************************/
14. #define DIGIT 4 //一个数组元素存放的位数个数
15. #define DEPTH 10000 //DEPTH = 10DIGIT
16. #define MAX 500 //数组的长度
18. typedef int bignum_t[MAX+1];
20. //赋值
21. void value(bignum_t a,const char val[]){    //val为数字字符串 
22.     char buf[MAX*DIGIT+1],ch;
23.     int i,j;
24.     memset((void*)a,0,sizeof(bignum_t));
25.     memcpy(buf,val,sizeof(val));
26.     
27.     for (a[0]=strlen(buf),i=a[0]/2-1;i>=0;i--)
28.         ch=buf[i],buf[i]=buf[a[0]-1-i],buf[a[0]-1-i]=ch;
29.     for (a[0]=(a[0]+DIGIT-1)/DIGIT,j=strlen(buf);j<a[0]*DIGIT;buf[j++]='0');
30.     for (i=1;i<=a[0];i++)
31.     for (a[i]=0,j=0;j<DIGIT;j++)
32.         a[i]=a[i]*10+buf[i*DIGIT-1-j]-'0';
33.     for (;!a[a[0]]&&a[0]>1;a[0]--);
34. }
36. //输出
37. void Write(const bignum_t a){
38.     int i,j;
39.     for (printf("%d",a[i=a[0]]),i--;i;i--)
40.     for (j=DEPTH/10;j;j/=10)
41.         printf("%d",a[i]/j%10);
42.     puts("");
43. }
45. //加法
46. void add(bignum_t a,const bignum_t b){//高精度数a+b,结果在a中
47.     int i;
48.     for (i=1;i<=b[0];i++)
49.     if ((a[i]+=b[i])>=DEPTH) a[i]-=DEPTH,a[i+1]++;
50.     if (b[0]>=a[0]) a[0]=b[0];
51.     else
52.     for (;a[i]>=DEPTH&&i<a[0];a[i]-=DEPTH,i++,a[i]++);
53.     a[0]+=(a[a[0]+1]>0);
54. }
56. //大数乘大数 
57. void mul(bignum_t c,const bignum_t a,const bignum_t b){//高精度数a*b,结果在c中
58.     int i,j;
59.     memset((void*)c,0,sizeof(bignum_t));
60.     for (c[0]=a[0]+b[0]-1,i=1;i<=a[0];i++)
61.     for (j=1;j<=b[0];j++)
62.     if ((c[i+j-1]+=a[i]*b[j])>=DEPTH)
63.     c[i+j]+=c[i+j-1]/DEPTH,c[i+j-1]%=DEPTH;
64.     for (c[0]+=(c[c[0]+1]>0);!c[c[0]]&&c[0]>1;c[0]--);
65. }
67. /**************************************************/
69. bignum_t f[N][N];
70. bignum_t g[N][N];
71. bignum_t x,y;
73. #define G(a,b) ((a)<(b)?(g[a][a]):(g[a][b]))
75. int main()
76. {
77.     int i,j,k,n,m,T=0;
78.     for(i=0;i<N;i++) for(j=0;j<N;j++){
79.         value(f[i][j],"0"); value(g[i][j],"0");
80.     }
82.     //DP
83.     value(f[0][0],"1"); value(f[1][1],"1");
84.     value(g[0][0],"1"); value(g[1][1],"1");
85.     for(i=2;i<N;i++) for(j=1;j<=i;j++){
86.         for(k=0;k<i;k++){
87.             mul(x,f[k][j-1],G(i-k-1,j-1));
88.             mul(y,G(k,j-1),f[i-k-1][j]);
89.             add(x,y);
90.             add(f[i][j],x);
91.         }
92.         cpy(g[i][j],g[i][j-1]);
93.         add(g[i][j],f[i][j]);
94.     }
95.     
96.     //input
97.     while(scanf("%d%d",&n,&m),n+m){
98.         if(T) puts("");
99.         printf("Case %d: ",++T);
100.         Write(f[n][m]);
101.     }
102.     
103.     return 0;
104. }