稀疏矩阵三元组乘法

#include
#include <stdlib.h>
using namespace std;
#define MaxSize 100
int M,N;
typedef struct
{
int r;
int c;
int d; //元素值
}TupNode; //三元组定义
typedef struct
{
int rows;
int cols;
int nums; //非零元素个数
TupNode data[MaxSize];
}TSMatrix; //三元组顺序表定义
/void CreatMat(TSMatrix &t,int a[100][100],int M,int N)
{
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++;
}
}
}
}
bool Value(TSMatrix &t,int x,int i,int j) //将i行j列元素替换为x 分以前就是0或者不是0
{
int k=0,k1;
if(i>=t.rows||j>t.cols)
return false;
while (k<t.nums && i>t.data[k].r) k++; //因为data[]里边存放的顺序就是i从小到大，根本不用考虑小于的情况
while (k<t.nums && it.data[k].r&&j>t.data[k].c) k++;
if(t.data[k].ri &&t.data[k].cj) //存在这样的元素 如果为零是不会存储到data中的；
t.data[k].d=x;
else
{
for(k1=t.nums-1;k1>=k;k1–)
{
t.data[k1+1].r=t.data[k1].r;
t.data[k1+1].c=t.data[k1].c;
t.data[k1+1].d=t.data[k1].d;
}
t.data[k].r=i;
t.data[k].c=j;
t.data[k].d=x;
t.nums++;
}
} /
void DispMat(TSMatrix t) //打印
{
int i;
if(t.nums<=0) return ;
for(i=0;i<t.nums;i++)
cout<<t.data[i].r<<’ ‘<<t.data[i].c<<’ '<<t.data[i].d<<endl;
}
/
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++) //按照0，1，2，3……n-1列进行转换
for(p=0;p<t.nums;p++)
if(t.data[p].cv) //v是列 从0-n-1
{
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++;
}
}
}
/
bool MutiPlyMatrix(TSMatrix t1,TSMatrix t2,TSMatrix &t3)
{
int count1=0,v=0;
int count2=0,t=0;;
if(t1.cols!=t2.rows)
{
cout<<“ERROR”<<endl;
return false;
}
for(int i=0;i<t1.nums;i++)
{
for(int j=0;j<t2.nums;j++)
if(t1.data[i].ct2.data[j].r)
{
int num=0;
for(int k=0;k<t3.nums;k++)
{
if(t3.data[k].rt1.data[i].r&&t3.data[k].c==t2.data[j].c) break;
num++; //如果没匹配到则 num=t3.nums;
}
if(num!=t3.nums)
{
t3.data[num].d=t3.data[num].d+t1.data[i].dt2.data[j].d;
}
if(num==t3.nums)
{
t3.data[count1].d=t1.data[i].dt2.data[j].d;
if(t3.data!=0)
{
t3.data[count1].r=t1.data[i].r;
t3.data[count1].c=t2.data[j].c;
t3.nums++;
count1++;
count2++;
}
}
}
}
if(count20)
{
cout<<“The answer is a Zero Matrix”<<endl;
t=1;
}
if(t0)
{DispMat(t3);}
return true;
}
int main()
{
TSMatrix t1,t2,t3;
int M1,N1,n1,M2,N2,n2;

cin>>M1>>N1>>n1;
t1.rows=M1;
t1.cols=N1;
t1.nums=0;
for(int i=0;i<n1;i++)
{
    int m1,n1,d1;
    cin>>m1>>n1>>d1;
    t1.data[t1.nums].r=m1;
    t1.data[t1.nums].c=n1;
    t1.data[t1.nums].d=d1;
    t1.nums++;
}
cin>>M2>>N2>>n2;
t2.rows=M2;
t2.cols=N2;
t2.nums=0;
for(int k=0;k<n2;k++)
{
    int m2,n2,d2;
    cin>>m2>>n2>>d2;
    t2.data[t2.nums].r=m2;
    t2.data[t2.nums].c=n2;
    t2.data[t2.nums].d=d2;
    t2.nums++;
   // cout<<t2.nums;
}
t3.nums=0;
//cout<<"the result is"<<endl;
MutiPlyMatrix(t1,t2,t3);
return 0;

}