假设稀疏矩阵A和B均以三元组表作为存储结构。试写出矩阵相加的算法，另设三元组表C存放结果矩阵。

稀疏矩阵的三元组顺序表类型TSMatrix的定义：

#define MAXSIZE 20 // 非零元个数的最大值

typedef struct {

int i,j; // 行下标,列下标

ElemType e; // 非零元素值

}Triple;

typedef struct {

Triple data[MAXSIZE+1]; // 非零元三元组表,data[0]未用

int mu,nu,tu; // 矩阵的行数、列数和非零元个数

}TSMatrix;实现函数如下：

Status AddTSM(TSMatrix A,TSMatrix B,TSMatrix &C)

/* 三元组表示的稀疏矩阵加法: C=A+B */

{

int ai,bi,ci,aj,bj,cj,ap,bp,cp;

ap = bp = cp = 1;

if(A.mu != B.mu || A.nu != B.nu){

return ERROR;

}

C.mu = A.mu;

C.nu = A.nu;

while(ap <= A.tu && bp <= B.tu){

ai = A.data[ap].i;

bi = B.data[bp].i;

if(ai > bi){

ci = bi;

while(ci == B.data[bp].i){

C.data[cp].i = ci;

C.data[cp].j = B.data[bp].j;

C.data[cp].e = B.data[bp].e;

++bp;

++cp;

}

}else if(ai < bi){

ci = ai;

while(ci == A.data[ap].i){

C.data[cp].i = ci;

C.data[cp].j = A.data[ap].j;

C.data[cp].e = A.data[ap].e;

++ap;

++cp;

}

}else if(ai == bi){

ci = ai;

aj = A.data[ap].j;

bj = B.data[bp].j;

if(aj > bj){

C.data[cp].i = ci;

C.data[cp].j = bj;

C.data[cp].e = B.data[bp].e;

++cp;

++bp;

}else if(aj < bj){

C.data[cp].i = ci;

C.data[cp].j = aj;

C.data[cp].e = A.data[ap].e;

++cp;

++ap;

}else if(aj == bj){

if(A.data[ap].e + B.data[bp].e != 0){

C.data[cp].i = ci;

C.data[cp].j = aj;

C.data[cp].e = A.data[ap].e + B.data[bp].e;

++cp;

}

++ap;

++bp;

}

}

}

//以上为稀疏矩阵A或B中的元素完全加完的情况

//以下为稀疏矩阵A或B中的元素部分剩余的情况

while(ap <= A.tu){

C.data[cp].i = A.data[ap].i;

C.data[cp].j = A.data[ap].j;

C.data[cp].e = A.data[ap].e;

++cp;

++ap;

}

while(bp <= B.tu){

C.data[cp].i = B.data[bp].i;

C.data[cp].j = B.data[bp].j;

C.data[cp].e = B.data[bp].e;

++cp;

++bp;

}

C.tu = --cp;

return OK;

}