hdu 4568 Hunter

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

题意：给定一个矩阵图，每走一个格子有一个花费，可以从任意边进入，任意边出去，并且图中有k个点要求全部走到，问最少的花费是多少，如果不能够全部走到，输出0。

思路：先用spfa对图做预处理，求出所有必须走的点和地图外构成的图的任意两点的最短路，然后就是枚举必须走的点的访问顺序问题，由于k很小，可以用状态压缩DP，dp[i][j]表示最后选的点是i,j表示已经选的点。

#include<cstdio>#include<cstring>#include<algorithm>#include<queue>using namespace std;const int inf=2099999999;int s[205][205],ma[205][205];int dis[15][15],c[15];int dir[4][2]={{0,1},{0,-1},{1,0},{-1,0}};int n,m;int f[15];int dp[15][1<<14];int d[40005],vis[40005];struct node{    int di,u;    bool operator <(const node& rhs) const{        return di>rhs.di;    }};void spfa(int st,int num){    int i;    priority_queue<node> q;    for(i=0;i<n*m;i++) d[i]=inf;    d[st]=0;    memset(vis,0,sizeof(vis));    q.push((node){0,st});    while(!q.empty()){        node p=q.top();        q.pop();        int u=p.u;        if(vis[u]) continue;        vis[u]=1;        int x,y,xx,yy,xu;        x=u/m,y=u%m;        if(ma[x][y]>-1) dis[num][ma[x][y]]=d[u];        if(x==0||x==n-1||y==0||y==m-1) c[num]=min(c[num],d[u]);        for(i=0;i<4;i++){            xx=x+dir[i][0];            yy=y+dir[i][1];            if(s[xx][yy]==-1) continue;            if(xx>n-1||xx<0||yy>m-1||yy<0) continue;            xu=xx*m+yy;            if(d[xu]>d[u]+s[xx][yy]){                d[xu]=d[u]+s[xx][yy];                q.push((node){d[xu],xu});            }        }    }}int main(){    int i,j,k,K,x,y,state;    int t;    scanf("%d",&t);    while(t--){        scanf("%d%d",&n,&m);        for(i=0;i<n;i++)         for(j=0;j<m;j++)           scanf("%d",&s[i][j]);        memset(ma,-1,sizeof(ma));        memset(dp,0x7f,sizeof(dp));        memset(dis,0,sizeof(dis));        scanf("%d",&K);        for(i=0;i<K;i++){            int a,b;            scanf("%d%d",&a,&b);            ma[a][b]=i;            f[i]=a*m+b;            c[i]=inf;        }        for(i=0;i<K;i++) spfa(f[i],i);        for(i=0;i<K;i++){            x=f[i]/m,y=f[i]%m;            dp[i][1<<i]=c[i]+s[x][y];        }        for(i=1;i<(1<<K);i++){        for(j=0;j<K;j++){            for(k=0;k<K;k++){                if((i&(1<<j))&&(i&(1<<k))&&j!=k){                    if(dp[j][i^(1<<k)]!=inf) dp[k][i] = min(dp[k][i], dp[j][i^(1<<k)]+dis[j][k]);                }            }        }    }        int ans=inf;        for(i=0;i<K;i++)          ans=min(ans,dp[i][(1<<K)-1]+c[i]);        printf("%d\n",ans);    }    return 0;}