动态规划之矩阵连乘

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

* p是一系列矩阵的系数
* n是矩阵的个数
* m[x][y]是当前第x个矩阵到第y个矩阵的最优连乘
* s[x][y]是记录当前第x个矩阵到第y个矩阵的最优断开位置(k)
*/


<!-- more -->
void (int *p,int n,int **m,int **s) {
for(int i=1;i<=n;i++) {
m[i][i] = 0;
}
for(int r=2;r<=n;r++) { //自底向上的动规
for(int i=1;i<=n-r+1;i++) {
int j=i+r-1;
m[i][j] = m[i+1][j]+p[i-1]*p[i]*p[j];
s[i][j] = i;
for(int k=i+1;k<j;k++) {
int t = m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
if(t<m[i][j]){
m[i][j] = t;
s[i][j] = k;
}
}
}
}
}