【高斯消元】【poj 1830】开关问题

最新推荐文章于 2019-12-09 19:29:44 发布

原创最新推荐文章于 2019-12-09 19:29:44 发布 · 600 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

代码专栏收录该内容

22 篇文章

订阅专栏

本文提供了一道来自POJ编号为1222的问题的解题思路及C++实现代码，主要使用高斯消元法解决线性方程组问题。

http://poj.org/problem?id=1222

最近智商出了点问题，被这道题卡了，不和上一道一样么？

//#define _TEST _TEST
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
/************************************************
Code By willinglive    Blog:http://willinglive.cf
************************************************/
#define rep(i,l,r) for(int i=l,___t=(r);i<=___t;i++)
#define per(i,r,l) for(int i=r,___t=(l);i>=___t;i--)
#define MS(arr,x) memset(arr,x,sizeof(arr))
#define LL long long
#define INE(i,u,e) for(int i=head[u];~i;i=e[i].next)
inline const int read()
{int r=0,k=1;char c=getchar();for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1;
for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0';return k*r;}
/////////////////////////////////////////////////
int n;
int s[40],e[40];
int a[40][40];
/////////////////////////////////////////////////
int gauss()
{
	int i,j;
	for(i=1;i<=n&&j<=n;j++)
	{
		int k=i;
		for(;k<=n;k++) if(a[k][i]) break;
		if(a[k][j])
		{
			if(i^k) rep(j,i,n) swap(a[i][j],a[k][j]);
			rep(k,i+1,n) if(a[k][i])
			    rep(l,1,n+1) a[k][l]^=a[i][l];
			i++;
		}
	}
	for(j=i;j<=n;j++) if(a[j][n+1]) return -1;
	return 1<<(n-i+1);
}
/////////////////////////////////////////////////
void input()
{
	MS(a,0);
    n=read();
    rep(i,1,n) s[i]=read();
    rep(i,1,n) e[i]=read();
    rep(i,1,n) a[i][n+1]=s[i]^e[i],a[i][i]=1;
    int i,j;
    while(i=read(),j=read(),i+j) a[j][i]=1;
}
void solve()
{
    int res=gauss();
    if(res==-1) puts("Oh,it's impossible~!!");
    else printf("%d\n",res);
}
/////////////////////////////////////////////////
int main()
{
    #ifndef _TEST
    freopen("std.in","r",stdin); freopen("std.out","w",stdout);
    #endif
    per(i,read(),1)
    input(),solve();
    return 0;
}