第九周 项目二【对称矩阵压缩存储的实现与应用2】

/*   
Copyright (c)2016,烟台大学计算机与控制工程学院   
All rights reserved.   
文件名称：第9周项目2 - 对称矩阵压缩存储的实现与应用（2）.cpp   
作    者：李潇   
完成日期：2016年10月31日   
版 本 号：v1.0      
问题描述: （2）设计算法，实现两个用压缩形式存储的对称矩阵A和B的加法和乘法。实现中请使用好前面设计的基本运算。    
输入描述：a,b对称矩阵的全部数据。   
程序输出：矩阵a,b,a+b与a*b的结果。   
*/    

#include <stdio.h>
#define N 4
#define M 10
int Value(int b[],int i,int j); //返回存储在b[M]中，对应二维数组A[i][j]的值
void madd(int a[],int b[],int c[][N]);
void mult(int a[],int b[],int c[][N]);
void Disp(int b[]);
void disp2(int c[][N]);

int Value(int b[],int i,int j)
{
	if(i>=j)
		return b[(i*(i+1))/2+j];
	else
		return b[(j*(j+1))/2+i];
}
void madd(int a[],int b[],int c[][N])
{
		int i,j;
		for(i=0;i<N;i++)
		{
			for(j=0;j<N;j++)
				c[i][j]=Value(a,i,j)+Value(b,i,j);
		}
}
void mult(int a[],int b[],int c[][N])
{
    int i,j,k,s;
    for (i=0; i<N; i++)
        for (j=0; j<N; j++)
        {
            s=0;
            for (k=0; k<N; k++)
                s=s+Value(a,i,k)*Value(b,k,j);
            c[i][j]=s;
        }
}
void Disp(int b[])
{
    int i,j;
	for(i=0;i<N;i++)
	{
		for(j=0;j<N;j++)
			printf("%4d",Value(b,i,j));
		printf("\n");
	}
}
void disp2(int c[][N])
{
   int i,j;
    for (i=0; i<N; i++)
    {
        for (j=0; j<N; j++)
            printf("%4d",c[i][j]);
        printf("\n");
    }
}
int main()
{
    int a[M]= {1,2,3,4,5,6,7,8,9,10};  //a表示压缩存储的对称矩阵
    int b[M]= {1,1,1,1,1,1,1,1,1,1};
    int c1[N][N],c2[N][N];
    madd(a,b,c1);
    mult(a,b,c2);
    printf("a矩阵:\n");
    Disp(a);
    printf("b矩阵:\n");
    Disp(b);
    printf("a+b:\n");
    disp2(c1);
    printf("a×b:\n");
    disp2(c2);
    printf("\n");
    return 0;
}

运行结果：

知识点总结：

对称矩阵的压缩存储与相加与相乘算法

心得体会：

我感觉相乘的算法不对，可是改了好多次试了好多次都不行，谁能给我解答一下