蓝书（算法竞赛进阶指南）刷题记录——POJ1830 开关问题（高斯消元）

题目：POJ1830.
题目大意：给定一个长度为为$n$的初始$01$序列$a$，要求你把序列操作（给一些位置取反）成一个给定的$01$序列.给出一些关系$ij$，表示若$i$被操作则$j$也要被操作.
$1\le n<29$.

设每一个开关的初始状态为$b_i$，结束状态$e_i$，$a_{i,j}=1/0$表示第$i$个数会/不会被第$j$个数影响，$x_i$表示第$i$个点被操作几次.我们发现这个东西可以列成一个异或方程组：
$$\{\begin{array}{c}{a}_{1,1}{x}_{1}xor{a}_{1,2}{x}_{2}xor\cdots xor{a}_{1,n}x=b_{1}xor{e}_1\\ a_{2,1}{x}_{1}xor{a}_{2,2}{x}_{2}xor\cdots xor{a}_{2,n}x=b_{2}xor{e}_2\\ ⋮\\ a_{n,1}{x}_{1}xor{a}_{n,2}{x}_{2}xor\cdots xor{a}_{n,n}x=b_{n}xor{e}_n\end{array}$$

我们得到了一个异或方程组，发现异或方程组其实很容易解，把加减操作换成异或即可.而且我们发现这个方程中一行只有$30$个$0/1$，所以可以状态压缩.

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
  using namespace std;

#define Abigail inline void
typedef long long LL;

const int N=29;

int a[N+9],n,ans;

int Gauss(int *a,int n){
  for (int i=1;i<=n;++i){
  	for (int j=i+1;j<=n;++j)
  	  if (a[j]>a[i]) swap(a[j],a[i]);
  	if (a[i]==0) return n-i+1;
  	if (a[i]==1) return -1;
  	for (int k=n;k>=1;--k)
  	  if (a[i]>>k&1){
  	    for (int j=1;j<=n;++j)
  	      if (i^j&&a[j]>>k&1) a[j]^=a[i];
  	    break;
	  }
  }
  return 0;
}

Abigail into(){
  scanf("%d",&n);
  int x,y;
  for (int i=1;i<=n;++i)
  	scanf("%d",&a[i]);
  for (int i=1;i<=n;++i){
    scanf("%d",&x);
    a[i]^=x;
  }
  while (~scanf("%d%d",&x,&y)&&x+y)
  	a[y]|=1<<x;
}

Abigail work(){
  for (int i=1;i<=n;++i)
    a[i]|=1<<i;
  ans=Gauss(a,n);
}

Abigail outo(){
  if (ans==-1) puts("Oh,it's impossible~!!");
  else printf("%d\n",1<<ans);
}

int main(){
  int T;
  scanf("%d",&T);
  while (T--){
    into();
    work();
    outo();
  }
  return 0;
}