设 a i , j a_{i,j} ai,j 表示按j后i是否会发生变化(1是0否), x i x_i xi 表示第i个开关是否按(1按0不按)
则可以得到 一个开关变化的过程最终结果为 sa ^ en 开始的状态异或 结束的状态。即你不管怎么按,这个开关经过层层抵消,最终要按1或者不按0.
而这个过程为
a
1
,
1
x
1
a_{1,1}x_1
a1,1x1 xor
a
1
,
2
x
2
a_{1,2}x_2
a1,2x2 xor … xor
a
1
,
n
x
n
a_{1,n}x_n
a1,nxn
… . . . . . . . . . . . . . . . . … . . . . . . . . .
我们可以构建增广矩阵,(1~n)位置对应1 - n 列
再用高斯消元,求出该系数矩阵的秩,即可以知道有多少个(cnt) x不确定,就是随便从0,1选,则方案就为 1<<cnt.
可以直接用二进制来压缩状态。
某些部分讲解不是很清楚。
代码:
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<iomanip>
#define int long long
#define sc scanf
#define pf printf
using namespace std;
typedef pair<int,int> pii;
typedef long long ll;
typedef unsigned long long LL;
const int INF = 0x3f3f3f3f;
const double eps = 1e-4;
const int mod = 10000;
const int N = 210;
int a[N];
signed main(){
// IOS;
#ifdef ddgo
freopen("C:/Users/asus/Desktop/ddgoin.txt","r",stdin);
#endif
int tt; cin>>tt;
while(tt --){
int n; cin>>n;
for(int i=1;i<=n;i++) cin>>a[i];
for(int i=1,j;i<=n;i++){
cin>>j;
a[i] ^= j; // 表示常数项的结果
a[i] |= (1<<i);//把a[i][i] 标记为1
}
int x,y;
while(cin>>x>>y && x &&y) a[y] |= (1<<x);//把a[y][x] 标记为1
int ans = 1;//特殊,满秩的情况
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){
ans = 1 << (n-i+1);
break;
}
if(a[i] == 1){//无解
ans = 0;
break;
}
for(int k=n;k;k--){// 把那列的1全部消去,异或消元
if(a[i] >> k & 1){
for(int j=1;j<=n;j++)
if(i != j && (a[j] >> k & 1)) a[j] ^= a[i];
break;
}
}
}
if(ans == 0) cout<<"Oh,it's impossible~!!"<<endl;
else cout<<ans<<endl;
}
return 0;
}