01矩阵求解的个数

Posted 2021-12-03 thusloop

int Gauss()//0 1 矩阵求解的个数 
{
	int i, j, r, x, y;
	for(i=1,j=1;i<=n&&j<=n;++j)
	{
		r=i;
		while(r<=n&&!tmp[r][j])++r;//找到第j列第一个非零元素所在行
		if(tmp[r][j]) //若该列存在非零元素
		{
			swap(tmp[i], tmp[r]);//将非零元素所在的行交换到当前行第i行
			for(x=1;x<=n;++x) //除当前行第i行外，若其他n-1行中有第j列元素非零的行，利用第i行对该行进行消元（此题中即进行异或操作）
			{
				if(x!=i&&tmp[x][j])
				{
					for(y=i;y<=n+1;++y) tmp[x][y]^=tmp[i][y];
				}
			}
			++i;//继续处理下一行
		}
	}
	for(x=i;x<=n;++x) //若最后增广列没有全部变成0，则说明该增广矩阵无解
	{
		if(tmp[x][n+1]) return -1;
	}
	return n-i+1;//自由变量个数为n - i + 1个，解的个数就是1<<(n - i - 1)
}