第八周(2) 项目二(2) 压缩存储的对称矩阵的运算

1. <pre name="code" class="cpp">/*  
2. Copyright (c)2017,烟台大学计算机与控制工程学院  
3. All rights reserved.  
4. 文件名称：第8周(2)项目2 - 对称矩阵压缩存储的实现与应用（2）.cpp  
5. 作    者：姬长水
6. 完成日期：2017年11月9日  
7. 版 本 号：v1.0  
8.   
9. 问题描述: （2）设计算法，实现两个用压缩形式存储的对称矩阵A和B的加法和乘法。实现中请使用好前面设计的基本运算。   
10. 输入描述：a,b对称矩阵的全部数据。  
11. 程序输出：矩阵a,b,a+b与a*b的结果。  
12. */   
13. </pre><br>  
14. <br>  
15. <pre></pre>  
16. <pre name="code" class="cpp">#include <iostream>  
17. #include <cstdio>  
18. using namespace std;  
19. #define N 4  
20. #define M 10  
21. int value(int a[],int i,int j)  
22. {  
23.     if(i>=j)  
24.         return a[(i*(i+1))/2+j];  
25.     else  
26.         return a[(j*(j+1))/2+i];  
27. }  
28. void madd(int a[],int b[],int c[][N])  
29. {  
30.     int i,j;  
31.     for(i=0;i<N;i++)  
32.     for(j=0;j<N;j++)  
33.         c[i][j]=value(a,i,j)+value(b,i,j);  
34.   
35.     }  
36. void mult(int a[],int b[],int c[][N])  
37. {  
38.     int i,j,k,s;  
39.     for(i=0;i<N;i++)  
40.         for(j=0;j<N;j++)  
41.     {  
42.         s=0;  
43.         for(k=0;k<N;k++)  
44.             s=s+value(a,i,k)*value(b,k,j);  
45.         c[i][j]=s;  
46.     }  
47. }  
48. void disp1(int a[])  
49. {  
50.     int i,j;  
51.     for(i=0;i<N;i++)  
52.     {  
53.         for(j=0;j<N;j++)  
54.             printf("%4d",value(a,i,j));  
55.         printf("\n");  
56.     }  
57. }  
58. void disp2(int c[][N])  
59. {  
60.     int i,j;  
61.     for(i=0;i<N;i++)  
62.     {  
63.         for(j=0;j<N;j++)  
64.             printf("%4d",c[i][j]);  
65.         printf("\n");  
66.     }  
67. }  
68. int main()  
69. {  
70.     int a[M]={1,2,3,4,5,6,7,8,9,0};  
71.     int b[M]={0,1,2,3,4,5,6,7,8,9};  
72.     int c1[N][N],c2[N][N];  
73.     madd(a,b,c1);  
74.     mult(a,b,c2);  
75.   
76. cout<<"a矩阵"<<endl;  
77. disp1(a);  
78. cout<<"b矩阵"<<endl;  
79. disp1(b);  
80. cout<<"a+b："<<endl;  
81. disp2(c1);  
82. cout<<"a*b："<<endl;  
83. disp2(c2);  
84. cout<<endl;  
85.   
86.     return 0;  
87. }  
88. </pre><pre name="code" class="cpp"><img src="https://img-blog.csdn.net/20171109111955514?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQvcmVueXVhbnN1bg==/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center" alt="">  
89. </pre><br>