【TSP】 HDU 4568 Hunter

点击打开链接

题意：

必须经过 从图外进入图中，必须经过 K个格子，再走出格子，求最小的值。

先对全图BFS，计算出最短路和格子 走出图的最短路

对于K个格子用TSP 处理

//#include <bits/stdc++.h>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cctype>
#include <climits>
#include <iostream>
#include <algorithm>
#include <queue>
using namespace std;
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define LL long long
int n,m,mp[333][333],x[15],y[15];
int len[333][333],vis[333][333];
int zdl[15][15],dian[15];
int dd[4][2]={1,0,-1,0,0,1,0,-1};
struct node
{
    int x,y,len;
    bool operator < (const node &b) const {
        return len>b.len;
    }
};
priority_queue<node> q;
int gl(int x,int y)
{
    if(x==0||x==n-1||y==0||y==m-1) return 1;
    return 0;
}
void bfs(int x,int y,int iii)
{
    memset(vis,0,sizeof(vis));
    while(!q.empty()) q.pop();
    node front , rear;
    front.x=x,front.y=y,front.len=len[x][y]=mp[x][y];
    q.push(front);
    vis[x][y]=1;
    if(gl(x,y)) dian[iii]=min(dian[iii],front.len-mp[x][y]);
    while(!q.empty())
    {
        front=q.top();
        q.pop();
        for(int i=0;i<4;i++)
        {
            int dx=front.x+dd[i][0];
            int dy=front.y+dd[i][1];
            if(dx<0||dy<0||dx>=n||dy>=m||vis[dx][dy]||mp[dx][dy]==-1)
                continue;
            vis[dx][dy]=1;
            len[dx][dy]=front.len;
            rear.x=dx,rear.y=dy,rear.len=front.len+mp[dx][dy];
            if(gl(dx,dy))
                dian[iii]=min(dian[iii],rear.len-mp[x][y]);
            q.push(rear);
        }
    }
}
/***************/
int dp[14][1<<13][14];
int main()
{
    int t,k;
    cin>>t;
    while(t--)
    {
        memset(zdl,-1,sizeof(zdl));
        cin>>n>>m;
        for(int i=0;i<n;i++)
        for(int j=0;j<m;j++)
        {
            scanf("%d",&mp[i][j]);
            if(mp[i][j]<0) mp[i][j]=-1;
        }
        scanf("%d",&k);
        for(int i=0;i<k;i++)
        {
            scanf("%d%d",&x[i],&y[i]);
        }
        for(int i=0;i<k;i++)
        {
            zdl[i][i]=0;
            memset(len,0,sizeof(len));
            dian[i]=999999;
            bfs(x[i],y[i],i);
            for(int j=0;j<k;j++)
            {
                if(i==j) continue;
                zdl[i][j]=len[x[j]][y[j]];
            }

        }
        /********/
        memset(dp,-1,sizeof(dp));
        for(int s=0;s<k;s++)//star
        {
            dp[s][1<<s][s]=0;
            for(int i=0;i<(1<<k);++i)
            {
                for(int j=0;j<k;++j)
                {
                    if(dp[s][i][j]==-1) continue;
                    for(int kk=0;kk<k;++kk)
                    {
                        if(i&(1<<kk)) continue;
                        if(dp[s][i|(1<<kk)][kk]==-1)
                            dp[s][i|(1<<kk)][kk]=dp[s][i][j]+zdl[j][kk];
                        else dp[s][i|(1<<kk)][kk]=min(dp[s][i][j]+zdl[j][kk],dp[s][i|(1<<kk)][kk]);
                    }
                }
            }
            for(int e=0;e<k;e++)
            {
                if(s==e) continue;
                dp[s][(1<<k)-1][e]+=mp[x[e]][y[e]];
            }
            dp[s][(1<<k)-1][s]=999999;
            for(int e=0;e<k;e++)
            {
                dp[s][(1<<k)-1][s]=min(dp[s][(1<<k)-1][s],dp[s][(1<<k)-1][e]+zdl[s][e]);
            }
        }
        if(k==1) dp[0][1][0]=mp[x[0]][y[0]];
        int ans=9999999;
        for(int s=0;s<k;s++)
        {
            for(int e=0;e<k;e++)
            {
                //cout<<s<<" "<<dian[s]<<" "<<e<<" "<<dian[e]<<" "<<dp[s][(1<<k)-1][e]<<endl;
                if(dp[s][(1<<k)-1][e]!=-1)
                    ans=min(ans,dp[s][(1<<k)-1][e]+dian[s]+dian[e]);
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}
/*
3 6
-1 -1 -1 -1 -1 -1
1 1 1 1 1 -1
-1 -1 -1 -1 -1 -1
5
1 0
1 1
1 2
1 3
1 4
*/