bzoj3669/lg2387:[NOI2014]魔法森林(LCT)

bzoj3669/lg2387:[NOI2014]魔法森林
题意简单来说,就是找到一条1到n的路径,使得这条路径上最大的A妖怪数+B妖怪数最小

显然,这是一道维护边权的题,那我们就把边权拆成点权
这道题要维护两种怪物数,因此我们可以按A从小到大排序,然后去比较B的数
对于每一条边 x->y,若之前不连通,连上即可
若联通,那就说明,如果这条边再加进去,x->y部分就会出现环,因此我们要选择一条边把他删去。找到x->y中最大的那条边,如果新边B更大,那就不加进去,否则就把原边删去假加入新边

#include <cstdio>
#include <algorithm>
using namespace std;
#define M 150010
#define inf 0x7fffffff 
struct node{int x,y,w1,w2;}a[M];
int n,m,ans,FA[M],ch[M][2],v[M],fa[M],rev[M],mx[M];
inline char gc(){
    static char now[1<<16],*S,*T;
    if(S==T){T=(S=now)+fread(now,1,1<<16,stdin);if(S==T) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0' || ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while('0'<=ch && ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
inline bool cmp(node x,node y){return x.w1<y.w1;}
inline int find(int x){
    if(FA[x]!=x) FA[x]=find(FA[x]);
    return FA[x];
}
inline bool isroot(int x){return x!=ch[fa[x]][0] && x!=ch[fa[x]][1];}
inline void pushdown(int x){
    if(rev[x]){
        swap(ch[x][0],ch[x][1]);
        rev[ch[x][0]]^=1;rev[ch[x][1]]^=1;
        rev[x]^=1;
    }
}
inline void push(int x){
    if(!isroot(x)) push(fa[x]);
    pushdown(x);
}
inline void update(int x){
    mx[x]=x;int l=ch[x][0],r=ch[x][1];
    if(v[mx[l]]>v[mx[x]]) mx[x]=mx[l];
    if(v[mx[r]]>v[mx[x]]) mx[x]=mx[r];
}
inline void rotate(int x){
    int y=fa[x],z=fa[y],t=ch[y][0]==x;
    if(!isroot(y)) ch[z][ch[z][1]==y]=x;
    fa[y]=x;fa[x]=z;fa[ch[x][t]]=y;
    ch[y][t^1]=ch[x][t];ch[x][t]=y;
    update(y);update(x);
}
inline void splay(int x){
    push(x);
    while(!isroot(x)){
        int y=fa[x];
        if(isroot(y)){rotate(x);return;}
        if(ch[y][0]==x^ch[fa[y]][0]==y) rotate(x);
        else rotate(y);rotate(x);
    }
}
inline void access(int x){
    int y=0;
    while(x){
        splay(x);ch[x][1]=y;update(x);
        y=x;x=fa[x];
    }
}
inline void makeroot(int x){
    access(x);splay(x);rev[x]^=1;
}
inline void link(int x,int y){
    makeroot(x);fa[x]=y;
} 
inline int query(int x,int y){
    makeroot(x);access(y);splay(y);return mx[y];
}
inline void cut(int x,int y){
    makeroot(x);access(y);splay(y);fa[x]=ch[y][0]=0;update(y);update(x);
}
int main(){
    n=read();m=read();
    for(int i=1;i<=n;i++) FA[i]=i;
    for(int i=1;i<=m;i++) a[i].x=read(),a[i].y=read(),a[i].w1=read(),a[i].w2=read();
    sort(a+1,a+m+1,cmp);ans=inf;
    for(int i=1;i<=m;i++){
        int x=find(a[i].x),y=find(a[i].y);v[i+n]=a[i].w2; 
        if(x!=y) FA[x]=y,link(a[i].x,i+n),link(i+n,a[i].y);
        else{
            int pos=query(a[i].x,a[i].y);
            if(v[pos]>a[i].w2){
                cut(a[pos-n].x,pos);cut(pos,a[pos-n].y);
                link(a[i].x,i+n),link(i+n,a[i].y);
            }
        }if(find(1)==find(n)) ans=min(ans,a[i].w1+v[query(1,n)]);
    }
    if(ans!=inf) printf("%d\n",ans);
    else printf("-1");
    return 0;
} 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值