第九周 项目三【系数矩阵的三元组表示法的实现与应用2】

/*问题及代码 
 *Copyright(c)2016,烟台大学计算机学院 
 *All right reserved. 
 *文件名称：稀疏矩阵的相加 .cpp 
 *作者：李潇 
 *完成日期；2016年11月3日 
 *版本号；v1.0 
 * 
 *问题描述：　采用三元组存储稀疏矩阵，设计两个稀疏矩阵相加的运算算法  
提示1：两个行数、列数相同的矩阵可以相加  
提示2：充分利用已经建立好的算法库解决问题 
 
 
 *输入描述：两个稀疏矩阵 
 *程序输出：稀疏矩阵的三元组表示结果以及相加的结果 
*/

//tup.h
#include "stdio.h" 
#define M 6
#define N 7
#define MaxSize 100
typedef int ElemType;

typedef struct
{
    int r;
	int c;
	ElemType d;
}TupNode;
typedef struct
{
	int rows;
	int cols;
	int nums;
	TupNode data[MaxSize];
}TSMatrix;
void CreatMat(TSMatrix &t,ElemType A[M][N]);
bool Value(TSMatrix &t,ElemType x,int i,int j);
bool Assign(TSMatrix t,ElemType &x,int i,int j);
void DispMat(TSMatrix t);
void TranTat(TSMatrix t,TSMatrix &tb);
bool MatAdd(TSMatrix a,TSMatrix b,TSMatrix &c);

//tup.cpp
#include "tup.h"
#include <stdio.h>
void CreatMat(TSMatrix &t,ElemType A[M][N])
{
	int i,j;
	t.rows=M;t.cols=N;t.nums=0;
	for(i=0;i<M;i++)
	{
		for(j=0;j<N;j++)
		{
			if(A[i][j]!=0)
			{
				t.data[t.nums].r=i;t.data[t.nums].c=j;
				t.data[t.nums].d=A[i][j];t.nums++;
			}
		}
	}
}
bool Value(TSMatrix &t,ElemType x,int i,int j)
{
	int k=0,k1;
	if(i>=t.rows || j>=t.cols)
		return false;
	while (k<t.nums && i>t.data[k].r) k++;
	while (k<t.nums && i== t.data[k].r && j>t.data[k].c) k++;
	if(t.data[k].r==i && t.data[k].c==j)
		t.data[k].d=x;
	else
	{
		for(k1=t.nums-1;k1>=k;k1--)
		{
			t.data[k1+1].r=t.data[k1].r;
			t.data[k1+1].c=t.data[k1].c;
			t.data[k1+1].d=t.data[k1].d;
		}
		t.data[k].r=i;t.data[k].c=j;t.data[k].d=x;
		t.nums++;
	}
	return true;
}
bool Assign(TSMatrix t,ElemType &x,int i,int j)
{
   int k=0;
   if(i>=t.rows || j>=t.cols)
	   return false;
   while (k<t.nums && i>t.data[k].r) k++;
   while (k<t.nums && i==t.data[k].r && j>t.data[k].c) k++;
   if (t.data[k].r==i && t.data[k].c==j)
	   x=t.data[k].d;
   else
	   x=0;
   return true;
}
void DispMat(TSMatrix t)
{
	 int i;
	 if (t.nums<=0)
		 return ;
	 printf("\t%d\t%d\t%d\n",t.rows,t.cols,t.nums);
	 printf("\t--------------------\n");
	 for(i=0;i<t.nums;i++)
	    printf("\t%d\t%d\t%d\n",t.data[i].r,t.data[i].c,t.data[i].d);
}
void TranTat(TSMatrix t,TSMatrix &tb)
{
	int p,q=0,v;
	tb.rows=t.cols;tb.cols=t.rows;tb.nums=t.nums;
	if(t.nums!=0)
	{
		for (v=0;v<t.cols;v++)
		   for(p=0;p<t.nums;p++)
			   if(t.data[p].c==v)
			   {   
			      tb.data[q].r=t.data[p].c;
				  tb.data[q].c=t.data[p].r;
				  tb.data[q].d=t.data[p].d;
				  q++;
			   }
	}
}
bool MatAdd(TSMatrix a,TSMatrix b,TSMatrix &c)
{
    int i,j;
    ElemType va,vb,vc;
    if (a.rows!=b.rows || a.cols!=b.cols)
        return false;                        //行数或列数不等时不能进行相加运算
    c.rows=a.rows;
    c.cols=a.cols;       //c的行列数与a的相同
    c.nums=0;
    for(i=0; i<M; i++)
        for(j=0; j<N; j++)
        {
            Assign(a,va,i,j);
            Assign(b,vb,i,j);
            vc=va+vb;
            if(vc)
                Value(c,vc,i,j);
        }
    return true;
}

//main.cpp
#include <stdio.h>
#include "tup.h"
int main()
{
   TSMatrix ta,tb,tc;
    int A[M][N]=
    {
        {0,0,1,0,0,0,0},
        {0,2,0,0,0,0,0},
        {3,0,0,0,0,0,0},
        {0,0,0,5,0,0,0},
        {0,0,0,0,6,0,0},
        {0,0,0,0,0,7,4}
    };
    int B[M][N]=
    {
        {0,0,10,0,0,0,0},
        {0,0,0,20,0,0,0},
        {0,0,0,0,0,0,0},
        {0,0,0,50,0,0,0},
        {0,0,20,0,0,0,0},
        {0,0,0,10,0,0,4}
    };
    CreatMat(ta,A);
    CreatMat(tb,B);
    printf("A:\n");
    DispMat(ta);
    printf("B:\n");
    DispMat(tb);
    if(MatAdd(ta, tb, tc))
    {
        printf("A+B:\n");
        DispMat(tc);
    }
    else
    {
        printf("相加失败\n");
    }
    return 0;
}

运行结果：

知识点总结：

 稀疏矩阵三元组的相加运算

心得体会：

 思路和二元组的是一样的，但是得在二元组的基础上做变通