離散实验报告(二).docVIP

“离散数学”实验报告 专业：计算机嵌入式 班级学号： 1045534125 姓名：孟腾腾 指导老师：闫仁武 日期：2011/11/15 一 、实验内容 求有限集上给定关系的自反、对称和传递闭包 二、实验思路 设计思路 在三种闭包中自反和对称闭包的求解很容易，对矩阵表示的关系，其自反闭包只要将矩阵的主对角线全部置为1就可；对称闭包则加上关系的转置矩阵，并将相加后大于1的元素值设置为1就可得到；传递闭包则有两种算法： 算法 R是集合X上的二元关系，则其自反闭包为 r（R）=RUI (x)( I (x)为恒等关系) R是集合X上的二元关系，则其对称闭包为 s（R）=RUR^c R是集合X上的二元关系，则其传递闭包为 直接根据计算，过程略 方法二：Warshall算法(1962) 设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)，否则结束． 三、流程图 主函数 四、程序清单 // Relation.cpp : Defines the entry point for the console application. // #include iostream.h int const Max=20; class Relation { protected: int Len,A[Max],R[Max][Max]; int Index(int x); public: Relation(){Len=0;}; Relation(int *p,int Q[Max][Max],int n); int IsofA(int x); int IsofR(int x,int y); Relation UnionA(Relation R1); Relation CrossR(Relation R1); Relation C_r(); Relation C_s(); Relation C_t(); Relation C_t(int); int Reflexive(); int Symmetric(); 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); return(-1); }; Relation::Relation(int *p,int Q[Max][Max],int n) { int i,j; for(i=0;in;i++) { A[i]=p[i]; for(j=0;jn;j++) R[i][j]=Q[i][j]; }; Len=i; }; int Relation::Reflexive() { for(int i=0;iLen;i++) if(!R[i][i])return(0); return(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); }; int Relation::IsofA(int x) { for(int i=0;iLen;i++) if(A[i]==x)return(1); return(0); }; int Relation::IsofR(int x,int y) { int i,j; i=Index(x); j=Index(y); if(i==-1||j==-1||R[i][j]==0)return(0); return(1); }; Relation Relation::UnionA(Relation R1) { int i,j; Relation t; for(i=0;iLen;i++)t.A[i]=A[i]; t.Len=Len; for(j=0;jR1.Len;j++) { for(i=0;iLen;i++) if(t.A[i]==R1.A[j]