现在的位置: 首页 > 综合 > 正文

斯坦纳树 hdu4085

2019年04月14日 ⁄ 综合 ⁄ 共 1636字 ⁄ 字号 评论关闭
//斯坦纳树
//hdu 4085
/*复杂度:O(n^2*2^k+n*3^k)
枚举子树的形态:f[ i ][ j ]=min{ f[ i ][ j ],f[ k ][ j ]+f[ l ][ j ] },其中k和l是对i的一个划分。k^l==i子集枚举
按照边进行松弛:f[ i ][ j ]=min{ f[ i ][ j ],f[ i ][ j' ]+w[ j ][ j' ] },其中j和j'之间有边相连。
f[i][j]表示包含集合i对应的点,并包含j点的最小生成树
dp[i]=min{dp[k]+dp[k^i]} 表示包含集合i并满足的最小森林
*/
#include<iostream>
#include<stdio.h>
#include<set>
#include<string.h>
#include<math.h>
#include<string>
#include<queue>
#include<algorithm>
#define N 55
using namespace std;
const int inf=(int)1e8;
int f[1<<11][N],dp[1<<11],n,m,h,g[N][N],vis[N];
bool ok(int state)
{
    int num=0;
    for(int i=0;i<h;i++) if((state>>i)&1) num++;
    for(int i=0;i<h;i++) if((state>>i+h)&1) num--;
//    cout<<state<<" num="<<num<<endl;
    return num==0;
}
void calc()//计算斯坦纳树
{
    memset(vis,0,sizeof vis);
    for(int i=1;i<(1<<h+h);i++)
    {
        for(int j=1;j<=n;j++)
            for(int k=(i-1)&i;k;k=(k-1)&i)
                f[i][j]=min(f[k][j]+f[i^k][j],f[i][j]);

        for(int k=1;k<=n;k++)
        {
            int best=-1;
            for(int j=1;j<=n;j++) if(vis[j]!=i&&(best==-1||f[i][best]>f[i][j])) best=j;
            if(k==1) dp[i]=f[i][best];
            vis[best]=i;
//            cout<<best<<endl;
            for(int j=1;j<=n;j++) f[i][j]=min(f[i][j],f[i][best]+g[best][j]);
        }
//        cout<<"i="<<i<<" h="<<h<<endl;
//        for(int j=1;j<=n;j++) cout<<f[i][j]<<' '; puts("");
        if(ok(i))
            for(int k=i&(i-1);k;k=(k-1)&i) dp[i]=min(dp[k]+dp[k^i],dp[i]);
        else
            dp[i]=inf;
//        cout<<"h="<<h<<endl;
    }
    h=(1<<h*2)-1;
    if(dp[h]>=inf) puts("No solution");
    else printf("%d\n",dp[h]);
}
void solve()
{
    scanf("%d%d%d",&n,&m,&h);
    for(int i=0;i<N;i++) for(int j=0;j<N;j++) g[i][j]=(i==j?0:inf);
    int u,v,w;
    for(int i=0;i<m;i++)
    {
        scanf("%d%d%d",&u,&v,&w);
        if(w<g[u][v]) g[u][v]=g[v][u]=w;
    }
    for(int i=0;i<(1<<h*2);i++)
        for(int j=0;j<=n;j++) f[i][j]=inf;
    for(int i=0;i<h;i++) f[1<<i][i+1]=0;
    for(int i=0;i<h;i++) f[1<<i+h][n-i]=0;
    calc();
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--) solve();
    return 0;
}

抱歉!评论已关闭.