ACM--LCT

争取刷完LCT

魔法森林

贪心,从小到大枚举构成答案的a
维护的是b的边权,不断修改最小生成树
记录最大边权下标
如果当前生成树最大边权小于新边权b,说明在a联通的情况下删除不会影响最优解

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+10;
#define int long long
int n,m,p[N],stk[N];
int find(int x){
   
    if(p[x]!=x)p[x]=find(p[x]);
    return p[x];
}
struct Edge{
   
    int x,y,a,b;
    bool operator<(const Edge&t)const{
   
        return a<t.a;
    }
}e[N];
int val[N];
struct node{
   
    int s[2],p,v;
    int max,min,sum,rev,lazy;
    int ma;//记录最大权值的下标
    #define ls tr[u].s[0]
    #define rs tr[u].s[1]
    #define lx tr[x].s[0]
    #define rx tr[x].s[1]
}tr[N];

void pushup(int u){
   
    int &L=tr[ls].ma;
    int &R=tr[rs].ma;
    int &U=tr[u].ma;
    U=u;//LCT的修改都是在pushdown和pushup里面的,所以要在这里处理好赋值
    //这里ma表示最大权值下标,U=u表示先赋予自己
    if(val[L]>val[U])U=L;
    if(val[R]>val[U])U=R;
//     tr[u].sum=tr[ls].sum+tr[u].v+tr[rs].sum;
//     if(u>n)tr[u].max=tr[u].min=tr[u].v;
//     else tr[u].max=-1e9,tr[u].min=1e9;
//     if(ls){
   
//         tr[u].max=max(tr[u].max,tr[ls].max);
//         tr[u].min=min(tr[u].min,tr[ls].min);
//     }
//     if(rs){
   
//         tr[u].max=max(tr[u].max,tr[rs].max);
//         tr[u].min=min(tr[u].min,tr[rs].min);
//     }
}
void pushrev(int u){
   
    swap(ls,rs);
    tr[u].rev^=1;
}
void pushlazy(int u){
   
    tr[u].sum=-tr[u].sum;
    tr[u].max=-tr[u].max;
    tr[u].min=-tr[u].min;
    tr[u].v=-tr[u].v;
    swap(tr[u].max,tr[u].min);
    tr[u].lazy^=1;
}
void pushdown(int u){
   
    if(tr[u].rev){
   
        pushrev(ls);
        pushrev(rs);
        tr[u].rev=0;
    }
//     if(tr[u].lazy){
   
//         pushlazy(ls);
//         pushlazy(rs);
//         tr[u].lazy=0;
//     }
}
bool isroot(int u){
   
    return tr[tr[u].p].s[0]!=u&&tr[tr[u].p].s[1]!=u;
}
void rotate(int x){
   
    int y=tr[x].p,z=tr[y].p;
    int k=tr[y].s[1]==x;
    if(!isroot(y))tr[z].s[tr[z].s[1]==y]=x;
    tr[x].p=z;
    tr[y].s[k]=tr[x].s[k^1],tr[tr[x].s[k^1]].p=y;
    tr[x].s[k^1]=y,tr[y].p=x;
    pushup(y),pushup(x);
}
void splay(int x){
   
    int top=0,r=x;
    stk[++top]=r;
    while(!isroot(r))stk[++top]=r=tr[r].p;
    while(top)pushdown(stk[top--]);
    while(!isroot(x)){
   
        int y=tr[x].p,z=tr[y].p;
        if(!isroot(y)){
   
            if((tr[z].s[1]==y)^(tr[y].s[1]==x))rotate(x);
            else rotate(y);
        }
        rotate(x);
    }
}
void access(int x){
   
    int z=x;
    for(int y=0;x;y=x,x=tr[y].p){
   
        splay(x);
        tr[x].s[1]=y,pushup(x);
    }
    splay(z);
}
void makeroot(int x){
   
    access(x);
    pushrev(x);
}
int findroot(int x){
   
    access(x);
    while(lx)pushdown(x),x=lx;
    splay(x);
    return x;
}
void split(int x,int y){
   
    makeroot(x);
    access(y);
}
void link(int x,int y){
   
    makeroot(x);
    if(findroot(y)!=x)tr[x].p=y;
}
void cut(int x,int y){
   
    makeroot(x);
    if(findroot(y)==x&&rx==y&&!tr[y].s[0]){
   
        tr[y].p=rx=0;
        pushup(x);
    }
}
void change(int x,int w){
   
    splay(x);
    tr[x].v=w;
    pushup(x);
}
signed main(){
   
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    cin>>n>>m;
    for(int i=1;i<=m;i++){
   
        int u,v,a,b;
        cin>>u>>v>>a>>b;
        e[i]={
   u,v,a,b};
    }
    for(int i=1;i<=n+m;i++)p[i]=i;
    sort(e+1,e+1+m);
    int ans=1e9;
    for(int i=1;i<=m;i++){
   
        int j=i+n;
        int x=e[i].x;
        int y=e[i].y;
        int a=e[i].a;
        int b=e[i].b;
        if(find(x)!=find(y)){
   
            p[find(x)]=find(y);
            tr[j].ma=j;
            val[j]=b;
            link(x,j),link(j,y);
        }
        else{
   
            split(x,y);
            int id=tr[y].ma;
            
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值