第八周(2) 项目二(2) 压缩存储的对称矩阵的运算

/*
Copyright (c)2017,烟台大学计算机与控制工程学院
All rights reserved.
文件名称：第8周(2)项目2 - 对称矩阵压缩存储的实现与应用（2）.cpp
作    者：孙仁圆
完成日期：2017年11月9日
版 本 号：v1.0

问题描述: （2）设计算法，实现两个用压缩形式存储的对称矩阵A和B的加法和乘法。实现中请使用好前面设计的基本运算。
输入描述：a,b对称矩阵的全部数据。
程序输出：矩阵a,b,a+b与a*b的结果。
*/

#include <iostream>
#include <cstdio>
using namespace std;
#define N 4
#define M 10
int value(int a[],int i,int j)
{
if(i>=j)
return a[(i*(i+1))/2+j];
else
return a[(j*(j+1))/2+i];
}
void madd(int a[],int b[],int c[]
)
{
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[]
)
{
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 disp1(int a[])
{
int i,j;
for(i=0;i<N;i++)
{
for(j=0;j<N;j++)
printf("%4d",value(a,i,j));
printf("\n");
}
}
void disp2(int c[]
)
{
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,0};
int b[M]={0,1,2,3,4,5,6,7,8,9};
int c1

,c2

;
madd(a,b,c1);
mult(a,b,c2);

cout<<"a矩阵"<<endl;
disp1(a);
cout<<"b矩阵"<<endl;
disp1(b);
cout<<"a+b："<<endl;
disp2(c1);
cout<<"a*b："<<endl;
disp2(c2);
cout<<endl;

return 0;
}