一看就知道是dp,再一看,才知道是矩阵链乘问题(不过刚开始时的状态方程还不是这样的,囧。。。)。在算导那里套了一下代码,外加解题报告的边界处理,水果了~
以下是代码:
- #include "stdio.h"
const int M=200;
const int inf=1<<29;- int main()
{
int dp[M][M],n,a[M];
int i,j,k,t,l;
scanf("%d",&n);
for(i=0;i<n;i++)
scanf("%d",&a[i]);
n--;
for(i=1;i<=n;i++)
dp[i][i]=0;
for(l=2;l<=n;l++)
{
for(i=1;i<=n-l+1;i++)
{
j=i+l-1;
dp[i][j]=inf;
for(k=i;k<j;k++)
{
t=dp[i][k]+dp[k+1][j]+a[i-1]*a[k]*a[j];
if(t<dp[i][j])
dp[i][j]=t;
}
}
}
printf("%d/n",dp[1][n]);
return 0;
}