稀疏矩阵的基本操作-顺序三元组实现

#include <stdio.h>
#include <stdlib.h>

#define MAXSIZE 100

typedef int DataType;

typedef struct         //三元组类型定义
{
    int i;
    int j;
    DataType e;
}Triple;

typedef struct
{
    Triple data[MAXSIZE];
    int m;            //矩阵行数
    int n;            //列数
    int len;          //非零元素个数
}TriSeqMatrix;

int CreateMatrix(TriSeqMatrix *M)
{
    int i,m,n;
    DataType e;
    int flag;
    printf("输入矩阵行数，列数，非零元素个数：");
    scanf("%d%d%d",&M->m,&M->n,&M->len);
    if(M->len > MAXSIZE)
        return 0;
    for(i=0;i<M->len;i++)
    {
        do                /*这个do-while 循环很诡异,可要可不要，除了flag有点价值外并卵用*/
        {
            printf("按行顺序输入第%d个非零元素所在的行(1~%d)，列(1~%d)，元素值：",i,M->m,M->n);
            scanf("%d%d%d",&m,&n,&e);
            flag = 0;
            if(m < 0||m>M->n||n<0||n>M->n)
                flag=1;
            if(i>0&&m<M->data[i-1].i||m == M->data[i-1].i&&n<=M->data[i-1].j)
                flag = 1;
        }while(flag);
        M->data[i].i=m;
        M->data[i].j=n;
        M->data[i].e=e;
    }
    return 1;
}

void DestroyMatrix(TriSeqMatrix *M)  /*静态存储结构 */
{
    M->m=M->n=M->len = 0;
}

void CopyMatrix(TriSeqMatrix M,TriSeqMatrix *N)
{
    int i;
    N->len=M.len;
    N->m=M.m;
    N->n=M.n;
    for(i = 0;i < M.len;i++)
    {
        N->data[i].i = M.data[i].i;
        N->data[i].j = M.data[i].j;
        N->data[i].e = M.data[i].e;
    }
}

int CompareElement(int a,int b)
{
    if(a<b)
        return -1;
    else if(a == b)
        return 0;
    else
        return 1;
}

/*稀疏矩阵的相加*/
int AddMatrix(TriSeqMatrix M,TriSeqMatrix N,TriSeqMatrix *Q)
{
    int m = 0,n = 0,k = 0;
    if(M.m!=N.m||M.n!=N.n)
        return 0;
    Q->m=M.m;
    Q->n=M.n;
    while(m<M.len && n<N.len)
    {
        switch(CompareElement(M.data[m].j,N.data[n].j))
        {
        case -1:
            Q->data[k++]=M.data[m++];
            break;
        case 0:
            Q->data[k++]=M.data[m++];         //Q->data[k++].e已加过
            Q->data[k].e += N.data[n++].e;   //k已加过
            if(Q->data[k].e == 0)            //若两数和为0，则不保存
                k--;
            break;
        case 1:
            Q->data[k++]=N.data[n++];
        }
    }
    /*非零元素处理完毕，接下来处理稀疏矩阵的0元素*/
    while(m<M.m)                        //若矩阵M元素未处理完毕，则直接赋给Q
        Q->data[k++]=M.data[m++];
    while(n<N.n)
        Q->data[k++]=N.data[n++];       //同上
    Q->len = k;                         //修改非零元素个数
    if(k > MAXSIZE)
        return 0;
    return 1;
}

/*稀疏矩阵相减·一个机智的做法*/
int SubMatrix(TriSeqMatrix M,TriSeqMatrix N,TriSeqMatrix *Q)
{
    int i;
    for(i = 0;i<N.len;i++)              //相减等于被加数*(-1)
        N.data[i].e *= -1;
    return AddMatrix(M,N,Q);
}

/*转置！！！第一种做法，第二种目前略*/
void TranMatrix(TriSeqMatrix M,TriSeqMatrix *N)
{
    int i,k,col;
    N->m=M.n;                           //对矩阵行，列进行操作
    N->n=M.m;
    N->len=M.len;
    if(N->len)                          //矩阵三元组进行操作！
    {
        k = 0;
        /*注意M.n的使用以及M.data[k].j的使用*/
        for(col = 0;col<M.n;col++)      //按照列号扫描
            for(i=0;i<M.len;i++)
                if(M.data[i].j == col)
        {
            N->data[k].i = M.data[i].j;
            N->data[k].j = M.data[i].i;
            N->data[k].e = M.data[i].e;
            k++;
        }
    }
}

void PrintMatrix(TriSeqMatrix M)
{
    int i;
    printf("稀疏矩阵是%d行×%d列，共%d个非零元素.\n",M.m,M.n,M.len);
    printf("行，列，元素值\n");
    for(i = 0;i<M.len;i++)
        printf("%2d,%6d,%8d\n",M.data[i].i,M.data[i].j,M.data[i].e);
}

int main()
{
    TriSeqMatrix M,N,Q;
    CreateMatrix(&M);
    TranMatrix(M,&N);
    AddMatrix(M,N,&Q);
    return 0;
}