矩阵连乘积 c语言程序,（用括号表示出结果）矩阵连乘C语言描述

1.矩阵连乘

矩阵相乘：

矩阵A和B相乘，p、q分别为A的行和列，q、r分别为B的行和列：

所以A和B相乘，两个矩阵相乘的计算量为pqr。

1.1 完全加括号的矩阵连乘积

完全加括号的矩阵连乘积可递归地定义为：

单个矩阵是完全加括号的；

矩阵连乘积 A是完全加括号的，则 A 可表示为2个完全加括号的矩阵连乘积 B 和 C 的乘积并加括号，即 A = (BC) 。

由于矩阵乘法满足结合律，所以计算矩阵的连乘可以有许多不同的计算次序。这种计算次序可以用加括号的方式来确定。

若一个矩阵连乘积的计算次序完全确定，也就是说该连乘积已完全加括号，则可以依此次序反复调用2个矩阵相乘的标准算法计算出矩阵连乘积。

给定n个矩阵｛A1,A2,…,An｝，其中Ai与Ai+1是可乘的，i=1，2…，n-1。如何确定计算矩阵连乘积的计算次序，使得依此次序计算矩阵连乘积需要的数乘次数最少。

代码如下：```

#include

#include

//---------------------------------------------------------------------------

//动态创建二维数组

template

T **creta2Darray(int n) {

T **p = new T *[n];

for(int i=0;i

p[i]=new int[n];

for(int j=0;j

}

return p;

}

//---------------------------------------------------------------------------

//输出二维数组，屏蔽第一维

void disp2Darray(int **p,int n) {

for(int i=1;i

for(int j=1;j

putchar('\n');

}

}

//---------------------------------------------------------------------------

//动态规划法计算矩阵连乘的最优解

void MatrixChain(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++){ //r为当前计算的链长(子问题规模)

for(int i=1;i<=n-r+1;i++){ //n-r+1为最后一个r链的前边界

int j = i+r-1; //计算前边界为r，链长为r的链的后边界

//将链ij划分为A(i) * ( A[i+1:j] )

m[i][j] = m[i+1][j] + p[i-1]*p[i]*p[j];

s[i][j] = i; //记录断开点的索引

for(int k = i+1 ; k

//将链ij划分为( A[i:k] )* (A[k+1:j])

int t = m[i][k] + m[k+1][j] + p[i-1] *p[k]*p[j];

if(t

m[i][j] = t; s[i][j] = k;

}

}

}

}

}

//---------------------------------------------------------------------------

//构造最优解

void Traceback(int i,int j,int **s) {

if(i==j)

return ;

if(i==s[i][j]&&(s[i][j]+1)==j)

{

printf("(");

printf("%d%d",i,j);

printf(")");

}

else if(i==s[i][j]&&(s[i][j]+1)!=j)

{

printf("(");

printf("%d",i);

Traceback(s[i][j]+1,j,s);

printf(")");

}

else if(i!=s[i][j]&&(s[i][j]+1)==j)

{

printf("(");

Traceback(i,s[i][j],s);

printf("%d",j);

printf(")");

}

else

{

Traceback(i,s[i][j],s);

Traceback(s[i][j]+1,j,s);

}

}

//---------------------------------------------------------------------------

//动态规划法

//---------------------------------------------------------------------------

void dyProg(int *p,int n,int **m,int **s) {

MatrixChain(p,n,m,s);

disp2Darray(m,n+1);

printf("\n");

disp2Darray(s,n+1);

Traceback(1,n,s);

}

//---------------------------------------------------------------------------

int main(int argc, char* argv[])

{

int p[] = {30,35,15,5,10,20,25};

int n = sizeof(p)/sizeof(int)-1;

int **m = creta2Darray(n+1);

int **s = creta2Darray(n+1);

//动态规划法

dyProg(p,n,m,s);

system("pause");

return 0;

}

//---------------------------------------------------------------------------

在构造最优解当中，利用了递归去打印最终结果的括号，结束条件是i=j，就是自己和自己打印的意思，然后我们分情况讨论，分四种情况，第一个if情况是只有两个矩阵打印，那直接两个相乘外面加括号即可，第二个else if情况是左边是就只有自己的就直接打印，第三个else if 是右边是就自己的直接打印，第四个else是左右两边都不只是能直接打印的，因此直接进入递归，从而我们能打印出来了最终结果。