ZOJ 3613 Wormhole Transport

斯坦纳树，$dp$。

先求出每个状态下连通的最小花费，因为可以是森林，所以$dp$一下。

#include<bits/stdc++.h>
using namespace std;

int n;
const int INF = 0x7FFFFFFF;
struct X
{
    int p,s,q;
}e[];
int id[];
vector<int>g[];
int val[][],d[][],dp[];
int f[*],k;
queue<int>Q;

int num[],cost[];

void spfa()
{
    while(!Q.empty())
    {
        int h = Q.front(); Q.pop(); f[h]=;
        int x=h/,y=h%;
        for(int i=;i<g[x].size();i++)
        {
            int to = g[x][i];

            if(to<k)
            {
                if(((<<to)&y)==)
                {
                    if(d[x][y]+val[x][to]<d[to][y|(<<to)])
                        d[to][y|(<<to)]=d[x][y]+val[x][to];
                }
            }

            else
            {
                if(d[x][y]+val[x][to]<d[to][y])
                {
                    d[to][y] = d[x][y]+val[x][to];
                    if(f[to*+y]==)
                    {
                        f[to*+y]=;
                        Q.push(to*+y);
                    }
                }
            }
        }
    }
}

bool cmp(X a,X b)
{
    if(a.p!=b.p) return a.p>b.p;
    return a.s>b.s;
}

void work(int x)
{
    int A=,B=;
    for(int i=;i<k;i++)
    {
        if(((<<i)&x)==) continue;
        A=A+e[i].p; B=B+e[i].s;
    }
    num[x] = min(A,B);
    cost[x] = dp[x];
}

int main()
{
    while(~scanf("%d",&n))
    {
        for(int i=;i<n;i++)
        {
            scanf("%d%d",&e[i].p,&e[i].s);
            e[i].p = min(,e[i].p);
            e[i].q = i;
        }
        sort(e,e+n,cmp);
        for(int i=;i<n;i++) id[e[i].q]=i;

        for(int i=;i<n;i++)
            for(int j=;j<n;j++) val[i][j]=INF;
        memset(f,,sizeof f);
        for(int i=;i<n;i++) g[i].clear();

        int m; scanf("%d",&m);
        for(int i=;i<=m;i++)
        {
            int a,b,c; scanf("%d%d%d",&a,&b,&c);
            a--; b--; a=id[a]; b=id[b];
            val[a][b]=min(val[a][b],c);
            val[b][a]=val[a][b];
            g[a].push_back(b);
            g[b].push_back(a);
        }

        k=; for(int i=;i<n;i++) if(e[i].p||e[i].s) k++;

        int st = <<k;
        for(int j=;j<st;j++)
            for(int i=;i<n;i++) d[i][j]=INF;

        for(int i=;i<n;i++)
        {
            if(i<k) d[i][<<i]=;
            else d[i][]=;
        }

        for(int j=;j<st;j++)
        {
            for(int i=;i<n;i++)
            {
                if(i<k)
                {
                    if(((<<i)&j)==) continue;
                    for (int x = j; x; x = (x-)&j)
                    {
                        int A=x ,B=j-A;
                        if(d[i][A|(<<i)]!=INF&&d[i][B|(<<i)]!=INF)
                            d[i][j] = min(d[i][j], d[i][A|(<<i)]+d[i][B|(<<i)]);
                    }
                }
                else
                {
                    for (int x = j; x; x = (x-)&j)
                    {
                        int A=x ,B=j-A;
                        if(d[i][A]!=INF&&d[i][B]!=INF)
                            d[i][j] = min(d[i][j], d[i][A]+d[i][B]);
                    }
                }
                if(d[i][j]!=INF) Q.push(i*+j);
            }
            spfa();
        }

        for(int j=;j<st;j++)
        {
            dp[j]=INF;
            for(int i=;i<n;i++) dp[j]=min(dp[j],d[i][j]);
        }

        memset(num,,sizeof num);
        memset(cost,,sizeof cost);

        for(int i=;i<st;i++)
        {
            if(dp[i]!=INF) work(i);
            for (int x = i; x; x = (x-)&i)
            {
                int A=x, B=i-x;
                if(num[A]+num[B]>num[i])
                {
                    num[i] = num[A]+num[B];
                    cost[i] = cost[A]+cost[B];
                }
                else if(num[A]+num[B]==num[i])
                {
                    if(cost[A]+cost[B]<cost[i])
                        cost[i] = cost[A]+cost[B];
                }
            }
        }

        int ans1=,ans2=;
        for(int i=;i<st;i++) ans1=max(ans1,num[i]);
        for(int i=;i<st;i++)
        {
            if(num[i]!=ans1) continue;
            ans2=min(ans2,cost[i]);
        }
        printf("%d %d\n",ans1,ans2);

    }
    return ;
}