代码:
#include
#include
#include
using namespace std;
#define M 4
#define N 4
#define MaxSize 100
typedef int ElemType;
typedef struct
{
int r;
int c;
ElemType d;///元素值
} TupNode; ///三元组定义
typedef struct
{
int rows;
int cols;
int nums;
TupNode data[MaxSize];
} TSMatrix; ///三元组顺序表定义
void CreatMat(TSMatrix &t,ElemType A[M][N])
{
t.rows=M;
t.cols=N;
t.nums=0;
for(int i=0; i
for(int j=0; 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,ElemType x,int i,int j)
{
int k=0,k1;
if(i>=t.rows||j>=t.cols)
return false;
while(kt.data[k].r)k++;
while(kt.data[k].c)k++;
if(t.data[k].r==i&&t.data[k].c==j)
t.data[k].d=x;
else
{
for(k1=t.nums-1; k1>=k; k1--)
{
t.data[k1+1].r=t.data[k].r;
t.data[k1+1].c=t.data[k].c;
t.data[k1+1].d=t.data[k].d;
}
t.data[k].r=i;
t.data[k].c=j;
t.data[k].d=x;
t.nums++;
}
return true;
}
bool Assign(TSMatrix t,ElemType &x,int i,int j)
{
int k=0;
if(i>=t.rows||j>=t.cols)
return false;
while(kt.data[k].r)k++;
while(kt.data[k].c)k++;
if(t.data[k].r==i&&t.data[k].c==j)
x=t.data[k].d;
else
x=0;
return true;
}
void DispMat(TSMatrix t)
{
if(t.nums<=0)
return ;
printf("\t%d\t%d\t%d\n",t.rows,t.cols,t.nums);
printf("\t-----------------\n");
for(int i=0; i
printf("\t%d\t%d\t%d\n",t.data[i].r,t.data[i].c,t.data[i].d);
}
void TranMat(TSMatrix t,TSMatrix &tb)
{
int i,j,k=0;
tb.rows=t.cols;
tb.cols=t.rows;
tb.nums=t.nums;
if(t.nums!=0)
{
for(i=0; i
for(j=0; j
if(t.data[j].c==i)
{
tb.data[k].r=t.data[j].c;
tb.data[k].c=t.data[j].r;
tb.data[k].d=t.data[j].d;
k++;
}
}
}
bool MatAdd(TSMatrix a,TSMatrix b,TSMatrix &c)
{
int i=0,j=0,k=0;
ElemType v;
if(a.rows!=b.rows||a.cols!=b.cols)
return false;
c.rows=a.rows;
c.cols=a.cols;
while(i
{
if(a.data[i].r==b.data[j].r)///先控制行相等
{
if(a.data[i].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].c>b.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[i].c;
c.data[k].d=v;
k++;
}
i++;
j++;
}
}
else if(a.data[i].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 true;
}
int getvalue(TSMatrix c,int i,int j)
{
int k=0;
while(k
k++;
if(k
return (c.data[k].d);
else
return (0);
}
bool MatMul(TSMatrix a,TSMatrix b,TSMatrix &c)
{
int i,j,k,p=0;
ElemType s;
if(a.cols!=b.rows)
return false;
for(i=0; i
for(j=0; j
{
s=0;
for(k=0; k
s+=getvalue(a,i,k)*getvalue(b,k,j);
if(s!=0)
{
c.data[p].r=i;
c.data[p].c=j;
c.data[p].d=s;
p++;
}
}
c.rows=a.rows;
c.cols=b.cols;
c.nums=p;
return true;
}
int main()
{
ElemType a1[N][N]={{1,0,3,0},{0,1,0,0},{0,0,1,0},{0,0,1,1}};
ElemType b1[M][M]={{3,0,0,0},{0,4,0,0},{0,0,1,0},{0,0,0,2}};
TSMatrix a,b,c;
CreatMat(a,a1);
CreatMat(b,b1);
printf("a的三元组:\n");
DispMat(a);
printf("b的三元组:\n");
DispMat(b);
printf("a转置为c\n");
TranMat(a,c);
printf("c的三元组\n");
DispMat(c);
printf("c=a+b\n");
MatAdd(a,b,c);
printf("c的三元组:\n");
DispMat(c);
printf("c=a*b\n");
MatMul(a,b,c);
printf("c的三元组:\n");
DispMat(c);
return 0;
}
总结
以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作具有一定的参考学习价值,谢谢大家对脚本之家的支持。如果你想了解更多相关内容请查看下面相关链接
您可能感兴趣的文章:Dijkstra算法最短路径的C++实现与输出路径
C++项目求Fibonacci数列的参考解答
C++实践IP地址类项目参考
C++实践数组作数据成员的参考
一张图总结C++中关于指针的那些事
C++实践数组类运算的实现参考
C++实践Time类中的运算符重载参考方法
C++实践分数类中运算符重载的方法参考
C++实践排序函数模板项目的参考方法
对称矩阵的压缩储存讲解