【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的时候会出错！！

  1 //BZOJ 2595
  2 #include<queue>
  3 #include<cmath>
  4 #include<cstdio>
  5 #include<cstring>
  6 #include<cstdlib>
  7 #include<iostream>
  8 #include<algorithm>
  9 #define rep(i,n) for(int i=0;i<n;++i)
 10 #define F(i,j,n) for(int i=j;i<=n;++i)
 11 #define D(i,j,n) for(int i=j;i>=n;--i)
 12 using namespace std;
 13 const int N=13,INF=~0u>>2;
 14 const int fx[]={1,0,-1,0},
 15           fy[]={0,1,0,-1};
 16 typedef long long LL;
 17 //#define debug
 18 int n,m,a[N][N],d[N][N],cnt=0;
 19 int f[N][N][1100],pre[N][N][1100];
 20 bool vis[N][N];
 21
 22 struct node{ int x,y; };
 23 queue<node>Q;
 24 inline int pack(int x,int y,int s){
 25     return x*1000000+y*10000+s;
 26 }
 27
 28 void spfa(int Set){
 29     while(!Q.empty()){
 30         int x=Q.front().x,y=Q.front().y; Q.pop();
 31         vis[x][y]=0;
 32         rep(i,4){
 33             int tx=x+fx[i],ty=y+fy[i];
 34             if (tx<1 || ty<1 || tx>n || ty>m) continue;
 35
 36             if (f[tx][ty][Set]>f[x][y][Set]+a[tx][ty]){
 37                 f[tx][ty][Set]=f[x][y][Set]+a[tx][ty];
 38                 pre[tx][ty][Set]=pack(x,y,Set);
 39                 if (!vis[tx][ty]) Q.push((node){tx,ty}),vis[tx][ty]=1;
 40             }
 41         }
 42     }
 43 }
 44
 45 void dfs(int i,int j,int Set){
 46     if (pre[i][j][Set]==INF || pre[i][j][Set]%5000==0)return;
 47     vis[i][j]=1;
 48     int x=pre[i][j][Set]/1000000,
 49     y=(pre[i][j][Set]-x*1000000)/10000,
 50     z=(pre[i][j][Set]-x*1000000-y*10000);
 51     #ifdef debug
 52     printf("%d,%d,%d----->%d,%d,%d\n",x,y,z,i,j,Set);
 53     #endif
 54     dfs(x,y,z);
 55     if (x==i && y==j) dfs(i,j,Set-z);
 56 }
 57
 58 void solve(){
 59     for(int Set=1;Set<(1<<cnt);Set++){
 60         F(x,1,n) F(y,1,m){
 61             for(int s=Set&(Set-1);s;s=(s-1)&Set){
 62                 if(f[x][y][Set]>f[x][y][s]+f[x][y][Set-s]-a[x][y]){
 63                     f[x][y][Set]=f[x][y][s]+f[x][y][Set-s]-a[x][y];
 64                     pre[x][y][Set]=pack(x,y,s);
 65                 }
 66             }
 67             if (f[x][y][Set]!=INF){Q.push((node){x,y}); vis[x][y]=1;}
 68         }
 69         spfa(Set);
 70     }
 71     int x,y;
 72     F(i,1,n) F(j,1,m) if (!a[i][j]) {x=i; y=j;break;}
 73     printf("%d\n",f[x][y][(1<<cnt)-1]);
 74
 75     #ifdef debug
 76     F(i,1,n) F(j,1,m)
 77         F(k,0,(1<<cnt)-1)
 78             if (pre[i][j][k]!=INF)
 79             printf("pre[%d][%d][%d]=%d\n",i,j,k,pre[i][j][k]);
 80     #endif
 81     dfs(x,y,(1<<cnt)-1);
 82     F(i,1,n) F(j,1,m){
 83         if (!a[i][j]) putchar(‘x‘);
 84         else if(vis[i][j]) putchar(‘o‘);
 85         else putchar(‘_‘);
 86         if (j==m) puts("");
 87     }
 88 }
 89
 90 int main(){
 91     #ifndef ONLINE_JUDGE
 92     freopen("file.in","r",stdin);
 93 //    freopen("file.out","w",stdout);
 94     #endif
 95     scanf("%d%d",&n,&m);
 96     F(i,0,10) F(j,0,10)
 97         F(k,0,(1<<10)) f[i][j][k]=pre[i][j][k]=INF;
 98     F(i,1,n)
 99         F(j,1,m){
100             scanf("%d",&a[i][j]);
101             if (!a[i][j])
102                 f[i][j][1<<(cnt++)]=0;
103         }
104     solve();
105     return 0;
106 }