离散数学 Warshall算法求传递闭包 C语言实现

 

求传递闭包有一种有效算法—Warshall算法，这种算法也便于计算机实现。 
(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加1；（i是行，j是列） 
(5)如果i≤n，则转到步骤3)，否则停止。 

例如：一个的矩阵M：

 

第一步：当i=1时;找到M[j][i]==1时的位置；即M[1][1]=1；

再将第 j 行加上第 i 行得到新的第 j 行，即将第1行加上第1行得到新的第1行；有以下新矩阵，i+1：

第二步：当i=2时；重复第一步的方法：有M[1][2]=1,M[3][2]=1;

再将第2行分别加到第1行和第3行上去得到新的第1、3行；有以下新矩阵，i+1：

第三步：当i=3时；重复上述步骤：有M[1][3]=1，M[2][3]=1,M[3][3]=1;

再将第3行分别加到第1行，第2行，第3行上去；得到新矩阵，i+1：

  由于i+1=4>3(行数)，停止，该矩阵为原矩阵的t(R)的关系矩阵；

代码：有每步过程

#include<bits/stdc++.h>
#include<stdio.h>
#include<iostream>
#include<string.h>
using namespace std;
int main()
{
    bool matrix[100][100];
    int m,n;
    memset(matrix,0,sizeof(matrix));
    while(scanf("%d%d",&m,&n)!=EOF)
    {
        int i,j,k;
        for(i=1;i<=m;i++)
        {
            for(j=1;j<=n;j++)
            {
                scanf("%d",&matrix[i][j]);
            }
        }
        cout<<"..............."<<endl;
        
        for(i=1;i<=m;i++)
        {
            for(j=1;j<=n;j++)
            {
                if (matrix[j][i]==1)
                {
                    for (k=0;k<=n;k++)
                    {
                        matrix[j][k]+=matrix[i][k];
                        continue;
                    }
                }
            }
                for(int a=1;a<=m;a++)
                {
                    for(int b=1;b<=n;b++)
                    {
                        cout<<matrix[a][b]<<" ";
                    }
                cout<<"\n";
                }
            cout<<"........."<<"\n";
        }

        cout<<"结果为:"<<"\n";
        for(i=1;i<=m;i++)
        {
            for(j=1;j<=n;j++)
            {
                cout<<matrix[i][j]<<" ";
            }
            cout<<endl;
        }
    }
return 0;
}