poj 1830 开关灯问题(高斯消元法)

题目分析:

注意:有个地方要反过来

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;

const int maxn=110;
int a[maxn][maxn];
int n,a1[maxn],a2[maxn];
inline int fab(int x)
{
	return x>0?x:-x;
}
inline void swap(int &x,int &y)
{
	int temp=x;
	x=y;
	y=temp;       
}
int Gauss()
{
	int i,j,k,max_r,t;
	for(i=0,j=0;i<n&&j<n;i++,j++)
	{  
		max_r=i;
		for(k=i+1;k<n;k++)
			if(fab(a[k][j])>fab(a[max_r][j]))
				max_r=k;
		if(max_r!=i)
		{
			for(t=j;t<=n;t++)
				swap(a[i][t],a[max_r][t]);
		}
		if(a[i][j]==0)
		{
			i--;
			continue;
		}
		for(k=i+1;k<n;k++)
		{
			/*
			if(a[k][j]!=0)
			{
			for(t=j;t<n+1;t++)
			a[k][t]=a[k][t]^a[i][t];
			}*/
			if(a[k][j]==0)
				continue;
			for(t=j;t<=n;t++)
				a[k][t]=a[k][t]^a[i][t];             }
	}
	for(k=i;k<n;k++)//无解
		if(a[k][n]!=0)
			return -1;
	if(i==n)//唯一解
		return 1;
	return 1<<(n-i);//有n-i个未知变量
}
int main()
{
	int k,i;//int k,n,i;
	scanf("%d",&k);
	while(k--)
	{
		scanf("%d",&n);
		memset(a,0,sizeof(a));
		for(i=0;i<n;i++)
			scanf("%d",&a1[i]);
		for(i=0;i<n;i++)
			scanf("%d",&a2[i]);
		int y,z;
		while(scanf("%d %d",&y,&z)&&(y||z))
		{
			a[z-1][y-1]=1;//a[y-1][z-1]=1;会错滴!!!
		}
		for(i=0;i<n;i++)
		{            
			a[i][i]=1;
			a[i][n]=a1[i]^a2[i];//why???
			/*if(a1[i]!=a2[i])
			a[i][n]=1;
			else
			a[i][n]=0;*/
		}
		int ans=Gauss();      
		if(ans==-1)
			printf("Oh,it's impossible~!!\n");
		else
			printf("%d\n",ans);  
	}
	system("pause");
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值