BZOJ2595 WC2008游览计划（斯坦纳树）

　　斯坦纳树板子题。

　　考虑状压dp，设f[i][j][S]表示当前在点(i,j)考虑转移，其所在的联通块包含的关键点集（至少）为S的答案。

　　转移时首先枚举子集，有f[i][j][S]=min{f[i][j][x]+f[i][j][y]-a[i][j]} (x&y=0,x|y=S)。

　　然后考虑从点(i,j)从哪拓展而来，有f[i][j][S]=min{f[x][y][S]}+a[i][j]，其中(x,y)为(i,j)的相邻点，使用spfa转移。

　　这里第二种转移仅在相同关键点集下进行，因为由更小点集转移而来的情况已在第一种转移中被考虑。

　　对于输出方案，记录如何转移而来即可。

　　其实并没有搞懂。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 12
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int n,m,s,a[N][N],id[N][N],f[N][N][<<N],from[N][N][<<N][],q[N*N*N*N][];
int wx[]={,,,-},wy[]={,,-,};
bool flag[N][N];
char b[N][N];
void inc(int &x){x++;if (x>n*m+) x-=n*m+;}
void spfa(int k)
{
    int head=,tail=;
    for (int i=;i<=n;i++)
        for (int j=;j<=m;j++)
        tail++,q[tail][]=i,q[tail][]=j,flag[i][j]=;
    do
    {
        inc(head);int x=q[head][],y=q[head][];flag[x][y]=;
        for (int i=;i<;i++)
        {
            int u=x+wx[i],v=y+wy[i];
            if (u&&u<=n&&v&&v<=m&&f[x][y][k]+a[u][v]<f[u][v][k])
            {
                f[u][v][k]=f[x][y][k]+a[u][v];
                from[u][v][k][]=x,from[u][v][k][]=y;
                if (!flag[u][v]) flag[u][v]=,inc(tail),q[tail][]=u,q[tail][]=v;
            }
        }
    }while (head!=tail);
}
void getans(int x,int y,int k)
{
    b[x][y]='o';
    while (from[x][y][k][]>)
    {
        int u=from[x][y][k][],v=from[x][y][k][];
        b[x=u][y=v]='o';
    }
    if (from[x][y][k][]<) getans(x,y,-from[x][y][k][]),getans(x,y,-from[x][y][k][]);
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj2595.in","r",stdin);
    freopen("bzoj2595.out","w",stdout);
#endif
    n=read(),m=read();
    for (int i=;i<=n;i++)
        for (int j=;j<=m;j++)
        {
            a[i][j]=read(),s+=(a[i][j]==);
            if (!a[i][j]) id[i][j]=s;
        }
    memset(f,,sizeof(f));
    for (int i=;i<=n;i++)
        for (int j=;j<=m;j++)
        {
            f[i][j][]=a[i][j];
            if (id[i][j]) f[i][j][<<id[i][j]-]=a[i][j];
        }
    for (int k=;k<(<<s);k++)
    {
        for (int i=;i<=n;i++)
            for (int j=;j<=m;j++)
                for (int x=k;x>=(k^x);x=x-&k)
                if (f[i][j][x]+f[i][j][k^x]-a[i][j]<f[i][j][k])
                {
                    f[i][j][k]=f[i][j][x]+f[i][j][k^x]-a[i][j];
                    from[i][j][k][]=-x,from[i][j][k][]=-(k^x);
                }
        spfa(k);
    }
    int ans=;
    for (int i=;i<=n;i++)
        for (int j=;j<=m;j++)
        ans=min(ans,f[i][j][(<<s)-]);
    cout<<ans<<endl;
    for (int i=;i<=n;i++)
        for (int j=;j<=m;j++)
        b[i][j]='_';
    for (int i=;i<=n;i++)
        for (int j=;j<=m;j++)
        if (f[i][j][(<<s)-]==ans)
        {
            getans(i,j,(<<s)-);
            for (int i=;i<=n;i++)
                for (int j=;j<=m;j++)
                if (!a[i][j]) b[i][j]='x';
            for (int i=;i<=n;i++)
            {
                for (int j=;j<=m;j++)
                putchar(b[i][j]);
                cout<<endl;
            }
            return ;
        }
}