HDU 4856 Tunnels(BFS+状压DP)

传送门：http://acm.hdu.edu.cn/showproblem.php?pid=4856
题意：先输入n，m，然后是一个n*n的图，然后输入m条隧道，每条隧道包括起点的x，y，终点的x，y。问你如何用最短的时间走遍所有隧道，每条隧道只能走一次，不能走遍则输出-1。
坑点：会有好几个隧道的起点是相同的，而且路过一个隧道时，可以选择不传送。
思路：最短时间走遍隧道，而且每个只能走一遍，这个很明显是一个状压dp，就像tsp那种。但是我们要处理出每条隧道之间的距离，所以我们就要对每个隧道的终点开始BFS，求出每个隧道的终点->其他隧道的起点的距离，处理完以后就直接状压dp一下就可以了。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <iostream>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
#define pb push_back
#define mp make_pair
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define calm (l+r)>>1

int n,m;
char pic[20][20];
pii st[20],ed[20];
int G[20][20];
int dp[1<<15][15];
#define F first
#define S second
void findID(int s,int x,int y,int step){
    for(int i=0;i<m;i++){
        if(st[i].F==x&&st[i].S==y)G[s][i]=step;
    }
}
bool vis[20][20];
int go[4][2]={{1,0},{-1,0},{0,1},{0,-1}};
inline bool judge(int x,int y){
    return x>=1&&x<=n&&y>=1&&y<=n;
}
void BFS(int s){
    queue< pair<pii,int> >Q;
    memset(vis,false,sizeof vis);
    Q.push(mp(ed[s],0));
    vis[ed[s].F][ed[s].S]=true;
    while(!Q.empty()){
        int x=Q.front().F.F;
        int y=Q.front().F.S;
        int step=Q.front().S;
        Q.pop();
        findID(s,x,y,step);
        for(int i=0;i<4;i++){
            int xx=x+go[i][0];
            int yy=y+go[i][1];
            if(judge(xx,yy)&&!vis[xx][yy]&&pic[xx][yy]!='#'){
                vis[xx][yy]=true;
                Q.push(mp(mp(xx,yy),step+1));
            }
        }
    }
}
int main(){
    //freopen("D://input.txt","r",stdin);
    while(scanf("%d%d",&n,&m)!=EOF){
        for(int i=1;i<=n;i++){
            scanf("%s",pic[i]+1);
        }
        for(int i=0;i<m;i++){
            scanf("%d%d%d%d",&st[i].F,&st[i].S,&ed[i].F,&ed[i].S);
        }
        memset(G,127,sizeof G);
        for(int i=0;i<m;i++){
            BFS(i);
        }
        memset(dp,127,sizeof dp);
        int INF=dp[0][0];
        for(int i=0;i<m;i++){
            dp[1<<i][i]=0;
        }
        for(int k=0;k<(1<<m);k++){
            for(int i=0;i<m;i++)if(dp[k][i]!=INF){
                for(int j=0;j<m;j++)if(!(k&(1<<j))&&(G[i][j]!=INF)){
                    dp[k|(1<<j)][j]=min(dp[k|(1<<j)][j],dp[k][i]+G[i][j]);
                }
            }
        }
        int ans=INF;
        for(int i=0;i<m;i++){
            ans=min(ans,dp[(1<<m)-1][i]);
        }
        if(ans==INF)ans=-1;
        printf("%d\n",ans);
    }
    return 0;
}