学习了一下黄学长的做法:http://hzwer.com/3845.html
边集按a值从小到大排序,依次加边,若构成环则删去环上最大边。
LCT维护的是点,由于是动态,所以不方便像树链剖分一样直接将边转到点上,一个很好用的技巧就是将边建新点,然后就好做了。只有边建成的点有权值,权值为b[i],LCT维护最大值的位置。当1和n连通时更新答案。(还是比较好实现的吧。。打完莫名其妙就胡对了)
(学长写的加了各种奇奇怪怪优化的最短路跑出了1100ms+。。我的LCT直接跑到了5700ms+。。)
【代码】
#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#define N 150005
#define M 2001005
#define INF 1e9
#define mod 51061
using namespace std;
typedef unsigned int ll;
typedef pair<ll,ll> pa;
int read()
{
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
int n,m,top,ans=INF;
int son[N][2],st[N],fa[N],mx[N],a[N];
bool rev[N];
class Edge{
public:
int x,y,A,B;
}e[100005];
bool cmp(Edge a,Edge b){
return a.A<b.A;
}
void Input_Init()
{
n=read();m=read();a[0]=-INF;
for(int i=1;i<=m;i++) e[i].x=read(),e[i].y=read(),e[i].A=read(),e[i].B=read();
sort(e+1,e+1+m,cmp);
}
bool isroot(int x){
return son[fa[x]][0]!=x&&son[fa[x]][1]!=x;
}
void pushup(int k)
{
int l=son[k][0],r=son[k][1];
mx[k]=k;
if(a[mx[k]]<a[mx[l]]) mx[k]=mx[l];
if(a[mx[k]]<a[mx[r]]) mx[k]=mx[r];
}
void pushdown(int k)
{
if(!rev[k]) return;
rev[k]=0,rev[son[k][0]]^=1,rev[son[k][1]]^=1;
swap(son[k][1],son[k][0]);
}
void Rotate(int x)
{
int y=fa[x],z=fa[y],l,r;
l=son[y][1]==x;r=l^1;
if(!isroot(y)) son[z][son[z][1]==y]=x;
fa[x]=z;fa[y]=x;fa[son[x][r]]=y;
son[y][l]=son[x][r];son[x][r]=y;
pushup(y);pushup(x);
}
void Splay(int x)
{
st[++top]=x;
for(int i=x;!isroot(i);i=fa[i]) st[++top]=fa[i];
while(top) pushdown(st[top--]);
while(!isroot(x))
{
int y=fa[x],z=fa[y];
if(!isroot(y)) {
if(son[z][0]==y^son[y][0]==x) Rotate(x);
else Rotate(y);
}
Rotate(x);
}
}
void Access(int x){
for(int t=0;x;t=x,x=fa[x])
Splay(x),son[x][1]=t,pushup(x);
}
void Make_Root(int x){
Access(x);Splay(x);rev[x]^=1;
}
int Find(int x)
{
Access(x);Splay(x);
while(son[x][0]) x=son[x][0];
return x;
}
void Link(int x,int y){
Make_Root(x);fa[x]=y;
}
void Cut(int x,int y){
Make_Root(x);Access(y);Splay(y);
son[y][0]=fa[x]=0;pushup(y);
}
int Query(int x,int y){
Make_Root(x);Access(y);Splay(y);
return mx[y];
}
void Solve()
{
for(int i=1;i<=m;i++)
{
static int x,y,A,B;
x=e[i].x,y=e[i].y,A=e[i].A,B=e[i].B;
if(Find(x)==Find(y))
{
int t=Query(x,y);
if(a[t]>B)
{
Cut(t,e[t-n].x);
Cut(t,e[t-n].y);
}
else {
if(Find(1)==Find(n))
ans=min(ans,A+a[Query(1,n)]);
continue;
}
}
a[n+i]=B;mx[n+i]=n+i;
Link(x,i+n);Link(y,i+n);
if(Find(1)==Find(n)) ans=min(ans,A+a[Query(1,n)]);
}
if(ans==INF) printf("-1\n");
else printf("%d\n",ans);
}
int main()
{
Input_Init();
Solve();
return 0;
}