离散数学 关系运算（c++）

离散数学 关系运算代码（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);
}