第九周--项目2对称矩阵压缩存储的实现与应用(1)
2015-10-29 18:57
302 查看
* * Copyright (c++) 2015 烟台大学计算机学院 * All right reserved. * 文件名称:project2.cpp * 作 者: 商文轲 * 完成日期:2015年10月29日 * 版 本 号:v1.9 * *问题描述:用压缩形式存储对称矩阵,实现下面的操作并测试 void Init(int *&b);//为N阶对称矩阵初始化存储数据的一维数组b int Value(int b[], int i, int j);//返回存储在b[M]中,对应二维数组A[i][j]的值 void Assign(int b[], int e, int i, int j);//将e赋值给对应二维数组元素A[i][j],要存储到b[M]中 void Disp(int b[]);//输出压缩存储在b中的对称矩阵 void Destroy(int b[]); //销毁存储空间 */
#include <stdio.h> #include <malloc.h> #define N 4 //为N阶对称矩阵初始化存储数据的一维数组B void Init(int *&b) { //压缩存贮只需要存储对角线及对角线以下元素 b = (int*)malloc(sizeof(int)*(N*(N+1)/2)); //数列求和N*(N+1)/2个元素 } //返回存储在b[M]中,对应二维数组A[i][j]的值 int Value(int b[], int i, int j) { if (i>=j) return b[(i*(i+1))/2+j]; //如果是对角线及对角线以下元素,从0到i-1有i行,共i*(i+1))/2个元素,再加上i行的j个元素 else return b[(j*(j+1))/2+i]; } //将e赋值给对应二维数组元素A[i][j],要存储到B[M]中 void Assign(int b[], int e, int i, int j) { if (i>=j) b[(i*(i-1))/2+j] = e; else b[(j*(j-1))/2+i] = e; return ; //回到调用地方 } //输出压缩存储在b中的对称矩阵 void Disp(int b[]) { int i,j; for (i=0; i<N; i++) { for (j=0; j<N; j++) printf("%4d",Value(b,i,j)); printf("\n"); } } //销毁存储空间 void Destroy(int b[]) { free(b); } int main() { int *b1; //指向整型的指针,待初始化 int i, j; int v; Init(b1); printf("请输入对称矩阵(只需要输入下三角部分即可)\n"); for(i=0;i<N;i++) { printf("输入第%d行的%d个数据元素: ", i+1, i+1); for(j=0; j<=i; j++) { scanf("%d", &v); Assign(b1, v, i, j); } } Disp(b1); Destroy(b1); return 0; }
总结:对称矩阵主要是利用对称的特点,储存矩阵时只储存对角线及上三角(或下三角)元素,实对称元素共享一个对称空间。思路不难想出,具体推算公式自己在至纸上试试更好理解记忆。
相关文章推荐
- C++primer学习:面象对象程序设计(5):容器与继承
- 数据结构之并查集
- javascript 时间处理问题
- appium for windows 环境搭建
- 第九周--项目1猴子选大王 (数组版)
- Dinic 算法求最大流(最小割) POJ 2536
- 关于UIInterfaceOritation 和 UIDeviceOritation
- export http://blog.csdn.net/hikaliv/article/details/4474835
- 机器学习2_逻辑回归和标准化的逻辑回归_AndrewNG
- 高通平台msm8909 LK 实现LCD 兼容
- Leetcode -- Best Time to Buy and Sell Stock IV
- 为Web程序员准备的10个最棒的jQuery视频插件
- 考研系列一-线性表类(顺序存储)
- 2015年9大优秀项目管理工具集锦
- 你的理想在哪
- 数据结构之后缀数组
- 2015年9大优秀项目管理工具集锦
- js移动div改进(兼容)
- php设计模式笔记--总结篇
- 在Eclipse中使用maven创建项目