[bzoj2595][WC2008]游览计划

2595: [Wc2008]游览计划

Time Limit: 10 Sec Memory Limit: 256 MBSec Special Judge
Submit: 849 Solved: 344
[Submit][Status][Discuss]
Description


Input

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

Output

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

Sample Input

4 4
0 1 1 0
2 5 5 1
1 5 5 1
0 1 1 0

Sample Output

6
xoox
___o
___o
xoox

HINT

对于100%的数据，N,M,K≤10，其中K为景点的数目。输入的所有整数均在[0,2^16]的范围内

Source

Ljcc930提供SPJ

这个东西好像叫斯坦纳树。
大体的思路是这样的：我们需要做两种dp
第一个是一个状压，我们将每个景点是否到达的压成二进制，然后dp一下。
p为k的子集：f[i][j][k]=f[i][j][p|b[i][j]]+f[i][j][(k-p)|b[x][y]]
我们可以这样去枚举子集：for(p=k;p;p=k&(p-1))
再对于每一个状态k，我们求完状压之后，再写一个类似于spfa的dp，用之前状压求出的数组去向各个方向更新。
设xx=x+xi[i],yy=y+yi[i];
f[xx][yy][k|b[xx][yy]]=f[x][y][k]+a[xx][yy]
其实这个东西的复杂度巨高无比，应该没有什么卵用吧。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define inf 210000000
struct S{int x,y,k;};
S use,pre[15][15][1<<11],q[1000000];
int n,m,tot=0,f[15][15][1<<11],a[15][15],b[15][15],ans[15][15],h,t;
int xi[4]={-1,0,0,1},yi[4]={0,-1,1,0};
bool check[15][15][1<<11];
void spfa(int k)
{
    int xx,yy,i;
    while(h<=t){
        S now=q[h];
        check[now.x][now.y][k]=true;
        for(i=0;i<4;i++){
            xx=now.x+xi[i];yy=now.y+yi[i];
            if(xx<=n&&xx>=1&&yy<=m&&yy>=1&&f[now.x][now.y][k]+a[xx][yy]<f[xx][yy][now.k|b[xx][yy]]){
                f[xx][yy][now.k|b[xx][yy]]=f[now.x][now.y][k]+a[xx][yy];
                use.x=now.x;use.y=now.y;use.k=k;
                pre[xx][yy][now.k|b[xx][yy]]=use;
                if(now.k|b[xx][yy]==k&&check[xx][yy][k]){
                    check[xx][yy][k]=false;
                    t+=1;q[t].x=xx;q[t].y=yy;q[t].k=k;
                }
            }
        }
        h+=1;
    }
}
void find(int x,int y,int k)
{
    if(!x||!y) return ;
    S now=pre[x][y][k];
    if(!a[x][y]) ans[x][y]=3;
    else ans[x][y]=2;
    find(now.x,now.y,now.k);
    if(now.x==x&&now.y==y)
      find(x,y,(k-now.k)|b[x][y]);
}
int main()
{
    int i,j,k,now,p,xx,yy;
    scanf("%d%d",&n,&m);
    memset(f,127/3,sizeof(f));
    memset(check,1,sizeof(check));  
    for(i=1;i<=n;i++)
      for(j=1;j<=m;j++){
        scanf("%d",&a[i][j]);
        if(!a[i][j]){
          tot+=1;
          f[i][j][1<<(tot-1)]=0;
          b[i][j]=1<<(tot-1);
        }
        f[i][j][0]=a[i][j];
        ans[i][j]=1;
      }
    for(k=1;k<=(1<<tot)-1;++k){
        h=1;t=0;
        for(i=1;i<=n;++i)
          for(j=1;j<=m;++j){
            if(!a[i][j]&&!(k&b[i][j])) continue;
            for(p=k;p;p=k&(p-1))
              if(f[i][j][k]>f[i][j][p|b[i][j]]+f[i][j][(k-p)|b[i][j]]-a[i][j]){
                f[i][j][k]=f[i][j][p|b[i][j]]+f[i][j][(k-p)|b[i][j]]-a[i][j];
                use.x=i;use.y=j;use.k=p|b[i][j];
                pre[i][j][k]=use;
              }
            if(f[i][j][k]<inf){
                check[i][j][k]=false;
                t+=1;q[t].x=i;q[t].y=j;q[t].k=k;
            }
          }
        spfa(k);
    }
    int minn=inf;
    for(i=1;i<=n;++i)
      for(j=1;j<=m;++j)
        if(!a[i][j]&&minn>f[i][j][(1<<tot)-1]){
            minn=f[i][j][(1<<tot)-1];
            xx=i;yy=j;
        }
    printf("%d\n",minn);
    find(xx,yy,(1<<tot)-1);
    for(i=1;i<=n;++i){
        for(j=1;j<=m;++j){
            if(ans[i][j]==1) printf("_");
            if(ans[i][j]==2) printf("o");
            if(ans[i][j]==3) printf("x");
        }
        printf("\n");
    }
}