离散数学实验报告分析.docVIP

离散数学闭包实验报告 专业 12计算机科学与技术 学号 1座机电话号码 姓名 周谦益 时间 2011—11--15 一、实验目的 通过上机程序，进一步加深对关系中自反闭包，对称闭包，传递闭包的理解。 掌握Warshall算法。 学会用程序解决离散数学中的问题。 增强我们编写程序的能力。 二、实验内容 求有限集上给定关系的自反、对称和传递闭包（用Warshall算法）。 三、实验环境 我的实验是在Vs2008实验环境下完成的，而所设计的程序也在这个环境下通过了编译，运行和测试。 四、实验原理和实现过程 设计思路 在三种闭包中自反和对称闭包的求解很容易，对矩阵表示的关系，其自反闭包只要将矩阵的主对角线全部置为1就可；对称闭包则只需要将矩阵中数值为1的元素关于主对角线对称的元素数值也设为1就可以了；而对于传递闭包，用Warshall算法可以很方便的计算出来。下面我就来具体分析一下每一种闭包运算的设计。 自反闭包的设计：我们只要把关系矩阵的对角线的元素全赋值为1就可以啦。具体程序如下： 求自反闭包的程序： int Relation::Reflexive //求自反闭包的函数 for int i 0;i Len;i++ if !R[i][i] return 0 ; return 1 ; ; 对称闭包的求法：对于对称闭包，我们只只需要将矩阵中数值为1的元素的对称位置的元素数值也设为1就可以了。具体程序如下： 对称闭包： int Relation::Symmetric //对称闭包的函数 int i,j,K Len-1; for i 0;i K;i++ for j i+1;j Len;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;i Len;i++ for int j 0;j Len;j++ if R[i][j]! t.R[i][j] return 0 ; return 1 ; ; 把上面的每个子函数串在一起，再加上主函数，就可以实现这次实验的要求。 程序的完整代码如下所示： #include stdafx.h #include 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;i Len;i++ if A[i] x return i ; return -1 ; ; Relation::Relation int *p,int Q[Max][Max],int n int i,j; for i 0;i n;i++ A[i] p[i]; for j 0;j n;j++ R[i][j] Q[i][j]; ; Len i; ; int Relation::Reflexive //求自反闭包的函数 for int i 0;i Len;i++ if !R[i][