矩阵的压缩存储
压缩存储:为多个值相同的元素只分配一个存储空间,对0不分配存储空间
特殊矩阵:具有许多相同元素或0元素,分布有一定规律的矩阵(对称矩阵、上下三角矩阵、对角矩阵)
对称矩阵: n阶方阵中,按主对角线对称,只存放主对角线和下三角元素(或上三角元素,又或按行或者按列存储)(共n(n+1)/2)
三角矩阵: 上或下三角矩阵,其另外一角都为同一常数,存储完对角线和下三角元素后,紧接着存储上三角元素一次(共n(n+1)/2+1)
三对角矩阵: 所有元素集中在三条主对角线上,其余均为0;(i,j)元素在数组中的下标为k=2i+j-3,反之,若元素在数组中第k个位置,则 i=(k+1)/3+1向下取整,j=k-2i+3i 取下限
稀疏矩阵
定义:非零元较零元少,且分布没有一定规律的矩阵
压缩存储原则:只存矩阵的行列维数和每个非零元的行列下标及其值
三元组:一个稀疏矩阵所有非零元素三元组的集合,再增加一个表示矩阵行数,列数和非零元个数的特殊三元组,就可以唯一地确定一个稀疏矩阵,把这些三元组的总体,称为该稀疏矩阵的“三元组表
例题
6-1 三元组顺序表表示的稀疏矩阵转置
本题要求实现一个函数,实现三元组顺序表表示的稀疏矩阵转置。
函数接口定义:
struct tripletable * trans(struct tripletable *t1);
其中 t1
是用户传入的参数。 函数须返回转置后的三元组顺序表指针。
裁判测试程序样例:
#include <stdio.h>
#include <stdlib.h>
#define M 100
struct node{
int i,j,v;//三元组单个节点的行列号与元素的值
};
struct tripletable
{
struct node S[M];//结点集合
int m,n,t;//行列非零元个数
};
struct tripletable * create()
{ int i;
struct tripletable *head=(struct tripletable *)malloc(sizeof(struct tripletable));
scanf("%d%d%d",&(head->m),&(head->n),&(head->t));
for(i=0;i<head->t;i++)
scanf("%d%d%d",&(head->S[i].i),&(head->S[i].j),&(head->S[i].v));
return head;
}
void print(struct tripletable * head)
{
int i;
for(i=0;i<head->t;i++)
printf("%d %d %d\n",(head->S[i].i),(head->S[i].j),(head->S[i].v));
}
struct tripletable * trans(struct tripletable *t1);
int main()
{
struct tripletable * head,*t2;
head=create();
t2=trans(head);
print(t2);
return 0;
}
/* 请在这里填写答案 */
输入样例:
输入第1行为矩阵行数m、列数n及非零元素个数t。
按行优先顺序依次输入t行,每行3个数,分别表示非零元素的行标、列标和值。
3 4 3
0 1 -5
1 0 1
2 2 2
输出样例:
输出转置后的三元组顺序表结果,每行输出非零元素的行标、列标和值,行标、列标和值之间用空格分隔,共t行。
0 1 1
1 0 -5
2 2 2
struct tripletable* trans(struct tripletable* t1)//转置,行列互换
{
int** a = new int* [t1->t], j = 0;
//创建二维数组a,动态分配创建
for (int i = 0; i < t1->t; i++)
{
int* q = new int[3];
a[i] = q;
}
for (int i = 0; i < t1->t; i++)//遍历tripletable中元素
{
a[i][j++] = t1->S[i].j;
a[i][j++] = t1->S[i].i;
a[i][j++] = t1->S[i].v;
t1->S[i].v = 0;
j = 0;
}
for (int i = 0; i < t1->t; i++)//对数组进行排序
{
for (int j = 0; j < t1->t - 1 - i; j++)
{
if (a[j][0] > a[j + 1][0])
{
int *T1 = a[j];
int* T2 = a[j + 1];
a[j] = T2; a[j + 1] = T1;
}
}
}
for (int i = 0; i < t1->t; i++)
{
t1->S[i].i = a[i][0];
t1->S[i].j = a[i][1];
t1->S[i].v = a[i][2];
}
return t1;
}