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； 
}