BZOJ2594: [Wc2006]水管局长数据加强版

题目链接

既然有边的删除操作。又询问路径最小值，就可以用LCT。
正着直接做的话，成环直接割掉环上最大边显然不可取。因为后面还有删除操作，当前删去的边可能会造成影响。（：废话吧。。但还是说一下比较好。。刚开始我也小小的纠结了一下。。）
所以将边离线下来，删边反过来就是倒序加边。先将原图中，后面不会删除的边做最小生成树。然后倒序加边，成环删去环上最大边，LCT维护区间最大值的位置。

果然是人丑自带常数大呢。。

【代码】

#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#define N 1100005
#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,Q;
int son[N][2],st[N],fa[N],mx[N],a[N];
bool rev[N];

class Offline_Query{
    public:
        int k,x,y,id,ans;
}q[100005];

class Edge{
    public:
        int x,y,z,id;
        bool flag;
}e[1000005];

bool cmp(Edge a,Edge b){
    if(a.x==b.x) return a.y<b.y;
    return a.x<b.x;
}

bool Cmp(Edge a,Edge b){
    return a.z<b.z;
}

bool CMP(Edge a,Edge b){
    return a.id<b.id;
}

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][0],son[k][1]);
}

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;
}

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];
}

int Find(int x){
    Access(x);Splay(x);
    while(son[x][0]) x=son[x][0];
    return x;
}

int Find_Pos(int x,int y)
{
    int l=1,r=m;
    while(l<=r)
    {
        int mid=l+r>>1;
        if(e[mid].x>x||(e[mid].x==x&&e[mid].y>y)) r=mid-1;
        else if(e[mid].x==x&&e[mid].y==y) return mid;
        else l=mid+1;
    }   
}

void Input_Init()
{
    n=read(),m=read();Q=read();
    for(int i=1;i<=m;i++) 
    {
        e[i].x=read(),e[i].y=read(),e[i].z=read();
        if(e[i].x>e[i].y)swap(e[i].x,e[i].y);
    }
    sort(e+1,e+1+m,Cmp);
    for(int i=1;i<=m;i++)
    {
        e[i].id=i;
        a[i+n]=e[i].z;
        mx[n+i]=n+i;
    }
    sort(e+1,e+1+m,cmp);
    for(int i=1;i<=Q;i++)
    {
        q[i].k=read(),q[i].x=read(),q[i].y=read();
        if(q[i].k==2)
        {
            if(q[i].x>q[i].y) swap(q[i].x,q[i].y);
            int t=Find_Pos(q[i].x,q[i].y);
            e[t].flag=1;q[i].id=e[t].id;
        }
    }
    sort(e+1,e+1+m,CMP);
    int tot=0;
    for(int i=1;i<=m;i++)
    {
        if(e[i].flag) continue;
        static int x,y;
        x=e[i].x,y=e[i].y;
        if(Find(x)!=Find(y))
        {
            Link(x,i+n);
            Link(y,i+n);
            tot++;if(tot==n-1) break;
        }
    }
}

void Solve()
{
    for(int i=Q;i;i--)
    {
        if(q[i].k==1) q[i].ans=a[Query(q[i].x,q[i].y)];
        else
        {
            int x=q[i].x,y=q[i].y,k=q[i].id;
            int t=Query(x,y);
            if(a[t]>e[k].z)
            {
                Cut(t,e[t-n].x);Cut(t,e[t-n].y);
                Link(x,k+n);Link(y,k+n);
            }
        }
    }
    for(int i=1;i<=Q;i++) if(q[i].k==1) printf("%d\n",q[i].ans);
}

int main()
{
    Input_Init();
    Solve(); 
    return 0;
}