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 ;
}
05-21 20:37