WC2008游览计划（BZOJ2595）

传送门

 

2595: [Wc2008]游览计划

Time Limit: 10 Sec  Memory Limit: 256 MBSec  Special Judge
[Submit][Status]

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

出题者XXOO的真丧病-_-#

写了好长时间了。。比较简单的插头dp。。保存路径也比较简单。。WA了N次是因为。。我把X、O写倒了。。

轮廓线以上的格子表示有没有选中，那么对于一个非景点格子，如果可以不选，那么要么是i-1,j不选 ，要么是i-1,j选了但是轮廓线上存在其他的格子和i-1,j是一个联通块。

Codes:

#include<set>
#include<queue>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N = 1<<15;
const int Hash = 10007;
#define nxt (cur^1)
#define For(i,n) for(int i=1;i<=n;i++)
#define Rep(i,l,r) for(int i=l;i<=r;i++)
#define Down(i,r,l) for(int i=r;i>=l;i--)
 
pair< pair<int,int> ,int > opt[N*110];
int pre[N*110],optot;
int loc,cur,n,m,maze[15][15],ans=1<<30,code[15];
char Ans[15][15];
 
struct statedp{
    int Loc[N],size,f[N],head[Hash],next[N],st[N];
    void clear(){size=0;memset(head,-1,sizeof(head));}
    void push(int state,int ans,int last,int x,int y,int kind){
        int Key = state % Hash;
        for(int p=head[Key];p!=-1;p=next[p])
            if(st[p]==state) {  
                if(ans<f[p]) {
                    f[p]=ans;pre[Loc[p]]=last;
                    opt[Loc[p]]=make_pair(make_pair(x,y),kind);
                }
                return;
            }
        f[size]=ans;st[size]=state;
        Loc[size]=++optot;
        pre[optot]=last;opt[optot]=make_pair(make_pair(x,y),kind);
        next[size]=head[Key];head[Key]=size++;
    }
}dp[2];
 
void init(){
    scanf("%d%d",&n,&m);
    For(i,n)
      For(j,m)
        scanf("%d",&maze[i][j]);
}
 
void decode(int st){
    For(i,m) code[i]=st&7,st>>=3;
}
 
int encode(){
    int st=0,cnt=1,h[10];
    memset(h,-1,sizeof(h));
    h[0]=0;
    Down(i,m,1){
        if(h[code[i]]==-1) h[code[i]]=cnt++;
        code[i]=h[code[i]];
        st = (st<<3)|code[i];
    }
    return st;
}
 
void trans(int Left,int Up){
    For(i,m)
      if(code[i]==Up) code[i] = Left;//合并联通快
}
 
void dpblank(int i,int j,int cur,int kind){
    Rep(k,0,dp[cur].size-1){
        decode(dp[cur].st[k]);
        int Left = (j==1)?(0):(code[j-1]) ,Up = (i==1)?(0):(code[j]),cnt=0;
        if(kind){
            int tot = 0;
            For(scan,m)  tot+=(code[scan]==Up);
            if(tot>1 || !code[j]){
                int tmp = code[j];code[j] = 0;
                dp[nxt].push(encode(),dp[cur].f[k],dp[cur].Loc[k],i,j,0);
                code[j]=tmp;
            }
        }
        if(Left&&Up) trans(Left,Up);else
        if(Left||Up) code[j] = Left + Up;
        else         code[j] = (5+(cnt++));
        dp[nxt].push(encode(),dp[cur].f[k]+kind,dp[cur].Loc[k],i,j,1);
    }
}
void DP(){
    int cur=0;dp[cur].clear();dp[cur].push(0,0,0,0,0,0);
    For(i,n)
        For(j,m){
            dp[nxt].clear();
            dpblank(i,j,cur,maze[i][j]);
            cur^=1;
        }
    Rep(k,0,dp[cur].size-1){
        decode(dp[cur].st[k]);
        int i;for(i=1;i<=m;i++) if(code[i]>1) break;
        if(i<m) continue;
        if(dp[cur].f[k]<ans){
            loc = dp[cur].Loc[k];
            ans=dp[cur].f[k];
        }
    }
    printf("%d\n",ans);
}
 
void PrintPath(int i){
    while(i){
        pair<pair<int,int>,int> T = opt[i];
        int x=T.first.first,y=T.first.second,c=T.second;
        if(!maze[x][y]) Ans[x][y]='x';
        else if(c==0)   Ans[x][y]='_';
        else            Ans[x][y]='o';
        i=pre[i];
    }
}
 
int main(){
    init();
    DP();
    PrintPath(loc);
    For(i,n)
      For(j,m)
        if(j==m) printf("%c\n",Ans[i][j]);
        else     printf("%c",Ans[i][j]);
    return 0;
}