【UOJ #108】【APIO 2013】TOLL

http://uoj.ac/problem/108
好神的一道题啊。
原图边权互不相同是重点！
如果有一个点集，有两组边集，要求这两组边集的并集的最小生成树，可以对两组边集分别求一下最小生成树构成新的两组边集，再合并两组新的边集。
还有，对于一个给定的最小生成树，要加入一条边，求加入一条边之后的最小生成树，那么就有这条边把原生成树上的一条边踢掉，或者这条边踢不动原生成树上的任何一条边（这条边加入后构成的环中的任意一条边的边权都小于这条边）。
这道题先考虑一个暴力的做法：先\(2^k\)枚举哪些新道路一定在最小生成树中，先在最小生成树中加入这些新道路，在这个基础上再加原图的边形成最小生成树，然后树的形态确定了，就要确定新道路的值了。对于原图中的边没能加入最小生成树的，肯定要想踢掉它连上后形成的环上的边，但环上原图中的边一定不会被踢掉，只可能踢掉环上的新道路。新道路为了防止自己被踢，就把自己的边权和这条想踢掉它的边的边权取min。
正解就是：改造原图！
先将所有新道路加入最小生成树，原图中最小生成树上没有通过新道路连通的点都缩成一个点（因为这些点之间的边不会对新道路的边权产生任何影响）。缩点之后的原图去掉新道路再求一遍最小生成树，设边集为\(M\)，在最外层\(2^k\)枚举哪些新道路一定在最小生成树中，然后拿\(M\)中的边继续加入，无法加入就用来更新新道路的权值。在这里用\(M\)中的边和用缩点后原图的边是等价的，因为\(M\)中的边能代表原图的连通性，而且原图中的边能加入最小生成树的边一定在\(M\)中（是\(M\)的子集），能更新新道路边权的边也一定能被\(M\)中的一条边代替。
时间复杂度\(O(m\log m+2^kk^2)\)。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;

const int N = 100003;
const int M = 300003;
const int K = 23;

struct Enode {
    int u, v, e;
    bool operator < (const Enode &AA) const {
        return e < AA.e;
    }
} G[M], A[K], New[M], MN[K];
int p[N], Gtot = 0, n, m, k, fa[N], fa2[N], Atot = 0, Newtot = 0, MNtot = 0, tot = 0, remark[N];
ll sp[N];

int find(int x) {return x == fa[x] ? x : fa[x] = find(fa[x]);}
int find2(int x) {return x == fa2[x] ? x : fa2[x] = find2(fa2[x]);}

ll sum[K];
bool used[K], mark[K];
struct node {int nxt, to, w;} E[K * K << 1];
int cnt, point[K], deep[K], fadis[K];
void ins(int u, int v, int w) {E[++cnt] = (node) {point[u], v, w}; point[u] = cnt;}

void dfs(int x) {
    sum[x] = sp[x];
    for (int i = point[x]; i; i = E[i].nxt)
        if (E[i].to != fa[x]) {
            fa[E[i].to] = x;
            fadis[E[i].to] = E[i].w;
            deep[E[i].to] = deep[x] + 1;
            if (E[i].w == 0x7fffffff) mark[E[i].to] = true;
            else mark[E[i].to] = false;
            dfs(E[i].to);
            sum[x] += sum[E[i].to];
        }
}

void minit(int u, int v, int e) {
    if (deep[u] < deep[v]) u ^= v ^= u ^= v;
    while (deep[u] > deep[v]) {
        if (mark[u] && fadis[u] > e)
            fadis[u] = e;
        u = fa[u];
    }
    if (u == v) return;
    while (u != v) {
        if (mark[u] && fadis[u] > e)
            fadis[u] = e;
        if (mark[v] && fadis[v] > e)
            fadis[v] = e;
        u = fa[u]; v = fa[v];
    }
}

int main() {
    scanf("%d%d%d", &n, &m, &k);
    int u, v, e, etot = 0;
    for (int i = 1; i <= m; ++i) {
        scanf("%d%d%d", &u, &v, &e);
        G[++Gtot] = (Enode) {u, v, e};
    }
    
    stable_sort(G + 1, G + Gtot + 1);
    for (int i = 1; i <= n; ++i) fa[i] = fa2[i] = i;
    for (int i = 1; i <= k; ++i) {
        scanf("%d%d", &u, &v);
        A[++Atot] = (Enode) {u, v, 0};
        u = find(u); v = find(v);
        if (u != v) fa[u] = v, ++etot;
    }
    
    for (int i = 1; i <= n; ++i) scanf("%d", p + i);
    
    for (int i = 1; i <= Gtot; ++i) {
        u = find(G[i].u); v = find(G[i].v);
        if (u != v) {
            ++etot;
            fa[u] = v;
            u = find2(G[i].u); v = find2(G[i].v);
            if (u != v) fa2[u] = v;
            if (etot == n - 1) break;
        }
    }
    
    for (int i = 1; i <= n; ++i) {
        if (!remark[find2(i)])
            remark[fa2[i]] = ++tot;
        remark[i] = remark[fa2[i]];
    }
    
    for (int i = 1; i <= n; ++i) sp[remark[i]] += p[i];
    
    for (int i = 1; i <= m; ++i) {
        u = remark[G[i].u]; v = remark[G[i].v];
        if (u != v)
            New[++Newtot] = (Enode) {u, v, G[i].e};
    }
    
    stable_sort(New + 1, New + Newtot + 1);
    for (int i = 1; i <= k; ++i) A[i].u = remark[A[i].u], A[i].v = remark[A[i].v];
    
    etot = 0;
    for (int i = 1; i <= tot; ++i) fa[i] = i;
    for (int i = 1; i <= Newtot; ++i) {
        u = find(New[i].u); v = find(New[i].v);
        if (u != v) {
            fa[u] = v;
            ++etot;
            MN[++MNtot] = New[i];
            if (etot == tot - 1) break;
        }
    }
    
    int S = (1 << k); ll ans = 0;
    for (int s = 1; s < S; ++s) {
        for (int i = 1; i <= tot; ++i) fa[i] = i, point[i] = 0;
        cnt = etot = 0;
        for (int i = 0; i < k; ++i)
            if ((s >> i) & 1) {
                u = A[i + 1].u; v = A[i + 1].v;
                if (find(u) != find(v)) {
                    fa[fa[u]] = fa[v];
                    ++etot;
                    ins(u, v, 0x7fffffff);
                    ins(v, u, 0x7fffffff);
                }
            }
        
        for (int i = 1; i <= MNtot; ++i) used[i] = false;
        if (etot < tot - 1) {
            for (int i = 1; i <= MNtot; ++i) {
                u = find(MN[i].u); v = find(MN[i].v);
                if (u != v) {
                    used[i] = true;
                    ++etot;
                    fa[u] = v;
                    ins(MN[i].u, MN[i].v, MN[i].e);
                    ins(MN[i].v, MN[i].u, MN[i].e);
                    if (etot == tot - 1) break;
                }
            }
        }
        
        fa[remark[1]] = 0; dfs(remark[1]);
        for (int i = 1; i <= MNtot; ++i)
            if (!used[i])
                minit(MN[i].u, MN[i].v, MN[i].e);
        
        ll ret = 0;
        for (int i = 1; i <= tot; ++i)
            if (mark[i])
                ret += sum[i] * fadis[i];
        if (ret > ans) ans = ret;
    }
    
    printf("%lld\n", ans);
    
    return 0;
}

转载于:https://www.cnblogs.com/abclzr/p/6720239.html