矩阵连乘问题动态规划算法

#include"stdio.h"

int m[7][7],s[7][7];
int p[7];

	void MatrixChain(int n, int *p, int m[7][7], int s[7][7])
	{
		int i, j, k, r, tmp;
		for(i=1;i<=n;i++) m[i][i]=0;
		for(r=2;r<=n;r++)
		{
			for(i=1;i<=n-r+1;i++)
			{
				j=i+r-1;
				m[i][j]= m[i+1][j]+p[i-1]*p[i]*p[j];   
				s[i][j]=i;
				for(k=i+1;k<j;k++)
				{
					tmp=m[i][k]+m[k+1][j]+p[i-1]*p[k]*p[j];
					if(tmp<m[i][j])
					{
						m[i][j]=tmp;   
						s[i][j]=k;
					}
				}
			}
		}
	}

	void Traceback(int i, int j, int s[7][7]){
		if(i==j) return;
		Traceback(i, s[i][j], s);
		Traceback(s[i][j]+1, j, s);
		printf("Multiply A%d,%d",i,s[i][j]);
		printf("and A%d,%d\n",s[i][j]+1,j);
	}

int main()
{
	p[0]=30;p[1]=35;p[2]=15;p[3]=5;p[4]=10;p[5]=20;p[6]=25;
	MatrixChain(6, p, m, s);
	Traceback(1,6,s);
	return 0;
}