刷题记录（NC16527 货车运输）（树上倍增）

NC16527 货车运输

题目链接

关键点：

1、题目大意是给出图的点和边，要求求两点间经过的路径的最小权值最大值。

2、通过画图观察样例可以发现，要求两点间最小权值的最大值，首先这两点要连通，其次才是找一条有着尽可能大的权值分布的路径，再在这条路径上找权值的最小值，那么选的边就为最大生成图上的边。

3、对于最大生成图的计算，我们可以采用kruskal算法，那么判断询问输出的-1，就可以判断两点的父亲节点是否相同

4、然后是如何求该最大生成图上，两点间边的最小权值，利用树上倍增的思想，维护一个数组f，f[i][j]表示i节点向上2^j的节点，w[i][j]表示i节点向上2^j路径上的最小权值

则f[i][j] = f[f[i][j-1]][j-1];

w[i][j] = min(w[i][j-1], w[f[i][j-1]][j-1]);

5、初始值，f[i][0] = i节点的父亲，w[i][0] = （若i节点为根结点，那么其为一个最大值，若不是根结点，为与父亲相连的边的权值）

6、这样一来，对于两点间的最小权值的最大值，就类似为两点间的lca求法，在两点深度不断逼近的过程中，更新最小权值答案。

完整代码：

# include <bits/stdc++.h>
# include <algorithm>
using namespace std;
const int N = 100000*2+10;
const int INF = 100000+10;
int n, m, q, cnt;
int f[N][20], w[N][20];
int fa[N], head[N], dep[N], a[N];//f代表i点往上2^j的点，w代表i点往上2^j的最小权值
bool vis[N];
struct ty1{
    int x, y, len;
}edge1[N];//原来图的边
struct ty2{
    int to, len, next;
}edge2[N];//最大生成树的边
void addedge(int x, int y, int len)
{
    edge2[++cnt].to = y;
    edge2[cnt].len = len;
    edge2[cnt].next = head[x];
    head[x] = cnt;
};
bool cmp(ty1 t1, ty1 t2)
{
    return t1.len>t2.len;
}
int find(int x)
{
    return x == fa[x]? x: find(fa[x]);
}
void kruskal()
{
    sort(edge1+1, edge1+1+m, cmp);
    for (int i=1; i<=n; i++)
        fa[i] = i;
    for (int i=1; i<=m; i++)
    {
        int fx = find(edge1[i].x);
        int fy = find(edge1[i].y);
//         cout<<fx<<" "<<fy<<endl;
        if (fx != fy)
        {
            fa[fx] = fy;
            addedge(edge1[i].x, edge1[i].y, edge1[i].len);
            addedge(edge1[i].y, edge1[i].x, edge1[i].len);
        }
    }
}
void dfs(int x)
{
    vis[x] = true;
    for (int i=head[x]; i!=-1; i=edge2[i].next)
    {
        int y = edge2[i].to;
        if (vis[y]) continue;
        dep[y] = dep[x] + 1;
        f[y][0] = x;
        w[y][0] = edge2[i].len;
        dfs(y);
    }
}
int lca(int x, int y)
{
    int ans = INF;
//     cout<<find(x)<<" "<<find(y)<<endl;
    if (find(x) != find(y))
        return -1;
    if (dep[x] > dep[y]) swap(x, y);
    for (int i=19; i>=0; i--)
    {
        if (dep[f[y][i]] >= dep[x])
        {
            ans = min(ans, w[y][i]);
            y = f[y][i];
        }
    }
    if (y == x)
        return ans;
    for (int i=19; i>=0; i--)
    {
        if (f[x][i] != f[y][i])
        {
            ans = min(ans, min(w[x][i], w[y][i]));
            x = f[x][i];
            y = f[y][i];
        }
    }
    ans = min(ans, min(w[x][0], w[y][0]));
    return ans;
}
int main()
{
    memset(head, -1, sizeof(head));
    scanf("%d%d", &n, &m);
    for (int i=1; i<=m; i++)
    {
        scanf("%d%d%d", &edge1[i].x, &edge1[i].y, &edge1[i].len);
    }
    kruskal();
    for (int i=1; i<=n; i++)
    {
        if (!vis[i])
        {
            dep[i] = 1;
            dfs(i);
            f[i][0] = i;
            w[i][0] = INF;
        }
    }
    for (int i=1; i<=19; i++)
    {
        for (int j=1; j<=n; j++)
        {
            f[j][i] = f[f[j][i-1]][i-1];
            w[j][i] = min(w[j][i-1], w[f[j][i-1]][i-1]);
        }
    }
    scanf("%d", &q);
    for (int i=1; i<=q; i++)
    {
        int x, y;
        scanf("%d%d", &x, &y);
        cout<<lca(x, y)<<endl;
    }
    
    
    
    return 0;
}