實驗一   判斷關系R是否為自反關系及對稱關系

A）功能

已知關系R由關系矩陣M給出，要求判斷由M表示的這個關系是否為自反關系和對稱關系。

B）基本思想

   從給定的關系矩陣來判關系R是否為自反、對稱、是很容易的。若M（R的關系矩陣）的主對角線元素均為1，則R是自反關系；若M為對稱矩陣，則R是對稱關系。因為R為自反的、對稱的是等價關系的必要條件，所以，本算法可以作為判等價關系算法的子程序給出。因此，在程序中設置標志變量F，若R是自反的、對稱的。則F=1，否則F=0。

實驗二  判關系R是否為可傳遞關系

A）功能

給出關系R的關系矩陣M，判關系R是否為可傳遞的。

B)基本思想

一個關系R的可傳遞性定義告訴我們，若關系R是可傳遞的，則必有：m_ik=1∧m_kj=1  m_ij=1。這個式子也可改寫成為： m_ij =0  m_ik =0∨m_kj=0。我們就是根據后一個公式來完成判斷可傳遞性這一功能的。同算法3.1一樣，可傳遞性也是等價關系的必要條件。因此，在這里也設置標志變量F，F=1標志關系R是可傳遞的，F=0標志關系R是不可傳遞的。

頭文件
Matrix.h
class Matrix
{
public:
    Matrix(int n);
    ~Matrix();
    void revilt(); //自反
    void sym();    //對稱
    void printMatrix(); //輸出
    void transfer();//傳遞
    int F();        //記錄R的性質
private:
    int n;
    int **a;
    bool f1;
    bool f2;
    bool f3;
};
//Matrix.cpp
#include "matrix.h"
#include <stdlib.h>
#include <iostream>
#include <time.h>

using namespace std;

 Matrix::Matrix(int h)
{
    bool f2=false;
    bool f1=false;
    bool f3=false;
    n=h;
    a= new int *[n];
    for(int k=0;k<n;k++)
        a[k] = new int[n];
    srand((unsigned)time(NULL));
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            if(rand()%2==0)
                a[i][j]=0;
            else
                a[i][j]=1;
        }
    }
}


 void Matrix::printMatrix()
 {
     for(int m=0;m<n;m++)
     {
         for(int j=0;j<n;j++)
            cout<<" "<<a[m][j]<<flush;
         cout<<endl;
     }
 }


 Matrix::~Matrix()
 {
 }


void Matrix::revilt()
{
    bool flag1=true;
    bool flag2=true;
    
    for(int i=0;i<n;i++)
    {
        if(a[i][i]!=1)
            flag1=false;
        if(a[i][i]!=0)
            flag2=false;
    }
    if(flag1==true)
    {
        cout<<"R是自反的！"<<endl;
        f1=true;
    }
    else if(flag2==true)
        cout<<"R是反自反的！"<<endl;
    else
        cout<<"R既不是自反的也不是反自反的！"<<endl;
}


void Matrix::sym()
{
    int i,j;
    bool flag3=true;
    for(i=0;i<n;i++)
        for(j=0;j<n;j++)
            if(a[i][j]!=a[j][i])
                flag3=false;
    if(flag3)
    {
        cout<<"R是對稱的！！"<<endl;
        f2=true;
    }
    else
        cout<<"R非對稱"<<endl;

}


int Matrix::F()
{
    if(f1==true && f2==true && f3==true)
        return 1;
    else
        return 0;
}


void Matrix::transfer()
{
    bool flag=true;
    int a2[20][20];
    for(int i=0;i<n;i++){
                    
        for(int x=0;x<n;x++){
            int m=0;
            for(int y=0;y<n;y++)
                m+=a[i][y]*a[y][x];
            if(m>1)
                m=1;
            a2[i][x] = m;
            }
        }
                
    cout<<"the M*M is "<<endl;
    for(int i=0;i<n;i++) {
        for(int j=0;j<n;j++) {
            cout<<" "<<a2[i][j]<<flush;
            }
            cout<<endl;
        }
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                {
                    if(a2[i][j]==1 && a[i][j]!=1) flag=false; 
                }
                if(flag)
                {
                    cout<<"這個是傳遞矩陣"<<endl;
                    f3 =true;
                }
                else
                    cout<<"這個是不傳遞矩陣"<<endl;
}
//主函數 MatrixMain.cpp
#include <iostream>
#include "Matrix.h"
using namespace std;
int main()
{
    Matrix a(3);
    cout<<"關系R所對應矩陣為"<<endl;
    a.printMatrix();
    a.revilt();
    a.sym();
    a.transfer();
    cout<<"F="<<a.F()<<endl;
    return 0;
}