括号

源代码:

#include<cstdio>
int n,k;
double L[201][51],R[201][51],f[201][201][2];
int main()
{
    scanf("%d%d",&n,&k);
    for (int a=1;a<=n;a++)
      for (int b=1;b<=k;b++)
        scanf("%lf%lf",&L[a][b],&R[a][b]);
    for (int a=0;a<=n;a++) //其实没有意义,就是套进去进行更简单的预处理。
      f[a+1][a][0]=1;
    for (int Length=2;Length<=n;Length+=2)
      for (int a=1;a<=n-Length+1;a++)
      {
        for (int b=1;b<=k;b++)
          f[a][a+Length-1][0]+=(f[a+1][a+Length-2][0]+f[a+1][a+Length-2][1])*L[a][b]*R[a+Length-1][b];
        for (int b=a+1;b<=a+Length-2;b+=2)
          f[a][a+Length-1][1]+=f[a][b][0]*(f[b+1][a+Length-1][0]+f[b+1][a+Length-1][1]);
      }
    double Ans=f[1][n][0]+f[1][n][1];
    printf("%.5lf",Ans); //这个数据类型还是挺烦人的。
    return 0;
}

/*
    好神奇的区间型动态规划。
    设f[i][j][0]表示(A)的方案数,f[i][j][1]表示AB的方案数,则有如下状态转移方程:
        f[i][j][0]=L[i][k]*R[i][k]*(f[i+1][j-1][0]+f[i+1][j-1][1])
        f[i][j][1]=f[i][S][0]*(f[S+1][j][0]+f[S+1][j][1])
    第二个避免了断点处理的重复。
    那么最终答案为(f[1][n][0]+f[1][n][1])。
*/
posted @ 2016-11-13 14:27  前前前世。  阅读(195)  评论(0编辑  收藏  举报