稀疏矩阵运算器-数据结构课程设计|数据结构课程设计
/*******************************稀疏矩阵运算器*************************************/#include "stdio.h"#include "stdlib.h"#define MAXSIZE 20 /*非零元个数*/#define MAXRC 20#define MAXR 20#define MAXC 20typedef struct{ int i,j; /*非零元行下标和列下标*/ int e;} Triple ; typedef struct{ Triple data[MAXSIZE+1]; /*非零元三元组，data[0]未用*/ int rpos[MAXRC+1]; /*各行第一非零元的位置表*/ int mu,nu,tu; /*矩阵的行数，列数，非零元数*/} RLSMatrix ;
void Print(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T);void ADD(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*加法*/{ int x=1; int y=1; int z=1; T->tu=0; T->mu=A->mu; T->nu=A->nu; while(x<=A->tu) { if(A->data[x].i
data[y].i) { T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; T->data[z].e=A->data[x].e; x++;z++; } else if(A->data[x].i==B->data[y].i) { if(A->data[x].j>B->data[y].j) { T->data[z].i=B->data[y].i; T->data[z].j=B->data[y].j; T->data[z].e=B->data[y].e; y++;z++; } else if(A->data[x].j==B->data[y].j){ if((T->data[z].e=A->data[x].e+B->data[y].e)!=0) { T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; x++; y++; z++; } else { x++; y++; } } else if(A->data[x].jdata[y].j) { T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; T->data[z].e=A->data[x].e; x++;z++; } } else if(A->data[x].i>B->data[y].i) { T->data[z].i=B->data[y].i; T->data[z].j=B->data[y].j; T->data[z].e=B->data[y].e; y++;z++; } if(y>=B->tu) break; } if(x<=A->tu){ while(x<=A->tu){ T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; T->data[z].e=A->data[x].e; x++;z++; } } if(y<=B->tu){ while(y<=B->tu){ T->data[z].i=B->data[y].i; T->data[z].j=B->data[y].j; T->data[z].e=B->data[y].e; y++;z++; } } T->tu=z; Print(A,B,T);}
void SUB(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*减法*/{ int x=1; int y=1; int z=1; T->tu=0; T->mu=A->mu; T->nu=A->nu; while(x<=A->tu) { if(A->data[x].idata[y].i) { T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; T->data[z].e=A->data[x].e; x++;z++; } else if(A->data[x].i==B->data[y].i) { if(A->data[x].j>B->data[y].j) { T->data[z].i=B->data[y].i; T->data[z].j=B->data[y].j; T->data[z].e=0-B->data[y].e; y++;z++; } else if(A->data[x].j==B->data[y].j){ if((T->data[z].e=A->data[x].e-B->data[y].e)!=0) { T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; x++; y++; z++; } else { x++; y++; } } else if(A->data[x].jdata[y].j) { T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; T->data[z].e=A->data[x].e; x++;z++; } }
else if(A->data[x].i>B->data[y].i) { T->data[z].i=B->data[y].i; T->data[z].j=B->data[y].j; T->data[z].e=0-B->data[y].e; y++;z++; } if(y>=B->tu) break; } if(x<=A->tu) { while(x<=A->tu) { T->data[z].i=A->data[x].i; T->data[z].j=A->data[x].j; T->data[z].e=A->data[x].e; x++;z++; } } if(y<=B->tu) { while(y<=B->tu){ T->data[z].i=B->data[y].i; T->data[z].j=B->data[y].j; T->data[z].e=0-B->data[y].e; y++;z++; } } T->tu=z; Print(A,B,T);}
int MUL(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*乘法*/{ int arow,ccol,tp,p,q,blow,t,i; int ctemp[MAXRC]; if(A->nu!=B->mu) return 0; T->mu=A->mu;T->nu=B->nu;T->tu=0; if(A->tu*B->tu!=0) { for(arow=1;arow<=A->mu;++arow) { for(i=1;i<=T->nu;i++) ctemp[i]=0; T->rpos[arow]=T->tu+1; if(arowmu) tp=A->rpos[arow+1]; else{tp=A->tu+1;} for(p=A->rpos[arow];pdata[p].j; if(blowmu) t=B->rpos[blow+1]; else t=B->tu+1; for(q=B->rpos[blow];qdata[q].j; ctemp[ccol] += A->data[p].e*B->data[q].e; } /*for q*/ } /*求得Q中的crow( arow)行的非零元*/ for(ccol=1;ccolnu;++ccol) { if(ctemp[ccol]) { if(++T->tu>MAXSIZE) return 0; T->data[T->tu].i=arow; T->data[T->tu].j=ccol; T->data[T->tu].e=ctemp[ccol]; } } } Print(A,B,T); return 1; }}
void Printt(RLSMatrix *A) /*打印矩阵*/{ int row,col,i=1,k=0; for(row=1;row<=A->mu;row++) { for(col=1;col<=A->nu;col++) { if(A->data[i].i==row&&A->data[i].j==col&&i<=A->tu) printf("%-4d",A->data[i++].e); else printf("%-4d",k); } printf("\n"); }}
void Print(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*打印函数*/{ clrscr(); if(A->mu!=0) {printf("Print Matrix A:\n");Printt(A); printf("\n");} else printf("Matrix A is blank!\n"); if(B->tu!=0) {printf("Print Matrix B:\n");Printt(B); printf("\n");} else printf("Matrix B is blank!\n"); if(T->tu!=0) {printf("Print Matrix T:\n");Printt(T); printf("\n");} else {printf("Matrix T is blank!\n");getch();} getch();}
void Enter(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*输入函数*/{ int i,k,num[MAXRC]; char cmd; clrscr(); printf("Enter the mu and nu and tu of Matrix A(mu,nu,tu): \n"); /*输入矩阵A*/ scanf("%d,%d,%d",&A->mu,&A->nu,&A->tu); while(A->mu<=0||A->mu>MAXR||A->nu<=0||A->nu>MAXC) { printf("Enter the mu and nu and tu of Matrix A(mu,nu,tu): \n"); scanf("%d,%d,%d",&A->mu,&A->nu,&A->tu); } printf("Enter the Matrix A(i-j:e):\n"); for(i=1;i<=A->tu;i++) scanf("%d-%d:%d",&A->data[i].i,&A->data[i].j,&A->data[i].e); for(i=1;i<=A->mu;i++) num[i]=0; for(i=1;i<=A->tu;i++) ++num[A->data[i].i]; A->rpos[1]=1; for(i=2;i<=A->mu;i++) A->rpos[i]=A->rpos[i-1]+num[i-1];
printf("Enter the mu and nu and tu of Matrix B(mu,nu,tu): \n"); /*输入矩阵B*/ scanf("%d,%d,%d",&B->mu,&B->nu,&B->tu); while(A->mu<=0||A->mu>MAXR||A->nu<=0||A->nu>MAXC) { printf("Enter the mu and nu and tu of Matrix B(mu,nu,tu): \n"); scanf("%d,%d,%d",&B->mu,&B->nu,&B->tu); } printf("Enter the Matrix B(i-j:e):\n"); for(i=1;i<=B->tu;i++) scanf("%d-%d:%d",&B->data[i].i,&B->data[i].j,&B->data[i].e); for(i=1;i<=B->mu;i++) num[i]=0; for(i=1;i<=B->tu;i++) ++num[B->data[i].i]; B->rpos[1]=1; for(i=2;i<=B->mu;i++) B->rpos[i]=B->rpos[i-1]+num[i-1];
do{ printf("Enter the operation you want:"); /*选择需要的运算*/ cmd=getche(); }while(cmd!='+'&&cmd!='-'&&cmd!='*'); getch(); switch(cmd) { case '+' : if(A->mu==B->mu&&A->nu==B->nu) ADD(A,B,T); break; case '-' : if(A->mu==B->mu&&A->nu==B->nu) SUB(A,B,T); break; case '*' : if(A->nu==B->mu) MUL(A,B,T); break; default : printf("No thia operation!"); break; }}void InitializationMatrix(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T) /*初始化矩阵*/{ A->mu=A->nu=A->tu=0; B->mu=B->nu=B->tu=0; T->mu=T->nu=T->tu=0;}/*************************************主函数************************************/void InitializationMatrix(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T);void Enter(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T);void Print(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T);void ADD(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T);void SUB(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T);int MUL(RLSMatrix *A,RLSMatrix *B,RLSMatrix *T);void Initialization();void Myinformation();void Exit();void main(){
int i,k,num[MAXRC]; char cmd; RLSMatrix *A,*B,*T; A=( RLSMatrix *)malloc(sizeof( RLSMatrix)); /*分配空间*/ B=( RLSMatrix *)malloc(sizeof( RLSMatrix)); T=( RLSMatrix *)malloc(sizeof( RLSMatrix));
Myinformation(); InitializationMatrix(A,B,T); while(1) { Initialization(); do{ gotoxy(10,10); printf("Enter E/e ,Q/q ,P/p to continue:"); cmd=getche(); printf("\n"); }while(cmd!='q'&&cmd!='Q'&&cmd!='e'&&cmd!='E'&&cmd!='p'&&cmd!='P'); switch(cmd) { case 'q' : Exit(); break; case 'Q' : Exit(); break; case 'e' : Enter(A,B,T); break; case 'E' : Enter(A,B,T); break; case 'p' : Print(A,B,T); break; case 'P' : Print(A,B,T); break; } }}
void Initialization() /*初始化函数*/{ clrscr(); printf("*********************************************************************************"); gotoxy(7,2); printf("Enter the RLSMatrix--E"); gotoxy(38,2); printf("Print--P"); gotoxy(55,2); printf("Quit--Q"); gotoxy(80,2); printf("*"); printf("********************************************************************************"); gotoxy(1,22); printf("*********************************************************************************"); gotoxy(10,23); printf("Enter a operation code: E/e ,Q/q ,P/p to continue: "); gotoxy(80,23); printf("*"); printf("********************************************************************************");
}void Myinformation() /*我的信息函数*/{ int i; clrscr(); for(i=1;i<=240;i++) printf("\1"); gotoxy(15,8); printf("This is a Matrix arithmetic calculator!"); gotoxy(15,12); printf("Name: liyongjun."); gotoxy(15,13); printf("Num: 3104006893."); gotoxy(15,14); printf("Grade: 2004."); gotoxy(15,15); printf("Class: cumputer science and technology 10"); gotoxy(1,22); for(i=1;i<=240;i++) printf("\1"); gotoxy(1,18); printf("Print any key to continue..."); getch();}
void Exit() /*退出函数*/ { char cmd; do{ gotoxy(10,11); printf("Do you want to out?Y/N:"); cmd=getche(); }while(cmd!='y'&&cmd!='Y'&&cmd!='n'&&cmd!='N'); if(cmd=='y'||cmd=='Y') exit(0);
}
