题目:http://poj.org/problem?id=1830
学完高斯消元做的第一题, 高斯消元的难点在于怎样建立方程组,来建立系数矩阵。如何建方程组:http://blog.csdn.net/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);
}
}