hdu3681 状态压缩最优值

写起来有点麻烦，没什么新思维，注意好细节就行。

ACcode：

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;

const int NS=16,INF=1<<30;

struct node
{
    int x,y,id;
    bool operator < (const node &cmp) const
    {
        return id>cmp.id;
    }
} p[NS],start;

int n,m,top,cnt,ans;
int mp[NS][NS];
int dis[NS][NS];
int vis[NS][NS];
int val[1<<NS][NS];
int dp[1<<NS][NS];
int dx[4]= {1,-1,0,0};
int dy[4]= {0,0,1,-1};
char ch[NS];

int judge(int x,int y)
{
    if(x>=0&&x<n&&y>=0&&y<m&&mp[x][y]>=0&&!vis[x][y])
    return 1; return 0;
}

void bfs(int k)
{
    for (int i=0; i<=top; i++)
    dis[k][i]=INF;
    memset(vis,0,sizeof(vis));
    queue<node> q;

    node t,cur=p[k],next;
    cur.id=0,vis[cur.x][cur.y]=1;
    q.push(cur);
    while (!q.empty())
    {
        cur=q.front(),q.pop();
        for (int i=0; i<=top; i++)
        {
            if (cur.x==p[i].x&&cur.y==p[i].y)
                dis[k][i]=cur.id;
        }

        for (int i=0; i<4; i++)
        {
            next.x=cur.x+dx[i];
            next.y=cur.y+dy[i];
            if (judge(next.x,next.y))
            {
                vis[next.x][next.y]=1;
                next.id=cur.id+1;
                q.push(next);
            }
        }
    }
}

int Max(int a1,int b1)
{
    return a1>b1?a1:b1;
}

int Min(int a1,int b1)
{
    return a1<b1?a1:b1;
}

int main()
{
    while (~scanf("%d%d",&n,&m)&&n+m)
    {
        top=cnt=0,ans=INF;
        for (int i=0; i<n; i++)
        {
            scanf("%s",ch);
            for (int j=0; j<m; j++)
            {
                if (ch[j]=='S') mp[i][j]=1;
                if (ch[j]=='D') mp[i][j]=-1;
                if (ch[j]=='F') mp[i][j]=0,start.x=i,start.y=j,start.id=0;
                if (ch[j]=='G') mp[i][j]=2,p[top].x=i,p[top].y=j,p[top++].id=2;
                if (ch[j]=='Y') mp[i][j]=3,p[top].x=i,p[top].y=j,p[top++].id=3,cnt++;
            }
        }
        sort(p,p+top),p[top]=start;
        for (int i=0; i<=top; i++) bfs(i);

        int t,lim=1<<top,st=(1<<cnt)-1;
        for (int i=0;i<lim;i++)
        for (int j=0;j<top;j++)
        dp[i][j]=val[i][j]=INF;
        for (int i=0;i<top;i++)
        dp[1<<i][i]=dis[top][i],
        val[1<<i][i]=i<cnt?dis[top][i]:0;
        if (!cnt) ans=0;
        for (int i=1;i<lim;i++)
        {
            for (int j=0;j<top;j++)
            {
                if (i&(1<<j))
                {
                    t=i^(1<<j);
                    for (int k=0;k<top;k++)
                    if (t&(1<<k))
                    {
                        if (j<cnt)
                        {
                            if (dp[i][j]>Max(dp[t][k],dis[k][j]+val[t][k]))
                            dp[i][j]=Max(dp[t][k],dis[k][j]+val[t][k]),
                            val[i][j]=dis[k][j]+val[t][k];
                            else if (dp[i][j]==Max(dp[t][k],dis[k][j]+val[t][k]))
                            val[i][j]=Min(dis[k][j]+val[t][k],val[i][j]);
                        }
                        else
                        {
                            if (dp[i][j]>Max(dp[t][k],dis[k][j]+val[t][k]))
                            dp[i][j]=Max(dp[t][k],dis[k][j]+val[t][k]);
                            val[i][j]=0;
                        }
                    }
                }
                if ((st&i)==st) ans=Min(ans,dp[i][j]);
            }
        }
        if (ans>=INF) ans=-1;
        printf("%d\n",ans);
    }
    return 0;
}