假设nn的稀疏矩阵A采用三元组表示，设计一个程序exp6-4.cpp实现如下功能：

（1）生成如下两个稀疏矩阵的三元组a和b；

（2）输出a转置矩阵的三元组；
（3）输出a+b的三元组；
（4）输出ab的三元组。
#include <iostream>
#include<cstdio>
using namespace std;
#define M 4
#define N 4
typedef struct
{
int r;
int c;
int d;
} TupNode;
typedef struct
{
int rows;
int cols;
int nums;
TupNode data[1000];
} TSMatrix;
void CreatMat(TSMatrix &t,int A[M]
)
{
int i,j;
t.rows=M;
t.cols=N;
t.nums=0;
for(i=0; i<M; i++)
for(j=0; j<N; 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 ;
printf("\t%d\t%d\t%d\n",t.rows,t.cols,t.nums);
printf("\t--------------------\n");
for(i=0; i<t.nums; i++)
printf("\t%d\t%d\t%d\n",t.data[i].r,t.data[i].c,t.data[i].d);
}
void TranTat(TSMatrix t,TSMatrix &tb)
{
int p,q=0,v;
tb.rows=t.cols;
tb.cols=t.rows;
tb.nums=t.nums;
if(t.nums!=0)
{
for(v=0; v<t.cols; v++)
for(p=0; p<t.nums; p++)
{
tb.data[q].r=t.data[p].c;
tb.data[q].c=t.data[p].r;
tb.data[q].d=t.data[p].d;
q++;
}
}
}
int main()
{
TSMatrix t,t1,tb;
int c[4][4];
int a[M]
= {{1,0,3,0},{0,1,0,0},{0,0,1,0},{0,0,1,1}};
CreatMat(t,a);
TranTat(t,tb);
printf("a的三元组：\n");
DispMat(t);
printf("a的转置三元组：\n");
DispMat(tb);
int b[M]
= {{3,0,0,0},{0,4,0,0},{0,0,1,0},{0,0,0,2}};
CreatMat(t1,b);
printf("b的三元组：\n");
DispMat(t1);
for (int i = 0; i < M; i++)
for (int j = 0; j < N; j++)

4000
c[i][j] = a[i][j] + b[i][j];
CreatMat(t,c);
printf("a+b的和：\n");
DispMat(t);
int sum;
for (int i = 0; i < M; i++)
{
for (int j = 0; j < N; j++)
{
sum = 0;
for (int k = 0; k < N; k++)
{
sum += a[i][k] * b[k][j];
}
c[i][j]=sum;
}
}
CreatMat(t,c);
printf("a*b的三元组：\n");
DispMat(t);
return 0;
}