i++; } else /*若triple_a的当前项的行号大于triple_b的当前项的行号*/ /*则将triple_b的项存入triple_c中*/ { triple_c.data[k].i = triple_b.data[j].i; triple_c.data[k].j = triple_b.data[j].j; triple_c.data[k].v = triple_b.data[j].v; k++; j++; } } triple_c.mu=triple_a.mu; /*产生第0行的结果*/ triple_c.nu=triple_a.nu; triple_c.tu=k-1; }
}/*TripleTableAdd*/
5.有两个稀疏矩阵的三元组triple_a和triple_b,请写出利用三元组这种数据结构,实现将两个稀疏矩阵相加,最后的和存入三元组表triple_c之中的算法。
#define MAXSIZE 1000 /*假设非零元个数的最大值是1000*/
typedef struct { int i, j; elemtype v; }triple;
typedef struct { triple data[MAXSIZE+1]; /*data[0]用于存放稀疏矩阵行,列和非零元个数*/ int mu, nu, tu; /*稀疏矩阵行、列和非零元的个数*/
} spmatrix;
int TripleValue(spmatrix c, int i, int j) /*该函数用在TripleMul函数之中,计算三元组c相应的值*/ { int k=1; while (k<=c.tu && (c.data[k].i!=i || c.data[k].j!=j)) k++; if (k<=c.tu) return(c.data[k].v); /*找到了返回该位置的值*/ else return(0); /*未找到说明该元素为0*/ }/*TripleValue*/
void TripleMul(int m,int n,int k,spmatrix a,spmatrix b,spmatrix c) { int i,j,h,p=1,s; for (i=0;i }/*TripleMul*/ 6.有两个稀疏矩阵的三元组triple_a和triple_b,请写出利用三元组这种数据结构,实现将两个稀疏矩阵相乘,并将最后的乘积存人二元组表triple_c之中的算法。