zoj 1361

又是一道用到位运算的题，继续熟悉了下位运算。。

/*
zoj_1361    搜索
不错的bfs题。
转一下一些acmer的解释：
1、宽搜是必要的，那么如何记录蛇当前的状态，以避免以后重复的访问变成了关键，
在这里我们将蛇的状态描述为如下三元组（x,y,state）,其中 (x,y)是蛇头的坐标，
state 记录的是身体的状态，由于身体最长为七段，每一段相对于前一段只有上下左右
四种状态，仅需要两位表示，则身体所有块的状态最多需要 7×2=14 位二进制数表示，
则我们可以构建矩阵 pMatrix[21][21][16384] 来保存访问过的状态，以避免重复访问
相同的状态，后面的工作就是传统的宽搜过程了；
2、如果仅按上面的方法来做的话会超时，于是加入了下面的减枝优化。即先不考虑身
体，只考虑蛇头，使用宽搜算出仅蛇头一段从初始点走到结束点所需要的最少步数
minstep，然后将除了头以外其它部分的初始位置看作石块，再用宽搜算出仅蛇头从初始
点走到结束点所需要的最少步数 maxstep，这里可以证明所要求的最小步数就在区间
 [minstep,maxstep] 内，如果minstep == maxstep，则该值就是所求值，如果 minstep
 不存在（不可达），则所求蛇的路径也不存在，（这里有一点要说明的是，即使 maxstep
  不存在，也不能说明所求最短路径不存在，因为求 maxstep 时将蛇尾看作固定的障碍物，
  而实际情况下，蛇尾会随蛇头运动，会使求 maxstep 过程中不可行的路径变为可行。）
  否则，我们就执行 1 里定义的宽搜过程，在宽搜过程中，当已走的步数加上从当前点
  到终点的最短距离（x 坐标 + y 坐标）大于 maxstep 时，就不再考虑此路径，如此达
  到减枝提速的目的。
3、maxstep 是上限可以这么理解：当将身体看作固定的障碍时，从头到终点要绕过尾巴，
而实际情况下，尾巴随头移动，故要绕开尾巴所要走的步数 <= 固定尾巴时的情况，又只
有在初始情况下，尾巴可能会出现在头所要走的路径的前方，需要绕开，在第一次绕开之
后，尾巴就一直跟在头后面了，不会再成为障碍物，与没尾巴时的情况一样，因此 maxstep
即是步数的上限。

小技巧，如果多CASE输入的话，初始化vis操作可能会占用比较多的时间，我们可以将
vis[][]==cas这个条件来判断是否是本CASE中所访问的点。
*/
#include <iostream>
#include <cstdio>
#include <string.h>
#include <queue>
#define N 21
using namespace std;
bool stone[N][N];
bool flag[N][N][1<<14];
int way[4][2]={0,1,0,-1,1,0,-1,0};
struct node
{
    int x,y,step;
    int body;
    node( int a=0,int b=0,int s=0,int bo=0 )
    {
        x=a , y=b , step=s , body=bo ;
    }
};

bool legal( node t,node temp,int L )
{
    int i,x,y,tt;
    int re[10];
    x=temp.x , y=temp.y ;
    for( i=L-1;i>0;i-- )
    {
        re[i]=temp.body&3;
        temp.body>>=2;
    }
    for( i=1;i<L;i++ )
    {
        if( t.x==x && t.y==y )  return false;
        if( re[i]==0 )  y++;
        else if( re[i]==1 ) y--;
        else if( re[i]==2 ) x++;
        else if( re[i]==3 ) x--;
    }
    if( t.x==x && t.y==y )  return false;
    return true;
}

int main()
{
    int n,m,L,i,j,k;
    int a,b,aa,bb,count;
    bool mark;
    node head,temp,t;
    queue <node>q;
    count=1;
    while( scanf( "%d%d%d",&n,&m,&L ) && (n||m||L) )
    {
        memset( stone,0,sizeof(stone) );
        memset( flag,0,sizeof(flag) );
        head.body=0;
        for( i=0;i<L;i++ )
        {
            scanf( "%d%d",&a,&b );
            if( i==0 )  head.x=a , head.y=b , head.step=0 ;
            else
            {
                if( a==aa && b==bb+1 )   head.body=(head.body<<2)+0;
                else if( a==aa && b==bb-1 ) head.body=(head.body<<2)+1;
                else if( a==aa+1 && b==bb ) head.body=(head.body<<2)+2;
                else if( a==aa-1 && b==bb ) head.body=(head.body<<2)+3;
            }
            aa=a , bb=b;
        }
        scanf( "%d",&k );
        for( i=0;i<k;i++ )
        {
            scanf( "%d%d",&a,&b );
            stone[a][b]=1;
        }
        q.push( head ); flag[head.x][head.y][head.body]=1;
        mark=false;
        while( !q.empty() )
        {
            temp=q.front(); q.pop();
            if( temp.x==1 && temp.y==1 ) { mark=true; break; }
            for( i=0;i<4;i++ )
            {
                t.x=temp.x+way[i][0];
                t.y=temp.y+way[i][1];
                t.step=temp.step+1;
                t.body=temp.body;
                if( t.x>0 && t.x<=n && t.y>0 && t.y<=m && !stone[t.x][t.y] &&
                   !flag[t.x][t.y][t.body] && legal(t,temp,L) )
                {
                    flag[t.x][t.y][t.body]=1;
                    t.body>>=2;
                    if( i==0 )  t.body+=( 1<<(2*L-4) );
                    else if( i==1 )  t.body+=( 0<<(2*L-4) );
                    else if( i==2 )  t.body+=( 3<<(2*L-4) );
                    else  t.body+=( 2<<(2*L-4) );
                    q.push(t);
                }
            }
        }
        if( mark )  printf( "Case %d: %d\n",count,temp.step );
        else printf( "Case %d: -1\n",count );
        while( !q.empty() )   q.pop();
        count++;
    }
    return 0;
}