动态规划求解矩阵连乘问题
题目:
相关分析见参考链接,大神讲的很好
源代码
#include<stdio.h>
#define N 6
#define INF 0x3f3f3f
int r[N+1] = { 5,10,3,12,5,50,6};
int m[N + 1][N + 1];//存储矩阵Ai到Aj连乘的最少次数
int s[N + 1][N + 1];//存储矩阵连乘最少次数时对应的划分k值
void map()//迭代不重复的求出m[i][j]
{
int i,j,k,x;
for (i = 1; i <= N; i++)//单个矩阵m[i][i]=0
m[i][i] = 0;
for(x=1;x<N;x++)//控制循环次数
for (i = 1, j = x + 1; j <= N; i++, j++)
{
for (k = i; k < j; k++)//遍历k值
{
int x = m[i][k] + m[k + 1][j] + r[i - 1] * r[k] * r[j];
if (x < m[i][j])
{
m[i][j] = x;
s[i][j] = k;
}
}
//printf("m[%d][%d]:%d\n",i,j,m[i][j]);//查看m[i][j]值
}
}
void traceback(int i, int j)//输出次数
{
if (i == j)
{
printf("A%d", i);
return;
}
printf("(");
traceback(i, s[i][j]);//k值以前
printf("*");
traceback(s[i][j]+1, j);//k值以后
printf(")");
}
int main()
{
for (int i = 1; i <= N; i++)
for (int j = 1; j <= N; j++)
m[i][j] = INF;
map();
printf("最少乘法次数:%d\n", m[1][6]);
printf("最佳连乘次序为:");
traceback(1, 6);
return 0;
}