#include <stdio.h>
#include <stdlib.h>
#define MAXSIZE 100
typedef int DataType;
typedef struct //三元组类型定义
{
int i;
int j;
DataType e;
}Triple;
typedef struct
{
Triple data[MAXSIZE];
int m; //矩阵行数
int n; //列数
int len; //非零元素个数
}TriSeqMatrix;
int CreateMatrix(TriSeqMatrix *M)
{
int i,m,n;
DataType e;
int flag;
printf("输入矩阵行数,列数,非零元素个数:");
scanf("%d%d%d",&M->m,&M->n,&M->len);
if(M->len > MAXSIZE)
return 0;
for(i=0;i<M->len;i++)
{
do /*这个do-while 循环很诡异,可要可不要,除了flag有点价值外并卵用*/
{
printf("按行顺序输入第%d个非零元素所在的行(1~%d),列(1~%d),元素值:",i,M->m,M->n);
scanf("%d%d%d",&m,&n,&e);
flag = 0;
if(m < 0||m>M->n||n<0||n>M->n)
flag=1;
if(i>0&&m<M->data[i-1].i||m == M->data[i-1].i&&n<=M->data[i-1].j)
flag = 1;
}while(flag);
M->data[i].i=m;
M->data[i].j=n;
M->data[i].e=e;
}
return 1;
}
void DestroyMatrix(TriSeqMatrix *M) /*静态存储结构 */
{
M->m=M->n=M->len = 0;
}
void CopyMatrix(TriSeqMatrix M,TriSeqMatrix *N)
{
int i;
N->len=M.len;
N->m=M.m;
N->n=M.n;
for(i = 0;i < M.len;i++)
{
N->data[i].i = M.data[i].i;
N->data[i].j = M.data[i].j;
N->data[i].e = M.data[i].e;
}
}
int CompareElement(int a,int b)
{
if(a<b)
return -1;
else if(a == b)
return 0;
else
return 1;
}
/*稀疏矩阵的相加*/
int AddMatrix(TriSeqMatrix M,TriSeqMatrix N,TriSeqMatrix *Q)
{
int m = 0,n = 0,k = 0;
if(M.m!=N.m||M.n!=N.n)
return 0;
Q->m=M.m;
Q->n=M.n;
while(m<M.len && n<N.len)
{
switch(CompareElement(M.data[m].j,N.data[n].j))
{
case -1:
Q->data[k++]=M.data[m++];
break;
case 0:
Q->data[k++]=M.data[m++]; //Q->data[k++].e已加过
Q->data[k].e += N.data[n++].e; //k已加过
if(Q->data[k].e == 0) //若两数和为0,则不保存
k--;
break;
case 1:
Q->data[k++]=N.data[n++];
}
}
/*非零元素处理完毕,接下来处理稀疏矩阵的0元素*/
while(m<M.m) //若矩阵M元素未处理完毕,则直接赋给Q
Q->data[k++]=M.data[m++];
while(n<N.n)
Q->data[k++]=N.data[n++]; //同上
Q->len = k; //修改非零元素个数
if(k > MAXSIZE)
return 0;
return 1;
}
/*稀疏矩阵相减·一个机智的做法*/
int SubMatrix(TriSeqMatrix M,TriSeqMatrix N,TriSeqMatrix *Q)
{
int i;
for(i = 0;i<N.len;i++) //相减等于被加数*(-1)
N.data[i].e *= -1;
return AddMatrix(M,N,Q);
}
/*转置!!!第一种做法,第二种目前略*/
void TranMatrix(TriSeqMatrix M,TriSeqMatrix *N)
{
int i,k,col;
N->m=M.n; //对矩阵行,列进行操作
N->n=M.m;
N->len=M.len;
if(N->len) //矩阵三元组进行操作!
{
k = 0;
/*注意M.n的使用以及M.data[k].j的使用*/
for(col = 0;col<M.n;col++) //按照列号扫描
for(i=0;i<M.len;i++)
if(M.data[i].j == col)
{
N->data[k].i = M.data[i].j;
N->data[k].j = M.data[i].i;
N->data[k].e = M.data[i].e;
k++;
}
}
}
void PrintMatrix(TriSeqMatrix M)
{
int i;
printf("稀疏矩阵是%d行×%d列,共%d个非零元素.\n",M.m,M.n,M.len);
printf("行,列,元素值\n");
for(i = 0;i<M.len;i++)
printf("%2d,%6d,%8d\n",M.data[i].i,M.data[i].j,M.data[i].e);
}
int main()
{
TriSeqMatrix M,N,Q;
CreateMatrix(&M);
TranMatrix(M,&N);
AddMatrix(M,N,&Q);
return 0;
}
稀疏矩阵的基本操作-顺序三元组实现
最新推荐文章于 2024-01-24 05:45:00 发布