稀疏矩陣之三元組壓縮儲存並且轉置

---

程式碼演示：

// δ=t/(n?m) <= 0.05 則證明是稀疏矩陣
/*  稀疏矩陣的型別說明及轉置演算法  */
#include<stdio.h>
#include<stdlib.h>

#define MAX 12

typedef int datatype;

typedef struct{
    int i, j;          // 行號、列號
    datatype  v;       // 元素值
} node;

typedef struct{
    int m, n, t;        // 行數，列數，非零元素個數
    node data[MAX];     // 三元組表
}spMatrix;

spMatrix *Create();   // 建立一個n行m列的稀疏矩陣
spMatrix *TransMat(spMatrix *a);    //稀疏矩陣的轉置
void Output(spMatrix *a);     //在螢幕上以行列的形式輸出矩陣

int main(){
    printf("輸入行列數:\n");
    spMatrix *a = Create();
    spMatrix *b = TransMat(a);
    printf("原矩陣:\n");
    Output(a);
    printf("轉置矩陣:\n");
    Output(b);
    return 0;
}

spMatrix *Create(){ // 建立一個3行4列的稀疏矩陣
    int m , n , k = 0, t = 0;
    scanf("%d%d",&m,&n);
    datatype element;
    spMatrix *matrix;
    matrix = (spMatrix *)malloc(sizeof(spMatrix));   //建立一個稀疏矩陣a
    matrix->m = m;
    matrix->n = n;
    printf("輸入%d行，%d列元素:\n",m,n);
    while (k < m*n){
        scanf("%d", &element);
        if (element != 0){
            matrix->data[t].i = k / n;
            matrix->data[t].j = k % n;
            matrix->data[t].v = element;
            t++;
        }
        k++;
    }
    matrix->t = t;
    return matrix;
}

spMatrix *TransMat(spMatrix *a){   // 稀疏矩陣的轉置
    int p, q, bno = 0;
    spMatrix *b;
    b = (spMatrix *)malloc(sizeof(spMatrix));   //為矩陣b分配記憶體空間
    b->m = a->n;
    b->n = a->m;
    b->t = 0;
    if (a->t == 0)  // 若b中元素全為零，則將b返回
        return b;
    for (p = 0; p < a->n; p++)
        for (q = 0; q < a->t; q++)
            if (a->data[q].j == p){
                b->data[bno].i = a->data[q].j;
                b->data[bno].j = a->data[q].i;
                b->data[bno].v = a->data[q].v;
                bno++;
            }
    b->t = bno;
    return b;
}

void Output(spMatrix *a){
    //輸出:在螢幕上以行列的形式輸出矩陣
    int i = 0, j = 0, k = 0;
    for (i = 0; i < a->m; i++){
        for (j = 0; j < a->n; j++){
            if (i == a->data[k].i && j == a->data[k].j){
                printf("%d ", a->data[k].v);
                k++;
            }
            else
                printf("%d ", 0);
        }
        printf("\n");
    }
}

效果展示：