LCT 最小生成树
洛谷题目传送门
BZOJ题目传送门
改了好久。。。都跨月了 最后发现x>y的时候没有swap然后排序就炸了。。。
同这道题。
维护边权套路加一个点。
删边的话就把过程倒过来,从最终状态加边就好了。具体实现详见代码(记个id然后sort一下)
注意数据加强版要先在最终状态的那几条边跑MST,完成了就break,不能像那道题一样和加边一起算(别问我怎么知道的,1e6的边不是闹着玩的)。
代码:
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 100005
#define M 1000005
#define F inline
#define V F void
using namespace std;
struct edge{ int x,y,id,z; }ed[M<<1];
struct tree{ int x,mx,fa,to[2],f; }t[M<<1];
int n,m,q,k,tp,num,stk[M<<1],fa[N],ans[N];
F char readc(){
static char buf[100000],*l=buf,*r=buf;
if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
if (l==r) return EOF; return *l++;
}
F int _read(){
int x=0,f=1; char ch=readc();
while (!isdigit(ch)) { if (ch=='-') f=-1; ch=readc(); }
while (isdigit(ch)) x=(x<<3)+(x<<1)+(ch^48),ch=readc();
return x*f;
}
V writec(int x){ if (x>9) writec(x/10); putchar(x%10+48); }
V _write(int x){ writec(x),putchar('\n'); }
F bool cmp1(edge a,edge b){
return a.x<b.x||(a.x==b.x&&a.y<b.y)||(a.x==b.x&&a.y==b.y&&a.z<b.z);
}
F bool cmp2(edge a,edge b){ return a.id<b.id||(a.id==b.id&&a.z<b.z); }
#define addedge(x,y,z,w) (ed[++k]=(edge){x,y,z,w})
int findfa(int x){ return x==fa[x]?x:fa[x]=findfa(fa[x]); }
#define rt(x) (t[t[x].fa].to[0]!=(x)&&t[t[x].fa].to[1]!=(x))
V pshp(int x){
int l=t[x].to[0],r=t[x].to[1]; t[x].x=x;
if (t[t[l].x].mx>t[t[x].x].mx) t[x].x=t[l].x;
if (t[t[r].x].mx>t[t[x].x].mx) t[x].x=t[r].x;
}
V pshd(int x){
if (!t[x].f) return; int &l=t[x].to[0],&r=t[x].to[1];
t[l].f^=1,t[r].f^=1,swap(l,r),t[x].f=0;
}
V rtt(int x){
int y=t[x].fa,z=t[y].fa,l=(t[y].to[0]==x);
if (!rt(y)) t[z].to[t[z].to[0]!=y]=x;
t[x].fa=z,t[y].fa=x,t[t[x].to[l]].fa=y;
t[y].to[l^1]=t[x].to[l],t[x].to[l]=y;
pshp(y),pshp(x);
}
V splay(int x){
for (int i=stk[tp=1]=x;!rt(i);i=t[i].fa) stk[++tp]=t[i].fa;
while (tp) pshd(stk[tp--]);
for (int y=t[x].fa,z=t[y].fa;!rt(x);rtt(x),y=t[x].fa,z=t[y].fa)
if (!rt(y)) rtt((t[z].to[0]==y^t[y].to[0]==x)?x:y);
}
V ccss(int x){ for (int i=0;x;i=x,x=t[x].fa) splay(x),t[x].to[1]=i,pshp(x); }
V mkrt(int x){ ccss(x),splay(x),t[x].f^=1; }
V lnk(int x,int y){ mkrt(x),t[x].fa=y; }
V sprt(int x,int y){ mkrt(x),ccss(y),splay(y); }
V cut(int x,int y){ sprt(x,y); if (t[y].to[0]==x) t[x].fa=t[y].to[0]=0; }
F int srch(int x,int y){ sprt(x,y); return t[y].x; }
F int ef(int x,int y){
int l=1,r=m,mid;
while (l<=r){
int p=ed[mid=l+r>>1].x,q=ed[mid].y;
if (p<x||p==x&&q<y) l=mid+1;
else if (p>x||p==x&&q>y) r=mid-1;
else return mid;
}
}
int main(){
n=_read(),m=_read(),q=_read();
for (int i=1,x,y,z;i<=m;i++){
x=_read(),y=_read(),z=_read();
if (x>y) swap(x,y); addedge(x,y,0,z);
}
sort(ed+1,ed+k+1,cmp1);
for (int i=1,f,x,y;i<=q;i++){
f=_read(),x=_read(),y=_read(); if (x>y) swap(x,y);
if (f==1) addedge(x,y,q-i+1,-1e9);
else ed[ef(x,y)].id=q-i+1;
}
sort(ed+1,ed+k+1,cmp2); for (int i=1;i<=n;i++) fa[i]=i;
for (int i=1,p=0,x,y,fx,fy;i<=k-q;i++)
if ((fx=findfa(x=ed[i].x))!=(fy=findfa(y=ed[i].y))){
p++,fa[fx]=fy,t[i+n].mx=ed[i].z,t[i+n].x=i+k;
lnk(x,i+n),lnk(y,i+n); if (p==n-1) break;
}
for (int i=k-q+1,x,y,e,z;i<=k;i++)
if ((z=ed[i].z)!=-1e9){
if (t[e=srch(x=ed[i].x,y=ed[i].y)].mx<=z) continue;
cut(ed[e-n].x,e),cut(ed[e-n].y,e);
t[i+n].mx=z,t[i+n].x=i+k,lnk(x,i+n),lnk(y,i+n);
}
else ans[++num]=t[srch(ed[i].x,ed[i].y)].mx;
for (int i=num;i;i--) _write(ans[i]); return 0;
}