Bzoj2595: [Wc2008]游览计划

Time Limit: 10 Sec  Memory Limit: 256 MBSec  Special Judge
Submit: 1463  Solved: 679
[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

 

斯坦纳树。

记录状态简直恶心。

 

炸了好久好久，最后发现是第81行初始化的时候x=INF写成了w=INF

↑白弄了好久的标准代码比对

/*by SilverN*/
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
const int INF=0x3f3f3f3f;
const int mx[5]={0,1,0,-1,0};
const int my[5]={0,0,1,0,-1};
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0' || ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int mp[15][15];
int n,m,k;
int id[15];
int f[15][15][1100];
struct pt{
    int x,y,w;
}p[15][15][1100];
queue<pair<int,int> >q;
bool vis[30][30];
void SPFA(int s){
    while(!q.empty()){
        int x=q.front().first;
        int y=q.front().second;
        q.pop();
        vis[x][y]=0;
        for(int k=1;k<=4;k++){
            int nx=x+mx[k],ny=y+my[k];
            if(nx<1 || nx>n || ny<1 || ny>m)continue;
            if(f[nx][ny][s]>f[x][y][s]+mp[nx][ny]){
                f[nx][ny][s]=f[x][y][s]+mp[nx][ny];
                p[nx][ny][s]=(pt){x,y,s};
                if(!vis[nx][ny]){
                    q.push(make_pair(nx,ny));
                    vis[nx][ny]=1;
                }
            }
        }
    }
    return;
}
void DFS(int x,int y,int s){//标记路径 
    if(x>=INF || !p[x][y][s].w)return;
    vis[x][y]=1;
    DFS(p[x][y][s].x,p[x][y][s].y,p[x][y][s].w);
    if(p[x][y][s].x==x && p[x][y][s].y==y)DFS(p[x][y][s].x,p[x][y][s].y,s-p[x][y][s].w);
    return;
}
void Print(){
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if(!mp[i][j])printf("x");
            else if(!vis[i][j])printf("_");
            else printf("o");
        }
        printf("\n");
    }
    return;
}
int main(){
    int i,j;
    n=read();m=read();k=0;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++){
            mp[i][j]=read();
            if(!mp[i][j])k++;
        }
    //
    memset(f,0x3f,sizeof f);
    for(i=1;i<=k;i++)id[i]=1<<(i-1);
    int ed=(1<<k)-1;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
            for(int c=0;c<=ed;c++)
                p[i][j][c].x=INF;
    //init
    k=0;
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
            if(!mp[i][j]){
                f[i][j][id[++k]]=0;
            }
    for(int s=1;s<=ed;s++){
        for(i=1;i<=n;i++)
         for(j=1;j<=m;j++){
             for(int c=(s-1)&s;c;c=(c-1)&s){
                 if(f[i][j][s]>f[i][j][c]+f[i][j][s-c]-mp[i][j]){
                     f[i][j][s]=f[i][j][c]+f[i][j][s-c]-mp[i][j];
                     p[i][j][s]=(pt){i,j,c};
                 }
             }
             if(f[i][j][s]<INF){
                 q.push(make_pair(i,j));
                 vis[i][j]=1;
            }
         }
        SPFA(s);
    }
    memset(vis,0,sizeof vis);
    while(!q.empty())q.pop();
    for(i=n;i;i--)
        for(j=m;j;j--){
            if(!mp[i][j]){
                DFS(i,j,ed);
                printf("%d\n",f[i][j][ed]);
                Print();
                return 0;
            }
        }
    return 0;
}