码迷,mamicode.com
首页 > 其他好文 > 详细

【POJ1681】Painter's Problem 高斯消元,求最小∑系数的异或方程组

时间:2015-02-04 13:04:36      阅读:118      评论:0      收藏:0      [点我收藏+]

标签:poj1681   painters problem   高斯消元   求最小系数   异或方程组   

#include <stdio.h>
int main()
{
	puts("转载请注明出处[vmurder]谢谢");
	puts("网址:blog.csdn.net/vmurder/article/details/43483547");
}


-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+

题意:

多组数据、

有个n*n的正方形,然后你要对某些位置进行操作,使得最后灯的状态都变成y。

操作:这个灯位置的上下左右以及自己这五盏灯状态都取反。


然后求最小操作次数。

-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+

先说我自己WA了的原因:

给正方形上的点id[i][j]=++cnt赋hash,然后没清id,WA了几遍~

还自己拍极限(永远n=15),还一直过。


嗯,反正分析我代码就懂了。(ctrl+F搜索"id[")

-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+

然后说正经题解:

一些自由元乱取0/1 AC的都是因为数据弱,不要乱提。

首先每个位置都可以选择“操作”或者“不操作”


然后对于自由元我们可以随意选择它开还是不开,而非自由元我们则需要回代一遍出解。

然后深搜一遍就好了~


这里我们不能给所有自由元都取“不操作”,因为这样看起来少了好多次操作数,

但是实际上,这些自由元是参与非自由元求值的回代的,

所以一个自由元的不作为可能导致很多个非自由元的悲剧。

-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+

代码:

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 300
#define inf 0x3f3f3f3f
using namespace std;
const int dx[]={0,0,0,1,-1};
const int dy[]={0,1,-1,0,0};

bool a[N][N],x[N];
int crs[N],n,ans;

void dfs(int id,int which,int now)
{
	if(id<1)
	{
		ans=min(ans,now);
		return ;
	}
	if(now>=ans)return ;

	if(crs[id]==which)
	{
		bool ret=a[id][n+1];
		for(int i=which+1;i<=n;i++)ret^=(a[id][i]*x[i]);
		x[which]=ret;
		dfs(id-1,which-1,now+ret);
	}
	else {
		x[which]=true;
		dfs(id,which-1,now+1);
		x[which]=false;
		dfs(id,which-1,now);
	}
	return ;
}

int Gauss(int n,int m)
{
	int i,j,k,id;
	for(id=i=1;i<n;i++,id++)
	{
		for(j=id;j<=m&&!a[j][i];j++);
		if(j>m){id--;continue;}

		crs[id]=i;
		if(id!=j)for(k=i;k<=n;k++)swap(a[id][k],a[j][k]);
		for(j=id+1;j<=m;j++)if(a[j][i])for(k=i;k<=n;k++)
			a[j][k]^=a[id][k];
	}
	for(i=id;i<=m;i++)if(a[i][n])return -1;
	return id-1;
}

int id[N][N],cnt;
char src[N];

void init()
{
	cnt=0,ans=inf;
	memset(a,0,sizeof a);
	memset(x,0,sizeof x);
	memset(id,0,sizeof id);
	memset(crs,0,sizeof crs);
}

int main()
{
	freopen("test.in","r",stdin);

	int i,j,k,g,t;
	for(scanf("%d",&g);g--;)
	{
		init();

		scanf("%d",&n);
		for(i=1;i<=n;i++)for(j=1;j<=n;j++)id[i][j]=++cnt;
		for(i=1;i<=n;i++)
		{
			scanf("%s",src+1);
			for(j=1;j<=n;j++)if(src[j]=='w')
				a[id[i][j]][cnt+1]=true;
		}
		for(i=1;i<=n;i++)for(j=1;j<=n;j++)for(k=0;k<=4;k++)
			if(t=id[i+dx[k]][j+dy[k]])
				a[id[i][j]][t]=true;
		n*=n;
		t=Gauss(n+1,n);
		if(t==-1)puts("inf");
		else {
			dfs(t,n,0);
			printf("%d\n",ans);
		}

	}
	return 0;
}



【POJ1681】Painter's Problem 高斯消元,求最小∑系数的异或方程组

标签:poj1681   painters problem   高斯消元   求最小系数   异或方程组   

原文地址:http://blog.csdn.net/vmurder/article/details/43483547

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!