啊 斯坦纳树 好像很厉害啊 反正我之前不会。。。
其实不知道实用性怎么样 毕竟复杂度不小。。。
大概过程就是一个状压DP+spfa 枚举状态 从小状态更新大状态 再对当前状态做一次像spfa一样的松弛操作
下面这个人讲的不错 可以去看看
http://blog.csdn.net/gzh1992n/article/details/9119543
为什么我的代码又那么短 有点担心优美度了 有谁提一下建议吗2333 我觉得还过得去?
#include<bits/stdc++.h>
#define me(a,x) memset(a,x,sizeof a)
using namespace std;
const int N=11,inf=1e9+7;
const int dx[4]={0,0,1,-1};
const int dy[4]={1,-1,0,0};
inline int read()
{
char ch=getchar(); int x=0,f=1;
while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}
return x*f;
}
struct S{int x,y,k;}g[N][N][1030];
int f[N][N][1030],a[N][N],n,m,b[N][N],cnt;
bool v[N][N];
bool upd(int x,int y,int u,int i,int j,int k,int w){
if(f[x][y][u]>w)return f[x][y][u]=w,g[x][y][u]=(S){i,j,k},1;
return 0;
}
queue<S>q;
void spfa(int t){
while(!q.empty()){
S u=q.front(); q.pop();
v[u.x][u.y]=0;
for(int k=0;k<4;k++){
int x=u.x+dx[k],y=u.y+dy[k];
if(x<1 || y<1 || x>n || y>m)continue;
if(upd(x,y,t,u.x,u.y,t,f[u.x][u.y][t]+a[x][y]) && !v[x][y])
v[x][y]=1,q.push((S){x,y,t});
}
}
}
void dfs(int x,int y,int k){
if(!k)return; v[x][y]=1;
dfs(g[x][y][k].x,g[x][y][k].y,g[x][y][k].k);
if(x==g[x][y][k].x && y==g[x][y][k].y)dfs(x,y,k-g[x][y][k].k);
}
void op(){
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
if(!a[i][j])printf("x");
else if(v[i][j])printf("o");
else printf("_");
}
printf("\n");
}
}
int main()
{
n=read(),m=read(); me(f,63);
int i,j,k,s; cnt=0;
for(i=1;i<=n;i++)for(j=1;j<=m;j++){
a[i][j]=read();
if(!a[i][j]){
b[i][j]=1<<cnt,f[i][j][1<<cnt]=0;
cnt++; g[i][j][1<<cnt].k=0;
}
f[i][j][0]=a[i][j];
}
int p=(1<<cnt);
for(k=1;k<p;k++){
me(v,0);
for(i=1;i<=n;i++)for(j=1;j<=m;j++){
if(!a[i][j] && !(k&b[i][j]))continue;
for(s=k;s;s=(s-1)&k)
upd(i,j,k,i,j,s,f[i][j][s]+f[i][j][k-s]-a[i][j]);
if(f[i][j][k]<inf)q.push((S){i,j,k}),v[i][j]=1;
}
spfa(k);
}
for(i=1;i<=n;i++)for(j=1;j<=m;j++)if(!a[i][j]){
printf("%d\n",f[i][j][p-1]);
me(v,0); dfs(i,j,p-1);
op();
return 0;
}
}