bzoj 2595: [Wc2008]游览计划

wc！！！斯坦纳树！！！！！！！！！！！！！！！！！！！！

while (1) {iq--，膝盖++；}

还是不太理解的样子。。。。（过几天自己就会了，，不虚（大雾））

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<map>
#include<queue>
#define N 1000005
#define inf 1000000000
#define LL long long
using namespace std;
int a[15][15];
int bin[20];
int n,m,K;
int f[12][12][1024];
struct data{
    int fir,sed,thd;
}pre[15][15][1024];
int xx[4]={0,0,1,-1},yy[4]={1,-1,0,0};
queue<pair<int , int > >q;
bool inq[15][15],vis[15][15];
void spfa(int st)
{
    while (!q.empty())
    {
        int x=q.front().first,y=q.front().second;
        inq[x][y]=0; q.pop();
        for (int k=0; k<4; k++)
        {
            int nowx=x+xx[k],nowy=y+yy[k];
            if (nowx<1 || nowy<1 || nowx>n || nowy>m) continue;
            if (f[nowx][nowy][st]>f[x][y][st]+a[nowx][nowy])
            {
                f[nowx][nowy][st]=f[x][y][st]+a[nowx][nowy];
                pre[nowx][nowy][st]=(data){x,y,st};
                if (!inq[nowx][nowy])
                {
                    q.push(make_pair(nowx,nowy));
                    inq[nowx][nowy]=1;
                }
            }
        }
    }
}
void dfs(int x, int y, int st)
{
    if (x>inf || pre[x][y][st].thd==0) return;
    vis[x][y]=1;
    data t=pre[x][y][st];
    dfs(t.fir,t.sed,t.thd);
    if (t.fir==x && t.sed==y) dfs(x,y,st-t.thd);
}
int main()
{
    bin[0]=1; for (int i=1; i<20; i++) bin[i]=bin[i-1]<<1;
    scanf("%d%d",&n,&m);
    for (int i=1; i<=n; i++)    
        for (int j=1; j<=m; j++)
        {
            scanf("%d",&a[i][j]);
            if (!a[i][j]) K++;
        }
    for (int i=1; i<=n; i++)
        for (int j=1; j<=m; j++)
            for (int k=0; k<bin[K]; k++)
                f[i][j][k]=pre[i][j][k].fir=inf;
    K=0;
    for (int i=1; i<=n; i++)
        for (int j=1; j<=m; j++)
            if (!a[i][j]) f[i][j][bin[K++]]=0;
    for (int st=1; st<bin[K]; st++)
    {
        for (int i=1; i<=n; i++)
            for (int j=1; j<=m; j++)
            {
                for (int s=st&(st-1);s;s=st&(s-1))
                {
                    int t=f[i][j][s]+f[i][j][st-s]-a[i][j];      //?????
                    if (t<f[i][j][st])
                    {
                        printf("%d %d %d",i,j,s); while (1);
                        f[i][j][st]=t;
                        pre[i][j][st]=(data){i,j,s};
                    }
                }
                if (f[i][j][st]<inf) q.push(make_pair(i,j)),inq[i][j]=1;
            }
            spfa(st);
    }
    int x,y;
    for (int i=1; i<=n; i++)
        for (int j=1; j<=m; j++)
            if (!a[i][j])
            {
                x=i; y=j; break;
            }
    dfs(x,y,bin[K]-1);
    printf("%d
",f[x][y][bin[K]-1]);
    for (int i=1; i<=n; i++)
    {
        for (int j=1; j<=m; j++)
        {
            if (!a[i][j]) printf("x");
            else if (vis[i][j]) printf("o");
            else printf("_");
        }
        puts("");
    }
    return 0;
}