【BZOJ】3669 [Noi2014]魔法森林 kruskal+LCT

题目传送门

一句话题意:求一条路径,使得 max(ai)+max(bi) 最小。输出这个最小值。

还是ZZK最强了,一眼就秒掉了这道题。

首先我们把所有的边按 ai 排序,从前往后加入边,显然当前的边是最大的 ai ,我们只需要用LCT维护一个从1到n路径上的 max(bi) 就行了。

不过LCT好像没法处理边权,那么我们可以换一种建LCT的方法:假设当前边的两端为 x,y ,当前边的编号为 s ,那么可以建xs+n的边和 s+ny 的边,把这条边的边权赋值在 s+n 号点上即可。

附上AC代码:

#include <cstdio>
#include <cctype>
#include <algorithm>
using namespace std;

const int N=2e5+10;
struct note{
    int x,y,a,b;
    bool operator < (const note lyf) const {return a<lyf.a;}
}s[N];
int n,m,fa[N],f[N],ch[N][2],mk[N],w[N],wz[N],ans;

inline char nc(void){
    static char ch[100010],*p1=ch,*p2=ch;
    return p1==p2&&(p2=(p1=ch)+fread(ch,1,100010,stdin),p1==p2)?EOF:*p1++;
}

inline void read(int &a){
    static char c=nc();int f=1;
    for (;!isdigit(c);c=nc()) if (c=='-') f=-1;
    for (a=0;isdigit(c);a=(a<<3)+(a<<1)+c-'0',c=nc());
    return (void)(a*=f);
}

#define isrt(x) ((x)!=ch[f[x]][0]&&(x)!=ch[f[x]][1])
inline void updata(int x){
    wz[x]=x;
    if (w[wz[ch[x][0]]]>w[wz[x]]) wz[x]=wz[ch[x][0]];
    if (w[wz[ch[x][1]]]>w[wz[x]]) wz[x]=wz[ch[x][1]];
}
inline void rotate(int x){
    int y=f[x],op=(x==ch[y][1]);
    ch[y][op]=ch[x][op^1];
    if (ch[x][op^1]) f[ch[x][op^1]]=y;
    if (f[x]=f[y],!isrt(y)) ch[f[y]][ch[f[y]][1]==y]=x;
    return f[y]=x,ch[x][op^1]=y,updata(y),updata(x);
}
inline void push(int x){if (mk[x]) mk[x]=0,mk[ch[x][0]]^=1,mk[ch[x][1]]^=1,swap(ch[x][0],ch[x][1]);}
inline void splay(int x){
    static int top,sk[N];sk[top=1]=x;
    for (int i=x; !isrt(i); i=f[i]) sk[++top]=f[i];
    while (top) push(sk[top--]);
    for (int fa=f[x]; !isrt(x); rotate(x),fa=f[x])
        if (!isrt(fa)) rotate((x==ch[fa][0])==(fa==ch[f[fa]][0])?fa:x);
}
inline void access(int x){for (int lst=0; x; lst=x,x=f[x]) splay(x),ch[x][1]=lst,updata(x);}
inline void root(int x){access(x),splay(x),mk[x]^=1;}
inline void add(int x,int y){root(x),f[x]=y;}
inline void del(int x,int y){root(x),access(y),splay(y),f[x]=ch[y][0]=0,updata(y);}
inline int query(int x,int y){return root(y),access(x),splay(y),w[wz[y]];}

inline void clear(void){
    for (int i=0; i<=n+m; ++i){
        i>n?w[i]=s[i-n].b:w[i]=-2e9;
        wz[i]=i,ch[i][0]=ch[i][1]=f[i]=0,mk[i]=0;
    }
}

inline int gf(int x){return x==fa[x]?x:fa[x]=gf(fa[x]);}
inline void merge(int i){
    int x=s[i].x,y=s[i].y,fx=gf(x),fy=gf(y);
    if (fx!=fy) return fa[fx]=fy,add(x,n+i),add(n+i,y);
    if (query(x,y)>s[i].b){
        int t=wz[y];
        del(s[t-n].x,t),del(t,s[t-n].y);
        add(x,n+i),add(n+i,y);
    }
}

int main(void){
    read(n),read(m);
    for (int i=1; i<=m; ++i) read(s[i].x),read(s[i].y),read(s[i].a),read(s[i].b);
    sort(s+1,s+1+m),clear(),ans=2e9;
    for (int i=1; i<=n; ++i) fa[i]=i;
    for (int i=1; i<=m; ++i){
        merge(i);
        if (gf(1)!=gf(n)) continue;
        ans=min(ans,s[i].a+query(1,n));
    }
    return printf("%d\n",ans==2e9?-1:ans),0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值