3422. 【NOIP2013模拟】水叮当的舞步 (Standard IO)

 

Time Limits: 1000 ms  Memory Limits: 262144 KB  Detailed Limits  

Description

水叮当得到了一块五颜六色的格子形地毯作为生日礼物，更加特别的是，地毯上格子的颜色还能随着踩踏而改变。

为了讨好她的偶像虹猫，水叮当决定在地毯上跳一支轻盈的舞来卖萌~~~

地毯上的格子有N行N列，每个格子用一个0~5之间的数字代表它的颜色。

水叮当可以随意选择一个0~5之间的颜色，然后轻轻地跳动一步，地毯左上角的格子所在的联通块里的所有格子就会变成她选择的那种颜色。这里连通定义为：两个格子有公共边，并且颜色相同。

由于水叮当是施展轻功来跳舞的，为了不消耗过多的真气，她想知道最少要多少步才能把所有格子的颜色变成一样的。

Input

每个测试点包含多组数据。

每组数据的第一行是一个整数N，表示地摊上的格子有N行N列。

接下来一个N*N的矩阵，矩阵中的每个数都在0~5之间，描述了每个格子的颜色。

N=0代表输入的结束。

Output

对于每组数据，输出一个整数，表示最少步数。

Sample Input

2

0 0 

0 0

3

0 1 2

1 1 2

2 2 1

0

Sample Output

0

3

Data Constraint

对于30%的数据，N<=5

对于50%的数据，N<=6

对于70%的数据，N<=7

对于100%的数据，N<=8，每个测试点不多于20组数据。

Source / Author: floodit

题解：

dfs暴力 + A*(其实就是一个小剪枝) + 矩阵优化（扯远了 ， 其实就是记录边界）

 

矩阵优化：

我们引入一个N*N的v数组。

左上角的格子所在的联通块里的格子标记为1。

左上角联通块周围一圈格子标记为2，其它格子标记为0。

如果某次选择了颜色c，我们只需要找出标记为2并且颜色为c的格子，向四周扩展，并相应地修改v标记，就可以不断扩大标记为1的区域。

最终如果所有格子标记都是1，那么显然找到了答案。

 

#include<bits/stdc++.h>
#define mem(a,b) memset(a,b,sizeof(a))
#define mcy(a,b) memcpy(a,b,sizeof(a))
#define ll long long
#define re register
#define inf 2000000000
#define N 10
#define open(x) freopen(x".in","r",stdin);freopen(x".out","w",stdout)
using namespace std;

const int fx[4][2] = {{1,0},{0,1},{-1,0},{0,-1}};

int n,i,j,t,ans,a[N][N];
int f[N][N],l[N*N*10][2],c[6],v[N][N];

int sp(int now)
{
	int base = now;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++) if(a[i][j]!=base)return 0;
	return 1;
}


int find(int c[])
{
	int cnt=0;
	for(re int i=0;i<6;i++) if(c[i])++cnt;
	return cnt;
}

void updata(int x,int y,int A)
{	
	//	
	v[x][y] = 1;
	c[A] --;
	for(int k=0;k<4;k++)
	{
		int xx = x+fx[k][0] ,yy = y + fx[k][1];
		if(xx>0 && yy>0 && xx<=n && yy<=n)
		{
			if(v[xx][yy] == 0 && a[xx][yy] == A) 
				updata(xx,yy,A); 
			if(v[xx][yy] == 0) v[xx][yy] = 2;
		}
	}
	
}

int maxdep=0;
void dfs(int now,int dep)
{
	if(dep + find(c) > maxdep) return ;
	if(find(c)==0)
	{
		ans = dep;
		return ;
	}
	int C[6] , V[N][N];
	mcy(C,c);
	mcy(V,v);
	for(int i=0;i<6;i++)
	{
		int tag=0;
		for(int j=1;j<=n;j++)
			for(int k=1;k<=n;k++) if(v[j][k] == 2 && a[j][k]==i)
			
				 updata(j,k,i),v[j][k]=1,tag=1;
		if(tag)
		{
			dfs(i , dep+1 );
			mcy(c,C);
			mcy(v,V);
		}
	}
	return ;
}

void cinit()
{
	mem(c,0);
	for(int i=1;i<=n;i++)
	 	for(int j=1;j<=n;j++) ++c[a[i][j]];
}

void ID()
{
	maxdep = 0;
	while(ans == inf)
	{
		++maxdep;
		mem(v,0);
		cinit();
		updata(1,1,a[1][1] ); // v , l , c
		dfs(a[1][1] , 0 /*ans && dep*/ );
	}
	
	return ;
}


int main()
{
	while(scanf("%d",&n),n)
	{
		for(i=1;i<=n;i++)
			for(j=1;j<=n;j++) scanf("%d",&a[i][j]);
		if(sp(a[1][1]))
		{
			printf("0\n");
			continue;
		}
		ans  = inf;
		ID();
		
		printf("%d\n" , ans);
	}
	return 0;	
	
}

我后面会写迭代加深的博客。

时间复杂度()