「Poetize5」水叮当的舞步

Description

水叮当得到了一块五颜六色的格子形地毯作为生日礼物，更加特别的是，地毯上格子的颜色还能随着踩踏而改变。
为了讨好她的偶像虹猫，水叮当决定在地毯上跳一支轻盈的舞来卖萌~~~
地毯上的格子有N行N列，每个格子用一个0~5之间的数字代表它的颜色。
水叮当可以随意选择一个0~5之间的颜色，然后轻轻地跳动一步，左上角的格子所在的联通块里的所有格子就会变成她选择的那种颜色。这里连通定义为：两个格子有公共边，并且颜色相同。
由于水叮当是施展轻功来跳舞的，为了不消耗过多的真气，她想知道最少要多少步才能把所有格子的颜色变成一样的。

题解：

我来当搬运工。。。

类型：IDA* (迭代加深启发式搜索)

方法一：

枚举每次选取了哪种颜色，然后找出左上角的格子所在的联通块，改变颜色。

为了避免来回改变、搜索深度过大，采用迭代加深的dfs限制搜索步数。

迭代加深也就是，依次限制搜索深度为0、1、2、3……进行搜索，搜索过程中发现深度超过限制就马上退出。只要搜索成功就找到了答案，也可以立即退出。

期望得分：0~10分。

方法二：

加入一个小剪枝：如果改变颜色后，左上角格子所在的联通块大小没有改变，可以剪枝。这样可以避免来回往复地搜索。

期望得分：10~20分。

方法三：

采用IDA*算法，设计估价函数。可以发现如果当前矩阵中除了左上角的联通块之外，共有M种颜色，那么还需要的步数不小于M。因此如果当前搜索深度+估价函数的值>深度限制，可以回溯。

期望得分：50~70分。

方法四：

我们可以发现，每次寻找左上角的格子所在的联通块耗费的时间常数巨大。因此我们在这里寻求突破。

我们引入一个N*N的v数组。左上角的格子所在的联通块里的格子标记为1。左上角联通块周围一圈格子标记为2，其它格子标记为0。如果某次选择了颜色c，我们只需要找出标记为2并且颜色为c的格子，向四周扩展，并相应地修改v标记，就可以不断扩大标记为1的区域，最终如果所有格子标记都是1，那么显然找到了答案。

代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 500+100

 #define maxm 500+100

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define mod 1000000007

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }
 const int dx[]={-,,,},dy[]={,-,,};
 int a[][],v[][],f[],n,ID;
 inline int left()
 {
     int tmp=;
     memset(f,,sizeof(f));
     for1(i,n)for1(j,n)
     if(!f[a[i][j]]&&v[i][j]!=){f[a[i][j]]=;tmp++;}
     return tmp;
 }
 inline void dfs(int x,int y,int c)
 {
     v[x][y]=;
     for0(i,)
     {
         int xx=x+dx[i],yy=y+dy[i];
         if(xx<||xx>n||yy<||yy>n||v[xx][yy]==)continue;
         v[xx][yy]=;
         if(a[xx][yy]==c)dfs(xx,yy,c);
     }
 }
 inline int fill(int c)
 {
     int tmp=;
     for1(i,n)for1(j,n)
     if(a[i][j]==c&&v[i][j]==)
     {
         tmp++;
         dfs(i,j,c);
     }
     return tmp;
 }
 inline bool IDA(int dep)
 {
     int g=left();
     if(dep+g>ID)return ;
     if(!g)return ;
     int rec[][];
     for0(i,)
     {
         memcpy(rec,v,sizeof(v));
         if(fill(i)&&IDA(dep+))return ;
         memcpy(v,rec,sizeof(v));
     }
     return ;
 }

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);

     while(cin>>n&&n)
      {
          memset(a,,sizeof(a));
          memset(v,,sizeof(v));
          for1(i,n)for1(j,n)a[i][j]=read();
          dfs(,,a[][]);
          for(ID=;;ID++)
             if(IDA())break;
          printf("%d\n",ID);
     }

     return ;

 }