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

以下是使用C语言编写的关系矩阵及闭包算法代码，包括自反闭包、对称闭包和传递闭包的计算：

#include <stdio.h>
#include <stdlib.h>

#define MAX_SIZE 100

void print_matrix(int matrix[MAX_SIZE][MAX_SIZE], int size) {
    for (int i = 0; i < size; i++) {
        for (int j = 0; j < size; j++) {
            printf("%d ", matrix[i][j]);
        }
        printf("\n");
    }
}

void reflexive_closure(int matrix[MAX_SIZE][MAX_SIZE], int size) {
    for (int i = 0; i < size; i++) {
        matrix[i][i] = 1;
    }
}

void symmetric_closure(int matrix[MAX_SIZE][MAX_SIZE], int size) {
    int temp_matrix[MAX_SIZE][MAX_SIZE];
    for (int i = 0; i < size; i++) {
        for (int j = 0; j < size; j++) {
            temp_matrix[i][j] = matrix[i][j] || matrix[j][i];
        }
    }
    for (int i = 0; i < size; i++) {
        for (int j = 0; j < size; j++) {
            matrix[i][j] = temp_matrix[i][j];
        }
    }
}

void transitive_closure(int matrix[MAX_SIZE][MAX_SIZE], int size) {
    int temp_matrix[MAX_SIZE][MAX_SIZE];
    for (int i = 0; i < size; i++) {
        for (int j = 0; j < size; j++) {
            temp_matrix[i][j] = matrix[i][j];
        }
    }
    for (int k = 0; k < size; k++) {
        for (int i = 0; i < size; i++) {
            for (int j = 0; j < size; j++) {
                matrix[i][j] = matrix[i][j] || (temp_matrix[i][k] && temp_matrix[k][j]);
            }
        }
    }
}

int main() {
    int size;
    printf("Enter the size of the matrix: ");
    scanf("%d", &size);
    int matrix[MAX_SIZE][MAX_SIZE];
    printf("Enter the elements of the matrix: \n");
    for (int i = 0; i < size; i++) {
        for (int j = 0; j < size; j++) {
            scanf("%d", &matrix[i][j]);
        }
    }
    printf("Original matrix:\n");
    print_matrix(matrix, size);
    reflexive_closure(matrix, size);
    printf("Reflexive closure:\n");
    print_matrix(matrix, size);
    symmetric_closure(matrix, size);
    printf("Symmetric closure:\n");
    print_matrix(matrix, size);
    transitive_closure(matrix, size);
    printf("Transitive closure:\n");
    print_matrix(matrix, size);
    return 0;
}