Warshall算法求传递闭包

算法描述：
传递闭包的一种有效算法—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)，否则停止。
例如：



代码：

#include <stdio.h>
#define N 4    //宏定义 
int get_matrix(int a[N][N])
{
    int i = 0,j = 0;
    for (i = 0;i < N;i++) 
    {
        for (j = 0;j < N;j++)
         {
            scanf("%d",&a[i][j]);
            if (a[i][j] != 0 && a[i][j] != 1)
                return 1;            
         }
    }
    return 0;
}
int output_matrix(int a[N][N])
{
    int i = 0,j = 0;
    for (i = 0;i < N;i++) {
        for (j = 0;j < N;j++) {
            printf("%d  ",a[i][j]);
        }
        putchar('\n');
    }
}
int warshall(int a[][N])
{    
    //(1)i＝1；
    //(2)对所有j如果a[j，i]＝1，则对k＝0，1，…，n-1，a[j，k]＝a[j，k]∨a[i，k]；
    //(3)i加1；
    //(4)如果i<n，则转到步骤2，否则停止
    int i = 0;
    int j = 0;
    int k = 0;
    for (i = 0;i < N;i++) 
    {                  
        for (j = 0;j < N;j++) 
        {
            if (a[j][i]) 
            {
                for (k = 0;k < N;k++) 
                {
                    a[j][k] = a[j][k] | a[i][k];//逻辑加 
                }
            } 
        }
    }
}
int main()
{
    int a[N][N] = {0};
    printf("please input a matrix with %d * %d:\n",N,N);
    if (get_matrix(a))
     {  printf("Get matrix error!Only 0 or 1 in matrix!\n");
        return 1; //错误返回主函数，返回值为1； 
     }
    warshall(a);
    output_matrix(a);

    return 0;//正常返回主函数，返回值为0； 
}