矩阵三元组的运算.docVIP

矩阵三元组的运算 //矩阵三元组之矩阵相加 相乘 #include iostream using namespace std; typedef int Elemtype; #define MAXSIZE 12500 //最大非零元素 typedef struct Triple { Elemtype value; int row,col; }Triple; typedef struct TSMatrix { Triple data[MAXSIZE+1]; int mu,nu,tu; }TSMatrix; TSMatrix T; void InputMatrix(TSMatrix T) //输入t个非零元素 { cout请输入稀疏矩阵的信息,(行，列，非零元素个数)endl; cinT.muT.nuT.tu; int i; cout请输入非零元素的信息(行，列，值),提醒(下标从1开始)endl; for(i=1;i=T.tu;++i) { cinT.data[i].rowT.data[i].colT.data[i].value; } } void Output(TSMatrix T) { cout矩阵的三元组表示(ROW=)T.mu COL=T.nu非零个数=T.tuendl; int i; for(i=1;i=T.tu;++i) { coutROW(行):T.data[i].row COL(列):T.data[i].col Value(值)T.data[i].valueendl; } } void TransposeSMatrix(TSMatrix M,TSMatrix T) //矩阵的转置 { T.mu=M.nu;T.nu=M.mu;T.tu=M.tu; int i,j,k=1; for(i=1;i=M.nu;++i) { for(j=1;j=M.tu;++j) if(M.data[j].col==i) { T.data[k].row=i; T.data[k].col=M.data[j].row; T.data[k].value=M.data[j].value; ++k; } } } void AddMastrix(TSMatrix M,TSMatrix T,TSMatrix Q) //矩阵相加 { int index_a,index_b,i=1,j=1,k=1; Q.mu=M.mu; Q.nu=M.nu; while (i=M.tuj=T.tu) { index_a=(M.data[i].row)*(M.data[i].col)+M.data[i].col; index_b=(T.data[j].row)*(T.data[j].col)+T.data[j].col; if(index_aindex_b) { Q.data[k]=M.data[i]; i++; k++; } else if(index_aindex_b) { Q.data[k]=T.data[j]; j++; k++; } else if(index_a==index_b) { if((M.data[i].value+T.data[j].value)!=0) { Q.data[k]=M.data[i]; Q.data[k].value=M.data[i].value+T.data[j].value; k++; } ++i; ++j; } } //复制剩余元素 for(;i=M.tu;++i) { Q.data[k]=M.data[i]; k++; } for(;j=T.tu;++j) Q.data[k++]=T.data[j]; Q.tu=k-1; } void Multiply(TSMatrix M,TSMatrix T,TSMatrix Q) //相乘 { if(M.nu!=T.mu) { cerr两矩阵相乘不合法endl; return ; } int *rowSize=new int[T.mu+1]; //存放每行非零元素的个数 int *rowStart=new int[T.mu+2]; //矩阵每行在三元组开始位置 int *temp=new int[T.nu+1]; //存放结果矩阵中每行的计算结果 int i,Current,k,ROWM,COLM,COLB; for(i=1;i=T.mu;i++) rowSize[i]=0; for(i=