离散数学实验报告解析.doc

离散数学闭包实验报告 专业 12计算机科学与技术 学号 姓名 周谦益 时间 2011—11--15 一、实验目的 通过上机程序，进一步加深对关系中自反闭包，对称闭包，传递闭包的理解。 掌握Warshall算法。 学会用程序解决离散数学中的问题。 增强我们编写程序的能力。 二、实验内容 求有限集上给定关系的自反、对称和传递闭包（用Warshall算法）。 三、实验环境 我的实验是在Vs2008实验环境下完成的，而所设计的程序也在这个环境下通过了编译，运行和测试。 四、实验原理和实现过程 设计思路 在三种闭包中自反和对称闭包的求解很容易，对矩阵表示的关系，其自反闭包只要将矩阵的主对角线全部置为1就可；对称闭包则只需要将矩阵中数值为1的元素关于主对角线对称的元素数值也设为1就可以了；而对于传递闭包，用Warshall算法可以很方便的计算出来。下面我就来具体分析一下每一种闭包运算的设计。 自反闭包的设计：我们只要把关系矩阵的对角线的元素全赋值为1就可以啦。具体程序如下： 求自反闭包的程序： int Relation::Reflexive()//求自反闭包的函数 { for(int i=0;iLen;i++) if(!R[i][i])return(0); return(1); }; 对称闭包的求法：对于对称闭包，我们只只需要将矩阵中数值为1的元素的对称位置的元素数值也设为1就可以了。具体程序如下： 对称闭包： int Relation::Symmetric()//对称闭包的函数 { int i,j,K=Len-1; for(i=0;iK;i++) for(j=i+1;jLen;j++) if(R[i][j]!=R[j][i])return(0); return(1); }; 传递闭包设计：传递闭包我主要用Warshall算法来求。 书上的Walshall算法的伪代码如下： 设R的关系矩阵为M (1)令矩阵A=M (2)置i=1 (3)对所有的j，若A[j，i]=1，则对于 k=1，2，…，n，令A[j，k]=A[j，k]+A[i，k] (4) i=i+l． (5)若i≤n，则转到(3)，否则结束 根据Warshall算法，我设计求出了关系的传递闭包，具体程序如下： int Relation::Transitive()//传递闭包函数 { Relation t; t=C_t(1); for(int i=0;iLen;i++) for(int j=0;jLen;j++) if(R[i][j]!=t.R[i][j])return(0); return(1); }; 把上面的每个子函数串在一起，再加上主函数，就可以实现这次实验的要求。 程序的完整代码如下所示： #include stdafx.h #include iostream using namespace std; int const Max=20; class Relation { protected: int Len,A[Max],R[Max][Max]; int Index(int x); public: Relation(){Len=0;};//空集A，空关系R Relation(int *p,int Q[Max][Max],int n);//n个元素集合A及其关系 int IsofA(int x);//x是否属于A int IsofR(int x,int y); Relation UnionA(Relation R1);//AUR1 Relation CrossR(Relation R1); Relation C_r(); Relation C_s(); Relation C_t(); Relation C_t(int); int Reflexive(); int Symmetric();//判断R是否对称 int Transitive(); int Equivalent_Relations(); void PrintA(); void PrintR(); }; int Relation::Index(int x) { int i; for(i=0;iLen;i++) if(A[i]==x)return(i); r