稀疏矩阵的压缩存储与操作（实测程序）

第1关：基于三元表顺序存储的稀疏矩阵的转置与输出（测试）

 //矩阵的建立、转置、输出
 #include<iostream>
 using namespace std;

#define MAX_SIZE 101
typedef int Elemtype;
typedef struct
{   int i;     //  行下标
    int j;        //  列下标
    Elemtype e;      // 元素值
}Triple;

typedef struct
{
    int  mu;         //   行数
    int  nu;         //   列数
    int  tu;         //   非0元素个数
    Triple  data[MAX_SIZE+1];   //data[0]不用
}TSMatrix;

void InitMatrix(TSMatrix *M)
{
    int m,n,t;
    cin>>m>>n>>t;
    M->mu=m;
    M->nu=n;
    M->tu=t;
    for(int k=1;k<=t;k++)
    {
        cin>>M->data[k].i>>M->data[k].j>>M->data[k].e;
    }

}

void PrintMatrix(TSMatrix *M)
{
    int k=1,p,q;
    for(p=1;p<=M->mu;p++){
        for(q=1;q<=M->nu;q++){
            if(M->data[k].i==p&&M->data[k].j==q)
            {
                cout<<M->data[k].e<<" ";
                k++;
            }
            else
                cout<<0<<" ";
        }
        cout<<endl;
    }
}

void TransMatrix(TSMatrix M , TSMatrix &T)
{
    int p , q , col ;
    T.mu=M.nu ;
    T.nu=M.mu ;
    T.tu=M.tu ;
    if(T.tu==0)
        cout<<"The Matrix M=0\n";
    else
    {
        q=1;
        for(col=1;col<=M.nu;col++)
            for(p=1;p<=M.tu;p++)
                if(M.data[p].j==col)
                {
                    T.data[q].i=M.data[p].j;
                    T.data[q].j=M.data[p].i;
                    T.data[q].e=M.data[p].e;
                    q++;
                }
    }
}
int main()
{
    TSMatrix M,T;
    InitMatrix(&M);
    TransMatrix(M,T);
    PrintMatrix(&T);
    return 0;
}

第2关：基于三元组顺序表存储的稀疏矩阵相加（测试）

//矩阵的相加
 
#include<iostream>
using namespace std;
#define MAX_SIZE 101
typedef int Elemtype ;

typedef struct
{
    int i;
    int j;
    Elemtype e;   //行、列、值
}Triple;
typedef struct lnode
{
    Triple data[MAX_SIZE];   //非零元素信息
    int mu,nu,tu;     //矩阵总行数、列数、非零元素个数
}TSMatrix;

void Input(TSMatrix *M);   //矩阵输入
void Print(TSMatrix M);    //矩阵输出
int Add(TSMatrix *c,TSMatrix a,TSMatrix b);    //矩阵加法

int main()
{
    TSMatrix A,C,D;
    Input(&A);
    Input(&C);
    if (Add(&D,A,C))
    {
        Print(D);
    }
    else
        cout<<0;
    return 0;
}

void Input(TSMatrix *M)
{
    int i;
    cin>>M->mu>>M->nu>>M->tu;
    for (i = 0; i < M->tu; i++)
    {
        cin>>M->data[i].i>>M->data[i].j>>M->data[i].e;
    }

}
void Print(TSMatrix M)
{
    int k=0,p,q;
    for(p=1;p<=M.mu;p++){
        for(q=1;q<=M.nu;q++){
            if(M.data[k].i==p&&M.data[k].j==q)
            {
                cout<<M.data[k].e<<" ";
                k++;
            }
            else
                cout<<0<<" ";
        }
        cout<<endl;
    }
}

int Add(TSMatrix *c,TSMatrix a,TSMatrix b)
{
    if (a.mu!=b.mu||a.nu!=b.nu)
        return 0;                               /* A,B行列数不相等 */
    int p,q;
    p=0; q=0;                                   /* p,q分别指向两个矩阵三元表下标 */
    c->mu = a.mu;
    c->nu = a.nu;
    c->tu = 0;
    while (p<a.tu&&q<b.tu)
    {
        if (a.data[p].i<b.data[q].i)        /* 情况1 */
            c->data[c->tu++]=a.data[p++];
        else if (a.data[p].i>b.data[q].i)   /* 情况2 */
            c->data[c->tu++]=b.data[q++];
        else                                    /* 3：行相同位置 */
        {
            if (a.data[p].j<b.data[q].j)        /* 情况3.1 */
                c->data[c->tu++]=a.data[p++];
            else if (a.data[p].j>b.data[q].j)   /* 情况3.2 */
                c->data[c->tu++]=b.data[q++];
            else                                    /* 列相同位置 */
            {
                int sum = a.data[p].e + b.data[q].e;
                if (sum)
                {
                    c->data[c->tu] = a.data[p];
                    c->data[c->tu].e = sum;
                    c->tu++;
                }
                p++; q++;
            }
        }
    }
    while (p < a.tu)  c->data[c->tu++] = a.data[p++];
    while (q < b.tu)  c->data[c->tu++] = b.data[q++];
    return 1;
}

第3关：基于十字链表的稀疏矩阵的存储与操作

//十字链表 
#include<iostream>
using namespace std;

#define MAX_SIZE 101
typedef int Elemtype;
typedef struct
{   int i;    //  行下标
    int j;        //  列下标
    Elemtype e;      // 元素值
}Triple;

typedef struct
{
    int  mu;         //   行数
    int  nu;         //   列数
    int  tu;         //   非0元素个数
    Triple  data[MAX_SIZE+1];   //data[0]不用
}TSMatrix;

void InitMatrix(TSMatrix *M)
{
    int m,n,t;
    cin>>m>>n>>t;
    M->mu=m;
    M->nu=n;
    M->tu=t;
    for(int k=1;k<=t;k++)
    {
        cin>>M->data[k].i>>M->data[k].j>>M->data[k].e;
    }
}

void Print(TSMatrix M)
{
    int flag1,flag2,flag=0;
    cin>>flag1>>flag2;
    if(flag1==1)
    {
        for(int k=1;k<=M.tu;k++)
        {
            if(M.data[k].i==flag2){
                cout<<M.data[k].e<<" ";
                flag=1;
            }
        }
        if(!flag)cout<<0;
    }
    else if(flag1==2)
    {
        for(int k=1;k<=M.tu;k++)
        {
            if(M.data[k].j==flag2){
                cout<<M.data[k].e<<" ";
                flag=1;
            }
        }
        if(!flag)cout<<0;
    }
}

int main()
{
    TSMatrix M;
    InitMatrix(&M);
    Print(M);
    return 0;
}