离散数学 关系运算代码(c++)
包含以下功能:
1.求有限集上给定关系的自反闭包、对称闭包、传递闭包。
2.求有限集上等价关系的数目。
3.输入集合和等价关系,求对应的商集。
4.判断是否为等价关系。
Relation.h
#ifndef RELATION_H_INCLUDED
#define RELATION_H_INCLUDED
#include<string.h>
using namespace std;
struct Relation
{
bool **Mat; //存放关系矩矩阵
string *Arr; //存放元素集合
int n; //存放元素个数
};
void InitRelation(Relation &R,int n) //初始化关系矩阵
{
R.n=n;
R.Arr=new string[n];
R.Mat=new bool *[n];
for(int j=0;j<n;j++)
{
R.Mat[j]=new bool[n];
for(int k=0;k<n;k++)
R.Mat[j][k]=false;
}
}
void Show(Relation R) //输出关系矩阵
{
cout<<'{';
for(int i=0;i<R.n-1;i++)
cout<<R.Arr[i]<<", ";
cout<<R.Arr[R.n-1]<<'}'<<endl;
for(int i=0;i<R.n;i++)
{
for(int j=0;j<R.n;j++)
cout<<R.Mat[i][j];
cout<<endl;
}
}
int Locate(Relation R,string str) //寻找目标元素在集合中的位置
{
for(int i=0;i<R.n;i++)
{
if(str==R.Arr[i])
return i;
}
return -1;
}
void CreateRelation(Relation &R) //输入关系
{
string a,b;
cout<<"请输入各元素"<<endl;
for(int i=0;i<R.n;i++)
cin>>R.Arr[i];
cout<<"请给定各元素之间关系,以# #结束"<<endl;
cin>>a>>b;
while(a!="#"||b!="#")
{
int i,j;
i=Locate(R,a);
j=Locate(R,b);
if(i<0||j<0)
cout<<"元素不合法,请检查!"<<endl;
else
R.Mat[i][j]=true;
cin>>a>>b;
}
}
Relation ReflexiveClosure(Relation &R) //求自反闭包
{
Relation S;
InitRelation(S,R.n);
for(int i=0;i<R.n;i++)
{
S.Arr[i]=R.Arr[i];
for(int j=0;j<R.n;j++)
S.Mat[i][j]=R.Mat[i][j];
}
for(int i=0;i<S.n;i++)
S.Mat[i][i]=true;
return S;
}
Relation SymmetricClosure(Relation &R) //求对称闭包
{
Relation S;
InitRelation(S,R.n);
for(int i=0;i<R.n;i++)
{
S.Arr[i]=R.Arr[i];
for(int j=0;j<R.n;j++)
S.Mat[i][j]=R.Mat[i][j];
}
for(int i=0;i<S.n;i++)
{
for(int j=i;j<S.n;j++)
{
if(S.Mat[i][j]==true)
S.Mat[j][i]=true;
}
}
return S;
}
Relation TransitiveClosure(Relation &R) //求传递闭包
{
Relation S;
InitRelation(S,R.n);
for(int i=0;i<R.n;i++)
{
S.Arr[i]=R.Arr[i];
for(int j=0;j<R.n;j++)
S.Mat[i][j]=R.Mat[i][j];
}
for(int i=0;i<S.n;i++)
{
for(int j=0;j<S.n;j++)
{
if(S.Mat[i][j]==true)
{
for(int k=0;k<S.n;k++)
{
if(S.Mat[j][k]==true)
S.Mat[i][k]=true;
}
}
}
}
return S;
}
int NumOfEquivalent(int n) //求等价关系数
{
int a[n+1][n+1];
int sum=0;
for(int i=1;i<=n;i++)
{
a[1][i]=1;
a[i][i]=1;
}
for(int i=2;i<=n;i++)
{
for(int j=i+1;j<=n;j++)
a[i][j]=a[i-1][j-1]+i*a[i][j-1];
}
for(int i=1;i<=n;i++)
sum+=a[i][n];
return sum;
}
bool isReflexive(Relation R) //验证自反性
{
for(int i=0;i<R.n;i++)
{
if(R.Mat[i][i]!=true)
return false;
}
return true;
}
bool isSymmetric(Relation R) //验证对称性
{
for(int i=0;i<R.n;i++)
{
for(int j=i+1;i<R.n;j++)
{
if(R.Mat[i][j]!=R.Mat[j][i])
return false;
}
}
return true;
}
bool isTreansitive(Relation R) //验证传递性
{
for(int i=0;i<R.n;i++)
{
for(int j=0;j<R.n;j++)
{
if(R.Mat[i][j]==1)
{
for(int k=0;k<R.n;k++)
{
if((R.Mat[j][k]==1&&R.Mat[i][k]==1)==false)
return false;
}
}
}
}
return true;
}
void QuotientSet(Relation R) //根据集合与关系求商集
{
bool visited[R.n];
for(int i=0;i<R.n;i++)
visited[i]=false;
if(isReflexive(R)==false&&isSymmetric(R)==false&&isTreansitive(R)==false)
{
cout<<"您输入的不是等价关系!"<<endl;
return;
}
else
{
cout<<"商集为: "<<"{";
for(int i=0;i<R.n;i++)
{
if(visited[i]==false)
{
visited[i]=true;
cout<<"{";
cout<<R.Arr[i];
for(int j=0;j<R.n;j++)
{
if(R.Mat[i][j]==1&&visited[j]==false)
{
visited[j]=true;
cout<<","<<R.Arr[j];
}
}
cout<<"}";
}
}
cout<<"}";
}
}
#endif // RELATION_H_INCLUDED
Main.cpp
#include <iostream>
#include "Relation.h"
using namespace std;
int main()
{
Relation A;
cout<<"请输入元素个数"<<endl;
int n;
cin>>n;
InitRelation(A,n);
cout<<"该集合有"<<NumOfEquivalent(n)<<"种等价关系"<<endl;
CreateRelation(A);
Show(A);
Show(TransitiveClosure(A));
Show(ReflexiveClosure(A));
Show(SymmetricClosure(A));
QuotientSet(A);
}