vijos1843(货车运输)

每个货车可以运输的最大货物量就是起点到终点的某一条路径上的最小权值,为了使运输量最大,那就是选一条路使得这条路上的最小权值边最大。

于是想到最大生成树,然后就是树上查询两点之间路径上的最小值,熟练剖分可以解。

#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
#define Del(a,b) memset(a,b,sizeof(a))
const int N = 10005;
int dep[N],siz[N],fa[N],id[N],son[N],val[N],top[N];
int num;
vector<int> v[N];
struct Edge{
    int from;
    int to;
    int w;
    bool operator < (const Edge& T) const{
        return w>T.w;
    }
}ed[50000+50];
struct tree
{
    int x,y,val;
    void read(int a,int b,int c){
        x = a; y = b; val = c;
    }
};
tree e[N];
void dfs1(int u, int f, int d) {
    dep[u] = d;
    siz[u] = 1;
    son[u] = 0;
    fa[u] = f;
    for (int i = 0; i < v[u].size(); i++) {
        int ff = v[u][i];
        if (ff == f) continue;
        dfs1(ff, u, d + 1);
        siz[u] += siz[ff];
        if (siz[son[u]] < siz[ff])
            son[u] = ff;
    }
}
void dfs2(int u, int tp) {
    top[u] = tp;
    id[u] = ++num;
    if (son[u]) dfs2(son[u], tp);
    for (int i = 0; i < v[u].size(); i++) {
        int ff = v[u][i];
        if (ff == fa[u] || ff == son[u]) continue;
        dfs2(ff, ff);
    }
}
#define lson(x) ((x<<1))
#define rson(x) ((x<<1)+1)
struct Tree
{
    int l,r,val;
};
Tree tree[4*N];
void pushup(int x) {
    tree[x].val = min(tree[lson(x)].val, tree[rson(x)].val);
}
void build(int l,int r,int v)
{
    tree[v].l=l;
    tree[v].r=r;
    if(l==r)
    {
        tree[v].val = val[l];
        return ;
    }
    int mid=(l+r)>>1;
    build(l,mid,v*2);
    build(mid+1,r,v*2+1);
    pushup(v);
}
int query(int x,int l, int r)
{
    if (tree[x].l >= l && tree[x].r <= r) {
        return tree[x].val;
    }
    int mid = (tree[x].l + tree[x].r) / 2;
    int ans = 0x3f3f3f3f;
    if (l <= mid) ans = min(ans, query(lson(x),l,r));
    if (r > mid) ans = min(ans, query(rson(x),l,r));
    return ans;
}
int Yougth(int u, int v) {
    int tp1 = top[u], tp2 = top[v];
    int ans = 0x3f3f3f3f;
    while (tp1 != tp2) {
        if (dep[tp1] < dep[tp2]) {
            swap(tp1, tp2);
            swap(u, v);
        }
        ans = min(query(1,id[tp1], id[u]), ans);
        u = fa[tp1];
        tp1 = top[u];
    }
    if (u == v) return ans;
    if (dep[u] > dep[v]) swap(u, v);
    ans = min(query(1,id[son[u]], id[v]), ans);
    return ans;
}
int par[10000+10];
int find(int u)
{
    if(u==par[u]) return u;
    return par[u] = find(par[u]);
}
int main()
{
    int n,m;
    memset(val,0x3f3f3f3f,sizeof(val));
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++)
        scanf("%d%d%d",&ed[i].from,&ed[i].to,&ed[i].w);
    sort(ed,ed+m);
    for(int i=1;i<=n;i++) par[i] = i; int cn = 1;
    for(int i=0;i<m;i++)
    {
        int fx = find(ed[i].from);
        int fy = find(ed[i].to);
        if(fx!=fy){
            par[fx] = fy;
            e[cn].read(ed[i].from,ed[i].to,ed[i].w);
            v[e[cn].x].push_back(e[cn].y);
            v[e[cn].y].push_back(e[cn].x);
            cn++;
        }
    }
    num = 0;
    dfs1(1,0,1);
    dfs2(1,1);
    for (int i = 1; i < cn; i++) {
        if (dep[e[i].x] < dep[e[i].y]) swap(e[i].x, e[i].y);
        val[id[e[i].x]] = e[i].val;
    }
    build(1,num,1);
    int q,a,b; scanf("%d",&q);
    while(q--)
    {
        scanf("%d%d",&a,&b);
        int fa = find(a);
        int fb = find(b);
        if(fa==fb) printf("%d\n",Yougth(a,b));
        else printf("-1\n");
    }
    return 0;
}


评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值