【JLOI2015】通道连接

题目

Time Limits: 1000 ms Memory Limits: 262144 KB Detailed Limits

Description

Input

Output

Sample Input

输入1：
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
输入2：
5 8 4
1 2 2
1 3 2
1 5 1
2 4 2
2 5 2
3 4 2
3 5 1
4 5 1
1 1
1 2
2 3
2 4

Sample Output

输出1：
4
【样例说明1】
选择(1,5), (3, 5), (2, 5), (4, 5) 这4 对情报站连接。
输出2：
4
【样例说明二】
选择(1,2), (3, 4) 这2 对情报站连接。

Data Constraint

题解

多个点的联通容易联想到斯坦纳树
但是已经忘记了怎么打了
看了一下之后发现是松弛（spfa）加上自己与自己的配对
其实就是包含很多劣解以及最优解的一个做法

可是这题不用形成一棵树，只要每种颜色都在同一个联通块就可以了
那么再设dp[x]表示当前完成的状况为x，保证每一次新加的森林都是包含若干完整颜色的就好了

贴代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
#define min(x,y) ((x)<(y)?(x):(y))
using namespace std;

const int maxp=13,maxn=3005;

int fi[maxn],ne[maxn*2],dui[maxn*2],dui1[maxn*2],qc[maxn];
int f[1005][1<<maxp],a[maxn][3],cc[maxn],dp[1<<maxp];
int d[maxn],h[maxn*10],q[maxn*10];
bool zh[1<<maxp],bz[maxn];
int i,j,k,l,m,n,x,y,z,p,now,cz,be,ed;

void add(int x,int y,int z){
    if (fi[x]==0) fi[x]=++now; else ne[qc[x]]=++now;
    dui[now]=y; dui1[now]=z; qc[x]=now;
}

void spfa(int z){
    be=0;
    while (be<ed){
        be++;
        k=fi[h[be]]; x=h[be];
        while (k){
            y=dui[k];
            if (f[x][z]+dui1[k]<f[y][z]){
                f[y][z]=f[x][z]+dui1[k];
                if (bz[y]==false){
                    bz[y]=true; h[++ed]=y;
                }
            }
            k=ne[k];
        }
        bz[x]=false;
    }
}
int main(){
    freopen("4079.in","r",stdin);
    scanf("%d%d%d",&n,&m,&p);
    fo(i,1,m){
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z); add(y,x,z);
    }
    fo(i,1,p){
        scanf("%d%d",&a[i][1],&a[i][2]);
        d[a[i][1]]=i;
    }
    cc[0]=1;
    fo(i,1,10) cc[i]=cc[i-1]*2;
    fo(i,1,10){
        x=0;
        fo(j,1,p) if (a[j][1]==i) x+=cc[j-1];
        zh[x]=true;
    }
    memset(f,3,sizeof(f)); cz=f[1][1];
    fo(i,1,n) f[i][0]=0;
    fo(i,1,p) f[a[i][2]][cc[i-1]]=0;
    fo(i,1,cc[p]-1){
        ed=0; memset(bz,false,sizeof(bz));
        fo(j,1,n){
            x=(i-1)&i;
            while (x){
                f[j][i]=min(f[j][i],f[j][i-x]+f[j][x]);
                x=(x-1)&i;
            }
            if (f[j][i]<cz){
                h[++ed]=j; bz[j]=true;
            }
        }
        spfa(i);
    }
    fo(i,1,cc[p]-1){
        x=(i-1)&i;
        while (x){
            if (zh[x]==true && zh[i-x]==true) zh[i]=true;
            if (zh[i]==true) break;
            x=(x-1)&i;
        }
    }
    memset(dp,66,sizeof(dp)); dp[0]=0;
    fo(i,1,cc[p]-1)if (zh[i]){
        fo(j,1,n) dp[i]=min(dp[i],f[j][i]);
        x=(i-1)&i;
        while (x){
            if (zh[x]) dp[i]=min(dp[i],dp[x]+dp[i-x]);
            x=(x-1)&i;
        }
    }
    printf("%d",dp[cc[p]-1]);
    return 0;
}