离散数学自反闭包对称闭包传递闭包C语言代码实现

想要用程序实现，首先要明白自反闭包等的概念，自反闭包是主对角线上的数字为零，对称闭包是矩阵关于主对角线是对称的，传递闭包是找最短路径（传递闭包个人不是很理解，我是根据书上边的warshall算法写的代码）

下面直接上代码：

#include<stdio.h>
#define N 100//符号变量
int a[N][N],b[N][N],c[N][N];//初始化三个二维数组
int main()
{
	int n;//矩阵行和列都是n行
	int i,j,k;
	scanf("%d",&n);//输入n值
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			scanf("%d",&a[i][j]);//依次输入矩阵的每一个数
		}
	}
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			b[i][j]=a[i][j];//将矩阵a的值赋值给矩阵b 
			
		}
	}
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			if(b[i][j]==1){
				b[j][i]=1;//由对称性可得b[i][j]==b[j][i]==1
			}
		}
		}
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			c[i][j]=a[i][j];//将矩阵a的值赋值给b
			if(i==j){
				c[i][j]=1;//自反闭包主对角线上的数为1
			}
		}
	}
	for(k=0;k<n;k++)
	{
		for(i=0;i<n;i++)
		{
			for(j=0;j<n;j++)
			{
				a[i][j]=a[i][j]||(a[i][k]&&a[k][j]);//传递闭包处理第j行+第i行送回第j行
			}
		}
	}
	printf("r(R):\n");
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
		    printf("%d ",c[i][j]);//输出自反闭包
		}
		printf("\n");//每行输出完成换行
	}
	printf("s(R):\n");
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
		    printf("%d ",b[i][j]);//输出对称闭包
		}
		printf("\n");//每行输出完成换行
	}
	printf("t(R):\n");
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			printf("%d ",a[i][j]);
		}
		printf("\n");//每行输出完成换行
	}
	return 0;
}

运行结果：