#include<stdlib.h>
#include<stdio.h>#define MaxSize 200
typedef int DataType;
typedef struct{         //三元组类型定义int i,j;DataType e;
}Triple;
typedef struct{             //矩阵类型定义Triple data[MaxSize];int rpos[MaxSize];     //用于存储三元组中的每一行的第一非零元素的位置int m,n,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{printf("请以行序顺序输入第%d个非零元素所在的行(0～%d),列(0～%d)元素值:",i,M->m-1,M->n-1);scanf("%d,%d,%d",&m,&n,&e);flag=0;                            /*初始化标志位*/if(m<0||m>M->m||n<0||n>M->n)    /*如果行号或列号不正确，标志位为1*/flag=1;if(i>0&&m<M->data[i-1].i||m==M->data[i-1].i&&n<=M->data[i-1].j) /*如果输入的顺序不正确，标志位为1*/flag=1;}while(flag);M->data[i].i=m;M->data[i].j=n;M->data[i].e=e;}return 1;
}
//稀疏矩阵的输出
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);
}
//稀疏矩阵相乘
void MultMatrix(TriSeqMatrix A,TriSeqMatrix B,TriSeqMatrix *C){int i,k,t,p,q,arow,brow,ccol;int temp[MaxSize];  //累加器int num[MaxSize];if(A.n!=B.m)     //如果矩阵A的列与B的行不相等，则返回return;C->m=A.m;     //初始化C的行数、列数和非零元素的个数C->n=B.n;C->len=0;if(A.len*B.len==0)      //只要有一个矩阵的长度为0，则返回return;/*---------------求矩阵B中每一行第一个非零元素的位置------------*/for(i=0;i<B.m;i++)       //初始化numnum[i]=0;for(k=0;k<B.len;k++){       //num存放矩阵B中每一行非零元素的个数i=B.data[k].i;num[i]++;}B.rpos[0]=0;//0,0for(i=1;i<B.m;i++)      //rpos存放矩阵B中每一行第一个非零元素的位置，是在B.data中的位置B.rpos[i]=B.rpos[i-1]+num[i-1];/*---------------求矩阵A中每一行第一个非零元素的位置--------------*/for(i=0;i<A.m;i++)           /*初始化num*/num[i]=0;  for(k=0;k<A.len;k++){i=A.data[k].i;num[i]++;}A.rpos[0]=0;     for(i=1;i<A.m;i++)          /*rpos存放矩阵A中每一行第一个非零元素的位置*/A.rpos[i]=A.rpos[i-1]+num[i-1];/*----------------计算两个矩阵的乘积---------------------------------*/for(arow=0;arow<A.m;arow++){        //依次扫描矩阵A的每一行for(i=0;i<B.n;i++)       //初始化累加器temptemp[i]=0;C->rpos[arow]=C->len;//对每个非0元处理if(arow<A.m-1)t=A.rpos[arow+1];elset=A.len;for(p=A.rpos[arow];p<t;p++){brow=A.data[p].j;      /*取出A中元素的列号*/if(brow<B.m-1)t=B.rpos[brow+1];elset=B.len;for(q=B.rpos[brow];q<t;q++) /*依次取出B中的第brow行，与A中的元素相乘*/{ccol=B.data[q].j;        temp[ccol]+=A.data[p].e*B.data[q].e;/*把乘积存入temp中*/}}for(ccol=0;ccol<C->n;ccol++)  /*将temp中元素依次赋值给C*/if(temp[ccol]){if(++C->len>MaxSize) return;C->data[C->len-1].i=arow;C->data[C->len-1].j=ccol;C->data[C->len-1].e=temp[ccol];}}
}
void main(){//0,0不能为0TriSeqMatrix M,N,Q;CreateMatrix(&M);PrintMatrix(M);CreateMatrix(&N);PrintMatrix(N);printf("矩阵M和N的乘积为:\n");MultMatrix(M,N,&Q);PrintMatrix(Q);system("pause");
}