POJ 1222 高斯消元法

原创 2012年03月30日 15:13:25

CXLOVE的第一篇博文,先来一个水题

开灯关灯问题,5*6的灯阵,将每一个位置上的状态看做一个变元,30个变元,列出30个异或方程,高斯消元解方程即可

N年前就做过这题,当时是位运算枚举+递推

从此cxlove的博客就正式开通了,欢迎大牛拍砖

/*
ID:cxlove
PROB:POJ 1222
HINT:高斯消元
*/
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
using namespace std;
int a[31][31];  //系数矩阵
int ans[31];
void gauss(){
	for(int i=0;i<30;i++){
		int k=i;
		for(;k<30;k++)
			if(a[k][i]!=0)
				break;
		for(int j=0;j<=30;j++)
			swap(a[i][j],a[k][j]);
		for(int j=0;j<30;j++)
			if(i!=j&&a[j][i])
				for(int k=0;k<=30;k++)
					a[j][k]=a[i][k]^a[j][k];
	}
}
int main(){
	int t,tt=0;
	scanf("%d",&t);
	while(t--){
		memset(a,0,sizeof(a));
		for(int i=0;i<30;i++){
			scanf("%d",&a[i][30]);
			ans[i]=0;
		}
		for(int i=0;i<30;i++){
			a[i][i]=1;
			if(i%6!=0)
				a[i-1][i]=1;
			if(i%6!=5)
				a[i+1][i]=1;
			if(i>5)
				a[i-6][i]=1;
			if(i<24)
				a[i+6][i]=1;
		}
		gauss();
		/*for(int i=0;i<30;i++){
		for(int j=0;j<=30;j++)
		printf("%d ",a[i][j]);
		printf("\n");
		}*/
		for(int j=0;j<30;j++)  
		     ans[j]=a[j][30];	
	        printf("PUZZLE #%d\n",++tt);
	        for(int i=0;i<30;i++){
		    printf("%d",ans[i]);
		    if(i%6==5)
			printf("\n");
		    else
			printf(" ");
                }
	}
	return 0;
}



版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

NOI 2017 场外酱油记

NOI2017

此博客停止更新

感谢大家的支持,自从退役后,上CSDN的次数越来越i

小记

2013年11月17日,长沙区域赛结束之后,吃翔爱就成了退役狗了。。。。。        作为比较弱的混了几年,所以也不需要什么退役贴了。。。        回想最后半年多,只能说自作孽不可活,不好好...

NOIP 2008 解题报告(笨小猴,火柴棒等式,传纸条,双栈排序)

NOIP2008

博弈-翻硬币游戏

转自:http://blog.sina.com.cn/s/blog_8f06da99010125ol.html 翻硬币游戏     一般的翻硬币游戏的规则是这样的:       N&...

HDU 3535 AreYouBusy 综合背包

转载请注明出处,谢谢 http://blog.csdn.net/ACM_cxlove?viewmode=contents           by---cxlove 非常综合的背包。 必须得有较好...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)