6. 矩阵的快速转置算法

矩阵的快速转置算法

数据压缩是提高传输、存储效率一种技术。教材第5章介绍了两种简单的压缩存储方法。
本实验要求实现三元组顺序表表示下的矩阵快速转置算法。

输入：
稀疏矩阵的行数、列数、非零元个数(三个数都大于0)
以行为主序输入稀疏矩阵三元组表
输出：
辅助数组num[ ]
辅助数组cpot[ ]
以行为主序输出对应的转置矩阵三元组表

例
输入：
6 7 8
1 2 12
1 3 9
3 1 -3
3 6 14
4 3 24
5 2 18
6 1 15
6 4 -7
输出：
num:2,2,2,1,0,1,0,
cpot:1,3,5,7,8,8,9,
1,3,-3
1,6,15
2,1,12
2,5,18
3,1,9
3,4,24
4,6,-7
6,3,14

代码

#include<iostream>
#include<cstring>
using namespace std;

struct SparseMatrix //稀疏矩阵定义
{
    int Rows, Cols, Terms;//行数、列数、非零元素个数
    struct //三元组表
    {
        int row, col, value;//行号、列号、非零元素的值
    }elem[105];
};

int main()
{
    SparseMatrix Matrix_A, Matrix_T;
    scanf("%d %d %d", &Matrix_A.Rows, &Matrix_A.Cols, &Matrix_A.Terms);
    int num[105], cpot[105];//辅助数组
    memset(num, 0, sizeof(num));//先置零，用来统计矩阵中各列非零元素数
    for (int i = 1; i <= Matrix_A.Terms; i++)
        scanf("%d %d %d", &Matrix_A.elem[i].row, &Matrix_A.elem[i].col, &Matrix_A.elem[i].value);
    for (int i = 1; i <= Matrix_A.Terms; i++)
        num[Matrix_A.elem[i].col]++;
    cpot[1] = 1;//计算开始转置后各行开始位置
    for (int i = 2; i <= Matrix_A.Cols; i++)
        cpot[i] = cpot[i - 1] + num[i - 1];
    printf("num:");
    for (int i = 1; i <= Matrix_A.Cols; i++)
        printf("%d,", num[i]);
    printf("\ncpot:");
    for (int i = 1; i <= Matrix_A.Cols; i++)
        printf("%d,", cpot[i]);
    printf("\n");
    for (int i = 1; i <= Matrix_A.Terms; i++)//从Matrix_A向Matrix_T传送
    {
        int t = cpot[Matrix_A.elem[i].col];
        Matrix_T.elem[t].row = Matrix_A.elem[i].row;
        Matrix_T.elem[t].col = Matrix_A.elem[i].col;
        Matrix_T.elem[t].value = Matrix_A.elem[i].value;
        cpot[Matrix_A.elem[i].col]++;//修改转置后第i行存放位置
    }
    for (int i = 1; i <= Matrix_A.Terms; i++)
        printf("%d,%d,%d\n", Matrix_T.elem[i].col, Matrix_T.elem[i].row, Matrix_T.elem[i].value);
    //system("pause");
    return 0;
}