poj1651 矩阵链乘

最新推荐文章于 2019-07-09 17:05:30 发布

fisty

最新推荐文章于 2019-07-09 17:05:30 发布

阅读量780

点赞数

分类专栏： ------基础DP 文章标签： ACM dp

本文链接：https://blog.csdn.net/u013008291/article/details/43773433

版权

------基础DP 专栏收录该内容

76 篇文章 1 订阅

订阅专栏

以矩阵链ABCD为例

按照矩阵链长度递增计算最优值

矩阵链长度为1时，分别计算出矩阵链A、B、C、D的最优值

矩阵链长度为2时，分别计算出矩阵链AB、BC、CD的最优值

矩阵链长度为3时，分别计算出矩阵链ABC、BCD的最优值

矩阵链长度为4时，计算出矩阵链ABCD的最优值

动归方程：

k为矩阵链断开的位置

d数组存放矩阵链计算的最优值，d[i][j]是以第i个矩阵为首，第j个矩阵为尾的矩阵链的最优值，i > 0

m数组内存放矩阵链的行列信息，m[i-1]和m[i]分别为第i个矩阵的行和列（i = 1、2、3...）

代码实现时需要注意的问题：计算顺序！！！

因为你要保证在计算m[i][j]查找m[i][k]和m[k+1][j]的时候，m[i][k]和m[k+1][j]已经计算出来了。

观察坐标的关系如图：

#include <iostream>  
#include <cstdio>  
#include <cstring>  
using namespace std;  
#define SIZE 100  
#define INF 999999999  
int m[SIZE];        //存放矩阵链的行列信息，m[i-1]和m[i]分别为第i个矩阵的行和列（i = 1、2、3...）  
int d[SIZE][SIZE];  //存放矩阵链计算的最优值，d[i][j]为第i个矩阵到第j个矩阵的矩阵链的最优值，i > 0  

int Best_DP(int n)  
{  
    //把d[i][i]置为0，1 <= i < n  
    memset(d, 0, sizeof(d));  

    int len;  
    //递归计算矩阵链的连乘最优值  
    //len = 1，代表矩阵链由两个矩阵构成  
    for (len = 1; len < n; len++)  
    {  
        int i, j, k;  
        for (i = 1, j = i+len; j < n; i++, j++)  
        {  
            int min = INF; //无穷大  
            for (k = i; k < j; k++)
            {  
                count = d[i][k] + d[k+1][j] + m[i-1] * m[k] * m[j];
                if(count < min)  
                {  
                    min = count;  
                }  
            }  
            d[i][j] = min;  
        }  
    }  
    return d[1][n-1];  
}  

int main()  
{  
    int n;  
    while (scanf("%d", &n) != EOF)  
    {  
        int i;  
        for (i = 0; i < n; i++)  
        {  
            scanf("%d", &m[i]);  
        }  
        printf("%d\n", Best_DP(n));  
    }  
    return 0;  
}