bzoj 4006: [JLOI2015]管道连接

Description

小铭铭最近进入了某情报部门，该部门正在被如何建立安全的通道连接困扰。

该部门有 n 个情报站，用 1 到 n 的整数编号。给出 m 对情报站 ui;vi 和费用 wi，表示情

报站 ui 和 vi 之间可以花费 wi 单位资源建立通道。

如果一个情报站经过若干个建立好的通道可以到达另外一个情报站，那么这两个情报站就

建立了通道连接。形式化地，若 ui 和 vi 建立了通道，那么它们建立了通道连接；若 ui 和 vi 均

与 ti 建立了通道连接，那么 ui 和 vi 也建立了通道连接。

现在在所有的情报站中，有 p 个重要情报站，其中每个情报站有一个特定的频道。小铭铭

面临的问题是，需要花费最少的资源，使得任意相同频道的情报站之间都建立通道连接。

Input

第一行包含三个整数 n;m;p，表示情报站的数量，可以建立的通道数量和重要情报站的数

量。接下来 m 行，每行包含三个整数 ui;vi;wi，表示可以建立的通道。最后有 p 行，每行包含

两个整数 ci;di，表示重要情报站的频道和情报站的编号。

Output

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

Sample Input

5 8 4
1 2 3
1 3 2
1 5 1
2 4 2
2 5 1
3 4 3
3 5 1
4 5 1
1 1
1 2
2 3
2 4

Sample Output

4

HINT

 

选择 (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。

 

最小斯坦纳树。 

 

#include<bits/stdc++.h>
using namespace std;  
int const N=1000+10;  
int const inf=1e9;  
int const M=1024+10;   
struct edge{
  int to,nt,w;  
}e[N<<4];  
int vis[N],h[N],cnt,n,m,p,f[N][M],a[11],b[11],st[11],id[11],g[M];  
queue<int> q;   
void add(int a,int b,int c){
  e[++cnt].to=b;  
  e[cnt].w=c;  
  e[cnt].nt=h[a];  
  h[a]=cnt; 
} 
int main(){
  scanf("%d%d%d",&n,&m,&p);  
  while (m--){
    int x,y,z;  
    scanf("%d%d%d",&x,&y,&z);  
    add(x,y,z);  
    add(y,x,z);  
  } 
  for(int i=1;i<=n;i++)  
    for(int j=0;j<(1<<p);j++)  
      f[i][j]=inf;  
  for(int i=1;i<=p;i++){ 
    scanf("%d%d",&a[i],&b[i]);  
    f[b[i]][1<<(i-1)]=0;  
    st[a[i]]+=(1<<i-1);  
    id[a[i]]=b[i];     
  } 
  for(int k=1;k<(1<<p);k++){
    for(int i=1;i<=n;i++) {
      for(int x=k&(k-1);x;x=(x-1)&k){
        int tmp=f[i][x]+f[i][k-x];  
        f[i][k]=min(f[i][k],tmp);  
      }
      if(f[i][k]<inf){
        q.push(i); vis[i]=1;  
      }    
    }  
    while (!q.empty()){
      int x=q.front();  
      q.pop();  
      vis[x]=0; 
      for(int i=h[x];i;i=e[i].nt){
        int v=e[i].to;  
        if(f[v][k]>f[x][k]+e[i].w){
          f[v][k]=f[x][k]+e[i].w; 
          if(!vis[v]){
            q.push(v);  
            vis[v]=1;  
          }  
        }  
      }  
    }  
  }
  for(int k=1;k<(1<<p);k++){
    int s=0,t=0;  
    for(int i=0;i<p;i++)         
      if(k&(1<<i)){
        s+=st[i+1];
        if(!t)t=id[i+1]; 
      }  
    g[k]=f[t][s];  
    for(int x=k&(k-1);x;x=(x-1)&k) 
      g[k]=min(g[k],g[x]+g[k-x]);  
  }             
  cout<<g[(1<<p)-1]<<endl;  
  return 0; 
}