HDU-5067-Harry And Dig Machine（状压DP）

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

DP[i][j]表示在已经经过状态为i的格子的情况下，当前在j点的最短路径；

显然有 DP[i][j]=min( DP[i][j],   dp[i|k][p] + dis[p][j] )

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <iostream>
#include <queue>
#include <cmath>
#include <string>
#include <map>
using namespace std;
#define INF 0x3f3f3f3f
#define MAXN 105
#define MAXM 100005

const long long mod=1e9+7;
int x[MAXN],y[MAXN],cnt;
int dis[MAXN][MAXN];
int dp[1<<15][MAXN];
int a[53][53];
int main()
{
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        int d;
        x[0]=y[0]=0;
        x[1]=y[1]=0;
        cnt=2;
        for(int i=0; i<n; ++i)
            for(int j=0; j<m; ++j)
            {
                scanf("%d",&a[i][j]);
                if(a[i][j]>0)
                {
                    x[cnt]=i;
                    y[cnt]=j;
                    for(int k=0; k<cnt; ++k)
                        dis[cnt][k]=dis[k][cnt]=abs(i-x[k])+abs(j-y[k]);
                    cnt++;
                }
            }
        memset(dp,INF,sizeof(dp));
        dp[1][0]=0;
        for(int i=1; i<(1<<cnt); ++i)
            for(int p=0; p<cnt; ++p)
                if((1<<p)&i)
                    if(dp[i][p]!=INF)
                        for(int j=0; j<cnt; ++j)
                            if(!((1<<j)&i))
                                dp[i|(1<<j)][j]=min( dp[i|(1<<j)][j], dp[i][p]+dis[p][j] );
        cout<<dp[ (1<<cnt)-1 ][1]<<endl;
    }
    return 0;
}

/*
5
2 5

*/