第八周项目三 对矩阵的压缩存储2

1. 烟台大学计算机学院  
2. 
   
3. 作者:王雪行
4.   
5. 问题描述：设计算法，实现两个用压缩形式存储的对称矩阵A和B的加法和乘法。 
6. 实现中请使用好前面设计的基本运算。  
7.   
8. 输入描述：无 
9.   
10. 输出描述：矩阵元素相加相乘的结果 
11.   
12. */   
13.   
14.   
15.   
16. #include <stdio.h>  
17. #define N 4  
18. #define M 10  
19. int value(int a[],int i,int j)  
20. {  
21.     if (i>=j)  
22.         return a[(i*(i+1))/2+j];  
23.     else  
24.         return a[(j*(j+1))/2+i];  
25. }  
26. void madd(int a[],int b[],int c[][N])  
27. {  
28.     int i,j;  
29.     for (i=0; i<N; i++)  
30.         for (j=0; j<N; j++)  
31.             c[i][j]=value(a,i,j)+value(b,i,j);  
32. }  
33. void mult(int a[],int b[],int c[][N])  
34. {  
35.     int i,j,k,s;  
36.     for (i=0; i<N; i++)  
37.         for (j=0; j<N; j++)  
38.         {  
39.             s=0;  
40.             for (k=0; k<N; k++)  
41.                 s=s+value(a,i,k)*value(b,k,j);  
42.             c[i][j]=s;  
43.         }  
44. }  
45. void disp1(int a[])  
46. {  
47.     int i,j;  
48.     for (i=0; i<N; i++)  
49.     {  
50.         for (j=0; j<N; j++)  
51.             printf("%4d",value(a,i,j));  
52.         printf("\n");  
53.     }  
54. }  
55. void disp2(int c[][N])  
56. {  
57.     int i,j;  
58.     for (i=0; i<N; i++)  
59.     {  
60.         for (j=0; j<N; j++)  
61.             printf("%4d",c[i][j]);  
62.         printf("\n");  
63.     }  
64. }  
65. int main()  
66. {  
67.     int a[M]= {1,2,3,4,5,6,7,8,9,10};  //a表示压缩存储的对称矩阵  
68.     int b[M]= {1,1,1,1,1,1,1,1,1,1};  
69.     int c1[N][N],c2[N][N];  
70.     madd(a,b,c1);  
71.     mult(a,b,c2);  
72.     printf("a矩阵:\n");  
73.     disp1(a);  
74.     printf("b矩阵:\n");  
75.     disp1(b);  
76.     printf("a+b:\n");  
77.     disp2(c1);  
78.     printf("a×b:\n");  
79.     disp2(c2);  
80.     printf("\n");  
81.     return 0;  
82. }  

运行结果：