HDU 5067 (状态压缩DP+TSP)

时间：2023-03-08 22:07:15

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

题目大意：蓝翔挖掘机挖石子。把地图上所有石子都运回起点，问最少耗时。

解题思路：

首先得YY出来。

最少耗时肯定是从起点出发，把所有石子点走一遍且只走一遍，把石子装在车上，然后最后回到起点。

由于石子堆最多也就10个。不难看出这就是个裸的TSP。

首先BFS计算出每个石子间的最短路。当时打BC的时候SB了，这题图没有障碍，直接根据坐标算就行了。

然后进行TSP就行了。

起点的石子无所谓。所以TSP的起点就是（1,1)。

#include "cstdio"

#include "iostream"

#include "cstring"

#include "vector"

#include "queue"

using namespace std;

#define inf 0x3f3f3f3f

int area[][],dp[<<][],vis[][],dis[][],dir[][]={-,,,,,-,,},n,m,g;

struct status

{

    int x,y,dep;

    status(int x,int y,int dep):x(x),y(y),dep(dep) {}

};

int bfs(status s,status e)

{

    memset(vis,,sizeof(vis));

    vis[s.x][s.y]=true;

    queue<status> Q;

    Q.push(status(s.x,s.y,));

    while(!Q.empty())

    {

        status t=Q.front();Q.pop();

        for(int s=;s<;s++)

        {

            int X=t.x+dir[s][],Y=t.y+dir[s][];

            if(X<||X>n||Y<||Y>m||vis[X][Y]) continue;

            vis[X][Y]=true;

            if(X==e.x&&Y==e.y) return t.dep+;

            Q.push(status(X,Y,t.dep+));

        }

    }

}

int main()

{

    //freopen("in.txt","r",stdin);

    while(scanf("%d%d",&n,&m)!=EOF)

    {

        vector<status> stone;

        stone.push_back(status(,,));

        for(int i=;i<=n;i++)

          for(int j=;j<=m;j++)

            {

               scanf("%d",&g);

               if(i==&&j==) continue;

               if(g) stone.push_back(status(i,j,));

            }

        for(int i=;i<stone.size();i++)

        {

            for(int j=;j<stone.size();j++)

            {

                if(i==j) dis[i][j]=;

                else {if(!dis[i][j]) dis[i][j]=dis[j][i]=bfs(stone[i],stone[j]);}

            }

        }

        int s=stone.size();

        int cnt=<<s;

        for(int i=;i<cnt;i++)

        {

            for(int j=;j<s;j++)

            {

                if(!i&(<<j)) continue;

                if(i==(<<j)) dp[i][j]=dis[][j];

                else

                {

                    dp[i][j]=inf;

                    for(int k=;k<s;k++) if(i&(<<k)&&k!=j) dp[i][j]=min(dp[i][j],dp[i^(<<j)][k]+dis[k][j]);

                }

            }

        }

        int ans=inf;

        for(int i=;i<s;i++)

            ans=min(ans,dp[cnt-][i]+dis[i][]);

        printf("%d\n",ans);

        memset(dp,,sizeof(dp));

        memset(dis,,sizeof(dis));

    }

}
11906643 2014-10-18 22:12:42 Accepted 5067 31MS 2200K 2159 B C++ Physcal

相关文章