第五章稀疏矩阵三元组表示

}
//按矩阵形式输出M
void PrintTSMatrix1(TSMatrix M)
{
int i,j,k=1;
Triple *p=M.data;
p++; //p指向第一个非零元素
for(i=1;i<=M.mu;i++)
{
for(j=1;j<=M.nu;j++)
if(k<=M.tu&&p->i==i&&p->j==j) //p指向非零元，且p所指元素为当前处理元素
break;
case 0:
Q.data[++q]=M.data[m++]; //M、N矩阵当前非零元素的行列均相等
Q.data[q].e+=N.data[n++].e;//矩阵M、N的当前元素值求和并赋给矩阵Q
if(Q.data[q].e==0) //元素值为0，不存入压缩矩阵
q--;
break;
case 1: Q.data[++q]=N.data[n++];
for(j=1;j<=M.tu;j++) //对于M的每一个值
Tc[M.data[j].i]+=M.data[j].e*Nc[M.data[j].j];//Tc中存N的第i列与M相乘的结果
for(j=1;j<=M.mu;j++)
if(Tc[j]!=0)
{
T.data[++T.tu].e=Tc[j];
Tc=(ElemType*)malloc((M.nu+1)*sizeof(ElemType));//Tc为矩阵T一行的临时数组（非压缩，【0】不用）
if(!Nc||!Tc) //创建临时数组不成功
exit(ERROR);
for(i=1;i<=N.nu;i++) //对于N的每一列
{
for(j=1;j<=N.mu;j++)
}
//3.输出稀疏矩阵M
void PrintTSMatrix(TSMatrix M)
{
int i;
printf("%d行%d列%d个非零元素\n",M.mu,M.nu,M.tu);
printf("行列元素值\n");
for(i=1;i<=M.tu;i++)
printf("%2d%4d%8d\n",M.data[i].i,M.data[i].j,M.data[i].e);
int p, q, col;
T.mu = M.nu; T.nu = M.mu; T.tu = M.tu;
if (T.tu)
{
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;,*Tc;
TSMatrix T; //临时矩阵
if(M.nu!=N.mu)
return ERROR;
T.nu=M.mu; //临时矩阵T是Q的转置矩阵
T.mu=N.nu;
T.tu=0;
Nc=(ElemType*)malloc((N.mu+1)*sizeof(ElemType));//Nc为矩阵N一列的临时数组（非压缩，【0】不用）
printf("矩阵C3(A的转置)：\n");
PrintTSMatrix1(C);
printf("创建矩阵A2：");
CreateTSMatrix(A);
PrintTSMatrix1(A);
printf("创建矩阵B2：（行数应与矩阵A的列数相同=%d）\n",A.nu);
CreateTSMatrix(B);
T.data[T.tu].i=i;
T.data[T.tu].j=j;
}
}
if(T.tu>MAX_SIZE) //非零元素个数太多
return ERROR;
TransposeTSMatrix(T,Q); //将T的装置赋给Q
DestroyTSMatrix(T); //销毁临时矩阵T
free(Tc); //释放动态数组Tc和Nc
{
case -1:
Q.data[++q]=M.data[m++];//将矩阵M的当前元素值赋给矩阵Q
break;
case 0:
switch(comp(M.data[m].j,N.data[n].j)) //M、N矩阵当前元素的行相等，继续比较列
{
case -1:
Q.data[++q]=M.data[m++];
printf("请输入矩阵的行数，列数，非零元素数：");
scanf("%d%d%d",&M.mu,&M.nu,&M.tu);
if(M.tu>MAX_SIZE)
return ERROR;
M.data[0].i=0;//为以下比较顺序做准备
for(i=1;i<=M.tu;i++) //注意i从1开始的
Status TransposeTSMatrix(TSMatrix M, TSMatrix &T);
//主函数
void main()
{
TSMatrix A,B,C;
printf("创建矩阵A：\n");
CreateTSMatrix(A);
PrintTSMatrix(A);
printf("由矩阵A复制矩阵B：\n");
}
break;
case 1: Q.data[++q]=N.data[n++];//将矩阵N的当前元素值赋给矩阵Q
}
}
while(m<=M.tu)//矩阵N的元素全部处理完毕
Q.data[++q]=M.data[m++];
while(n<=N.tu)//矩阵M的元素全部处理完毕
Q.data[++q]=N.data[n++];
{
T=M;
}
//5.求稀疏矩阵的和Q=M+N，M和N的行数和列数相等
//AddTSMatrix()中要用到
int comp(int c1,int c2)
{
if(c1<c2)
return -1;
if(c1==c2)
return 0;
return 1;
}
Status AddTSMatrix(TSMatrix M,TSMatrix N,TSMatrix &Q)
}TSMatrix;
//函数声明
Status CreateTSMatrix(TSMatrix &M);
void DestroyTSMatrix(TSMatrix &M);
void PrintTSMatrix(TSMatrix M);
void PrintTSMatrix1(TSMatrix M);
void CopyTSMatrix(TSMatrix M,TSMatrix &T);
注意：
程序：加法、乘法不调用。
#include<stdio.h>
#include<stdlib.h> //稀疏矩阵的三元组顺序表存储表示
#define MAX_SIZE 12500 //假设非零元个数的最大值为12500
#define OK 1
#define ERROR 0
typedef int ElemType;
free(Nc);
return OK;
}
//8.求稀疏矩阵M的转置矩阵T TransposeTSMatrix(TSMatrix M,TSMatrix &T)
Status TransposeTSMatrix(TSMatrix M, TSMatrix &T)
{ //算法5.1
//采用三元组顺序表存储表示，求稀疏矩阵M的转置矩阵T
{
printf("%3d",p->e);//输出p所指元素的值
p++;//p指向下一个元素
k++;//计数器+1
}
else //p所指元素不是当前处理元素
printf("%3d",0);//输出0
printf("\n");
}
}
//4.由稀疏矩阵M复制得到T
void CopyTSMatrix(TSMatrix M,TSMatrix &T)
Q.tu=q;//矩阵Q的非零元素个数
if(q>MAX_SIZE)//非零元素个数太多
return ERROR;
return OK;
}
//6.求稀疏矩阵的差Q=M-N，M和N的行数和列数相等
Status SubtTSMatrix(TSMatrix M,TSMatrix N,TSMatrix &Q)
{
PrintTSMatrix1(B);
//MultTSMatrix(A,B,C);
printf("矩阵C5(A*B)：\n");