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

问题：

[cpp]  view plain  copy

 

1. /*  
2. * Copyright (c)2016,烟台大学计算机与控制工程学院  
3. * All rights reserved.  
4. * 文件名称：项目2.cpp  
5. * 作    者：梁凯  
6. * 完成日期：2016年10月27日  
7. * 版 本 号：v1.0   
8. *问题描述：设计算法，实现两个用压缩形式存储的对称矩阵A和B的加法和乘法。实现中请使用好前面设计的基本运算。  
9.               
10. *输入描述：无  
11. *程序输出：测试结果  
12. */   

代码：

[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. }  

运行结果：