【POJ 1681】 手机游戏 ---高斯消元

原创 2013年12月03日 17:51:32

Description

有一个有趣的手机游戏,有n*n个正方形的小按钮,有的按钮是黄色,有的按钮是白色。玩家的任务是通过点击按钮,让所有按钮都变成黄色,点按钮的次数越少,得分越高。 
但是按钮有个奇怪的特点,当你点击了坐标为(x,y)的按钮后,坐标为(x,y),(x+1,y),(x-1,y),(x,y-1),(x,y+1)的五个按钮会同时改变自身的颜色,是白色的变成黄色,黄色的变成白色。完成游戏最少需要点击多少次按钮呢?请找出答案。 
 

Input

第一行,一个整数n,表示有n*n个按钮。(1<=n<=40) 
接下来是一个由小写字母'y'和字符'w'构成的n*n的矩阵,描述了游戏开始时的情景。'y'表示该按钮式黄色,'w'表示该按钮式白色。

Output

如果能够完成游戏,输出一个整数,表示所需最小点击次数。 
如果无法完成游戏,输出"inf"

Sample Input

样例输入1:
2
yw
ww

样例输入2:

5
wwwww
wwwww
wwwww
wwwww
wwwww

Sample Output

样例输出1:
1

样例输出2:
15


【分析】

大概大部分的人看到这道题都认为是搜索。但仔细分析后我们知道,宽搜的话,状态无法存储。而深搜的话,状态数过多,会超时。仔细分析后得知,可以用高斯消元法解异或方程组

我们可以得到这样一个方程组:

v[1][1]*x[1] xor v[1][2]*x[2] xor v[1][3]*x[3] ...... xor v[1][n*n] =b[1]

v[2][1]*x[1] xor v[2][2]*x[2] xor v[2][3]*x[3] ...... xor v[2][n*n] =b[2]

......

v[n*n][1]*x[1] xor v[n*n][2]*x[2] xor v[n*n][3]*x[3] ...... xor v[n*n][n*n] =b[n*n]

其中,v[i][j]为1表示i,j这两个位置可以互相影响,为0表示不能互相影响。x[i]表示解。如果i这个位置需要被改变,b[i]为1,否则,b[i]为0;(我给坐标(x,y)标号为((x-1)*n+y))。

然后直接高斯消元法求出解的个数,即为答案。

现在来分析时间复杂度,原图为n×n,而我们高斯消元的增广矩阵为(n×n)×(n×n),我们记N=n×n。高斯消元的时间复杂度为O(N^3)。那么总的时间复杂度便为O(n*6)。

n的最大值为40,显然是会超时的。但是正确的做法就是这样,这又是为什么呢。

我们关注一下这道题的特点,与i位置关联的点最多只有4个,也就是说,在增广矩阵中,i这一列最多只有5个一(有一个是自己的一)。而高斯消元的主要代价在于,每次消元需要O(N)消去一行,要消O(N)次元。而现在我们只需要消最多4次,每次依然是O(N)。所以,之前的那个时间上界是相当松的。

新的时间复杂度为O(4*N^2),换为n便是O(4*n^4),是稳稳地能过的。


【代码】

/*
    ID:Ciocio
	LANG:C++
	DATE:2013-12-03
	TASK:NKOJ-1987
*/
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>

using namespace std;

#define MAXN 45
#define loc(x,y) ((x-1)*n+y)

int n,N;
int v[MAXN*MAXN][MAXN*MAXN],dir[4][2]={{-1,0},{1,0},{0,1},{0,-1}};

char _read(){char x;for(x=getchar();x!='y'&&x!='w';x=getchar());return x;}

void _init()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			v[loc(i,j)][n*n+1]=(_read()=='y')?0:1;
			v[loc(i,j)][loc(i,j)]=1;
		}
}

void _Guass()
{
	int i,j,k;
	for(i=j=1;i<=N&&j<=N;i++,j++)
	{
		for(k=i;v[k][j]==0&&k<=N;k++);
		if(k==N+1){i--;continue;}
		if(k!=i)
			for(int p=j;p<=N+1;p++)
				swap(v[i][p],v[k][p]);
		for(k=i+1;k<=N;k++)
			if(v[k][j])
			    for(int p=j;p<=N+1;p++)
				    v[k][p]^=v[i][p];
	}
	for(;i<=N;i++)
		if(v[i][N+1])
		{
			printf("inf\n");
			return;
		}
	int ans=0;
	for(i=N-1;i;i--)
		for(j=i+1;j<=N;j++)
			v[i][N+1]^=v[i][j]*v[j][N+1];
	for(int i=1;i<=N;i++)
		ans+=v[i][N+1];
	cout<<ans<<endl;
}

void _solve()
{
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			for(int k=0;k<=3;k++)
				if(1<=i+dir[k][0]&&i+dir[k][0]<=n)
					if(1<=j+dir[k][1]&&j+dir[k][1]<=n)
						v[loc(i,j)][loc(i+dir[k][0],j+dir[k][1])]=1;
	N=n*n;
	_Guass();
}

int main()
{
	_init();
	_solve();
	return 0;
}


相关文章推荐

【高斯消元】poj 1681 Painter's Problem

http://poj.org/problem?id=1681 高斯消元求最小步数,二进制枚举0-1变元,统计和的最小值 /* poj 1222 高斯消元-枚举自由变元 题意: 0-...

poj1681 高斯消元+dfs枚举

wa  N多次后,终于把高斯消元的模板给弄出来了,并且精简了许多。 用dfs枚举变元比二进制枚举变元简单多了。 (也终于明白 变元代表只得解的个数。) 上模板。 #include #include...

poj 1681 Painter's Problem(高斯消元)

http://poj.org/problem?id=1681 求最少经过的步数使得输入的矩阵全变为y。 思路:高斯消元求出自由变元,然后枚举自由变元,求出最优值。 注意根据自由变...

高斯消元 poj-1222、1681、1830、3185,sicily-1048

这几个题目都是翻转之后关联其他点的翻转,并且都是采用的01翻转,实属良心。 另外建立增广矩阵真的很费劲 poj 1222 #include #include #include  #in...

poj 1681 Painter's Problem(高斯消元 开关问题 )

Language: Default Painter's Problem Time Limit: 1000MS   Memory Limit: 10000K To...

poj 1681 Painter's Problem 高斯消元 枚举自由变元

Painter's Problem Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 5598 ...

NKOJ 1987 手机游戏(高斯消元)

手机游戏 问题描述有一个有趣的手机游戏,有n*n个正方形的小按钮,有的按钮是黄色,有的按钮是白色。玩家的任务是通过点击按钮,让所有按钮都变成黄色,点按钮的次数越少,得分越高。 但是按钮有个奇怪的特点,...

pku 1681 Painter's Problem(高斯消元)

Painter's Problem Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 3462 ...

poj1681 Painter's Problem 高斯消元

Language: Default Painter's Problem Time Limit: 1000MS   Memory Limit: 10000K To...

poj 1681 Painter's Problem 高斯消元

还是经典的高斯消元,和poj1753其实一模一样。要枚举自由变元得到最优解,不过这题不枚举自由变元使得自由变元都为0也能过,事实上必须要枚举的,数据太弱了~像1753不枚举自由变元的话就过不了。 我...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【POJ 1681】 手机游戏 ---高斯消元
举报原因:
原因补充:

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