第八周 项目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;  
}