求一个C++的稀疏矩阵类，要求能实现加法，减法和乘法，正常输出。还能转置。。。急！！！急！！！

#include
using namespace std;
const int MAXSIZE=100; // 定义非零元素的最多个数
const int MAXROW=10; // 定义数组的行数的最大值
typedef struct { // 定义三元组的元素
int i,j; //行数 列数
int e; //非零元的值
}Triple;
typedef struct { // 定义普通三元组对象
Triple data[MAXSIZE+1];
int a,b,c; //行数，列数，非零元个数
}TSMatrix;
typedef struct { // 定义带链接信息的三元组对象
Triple data[MAXSIZE+2];
int rpos[MAXROW+1];
int a,b,c;
}RLSMatrix;
template
int shuru(P & T,int y){ //输入矩阵，按三元组格式输入
cout<<"输入稀疏矩阵的行，列和非零元素个数:"<cin>>T.a>>T.b>>T.c;
cout<<"请输出非零元素的位置和值:"<int k=1;
for(;k<=T.c;k++)
cin>>T.data[k].i>>T.data[k].j>>T.data[k].e;
return 1;
}
template
void shuchu(P T){ // 输出矩阵，按标准格式输出
int m,n,k=1;
for(m=0;mfor(n=0;nif((T.data[k].i-1)==m&&(T.data[k].j-1)==n){
cout.width(4);
cout<else{
cout.width(4); cout<<"0"; }
}
cout<}
}
void zhuanzhi( ) // 求矩阵的转置矩阵
{
TSMatrix M,T; //定义预转置的矩阵
shuru(M, 0); //输入矩阵
int num[MAXROW+1];
int cpot[MAXROW+1]; // 构建辅助数组
int q,p,t;
T.c=M.c; T.a=M.b; T.b=M.a;
if(T.c){
for(int col=1;col<=M.c;col++) num[col]=0;
for(t=1;t<=M.c;t++) ++num[M.data[t].j];
cpot[1]=1;
for(int i=2;i<=M.b;i++) cpot[i]=cpot[i-1]+num[i-1]; // 求出每一列中非零元素在三元组中出现的位置
for(p=1;p<=M.c;p++){
col=M.data[p].j; q=cpot[col];
T.data[q].i=col; T.data[q].j=M.data[p].i;
T.data[q].e=M.data[p].e; ++cpot[col];
}
}
cout<<"输入稀疏矩阵的转置矩阵为"<shuchu(T);
}
void Count(RLSMatrix &T) // 求取每一行中非零元素在三元组中出现的位置
{
int num[MAXROW+1];
for(int col=1;col<=T.a;col++) num[col]=0;
for(col=1;col<=T.c;col++) ++num[T.data[col].i];
T.rpos[1]=1;
for(int i=2;i<=T.a;i++) T.rpos[i]=T.rpos[i-1]+num[i-1]; }
typedef struct OLNode // 定义十字链表元素
{
int i,j; //该非零元的行列下标
int e; //非零元值
struct OLNode *right,*down; // 该非零元所在行表和列表的后继元素
}OLNode,*OLink; typedef struct // 定义十字链表对象结构体
{
OLink *rhead,*chead; //
int a,b,c; // 系数矩阵的行数，列数，和非零元素个数
}CrossList; void CreateSMatrix_OL(CrossList & M) // 创建十字链表
{
int x,y,m;
cout<<"请输入稀疏矩阵的行，列，及非零元素个数"<cin>>M.a>>M.b>>M.c;
if(!(M.rhead=(OLink*)malloc((M.a+1)*sizeof(OLink)))) exit(0);
if(!(M.chead=(OLink*)malloc((M.b+1)*sizeof(OLink)))) exit(0);
for(x=0;x<=M.a;x++)
M.rhead[x]=NULL; // 初始化各行，列头指针，分别为NULL
for(x=0;x<=M.b;x++)
M.chead[x]=NULL;
cout<<"请按三元组的格式输入数组："<for(int i=1;i<=M.c;i++){
cin>>x>>y>>m; // 按任意顺序输入非零元，（普通三元组形式输入）
OLink p,q;
if(!(p=(OLink)malloc(sizeof(OLNode)))) exit(0); // 开辟新节点，用来存储输入的新元素
p->i=x; p->j=y; p->e=m;
if(M.rhead[x]==NULL||M.rhead[x]->j>y){
p->right=M.rhead[x]; M.rhead[x]=p;
}
else{
for(q=M.rhead[x];(q->right)&&(q->right->jright); // 查找节点在行表中的插入位置
p->right=q->right; q->right=p; // 完成行插入
}
if(M.chead[y]==NULL||M.chead[y]->i>x){
p->down=M.chead[y]; M.chead[y]=p;
}
else{
for(q=M.chead[y];(q->down)&&(q->down->idown); // 查找节点在列表中的插入位置
p->down=q->down; q->down=p; // 完成列插入 }
}
} void shuchu(CrossList T){ // 输出十字链表，用普通数组形式输出
for(int i=1;i<=T.a;i++){
OLink p=T.rhead[i];
for(int j=1;j<=T.b;j++){
if((p)&&(j==p->j)){
cout<e; p=p->right;
}
else
cout<<"0";
}
cout<}
}
void chengfa() //矩阵的乘法
{
CrossList M,N; // 创建两个十字链表对象，并初始化
CreateSMatrix_OL(M);
CreateSMatrix_OL(N);
cout<<"输入的两稀疏矩阵的积为："<OLink pa,pb,pre ,hl[MAXROW+1]; //定义辅助指针，pa，pb分别为M,N当前比较的元素，pre为pa的前驱元素
for(int x=1;x<=M.b;x++) hl[x]=M.chead[x];
for(int k=1;k<=M.a;k++){ // 对M的每一行进行操作
pa=M.rhead[k]; pb=N.rhead[k]; pre=NULL;
while(pb){ // 把N中此行的每个元素取出，
OLink p;
if(!(p=(OLink)malloc(sizeof(OLNode)))) exit(0); // 开辟新节点，存储N中取出的元素
p->e=pb->e; p->i=pb->i; p->j=pb->j;
if(NULL==pa||pa->j>pb->j){ // 当M此行已经检查完或者pb因该放在pa前面
if(NULL==pre)
M.rhead[p->i]=p;
else
pre->right=p;
p->right=pa; pre=p;
if(NULL==M.chead[p->j]){ // 进行列插入
M.chead[p->j]=p; p->down=NULL;
}
else{
p->down=hl[p->j]->down; hl[p->j]->down=p;
}
hl[p->j]=p;
pb=pb->right;
}
else
if((NULL!=pa)&&pa->jj){ // 如果此时的pb元素因该放在pa后面，则取以后的pa再来比较
pre=pa; pa=pa->right;
}
else
if(pa->j==pb->j){ // 如果pa，pb位于同一个位置上，则将值相加
pa->e *= pb->e;
if(!pa->e){ // 如果相加后的和为0，则删除此节点，同时改变此元素坐在行，列的前驱元素的相应值
if(NULL==pre) // 修改行前驱元素值
M.rhead[pa->i]=pa->right;
else
pre->right=pa->right;
p=pa; pa=pa->right;
if(M.chead[p->j]==p) M.chead[p->j]=hl[p->j]=p->down; // 修改列前驱元素值
else
hl[p->j]->down=p->down;
free(p); pb=pb->right;
}
else{
pa=pa->right; pb=pb->right;
}
}
}
}
shuchu(M);
}
void jianfa() //矩阵的减法
{
CrossList M,N; // 创建两个十字链表对象，并初始化
CreateSMatrix_OL(M);
CreateSMatrix_OL(N);
cout<<"输入的两稀疏矩阵的差矩阵为："<OLink pa,pb,pre ,hl[MAXROW+1]; //定义辅助指针，pa，pb分别为M,N当前比较的元素，pre为pa的前驱元素
for(int x=1;x<=M.b;x++) hl[x]=M.chead[x];
for(int k=1;k<=M.a;k++){ // 对M的每一行进行操作
pa=M.rhead[k]; pb=N.rhead[k]; pre=NULL;
while(pb){ // 把N中此行的每个元素取出，
OLink p;
if(!(p=(OLink)malloc(sizeof(OLNode)))) exit(0); // 开辟新节点，存储N中取出的元素
p->e=pb->e; p->i=pb->i; p->j=pb->j;
if(NULL==pa||pa->j>pb->j){ // 当M此行已经检查完或者pb因该放在pa前面
if(NULL==pre)
M.rhead[p->i]=p;
else
pre->right=p;
p->right=pa; pre=p;
if(NULL==M.chead[p->j]){ // 进行列插入
M.chead[p->j]=p; p->down=NULL;
}
else{
p->down=hl[p->j]->down; hl[p->j]->down=p;
}
hl[p->j]=p;
pb=pb->right;
}
else
if((NULL!=pa)&&pa->jj){ // 如果此时的pb元素因该放在pa后面，则取以后的pa再来比较
pre=pa; pa=pa->right;
}
else
if(pa->j==pb->j){ // 如果pa，pb位于同一个位置上，则将值相加
pa->e -= pb->e;
if(!pa->e){ // 如果相加后的和为0，则删除此节点，同时改变此元素坐在行，列的前驱元素的相应值
if(NULL==pre) // 修改行前驱元素值
M.rhead[pa->i]=pa->right;
else
pre->right=pa->right;
p=pa; pa=pa->right;
if(M.chead[p->j]==p) M.chead[p->j]=hl[p->j]=p->down; // 修改列前驱元素值
else
hl[p->j]->down=p->down;
free(p); pb=pb->right;
}
else{
pa=pa->right; pb=pb->right;
}
}
}
}
shuchu(M);
} void add() //矩阵的加法
{
CrossList M,N; // 创建两个十字链表对象，并初始化
CreateSMatrix_OL(M);
CreateSMatrix_OL(N);
cout<<"输入的两稀疏矩阵的和矩阵为："<OLink pa,pb,pre ,hl[MAXROW+1]; //定义辅助指针，pa，pb分别为M,N当前比较的元素，pre为pa的前驱元素
for(int x=1;x<=M.b;x++) hl[x]=M.chead[x];
for(int k=1;k<=M.a;k++){ // 对M的每一行进行操作
pa=M.rhead[k]; pb=N.rhead[k]; pre=NULL;
while(pb){ // 把N中此行的每个元素取出，
OLink p;
if(!(p=(OLink)malloc(sizeof(OLNode)))) exit(0); // 开辟新节点，存储N中取出的元素
p->e=pb->e; p->i=pb->i; p->j=pb->j;
if(NULL==pa||pa->j>pb->j){ // 当M此行已经检查完或者pb因该放在pa前面
if(NULL==pre)
M.rhead[p->i]=p;
else
pre->right=p;
p->right=pa; pre=p;
if(NULL==M.chead[p->j]){ // 进行列插入
M.chead[p->j]=p; p->down=NULL;
}
else{
p->down=hl[p->j]->down; hl[p->j]->down=p;
}
hl[p->j]=p;
pb=pb->right;
}
else
if((NULL!=pa)&&pa->jj){ // 如果此时的pb元素因该放在pa后面，则取以后的pa再来比较
pre=pa; pa=pa->right;
}
else
if(pa->j==pb->j){ // 如果pa，pb位于同一个位置上，则将值相加
pa->e += pb->e;
if(!pa->e){ // 如果相加后的和为0，则删除此节点，同时改变此元素坐在行，列的前驱元素的相应值
if(NULL==pre) // 修改行前驱元素值
M.rhead[pa->i]=pa->right;
else
pre->right=pa->right;
p=pa; pa=pa->right;
if(M.chead[p->j]==p) M.chead[p->j]=hl[p->j]=p->down; // 修改列前驱元素值
else
hl[p->j]->down=p->down;
free(p); pb=pb->right;
}
else{
pa=pa->right; pb=pb->right;
}
}
}
}
shuchu(M);
} int main() //主函数
{
cout<<"1:稀疏矩阵的加法。"<cout<<"2:稀疏矩阵的乘法。"<cout<<"3:稀疏矩阵的减法。"<cout<<"4:稀疏矩阵的转置."<cout<<"0:退出程序。"<char c=getchar();
if(c=='1')
add(); //调用矩阵相加函数
else
if(c=='2')
chengfa(); //调用矩阵相乘函数
else
if(c=='3')
jianfa(); //调用矩阵相减函数
else
if(c=='4')
zhuanzhi( ); //调用矩阵转置函数
else
exit(0); //退出
return 0;
}