第八周 项目3-对称矩阵压存储的实现与应用

目的：运用串解决问题。

代码：

#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[][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 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[][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,0};      int b[M]={0,1,2,3,4,5,6,7,8,9};      int c1[N][N],c2[N][N];      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;  }