[JLOI2015]管道连接

题目描述

小铭铭最近进入了某情报部门，该部门正在被如何建立安全的通道连接困扰。该部门有 n 个情报站，用 1 到 n 的整数编号。给出 m 对情报站 ui;vi 和费用 wi，表示情报站 ui 和 vi 之间可以花费 wi 单位资源建立通道。

如果一个情报站经过若干个建立好的通道可以到达另外一个情报站，那么这两个情报站就建立了通道连接。形式化地，若 ui 和 vi 建立了通道，那么它们建立了通道连接；若 ui 和 vi 均与 ti 建立了通道连接，那么 ui 和 vi 也建立了通道连接。

现在在所有的情报站中，有 p 个重要情报站，其中每个情报站有一个特定的频道。小铭铭面临的问题是，需要花费最少的资源，使得任意相同频道的情报站之间都建立通道连接。

输入输出格式

输入格式：

第一行包含三个整数 n;m;p，表示情报站的数量，可以建立的通道数量和重要情报站的数量。接下来 m 行，每行包含三个整数 ui;vi;wi，表示可以建立的通道。最后有 p 行，每行包含两个整数 ci;di，表示重要情报站的频道和情报站的编号。

输出格式：

输出一行一个整数，表示任意相同频道的情报站之间都建立通道连接所花费的最少资源总量。

输入输出样例

输入样例#1：

输出样例#1：

说明

选择 (1; 5); (3; 5); (2; 5); (4; 5) 这 4 对情报站连接。

对于 100% 的数据，0 <ci <= p <= 10; 0 <ui;vi;di <= n <= 1000; 0 <= m <= 3000; 0 <= wi <=

20000。

首先做一遍裸的斯坦纳树

f[i][sta]=min{f[i][sub]+f[i][sta-sub]}

因为题目可能出现斯坦纳森林，所以还要进行一次子集dp

将每一种频道作为二进制为储存，还要计算出ki[i]表示i频道需要覆盖的点(二进制数)

枚举1<=sta<(1<<p),k=sta频道状态需要覆盖的点

dp[sta]=min(f[j][k]) j=1~n

dp[sta]=min(dp[sta],dp[sub]+dp[sta-sub])

答案是dp[(1<<p)]

 #include<iostream>

 #include<cstdio>

 #include<cstring>

 #include<algorithm>

 #include<queue>

 using namespace std;

 struct Node

 {

   int next,to,dis;

 }edge[];

 int num,head[],f[][],dp[],ki[],n,m,p;

 bool vis[];

 void add(int u,int v,int dis)

 {

   num++;

   edge[num].next=head[u];

   head[u]=num;

   edge[num].to=v;

   edge[num].dis=dis;

 }

 void bfs(int s)

 {int i;

   memset(vis,,sizeof(vis));

   queue<int> Q;

   for (i=;i<=n;i++)

     vis[i]=,Q.push(i);

   while (Q.empty()==)

     {

       int u=Q.front();

       //cout<<u<<endl;

       Q.pop();

       vis[u]=;

       for (i=head[u];i;i=edge[i].next)

     {

       int v=edge[i].to;

       if (f[v][s]>f[u][s]+edge[i].dis)

         {

           f[v][s]=f[u][s]+edge[i].dis;

           if (vis[v]==)

         {

           vis[v]=;

           Q.push(v);

         }

         }

     }

     }

 }

 int main()

 {int i,u,v,w,x,y,j,k;

   cin>>n>>m>>p;

   for (i=;i<=m;i++)

     {

       scanf("%d%d%d",&u,&v,&w);

       add(u,v,w);add(v,u,w);

     }

   memset(f,/,sizeof(f));

   for (i=;i<=p;i++)

     {

       scanf("%d%d",&x,&y);

       ki[x]|=(<<i-);

       f[y][<<i-]=;

     }

   for (i=;i<(<<p);i++)

     {

       for (j=;j<=n;j++)

     {

       for (k=i;k;k=(k-)&i)

         if (f[j][i]>f[j][k]+f[j][i-k])

           f[j][i]=f[j][k]+f[j][i-k];

     }

       bfs(i);

     }

   memset(dp,/,sizeof(dp));

   for (i=;i<(<<p);i++)

     {int k=;

       for (j=;j<=p;j++)

     if (i&(<<j-)) k|=ki[j];

       for (j=;j<=n;j++)

     if (dp[i]>f[j][k])

     dp[i]=f[j][k];

       for (j=i;j;j=(j-)&i)

     if (dp[j]+dp[i-j]<dp[i])

     dp[i]=dp[j]+dp[i-j];

     }

   cout<<dp[(<<p)-];

 }