bzoj 1556 墓地秘密 —— 状压DP

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1556

预处理出两个障碍四个方向之间的距离(转弯次数)，就可以状压DP了；

但预处理很麻烦...参考了TJ...：https://blog.csdn.net/senyelicone/article/details/56668048

用 spfa ，记录当前位置带一个朝向，然后转移时判断一下如果朝向不同就+1；

最后再从起点出发同样预处理一下，作为初始状态即可；

注意读入的地图上的 '#' 不仅是机关石，还有墙...所以不能忽略。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
int const maxn=,maxm=,inf=0x3f3f3f3f;
int n,m,T,f[<<][maxm][],dis[maxm][][maxm][],d[maxn][maxn][];
int dx[]={,,,-},dy[]={,,-,},xx[maxm],yy[maxm],sx,sy,ans;
bool vis[maxn][maxn],roc[maxn][maxn];
queue<pair<int,int> >q;
bool ck(int x,int y){return x>&&y>&&x<=n&&y<=m;}
void spfa(int nw,int x,int y,int dr)
{
    if(roc[x][y]||!ck(x,y))return;
    while(q.size())q.pop();
    memset(d,0x3f,sizeof d);
    q.push(make_pair(x,y)); vis[x][y]=;
    for(int k=;k<;k++)d[x][y][k]=;
    while(q.size())
    {
        int nx=q.front().first,ny=q.front().second; q.pop(); vis[nx][ny]=;
        for(int i=;i<;i++)
        {
            int tx=nx+dx[i],ty=ny+dy[i];
            if(roc[tx][ty]||!ck(tx,ty))continue;//
            for(int j=;j<;j++)
                if(d[tx][ty][j]>d[nx][ny][i]+(i!=j))
                {
                    d[tx][ty][j]=d[nx][ny][i]+(i!=j);
                    if(!vis[tx][ty])vis[tx][ty]=,q.push(make_pair(tx,ty));
                }
        }
    }
    for(int i=;i<=T;i++)
        for(int j=;j<;j++)//从j撞击xi,yi
        {
            int tx=xx[i]+dx[j],ty=yy[i]+dy[j],tmp=inf;
            for(int k=;k<;k++)tmp=min(tmp,d[tx][ty][k]+(tx+dx[k]!=xx[i]||ty+dy[k]!=yy[i]));//反向
            dis[nw][dr][i][j]=tmp;
        }
}
int main()
{
    scanf("%d%d%d",&n,&m,&T);
    char ch[maxn];
    for(int i=;i<=n;i++)
    {
        cin>>ch;
        for(int j=;j<m;j++)
            if(ch[j]=='#')roc[i][j+]=;
    }
    for(int i=;i<=T;i++)
    {
        scanf("%d%d",&xx[i],&yy[i]);
//        roc[x][y]=1;//'#'表示墙，不一定是机关石!!!
    }
    for(int i=;i<=T;i++)
        for(int j=;j<;j++)
            spfa(i,xx[i]+dx[j],yy[i]+dy[j],j);
    scanf("%d%d",&sx,&sy);
    spfa(T+,sx,sy,);
    memset(f,0x3f,sizeof f); f[][T+][]=;
    int mx=(<<T); ans=inf;
    for(int s=;s<mx;s++)
      for(int i=;i<=T+;i++)
        for(int j=;j<=;j++)  if(f[s][i][j]!=inf)
          for(int k=;k<=T;k++)
            for(int l=;l<;l++)
              f[s|(<<(k-))][k][l]=min(f[s|(<<(k-))][k][l],f[s][i][j]+dis[i][j][k][l]+);
    for(int i=;i<=T;i++)
        for(int j=;j<;j++)
            ans=min(ans,f[mx-][i][j]);
    printf("%d\n",ans);
    return ;
}