第九周项目2-对称矩阵压缩存储的实现与应用（2）

/*  
 * Copyright (c) 2016, 烟台大学计算机与控制工程学院  
 * All rights reserved。  
 * 文件名称 ：1.cpp  
 * 作    者 ：杨俊杰
 * 完成日期 ：2016年 11月3日  
 * 版 本 号 ：v1.0  
 * 问题描述 ：设计算法，实现两个用压缩形式存储的对称矩阵A和B的加法和乘法。实现中请使用好前面设计的基本运算。
 * 输入描述 ：
 * 输出描述 ：

 */ 

[cpp] view plain copy

 

1. #include <stdio.h>    
2. #define N 4    
3. #define M 10    
4. int value(int a[],int i,int j)    
5. {    
6.     if (i>=j)    
7.         return a[(i*(i+1))/2+j];    
8.     else    
9.         return a[(j*(j+1))/2+i];    
10. }    
11. void madd(int a[],int b[],int c[][N])    
12. {    
13.     int i,j;    
14.     for (i=0; i<N; i++)    
15.         for (j=0; j<N; j++)    
16.             c[i][j]=value(a,i,j)+value(b,i,j);    
17. }    
18. void mult(int a[],int b[],int c[][N])    
19. {    
20.     int i,j,k,s;    
21.     for (i=0; i<N; i++)    
22.         for (j=0; j<N; j++)    
23.         {    
24.             s=0;    
25.             for (k=0; k<N; k++)    
26.                 s=s+value(a,i,k)*value(b,k,j);    
27.             c[i][j]=s;    
28.         }    
29. }    
30. void disp1(int a[])    
31. {    
32.     int i,j;    
33.     for (i=0; i<N; i++)    
34.     {    
35.         for (j=0; j<N; j++)    
36.             printf("%4d",value(a,i,j));    
37.         printf("\n");    
38.     }    
39. }    
40. void disp2(int c[][N])    
41. {    
42.     int i,j;    
43.     for (i=0; i<N; i++)    
44.     {    
45.         for (j=0; j<N; j++)    
46.             printf("%4d",c[i][j]);    
47.         printf("\n");    
48.     }    
49. }    
50. int main()    
51. {    
52.     int a[M]= {1,2,3,4,5,6,7,8,9,10};  //a表示压缩存储的对称矩阵    
53.     int b[M]= {1,1,1,1,1,1,1,1,1,1};    
54.     int c1[N][N],c2[N][N];    
55.     madd(a,b,c1);    
56.     mult(a,b,c2);    
57.     printf("a矩阵:\n");    
58.     disp1(a);    
59.     printf("b矩阵:\n");    
60.     disp1(b);    
61.     printf("a+b:\n");    
62.     disp2(c1);    
63.     printf("a×b:\n");    
64.     disp2(c2);    
65.     printf("\n");    
66.     return 0;    
67. }    

运行结果：