高斯消元解xor方程 poj1830

莫队长最喜欢的一类东西。。。

上次听赖神讲线性相关之前，莫队长说他讲的高斯消元最简单，一开始不信，后来赖神讲了之后。。。呃

其实赖神讲的也不是很难接受，但听他讲了之后现在却没什么印象了，果然跨度还是太大了。

poj1830 典型的xor方程组，每个开关根据与其他开关关系建立方程，其跟自己肯定相关，相关系数为1，否则为0

如共有4个开关 开关1与2、4相关，结果从0变为1，则建立方程 1*a1 xor 1*a2 xor 0*a3 xor 1*a4=1

系数矩阵（1,1,0,1,1）。建立总系数矩阵后进行消元，方案数=2^自由元数，因为自由元是不确定的，它们有0、1两种选择，

根据乘法原理，即可得出。由于不需输出方案，所以不需回带，若无解则是出现（0，0,0,0，.......，0,1）情况.

var a:array[0..29]of longint; t,n,ans:longint; procedure init; var i,j,k,l,x,y:longint; begin readln(n); for i:=1 to n do read(a[i]);readln; for i:=1 to n do begin read(x);a[i]:=(a[i] xor x)+1<<i end; while not(seekeof) do begin readln(x,y); if x+y=0 then break; a[y]:=a[y] or 1<<x end; ans:=0;j:=1; for i:=1 to n do begin if a[i]=1 then begin writeln('Oh,it''s impossible~!!');exit end; k:=j; while (a[k]>>i and 1=0)and(k<=n) do inc(k); if k<=n then begin a[0]:=a[k];a[k]:=a[j];a[j]:=a[0]; for l:=j+1 to n do if (a[l]>>i and 1=1) then a[l]:=a[l] xor a[j]; inc(j) end end; ans:=n-j+1; for i:=j to n do if a[j]=1 then begin writeln('Oh,it''s impossible~!!');exit end; writeln(1<<ans) end; begin assign(input,'1830.in');reset(input); readln(t); for t:=1 to t do init; close(input) end.