【WC2008】bzoj2595 游览计划

Description

Input

第一行有两个整数，N和 M，描述方块的数目。 接下来 N行， 每行有 M 个非负整数， 如果该整数为 0， 则该方块为一个景点；
否则表示控制该方块至少需要的志愿者数目。 相邻的整数用 （若干个） 空格隔开， 行首行末也可能有多余的空格。

Output

由 N + 1行组成。第一行为一个整数，表示你所给出的方案 中安排的志愿者总数目。 接下来 N行，每行M
个字符，描述方案中相应方块的情况： z ‘_’（下划线）表示该方块没有安排志愿者； z
‘o’（小写英文字母o）表示该方块安排了志愿者； z ‘x’（小写英文字母x）表示该方块是一个景点；
注：请注意输出格式要求，如果缺少某一行或者某一行的字符数目和要求不 一致（任何一行中，多余的空格都不允许出现） ，都可能导致该测试点不得分。

斯坦纳树。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define R() for (int i=1;i<=n;i++) for (int j=1;j<=m;j++)
const int oo=0x3f3f3f3f;
int a[15][15],dp[10000][15][15],num[15][15],vis[15][15],
pre1[10000][15][15],pre2[10000][15][15],use[15][15],
xx[]={1,-1,0,0},yy[]={0,0,1,-1},
n,m,s,tot;
bool ok(int x,int y)
{
    return x>=1&&x<=n&&y>=1&&y<=m;
}
struct point
{
    int x,y,w;
    bool operator < (const point &p) const
    {
        return w>p.w;
    }
}p1;
void init()
{
    scanf("%d%d",&n,&m);
    R()
      {
        scanf("%d",&a[i][j]);
        if (!a[i][j]) num[i][j]=tot++;
      }
}
void solve()
{
    int s1,x,y,x1,y1;
    memset(dp,0x3f,sizeof(dp));
    R()
    {
        dp[0][i][j]=a[i][j];
        if (!a[i][j]) dp[1<<num[i][j]][i][j]=0;
    }
    for (s=1;s<(1<<tot);s++)
    {
        R() for (s1=s;s1;s1=(s1-1)&s)
          if ((x=dp[s1][i][j]+dp[s^s1][i][j]-a[i][j])<dp[s][i][j])
          {
            dp[s][i][j]=x;
            pre1[s][i][j]=0;
            pre2[s][i][j]=s1;
          }
        priority_queue<point> q;
        R() if (dp[s][i][j]<oo) q.push((point){i,j,dp[s][i][j]});
        memset(vis,0,sizeof(vis));
        while (!q.empty())
        {
            p1=q.top();
            q.pop();
            x=p1.x;
            y=p1.y;
            if (vis[x][y]) continue;
            vis[x][y]=1;
            for (int k=0;k<4;k++)
              if (ok(x1=x+xx[k],y1=y+yy[k])&&dp[s][x][y]+a[x1][y1]<dp[s][x1][y1])
              {
                dp[s][x1][y1]=dp[s][x][y]+a[x1][y1];
                pre1[s][x1][y1]=x;
                pre2[s][x1][y1]=y;
                q.push((point){x1,y1,dp[s][x1][y1]});
              }
        }
    }
}
void mark(int s,int x,int y)
{
    int cnt=0;
    use[x][y]=1;
    if (pre1[s][x][y]) mark(s,pre1[s][x][y],pre2[s][x][y]);
    else if (pre2[s][x][y])
    {
        mark(pre2[s][x][y],x,y);
        mark(s^pre2[s][x][y],x,y);
    }
}
void out()
{
    int ans=oo,x,y;
    R() if (dp[(1<<tot)-1][i][j]<ans)
    {
        ans=dp[(1<<tot)-1][i][j];
        x=i;
        y=j;
    }
    printf("%d\n",ans);
    mark((1<<tot)-1,x,y);
    for (int i=1;i<=n;i++)
    {
        for (int j=1;j<=m;j++)
          if (!a[i][j]) putchar('x');
          else
          {
            if (use[i][j]) putchar('o');
            else putchar('_');
          }
        putchar('\n');
    }
}
int main()
{
    init();
    solve();
    out();
}