poj1830（高斯消元）

地址：http://poj.org/problem?id=1830

开关问题

Time Limit: 1000MS		Memory Limit: 30000K
Total Submissions: 5315		Accepted: 1971

Description

有N个相同的开关，每个开关都与某些开关有着联系，每当你打开或者关闭某个开关的时候，其他的与此开关相关联的开关也会相应地发生变化，即这些相联系的开关的状态如果原来为开就变为关，如果为关就变为开。你的目标是经过若干次开关操作后使得最后N个开关达到一个特定的状态。对于任意一个开关，最多只能进行一次开关操作。你的任务是，计算有多少种可以达到指定状态的方法。（不计开关操作的顺序）

Input

输入第一行有一个数K，表示以下有K组测试数据。 
每组测试数据的格式如下： 
第一行 一个数N（0 < N < 29） 
第二行 N个0或者1的数，表示开始时N个开关状态。 
第三行 N个0或者1的数，表示操作结束后N个开关的状态。 
接下来 每行两个数I J，表示如果操作第 I 个开关，第J个开关的状态也会变化。每组数据以 0 0 结束。 

Output

如果有可行方法，输出总数，否则输出“Oh,it's impossible~!!” 不包括引号

Sample Input

2
3
0 0 0
1 1 1
1 2
1 3
2 1
2 3
3 1
3 2
0 0
3
0 0 0
1 0 1
1 2
2 1
0 0

Sample Output

4
Oh,it's impossible~!!

Hint

第一组数据的说明： 
一共以下四种方法： 
操作开关1 
操作开关2 
操作开关3 
操作开关1、2、3 （不记顺序） 

思路：用高斯消元，组数通过自由变元求出，至于为什么是自由变元，请百度（查出来为什么给我说下）。

代码：

//根据我第一道高斯消元代码改过来的，有点乱

#include<iostream>
#include<queue>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
#define LL __int64
#define col 30
int num[40][40],ans[40],n;
int getans(){
    int s=0,i=0;
    for(int m=0;i<n,m<n;m++){
        int k=i;
        for(int j=i+1;j<n;j++){
            if(abs(num[j][m])>abs(num[k][m]))
                k=j;
        }
        if(k!=i){
            for(int j=m;j<=n;j++)
                swap(num[i][j],num[k][j]);
        }
        if(!num[i][m]) continue;
        for(int j=i+1;j<n;j++){
            if(!num[j][m]) continue;
            for(k=m;k<=n;k++)
                num[j][k]^=num[i][k];
        }
        i++;
    }
    for(int l=i;l<n;l++){
        if(num[l][n]) return -1;  //无解
    }
    if(i<n) return n-i;  //多组解
    return 0;  //唯一解
}
int main(){
    int t,a[50],b[50],ai,bi;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        memset(num,0,sizeof(num));
        for(int i=0;i<n;i++)
            scanf("%d",&a[i]);
        for(int i=0;i<n;i++)
            scanf("%d",&b[i]);
        for(int i=0;i<n;i++){
            num[i][n]=a[i]^b[i];  //异或的值是答案
            num[i][i]=1;  //注意本身
        }
        while(scanf("%d%d",&ai,&bi)>0,ai|bi)
            num[bi-1][ai-1]=1;
        int ss=getans();
        if(ss<0) puts("Oh,it's impossible~!!");
        else printf("%d\n",(1<<ss));
    }
    return 0;
}