C语言实现关系的闭包运算

问题描述：
利用矩阵求解有限集上给定关系的自反和对称闭包。
输入格式：
首先输入关系矩阵R的维数，回车之后输入矩阵每个元素，以空格或回车分开。 只能输入0或1。
输出格式：
输出自反闭包关系矩阵r®,对称闭包关系矩阵s®。
样例输入：
4
0 1 0 0
0 0 0 1
0 0 0 0
1 0 1 0
样例输出
r®:
1 1 0 0
0 1 0 1
0 0 1 0
1 0 1 1
s®:
0 1 0 1
1 0 0 1
0 0 0 1
1 1 1 0
关于关系的定义不再赘述，下面分别给出其求解公式：
自反闭包：设R是集合A上的二元关系。那么r®=R∪E (E是A上相等关系) ；
对称闭包：设R是集合A上的二元关系。那么s®=R∪R^-1;
不难发现，自反闭包本质上就是在原矩阵的基础上将对角线元素全改为1，而对称闭包其实就是将原矩阵中为1的元素
a[i][j],其对应的a[j][i]也变为1，这是很容易实现的

#include<stdio.h>
#define N 100
int main()
{
	int n,i,j,a[N][N],b[N][N],c[N][N]={0};
	scanf("%d",&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++)
		{
			if(i==j) b[i][j]=1;
			else b[i][j]=a[i][j];
			if(a[i][j]==1)
			{
				c[i][j]=1;
				c[j][i]=1;
			}
		}
	}
	printf("r(R):\n");
	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
		{
			printf("%d ",b[i][j]);
		}
		printf("\n");
	}
	printf("s(R):\n");
	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
		{
			printf("%d ",c[i][j]);
		}
		if(i<n-1) printf("\n");
	}
}

下面介绍传递闭包，


利用定义求解传递闭包的博客单独写了一篇，可进行参看