PKU 1830 开关问题

 题目：http://poj.org/problem?id=1830

学完高斯消元做的第一题， 高斯消元的难点在于怎样建立方程组，来建立系数矩阵。如何建方程组：http://blog.youkuaiyun.com/shiren_Bod/article/details/5766907，写的很好, 学习了～～。

此题需要注意的是，初始化矩阵时，开关I和J不要弄反了～～～ 如果有自由变元，返回1<<(n-row)。

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
using namespace  std;
const int maxn=30;
int ca, n, a, b;
int begin[maxn+10], end[maxn+10], map[maxn+10][maxn+10];
int gauss()
{
	int row, col;
	for (col=0, row=0; row<n, col<n; col++, row++)
	{
		int maxr=row;
		for (int i=row; i<n; i++) if (map[i][col]!=0) {maxr=i; break;}
		if (maxr!=row) 
			for (int i=0; i<=n; i++) swap(map[row][i], map[maxr][i]);
		if (map[row][col]==0) {row--; continue;}
		for (int i=row+1; i<n; i++)
		{
			if (!map[i][col]) continue;
			for (int j=col; j<=n; j++) map[i][j]=map[i][j]^map[row][j];
		}
	}
	for (int i=row; i<n; i++) if (map[i][n]!=0) return -1;
	return (1<<(n-row));
}
int main()
{
	//freopen("in.txt", "r", stdin);
	scanf("%d", &ca);
	while (ca--)
	{
		scanf("%d", &n);
		memset(map, 0, sizeof(map));
		for (int i=0; i<n; i++) scanf("%d", &begin[i]);
		for (int i=0; i<n; i++) scanf("%d", &end[i]), map[i][n]=begin[i]^end[i];
		for (int i=0; i<n; i++) map[i][i]=1;
		while (scanf("%d %d", &a, &b)==2)
		{
			if (a==0&&b==0) break;
			map[b-1][a-1]=1;
		}
		int flag=gauss();
		if (flag==-1) printf("Oh,it's impossible~!!\n");
		else printf("%d\n", flag);
	}
}