离散复习资料之一（Warshall算法）

我们先对R（n）进行解释一下，不然你们看不太懂。

就是R（n）表示，通过从 1...n 个点绕到到达所求点，注意原来的位置的点，发生改变的就是刚才的含义去解释。

 给一个例子解释一下。

 

其中的R（i）...R(N)；(0<=N<=4)

初始的R（0）就是可达矩阵。

R（1）表示在R(0)的基础上面变化

R（1）表示在R(0)的基础上面变化

R（2）表示在R(1)的基础上面变化

R（3）表示在R(2)的基础上面变化

R（4）表示在R(3)的基础上面变化

 

关键：

R(0)表示直接相连的关系，可以到达就记为1，否则就记为0.

R(1)表示从该点通过1可以到达目标点的记为1，否则就记为0.（多出来的点的关系）

R（2）表示从该点通过1,2可以到达目标点到达的记为1.否则就记为0.（与R（1）相比，多出来的点的关系。）

R（3）与R（4）表示的意思同上。

Warshall算法思想：

      R（0）就只是有向图的邻接矩阵。

      R（1）包含可以用第一个顶点作为中间点的路径信息。

      R（n）说明可以用所有的顶点作为中间点寻找有向路径，所以R（n）就是我们要求的传递闭包,即可达矩阵。

      这个算法的中心是我们可以通过R（k－1）来计算R（k）的所有元素。

算法的步骤如下：

R(k) <- R（k-1）

for(k<-1;k<=n;k++)

   for(i<-1;i<=n;i++)

     for(j<-1;j<=n;j++)

          R(k)[i,j] = R(k-1)[i,j] or R(k-1)[i,k] and R(k-1)[k,j]

return  R(n)

 

通过上面的算法步骤，我们也可以知道，推出R（n），只要知道R（n-1）的关系矩阵从而推出R（n）。

每一点的计算关系都会遵循一个计算公式。

  “R(k)[i,j] = R(k-1)[i,j] or R(k-1)[i,k] and R(k-1)[k,j]”

从而程序代码也可以按照上面的计算关系进行计算，改变。

看一个代码改变一下求R（N），已知的是R（0）。

代码如下：
 

#include<iostream>  
#include<cstdio>  
  
using namespace std;  
  
int temp[5][5],temp1[5][5];  
  
void exchange(int sum,int k)  
{  
    for(int i=1;i<=sum;i++)  
    {  
        for(int j=1;j<=sum;j++)  
        {  
            temp[i][j] = temp1[i][j] or temp1[i][k] and temp1[k][j];  
        }  
    }  
      
    for(int i=1;i<=sum;i++)  
    {  
        for(int j=1;j<=sum;j++)  
        {  
            temp1[i][j] = temp[i][j];  
        }  
    }  
}  
  
int main()  
{  
    int sum;  
    cout<<"请输入你的矩阵边数："<<endl;  
    cin>>sum;  
    cout<<"请输入各边的关系："<<endl;  
    for(int i=1;i<=sum;i++)  
    {  
        for(int j=1;j<=sum;j++)  
        {  
            cin>>temp1[i][j];  
        }  
    }  
    exchange(sum,1);  
    cout<<"R（1）的传递闭包关系图如下："<<endl;  
    for(int i=1;i<=sum;i++)  
    {  
        int ff=0;  
        for(int j=1;j<=sum;j++)  
        {  
            if(ff++!=0)  
              cout<<' ';  
            cout<<temp1[i][j];  
        }  
        cout<<endl;  
    }  
        exchange(sum,2);  
   cout<<"R（2）的传递闭包关系图如下："<<endl;  
    for(int i=1;i<=sum;i++)  
    {  
        int ff=0;  
        for(int j=1;j<=sum;j++)  
        {  
            if(ff++!=0)  
              cout<<' ';  
            cout<<temp1[i][j];  
        }  
        cout<<endl;  
    }  
        exchange(sum,3);  
   cout<<"R（3）的传递闭包关系图如下："<<endl;  
    for(int i=1;i<=sum;i++)  
    {  
        int ff=0;  
        for(int j=1;j<=sum;j++)  
        {  
            if(ff++!=0)  
              cout<<' ';  
            cout<<temp1[i][j];  
        }  
        cout<<endl;  
    }  
        exchange(sum,4);  
     cout<<"R（4）的传递闭包关系图如下："<<endl;  
    for(int i=1;i<=sum;i++)  
    {  
        int ff=0;  
        for(int j=1;j<=sum;j++)  
        {  
            if(ff++!=0)  
              cout<<' ';  
            cout<<temp1[i][j];  
        }  
        cout<<endl;  
    }  
    return 0;  
}  

 

 

 

菜鸟的总结：在计算机里面计算矩阵的话，or 和 | 和 || ，and和 & 和 &&是一样的。想知道为什么吗？ 哈哈，去看一下，位操作那一张呗！

 

运行结果如下：

请输入你的矩阵边数：
4
请输入各边的关系：
0 1 0 0
0 0 1 1
1 1 0 1
1 0 0 0
R（1）的传递闭包关系图如下：
0 1 0 0
0 0 1 1
1 1 0 1
1 1 0 0
R（2）的传递闭包关系图如下：
0 1 1 1
0 0 1 1
1 1 1 1
1 1 1 1
R（3）的传递闭包关系图如下：
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1
R（4）的传递闭包关系图如下：
1 1 1 1
1 1 1 1
1 1 1 1
1 1 1 1

--------------------------------
Process exited after 18.12 seconds with return value 0
请按任意键继续. . .