请用C语言写出以下算法的代码:首先将任意关系通过关系矩阵表示，关系的自反闭包只需将矩阵的主对角线元素全部置为1；关系的对称闭包即关系矩阵与其转置矩阵进行加法运算（逻辑加）；关系的传递闭包转化为矩阵的计算，计算公式为 Mt = M + M2 + M3 + …（逻辑加） 其中M为关系矩阵，Mt为传递闭包的矩阵。 关系的传递闭包也可用Floyd-Warshall算法计算，效率更高，Floyd-Warshall算法原理为： 设R的关系矩阵为n阶矩阵M，传递闭包矩阵为A = A[j, k]， (1)令矩阵A = M

这是一个将关系矩阵转化为传递闭包矩阵的代码实现，其中使用了 Floyd-Warshall 算法：

#include <stdio.h> 

#define MAX_NODES 100

int matrix[MAX_NODES][MAX_NODES];
int closure[MAX_NODES][MAX_NODES];
int num_nodes;

void find_transitive_closure() {
    int i, j, k;
    for (i = 0; i < num_nodes; i++) {
        for (j = 0; j < num_nodes; j++) {
            closure[i][j] = matrix[i][j];
        }
    }
    for (k = 0; k < num_nodes; k++) {
        for (i = 0; i < num_nodes; i++) {
            for (j = 0; j < num_nodes; j++) {
                closure[i][j] = closure[i][j] || (closure[i][k] && closure[k][j]);
            }
        }
    }
}

int main() {
    int i, j;
    // 读入关系矩阵
    printf("请输入节点个数：");
    scanf("%d", &num_nodes);
    printf("请输入关系矩阵：\n");
    for (i = 0; i < num_nodes; i++) {
        for (j = 0; j < num_nodes; j++) {
            scanf("%d", &matrix[i][j]);
        }
    }
    // 计算传递闭包矩阵
    find_transitive_closure();
    // 输出传递闭包矩阵
    printf("传递闭包矩阵：\n");
    for (i = 0; i < num_nodes; i++) {
        for (j = 0; j < num_nodes; j++) {
            printf("%d ", closure[i][j]);
        }
        printf("\n");
    }
    return 0;
}