【BZOJ】【2595】【WC2008】游览计划

Orz zky神犇http://blog.csdn.net/iamzky/article/details/42029921

spfa的灵活应用！（好像是求了一个叫做斯坦纳树的东西……）

o(︶︿︶)o 唉我就是太水了，离散化写跪了，x*1e5+y*1e4+k，但是这题里我x和y的范围是[1,10]所以在y==10的时候会出错！！

//BZOJ 2595
#include<queue>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define rep(i,n) for(int i=0;i<n;++i)
#define F(i,j,n) for(int i=j;i<=n;++i)
#define D(i,j,n) for(int i=j;i>=n;--i)
using namespace std;
const int N=13,INF=~0u>>2;
const int fx[]={1,0,-1,0},
          fy[]={0,1,0,-1};
typedef long long LL;
//#define debug
int n,m,a[N][N],d[N][N],cnt=0;
int f[N][N][1100],pre[N][N][1100];
bool vis[N][N];

struct node{ int x,y; };
queue<node>Q;
inline int pack(int x,int y,int s){
    return x*1000000+y*10000+s;
}

void spfa(int Set){
    while(!Q.empty()){
        int x=Q.front().x,y=Q.front().y; Q.pop();
        vis[x][y]=0;
        rep(i,4){
            int tx=x+fx[i],ty=y+fy[i];
            if (tx<1 || ty<1 || tx>n || ty>m) continue;
            
            if (f[tx][ty][Set]>f[x][y][Set]+a[tx][ty]){
                f[tx][ty][Set]=f[x][y][Set]+a[tx][ty];
                pre[tx][ty][Set]=pack(x,y,Set);
                if (!vis[tx][ty]) Q.push((node){tx,ty}),vis[tx][ty]=1;
            }
        }
    }
}

void dfs(int i,int j,int Set){
    if (pre[i][j][Set]==INF || pre[i][j][Set]%5000==0)return;
    vis[i][j]=1;
    int x=pre[i][j][Set]/1000000,
    y=(pre[i][j][Set]-x*1000000)/10000,
    z=(pre[i][j][Set]-x*1000000-y*10000);
    #ifdef debug
    printf("%d,%d,%d----->%d,%d,%d
",x,y,z,i,j,Set);
    #endif
    dfs(x,y,z);
    if (x==i && y==j) dfs(i,j,Set-z);
}

void solve(){
    for(int Set=1;Set<(1<<cnt);Set++){
        F(x,1,n) F(y,1,m){
            for(int s=Set&(Set-1);s;s=(s-1)&Set){
                if(f[x][y][Set]>f[x][y][s]+f[x][y][Set-s]-a[x][y]){
                    f[x][y][Set]=f[x][y][s]+f[x][y][Set-s]-a[x][y];
                    pre[x][y][Set]=pack(x,y,s);
                }
            }
            if (f[x][y][Set]!=INF){Q.push((node){x,y}); vis[x][y]=1;}
        }
        spfa(Set);
    }
    int x,y;
    F(i,1,n) F(j,1,m) if (!a[i][j]) {x=i; y=j;break;}
    printf("%d
",f[x][y][(1<<cnt)-1]);
    
    #ifdef debug
    F(i,1,n) F(j,1,m)
        F(k,0,(1<<cnt)-1) 
            if (pre[i][j][k]!=INF) 
            printf("pre[%d][%d][%d]=%d
",i,j,k,pre[i][j][k]);
    #endif 
    dfs(x,y,(1<<cnt)-1);
    F(i,1,n) F(j,1,m){
        if (!a[i][j]) putchar('x');
        else if(vis[i][j]) putchar('o');
        else putchar('_');
        if (j==m) puts("");
    }
}

int main(){
    #ifndef ONLINE_JUDGE
    freopen("file.in","r",stdin);
//    freopen("file.out","w",stdout);
    #endif
    scanf("%d%d",&n,&m);
    F(i,0,10) F(j,0,10)
        F(k,0,(1<<10)) f[i][j][k]=pre[i][j][k]=INF;
    F(i,1,n)
        F(j,1,m){
            scanf("%d",&a[i][j]);
            if (!a[i][j])
                f[i][j][1<<(cnt++)]=0;
        }
    solve();
    return 0;
}