稀疏矩阵A、B均采用三元组顺序表表示,验证实现矩阵A快速转置算法,并设计、验证矩阵A、B相加得到矩阵C的算法。
(1) 从键盘输入矩阵的行数和列数,随机生成稀疏矩阵。
(2) 设计算法将随机生成的稀疏矩阵转换成三元组顺序表形式存储。
(3) 设计算法将快速转置得到的与相加得到的三元组顺序表分别转换成矩阵形式。
(4) 输出随机生成的稀疏矩阵A、B及其三元组顺序表、快速转置得到的与相加得到的三元组顺序表及其矩阵形式。
#include"stdio.h"
#include"stdlib.h"
#include"string.h"
#include"math.h"
#include"time.h"
#define TURE 1
#define FLASE 0
#define ERROR 0
#define OK 1
#define INFLASEBLE -1
#define MAXSIZE 2500
#define Size 50
typedef int ElemType;
typedef int Status;
typedef struct {
int i,j;
ElemType e;
} Triple;
typedef struct {
Triple data[MAXSIZE+1];
int mu,nu,tu;
} TSMatrix;
Status CreateTSMatrix(TSMatrix *M,int row,int col) {
int t,m,n,p=1;
srand((int)time(NULL));
M->mu=row;
M->nu=col;
for(m=0; m<M->mu; m++) {
for(n=0; n<M->nu; n++) {
t=rand()%( (int)sqrt(M->mu*M->nu) );
if(!t) {
if(rand()%2)
M->data[p].e=rand()%20;
else
M->data[p].e=rand()%20-20;
M->data[p].i=m+1;
M->data[p].j=n+1;
p++;
}
}
}
M->tu=p-1; //p从1始,刚好等于 非零元个数-1
return OK;
}
Status PrintTSMatrix(TSMatrix M) {
int i;
if(M.tu==0)
printf("无非零元!\n");
else {
printf("该矩阵的行数为%d、列数为%d、非零元素个数为%d",M.mu,M.nu,M.tu);
printf("\n非零元的坐标及值:\n\n行\t列\t元素值\n");
for(i=1; i<=M.tu; i++) {
printf("%-4d\t%-4d\t%-4d\n",M.data[i].i,M.data[i].j,M.data[i].e);
}
printf("\n");
}
return OK;
}
Status PrintMatrix(TSMatrix M) {
int i,j,k=0;
Triple *p=M.data;
p++;
printf("阵列形式为:\n");
for(i=1; i<=M.mu; i++) {
for(j=1; j<=M.nu; j++)
if (k<M.tu&&p->i==i&&p->j==j) {
printf("%4d",p->e);
p++, k++;
} else
printf("%4d",0);
printf("\n");
}
printf("\n");
return OK;
}
Status JudgeTSMatrix(TSMatrix M,int row,int col) {
int order=1; //data[0]未用
while(order<=M.tu) {
if(M.data[order].i==row&&M.data[order].j==col)
return order;
order++;
}
return 0;
}
Status TSMatrixAdd(TSMatrix M,TSMatrix N,TSMatrix *Q) {
int row,col;
int order=1,order1,order2; //data[0]未用
if(M.mu==N.mu&&M.nu==N.nu) {
for(row=1; row<=M.mu; row++) {
for(col=1; col<=M.nu; col++) {
order1=JudgeTSMatrix(M,row,col);
order2=JudgeTSMatrix(N,row,col);
if(order1&&order2) {
Q->data[order].i=row;
Q->data[order].j=col;
Q->data[order].e=M.data[order1].e+N.data[order2].e;
order++;
} else if(order1&&(!order2)) {
Q->data[order].e=M.data[order1].e;
Q->data[order].i=M.data[order1].i;
Q->data[order].j=M.data[order1].j;
order++;
} else if((!order1)&&order2) {
Q->data[order].e=N.data[order2].e;
Q->data[order].i=N.data[order2].i;
Q->data[order].j=N.data[order2].j;
order++;
}
}
}
Q->mu=M.mu;
Q->nu=M.nu;
Q->tu=order-1;
printf("矩阵 A 加矩阵 B 得到矩阵 C\n");
return OK;
} else {
printf("\n不是同型矩阵不能进行相加!\n");
return ERROR;
}
}
/*Status AddSmatrix(TSMatrix M,TSMatrix N,TSMatrix *Q)
{
if(M.mu!=N.mu || M.nu!=N.nu)return ERROR;
Q->mu=M.mu;
Q->nu=N.nu;
Q->tu=0;
int atu=1,btu=1;
while(M.tu>0 && N.tu>0)
{
if(M.data[atu].i<N.data[btu].i)
Q->data[++Q->tu]=M.data[atu++];
else if(M.data[atu].i>B.data[btu].i)
C->data[++C->tu]=B.data[btu++];
else
{
if(A.data[atu].j<B.data[btu].j)
C->data[++C->tu]=A.data[atu++];
else if(A.data[atu].j>B.data[btu].j)
C->data[++C->tu]=B.data[btu++];
else
{
if(A.data[atu].e+B.data[btu].e==0)
{
atu++;btu++;
}
C->data[++C->tu].e=B.data[btu++].e+A.data[atu++].e;
}
}
}
for(;atu<A.tu;atu++)
{C->data[++C->tu]=A.data[atu];}
for(;btu<B.tu;btu++)
{C->data[++C->tu]=B.data[btu];}
return OK;
}*/
Status FastTransMatrix(TSMatrix M,TSMatrix *N) {
int i,t,p,q;
int num[100],cpot[100];
N->mu=M.nu;
N->nu=M.mu;
N->tu=M.tu;
if(N->tu) {
for(i=1; i<=M.nu; i++)
num[i]=0;
for(t=1; t<=M.tu; t++)
num[M.data[t].j]++;
cpot[1]=1;
for(i=2; i<=M.nu; i++)
cpot[i]=cpot[i-1]+num[i-1];
for(p=1; p<=M.tu; ++p) {
i=M.data[p].j;
q=cpot[i];
N->data[q].i=M.data[p].j;
N->data[q].j=M.data[p].i;
N->data[q].e=M.data[p].e;
cpot[i]++;
}
}
return OK;
}
int main() {
while(1) {
TSMatrix A,A1,B,C;
int row,col;
printf("输入A\n");
printf("请输入行列\n");
scanf("%d%d",&row,&col);
CreateTSMatrix(&A,row,col);
PrintTSMatrix(A);
PrintMatrix(A);
printf("由矩阵 A 转置得矩阵 A1\n");
FastTransMatrix(A,&A1);
PrintTSMatrix(A1);
PrintMatrix(A1);
printf("矩阵 B:\n");
printf("请输入行列\n");
scanf("%d%d",&row,&col);
CreateTSMatrix(&B,row,col);
PrintTSMatrix(B);
PrintMatrix(B);
if(TSMatrixAdd(A,B,&C)) {
PrintTSMatrix(C);
PrintMatrix(C);
}
printf("\n");
}
return 0;
}