poj 1830 开关问题（高斯消元）

最新推荐文章于 2022-08-10 10:15:09 发布

clover_hxy

最新推荐文章于 2022-08-10 10:15:09 发布

阅读量288

点赞数

分类专栏：高斯消元

本文链接：https://blog.csdn.net/clover_hxy/article/details/54141637

版权

高斯消元专栏收录该内容

23 篇文章 0 订阅

订阅专栏

开关问题

Time Limit: 1000MS		Memory Limit: 30000K
Total Submissions: 8070		Accepted: 3192

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

Sample Output

4
Oh,it's impossible~!!

Hint

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

Source

LIANGLIANG@POJ

[Submit] [Go Back] [Status] [Discuss]

题解：高斯消元

方案数就是2^(自由元的数量）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 1000
#define inf 1000000000
using namespace std;
int n,m;
int a[N][N],b[N],ans[N];
int s[N],t[N];
long long tot;
void gauss()
{
	tot=1;
	for (int i=1;i<=n;i++){
		int num=i;
		for (int j=i;j<=n;j++)
		 if (a[j][i]) {
		 	num=j;
		 	break;
		 }
		if (!a[num][i]) {
			tot*=2;
		    continue;
		}
		for (int j=1;j<=n;j++) swap(a[i][j],a[num][j]);
		swap(b[i],b[num]);
		for (int j=1;j<=n;j++)
		if (j!=i&&a[j][i]){
			for (int k=1;k<=n;k++)
			 a[j][k]^=a[i][k];
			b[j]^=b[i];
		}
    }
	for (int i=n;i>=1;i--){
		if (b[i]!=0&&a[i][i]==0){
			tot=inf;
			return;
		}
	}
}
int main()
{
	freopen("a.in","r",stdin);
//	freopen("my.out","w",stdout);
	int T; scanf("%d",&T);
	for (int t=1;t<=T;t++) {
		scanf("%d",&n); 
		memset(a,0,sizeof(a));
		memset(b,0,sizeof(b));
		for (int i=1;i<=n;i++) scanf("%d",&b[i]);
		for (int i=1;i<=n;i++){
			int x; scanf("%d",&x);
			b[i]^=x; a[i][i]=1;
		}
		//for (int i=1;i<=n;i++) cout<<b[i]<<" ";
	//	cout<<endl;
		int x,y;
		while (true) {
			scanf("%d%d",&x,&y);
			if (x==0&&y==0) break;
			a[y][x]=1;
		}
		gauss();
		if (tot==inf) printf("Oh,it's impossible~!!\n");
		else printf("%I64d\n",tot);
	}
}