稀疏矩阵加减乘以及三元组存储.docVIP

#includeiostream #includewindows.h using namespace std; #define N 4 int al[N][N],bl[N][N]; typedef int ElemType; #define MaxSize 100 struct TupNode { int r; int c; int d; }; struct TSMatrix { int rows; int cols; int nums; TupNode data[MaxSize]; }; void CreatMat(TSMatrix t,int A[N][N]) { int i,j; t.rows=N;t.cols=N;t.nums=0; for(i=0;iN;i++) { for(j=0;jN;j++) if(A[i][j]!=0) { t.data[t.nums].r=i; t.data[t.nums].c=j; t.data[t.nums].d=A[i][j]; t.nums++; } } } void DispMat(TSMatrix t) { int i; if(t.nums=0) return; cout行数，列数，非零元个数endl; coutt.rows t.cols t.nums; coutendl; cout三元组表为:endl; for(i=0;it.nums;i++) coutt.data[i].r t.data[i].c t.data[i].dendl; } int MatAdd(TSMatrix a,TSMatrix b,TSMatrix c) { int i=0,j=0,k=0; int v; if(a.rows!=b.rows||a.cols!=b.cols) return 0; c.rows=a.rows;c.cols=a.cols; while(ia.nums jb.nums) { if(a.data[i].r==b.data[j].r) { if(a.data[i].cb.data[j].c) { c.data[k].r=a.data[i].r; c.data[k].c=a.data[i].c; c.data[k].d=a.data[i].d; k++; i++; } else if (a.data[i].cb.data[j].c) { c.data[k].r=b.data[j].r; c.data[k].c=b.data[j].c; c.data[k].d=b.data[j].d; k++; j++; } else { v=a.data[i].d+b.data[j].d; if(v!=0) { c.data[k].r=a.data[i].r; c.data[k].c=a.data[j].c; c.data[k].d=v; k++; } i++; j++; } } else if(a.data[i].rb.data[j].r) { c.data[k].r=a.data[i].r; c.data[k].c=a.data[i].c; c.data[k].d=a.data[i].d; k++; i++; } else { c.data[k].r=b.data[j].r; c.data[k].c=b.data[j].c; c.data[k].d=b.data[j].d; k++; j++; } } c.nums=k; return 1; } int MatCut(TSMatrix a,TSMatrix b,TSMatrix c) { int i=0,j=0,k=0; int v; if(a.rows!=b.rows||a.cols!=b.cols) return 0; c.rows=a.rows;c.cols=a.cols; while(ia.nums jb.nums) { if(a.data[i].r==b.data[j].r) { if(a.data[i].cb.data[j].c) { c.data[k].r=a.data[i].r; c.data[k].c=a.data[i