• kruskal重构树


    解决中:

    1. 任意两节点(可以不连通)找到x<->y路径中边权的最小的最大值,反之亦然(也可以用树剖写)
    2. 给定起点,经过的路径边权有某限制下的(如小于等于某值)点权第k小(大),需要主席树。

    对于1:
    看着像二分。。
    对原图边权排序,生成树是直接并查集merge x,y两个节点,重构树的话会新生成一个节点,并把原来的x与y的边断开,用边权当做新节点的点权。

    //ff[]是并查集
    for(int i=0;i<m;++i){
            int fx=find(edge[i].x),fy=find(edge[i].y);
            if(fx^fy){
                ff[fx]=ff[fy]=++id;//id是从n+1开始的
                val[id]=edge[i].v;
                add(id,fx);//id是fx,fy的根,不需要加双向边了
                add(id,fy);
            }
        }
    

    这样做的话最多是2*n-1个节点,原图节点x在新图中还是x,新节点的编号从n+1开始,注意,新节点一定是其子树的根,这也保证了树的深度越小,点权越小。
    所以,树建完后,如果原图不连通,则是森林,加上一个虚的根节点,值为-1。
    对于需要查询的x,y两点,若他们在原图联通,在新图上一定是某真子树的节点。
    因为联通的话一定在同一颗树上,该树的根节点一定是一个新节点,由于根节点的值比子节点小,所以x,y最大路径上的最小边权就是xy的lca,如果排序的时候是从小到大排序的,则lca就是路径中的最大值最小。

    luogu P1967 货车运输
    每条路径都有限重,求x-y能通过的最重车辆,即求路径中权值最小值的最大值

    #include<bits/stdc++.h>
    using namespace std;
    #define gc getchar
    #define TT template<class T>inline
    TT bool read(T &x){
        x=0;register char c=gc();register bool f=0;
        while(c<48||c>57){if(c==EOF)return 0;f^=c=='-',c=gc();}
        while(47<c&&c<58)x=(x<<3)+(x<<1)+(c^48),c=gc();
        if(f)x=-x;return 1;
    }
    TT bool read(T&a,T&b){return read(a)&&read(b);}
    TT bool read(T&a,T&b,T&c){return read(a)&&read(b)&&read(c);}
    typedef long long ll;
    const ll MAXN=1e5+8;
    int head[MAXN<<1],cnt;
    int n,m;
    //原图的边
    struct Edge{int x,y,v;bool operator<(const Edge&o){return v>o.v;}}edge[MAXN*5];
    //新图的边
    struct Node{int y,nt;}node[MAXN<<1];
    //新图加边
    void add(int x,int y){
        node[++cnt].y=y;
        node[cnt].nt=head[x];
        head[x]=cnt;
    }
    //lca
    int deep[MAXN<<1],tot[MAXN<<1],son[MAXN<<1],fa[MAXN<<1];
    void dfs1(int x,int f){
        deep[x]=deep[f]+1,tot[x]=1,fa[x]=f;
        int max_son=-1;
        for(int i=head[x];i;i=node[i].nt){
            int y=node[i].y;
            dfs1(y,x);
            tot[x]+=tot[y];
            if(max_son<tot[y]){
                son[x]=y;
                max_son=tot[y];
            }
        }
    }
    int top[MAXN];
    void dfs2(int x,int tp){
        top[x]=tp;
        if(son[x])dfs2(son[x],tp);
        for(int i=head[x];i;i=node[i].nt){
            int y=node[i].y;
            if(top[y])continue;
            dfs2(y,y);
        }
    }
    //并查集
    int ff[MAXN<<1];
    int find(int x){return ff[x]^x?ff[x]=find(ff[x]):x;}
    //新图权值,编号为1~n的权值为0,因为没改动过。。
    int val[MAXN<<1],id;
    void kruskal(){
        id=n;//新节点从n+1开始
        sort(edge,edge+m);
        for(int i=0;i<=2*n+1;++i)ff[i]=i;
        for(int i=0;i<m;++i){
            int fx=find(edge[i].x),fy=find(edge[i].y);
            if(fx^fy){
                ff[fx]=ff[fy]=++id;
                val[id]=edge[i].v;
                add(id,fx);add(id,fy);
            }
        }
        val[++id]=-1;//原图可能不连通,添加根节点,如果lca是该节点,说明xy不连通。
        set<int>root;//森林的根节点
        for(int i=1;i<=n;++i){
            int f=find(i);
            if(!root.count(f))root.insert(f);
        }
        for(set<int>::iterator it=root.begin();it!=root.end();++it){
            add(id,*it);//根节点到森林加边
        }//求lca
        dfs1(id,0);
        dfs2(id,id);
    }
    int lca(int x,int y){
        while(top[x]^top[y]){
            if(deep[top[x]]>deep[top[y]])x=fa[top[x]];
            else y=fa[top[y]];
        }
        return deep[x]<deep[y]?x:y;
    }
    int main() {
        read(n,m);
        for(int i=0;i<m;++i)read(edge[i].x,edge[i].y,edge[i].v);
        kruskal();
        int q,x,y;
        read(q);
        while(q--){
            read(x,y);
            printf("%d
    ",val[lca(x,y)]);
        }
        return 0;
    }
    
  • 相关阅读:
    Maven项目避坑点 language Lever 自动变成5 配置文件无法自动生成
    File类使用
    枚举的使用
    window处理服务端口被占用的情况
    研修教学创新第一讲:什么是好的教学创新
    如何在大赛中做好教学创新与设计(第一讲)
    国家社科基金
    前端打印生成pdf
    print到文本文件
    windows批处理同步系统时间
  • 原文地址:https://www.cnblogs.com/foursmonth/p/14161763.html
Copyright © 2020-2023  润新知