9 对称矩阵压缩储存的实现与应用

[html]  view plain  copy

 

1. /*    
2. * Copyright (c)2016,烟台大学计算机与控制工程学院    
3. * All rights reserved.    
4. * 文件名称：项目2.cpp    
5. * 作    者：泮春宇  
6. * 完成日期：2016年11月17    
7. * 版 本 号：v1.0     
8. *问题描述：（1）用压缩形式存储对称矩阵，实现下面的操作并测试    
9.                 
10.             void Init(int *&b);//为N阶对称矩阵初始化存储数据的一维数组B    
11.             int Value(int b[], int i, int j);//返回存储在b[M]中，对应二维数组A[i][j]的值    
12.             void Assign(int b[], int e, int i, int j);//将e赋值给对应二维数组元素A[i][j]，要存储到B[M]中    
13.             void Disp(int b[]);//输出压缩存储在b中的对称矩阵    
14.             void Destroy(int b[]);//销毁存储空间    
15.                 
16. *输入描述：无    
17. *程序输出：测试结果    
18. */      

[cpp]  view plain  copy

 

1. #include <stdio.h>    
2. #include <malloc.h>    
3. #define N 4    
4.     
5. //为N阶对称矩阵初始化存储数据的一维数组B    
6. void Init(int *&b)    
7. {    
8.     b = (int*)malloc(sizeof(int)*(N*(N+1)/2));    
9. }    
10.     
11. //返回存储在b[M]中，对应二维数组A[i][j]的值    
12. int Value(int b[], int i, int j)    
13. {    
14.     if (i>=j)    
15.         return b[(i*(i+1))/2+j];    
16.     else    
17.         return b[(j*(j+1))/2+i];    
18. }    
19.     
20. //将e赋值给对应二维数组元素A[i][j]，要存储到B[M]中    
21. void Assign(int b[], int e, int i, int j)    
22. {    
23.     if (i>=j)    
24.         b[(i*(i+1))/2+j] = e;    
25.     else    
26.         b[(j*(j+1))/2+i] = e;    
27.     return;    
28. }    
29.     
30. //输出压缩存储在b中的对称矩阵    
31. void Disp(int b[])    
32. {    
33.     int i,j;    
34.     for (i=0; i<N; i++)    
35.     {    
36.         for (j=0; j<N; j++)    
37.             printf("%4d",Value(b,i,j));    
38.         printf("\n");    
39.     }    
40. }    
41.     
42. //销毁存储空间    
43. void Destroy(int b[])    
44. {    
45.     free(b);    
46. }    
47.     
48. int main()    
49. {    
50.     int *b1;  //指向整型的指针，待初始化    
51.     int i, j;    
52.     int v;    
53.     Init(b1);    
54.     printf("请输入对称矩阵(只需要输入下三角部分即可)\n");    
55.     for(i=0;i<N;i++)    
56.     {    
57.         printf("输入第%d行的%d个数据元素: ", i+1, i+1);    
58.         for(j=0; j<=i; j++)    
59.         {    
60.              scanf("%d", &v);    
61.              Assign(b1, v, i, j);    
62.         }    
63.     }    
64.     Disp(b1);    
65.     Destroy(b1);    
66.     return 0;    
67. }    

运行结果：

知识点总结：

对称矩阵的操作

学习心得：

以前没有接触过对称矩阵，感觉很陌生，通过这次学习对其有了比较深入的了解