开关问题[POJ1830]

欢迎大家访问我的老师的OJ———caioj.cn

题面描述

传送门

思路

设$x_i$表示第i个开关的操作情况，$x_i=1$表示按了这个开关，$x_i=0$表示没按。再统计$a_{i,j}$表示第i个开关和第j个开关的联系情况，$a_{i,j}=1$表示按下j会影响i的状态，$a_{i,j}=0$表示不会影响，特别地，令$a_{i,i}=1$。

一个开关最后的状态$ds{t}_i$,取决于它最初的状态$sr{c}_i$,以及所有与它有关系的开关的操作情况执行异或运算得到的结果，可列出异或方程组:
$$\{\begin{array}{l}{a}_{1,1}{x}_1\mathrm{xor}{a}_{1,2}{x}_2\mathrm{xor}...\mathrm{xor}{a}_{1,n}{x}_n=sr{c}_1\mathrm{xor}ds{t}_1\\ a_{2,1}{x}_1\mathrm{xor}{a}_{2,2}{x}_2\mathrm{xor}...\mathrm{xor}{a}_{2,n}{x}_n=sr{c}_2\mathrm{xor}ds{t}_2\\ ⋮\\ a_{n,1}{x}_1\mathrm{xor}{a}_{n,2}{x}_2\mathrm{xor}...\mathrm{xor}{a}_{n,n}{x}_n=sr{c}_n\mathrm{xor}ds{t}_n\end{array}$$

异或其实就是不进位加法，我们可以通过高斯消元来解决。

在执行高斯消元的过程中，把加、减法替换成异或，且不需要执行乘法。

最终我们可以得到该异或方程组对应的简化阶梯形矩阵。若存在形如$0=1$的方程，则方程组无解。

否则因为自由元可以取$0$或$1$，所以方程解的数量就是$2^{自由元数量}$

由于偷懒省事也省时，我们可以状态压缩一下$a_{i,j}$，即:

int a[100];//0位存 开始xor结束 
a[i]^=j;
a[i]|=1<<i;//a[i][i]=1;
a[y]|=1<<x;//按x影响y ,a[y][x]=1

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
using namespace std;
int a[100];//0位存 开始xor结束 
int main()
{
	int t;scanf("%d",&t);
	while(t--)
	{
		int n;scanf("%d",&n);
		for(int i=1;i<=n;i++)scanf("%d",&a[i]);
		for(int i=1;i<=n;i++)
		{
			int j;scanf("%d",&j);
			a[i]^=j;
			a[i]|=1<<i;//a[i][i]=1;
		}
		int x,y;
		while(~scanf("%d%d",&x,&y)&&x&&y)
			a[y]|=1<<x;//按x影响y 
		int ans=1;
		for(int i=1;i<=n;i++)
		{
			//寻找最高位
			int x=i;
			for(int j=i+1;j<=n;j++)if(a[j]>a[x])x=j;
			swap(a[x],a[i]);
			//消元成功 
			if(a[i]==0){ans=1<<(n-i+1);break;}
			//系数全为0，但结果却为1. 
			if(a[i]==1){ans=0;break;}
			for(int j=n-i+1;j>=1;j--)
				if((a[i]>>j&1))//a[i][j]==1
				{
					for(int k=i+1;k<=n;k++)
						if(i!=k&&(a[k]>>j&1))//a[k][j]==1
							a[k]^=a[i];
					break;
				}
		}
		if(ans==0)puts("Oh,it's impossible~!!");
		else printf("%d\n",ans);
	}
	return 0;
}