hdu 4568(SPFA预处理+TSP)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4568

思路：先用spfa预处理出宝藏与宝藏之间的最短距离，宝藏到边界的最短距离，然后就是经典的求TSP过程了。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 using namespace std;
 #define MAXN 222
 #define inf 1<<30

 struct Point{
     int x,y;
 }point[MAXN];

 int value[MAXN][MAXN];//原图
 int map[MAXN][MAXN];//宝藏间的最短距离
 int dist[MAXN][MAXN];//宝藏到边界的距离
 int dd[MAXN];//宝藏到达边界的最短距离
 bool mark[MAXN][MAXN];
 int dp[<<][];
 int n,m,k;
 int dir[][]={{-,},{,},{,-},{,}};

 void spfa(int num)
 {
     memset(mark,false,sizeof(mark));
     for(int i=;i<n;i++)
         for(int j=;j<m;j++)dist[i][j]=inf;
     queue<pair<int,int> >que;
     que.push(make_pair(point[num].x,point[num].y));
     if(dist[point[num].x][point[num].y]==-)return;
     dist[point[num].x][point[num].y]=;
     while(!que.empty()){
         int x=que.front().first,y=que.front().second;
         que.pop();
         mark[x][y]=false;
         if(x==||x==(n-)||y==||y==(m-)){
             dd[num]=min(dd[num],dist[x][y]);
         }
         for(int i=;i<;i++){
             int xx=x+dir[i][],yy=y+dir[i][];
             if(xx>=&&xx<n&&yy>=&&yy<m&&value[xx][yy]!=-){
                 if(dist[x][y]+value[xx][yy]<dist[xx][yy]){
                     dist[xx][yy]=dist[x][y]+value[xx][yy];
                     if(!mark[xx][yy]){
                         mark[xx][yy]=true;
                         que.push(make_pair(xx,yy));
                     }
                 }
             }
         }
     }
 }

 int main()
 {
     int _case;
     scanf("%d",&_case);
     while(_case--){
         scanf("%d%d",&n,&m);
         for(int i=;i<n;i++)
             for(int j=;j<m;j++)
                 scanf("%d",&value[i][j]);
         scanf("%d",&k);
         for(int i=;i<k;i++)scanf("%d%d",&point[i].x,&point[i].y);
         for(int i=;i<k;i++)
             for(int j=;j<k;j++)
                 map[i][j]=(i==j)?:inf;
         for(int i=;i<(<<k);i++)
             for(int j=;j<k;j++)dp[i][j]=inf;
         fill(dd,dd+MAXN,inf);
         for(int i=;i<k;i++){
             spfa(i);
             for(int j=;j<k;j++){
                 if(i==j)continue;
                 map[i][j]=min(map[i][j],dist[point[j].x][point[j].y]);//宝藏与宝藏之间的最近距离
             }
             dp[<<i][i]=dd[i]+value[point[i].x][point[i].y];//宝藏到边的最近距离
         }
         for(int s=;s<(<<k);s++){
             for(int i=;i<k;i++){
                 if(s&(<<i)==)continue;
                 if(dp[s][i]==inf)continue;
                 for(int j=;j<k;j++){
                     if(s&(<<j)==)continue;
                     dp[s|(<<j)][j]=min(dp[s|(<<j)][j],dp[s][i]+map[i][j]);
                 }
             }
         }
         int ans=inf;
         for(int i=;i<k;i++){
             ans=min(ans,dp[(<<k)-][i]+dd[i]);//拿到了所有的宝藏之后还要走出来
         }
         printf("%d\n",ans);
     }
     return ;
 }