矩阵运算器.docVIP

矩阵运算器 #includeiostream #includevector using namespace std; templateclass T struct Triple { int r,c; T elem; }; templateclass T class SparseMatrix //是顺序表的推广 { vector TripleT triList; int rows,cols,num; public: SparseMatrix(); SparseMatrix(vector TripleT tlist,int rs,int cs,int n); void Simple_trans(SparseMatrix B); void Quick_trans(SparseMatrix B); //转置 void Plus(SparseMatrixT A,SparseMatrixT B); //加法(C=A+B) void Mult(SparseMatrixT A,SparseMatrixT B); //乘法(C=A*B) void PrintSparseMatrix(); //输出矩阵 }; templateclass T SparseMatrixT::SparseMatrix() //无参构造函数 { rows=0,cols=0,num=0; } templateclass T SparseMatrixT::SparseMatrix(vector TripleT tlist,int rs,int cs,int n) //有参构造函数 { rows=rs;cols=cs;num=n; for(int i=0;in;i++) triList.push_back(tlist[i]); } templateclass T void SparseMatrixT::Simple_trans(SparseMatrixT B) //朴素转置，B=A { B.rows=cols;B.cols=rows;B.num=num; B.triList.resize(num); //先将B的triList初始化 if(num==0) return; int q=0; for(int col=0;colcols;col++) for(int p=0;pnum;p++) if(triList[p].c==col) { B.triList[q].r=triList[p].c; B.triList[q].c=triList[p].r; B.triList[q].elem=triList[p].elem; q++; } } templateclass T void SparseMatrixT::Quick_trans(SparseMatrixT B) { B.rows=cols; B.cols=rows; B.num=num; if(num==0) return; B.triList.resize(num); int* cnum=new int[cols]; int* cpot=new int[cols]; for(int col=0;colcols;col++) cnum[col]=0; for(int i=0;inum;i++) cnum[triList[i].c]++; //计算每一列非零元个数 cpot[0]=0; for(col=1;colcols;col++) //每列第一个非零元在triList数组中的位置 cpot[col]=cpot[col-1]+cnum[col-1]; int p,q; for(p=0;pnum;p++) { col=triList[p].c; q=cpot[col]; B.triList[q].r=triL