第3题答案
【参考答案】
#define MAXVSIZE 非零元个数的最大值
#define MAXBSIZE 矩阵行,列数最大值
typedef struct{
int V[MAXSIZE]; //用以保存非零元的一维数组分量
int B[MAXBSIZE][MAXBSIZE]; //用以标记零元素和非零元素的二维数组B
int mu; //矩阵行数
int nu; //矩阵列数
int tu; //非零元个数
}Matrix;
Status Add(Matrix A1,Matrix A2,Matrix &A){//A=A1+A2
if (A1.mu!=A2.mu||A1.nu!=A2.nu) return ERROR;//两矩阵行列数不同
A.mu=A1.mu; A.nu=A1.nu; //和矩阵的行列数与求和矩阵相同
for (i=0;ifor(j=0;j A.B[i][j]=0;
//初始化和矩阵的标记数组
k=0;k1=0;k2=0;
//k为A.V中存放新非零元的位置
//k1为A1中当前非零元下标,k2为A2中当前非零元下标
for (i=0;ifor(j=0;jif (A1.B[i][j]){ //A1的当前位置有非零元
if (A2.B[i][j]){ //A2的当前位置也有非零元,需求和
t=A1.V[k1++]+A2.V[k2++];
if (t){ //和非零,生成A中新非零元
A.V[k++]=t;
A.B[i][j]=1;
}
else{ //A2的当前位置无非零元,以A1的当前非零元为A中新非零元
A.V[k++]=A1.V[k1++];
A.B[i][j]=1;
}//else
}//if
else{ //A1的当前位置无非零元
if (A2.B[i][j]){ //A2的当前位置有非零元,以A2的当前非零元为A中新非零元
A.V[k++]=A2.V[k2++];
A.B[i][j]=1;
}//if
}//else
A.tu=k;
return OK;
}//add