源代码: #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])。 */