【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=,INF=~0u>>;
 const int fx[]={,,-,},
           fy[]={,,,-};
 typedef long long LL;
 //#define debug
 int n,m,a[N][N],d[N][N],cnt=;
 int f[N][N][],pre[N][N][];
 bool vis[N][N];

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

 void spfa(int Set){
     while(!Q.empty()){
         int x=Q.front().x,y=Q.front().y; Q.pop();
         vis[x][y]=;
         rep(i,){
             int tx=x+fx[i],ty=y+fy[i];
             if (tx< || ty< || 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]=;
             }
         }
     }
 }

 void dfs(int i,int j,int Set){
     if (pre[i][j][Set]==INF || pre[i][j][Set]%==)return;
     vis[i][j]=;
     int x=pre[i][j][Set]/,
     y=(pre[i][j][Set]-x*)/,
     z=(pre[i][j][Set]-x*-y*);
     #ifdef debug
     printf("%d,%d,%d----->%d,%d,%d\n",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=;Set<(<<cnt);Set++){
         F(x,,n) F(y,,m){
             for(int s=Set&(Set-);s;s=(s-)&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]=;}
         }
         spfa(Set);
     }
     int x,y;
     F(i,,n) F(j,,m) if (!a[i][j]) {x=i; y=j;break;}
     printf("%d\n",f[x][y][(<<cnt)-]);

     #ifdef debug
     F(i,,n) F(j,,m)
         F(k,,(<<cnt)-)
             if (pre[i][j][k]!=INF)
             printf("pre[%d][%d][%d]=%d\n",i,j,k,pre[i][j][k]);
     #endif
     dfs(x,y,(<<cnt)-);
     F(i,,n) F(j,,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,,) F(j,,)
         F(k,,(<<)) f[i][j][k]=pre[i][j][k]=INF;
     F(i,,n)
         F(j,,m){
             scanf("%d",&a[i][j]);
             if (!a[i][j])
                 f[i][j][<<(cnt++)]=;
         }
     solve();
     return ;
 }