1. 程式人生 > >BZOJ2594 [Wc2006]水管局長資料加強版【LCT】

BZOJ2594 [Wc2006]水管局長資料加強版【LCT】

Description

SC省MY市有著龐大的地下水管網路,嘟嘟是MY市的水管局長(就是管水管的啦),嘟嘟作為水管局長的工作就是:每天供水公司可能要將一定量的水從x處送往y處,嘟嘟需要為供水公司找到一條從A至B的水管的路徑,接著通過資訊化的控制中心通知路徑上的水管進入準備送水狀態,等到路徑上每一條水管都準備好了,供水公司就可以開始送水了。嘟嘟一次只能處理一項送水任務,等到當前的送水任務完成了,才能處理下一項。
在處理每項送水任務之前,路徑上的水管都要進行一系列的準備操作,如清洗、消毒等等。嘟嘟在控制中心一聲令下,這些水管的準備操作同時開始,但由於各條管道的長度、內徑不同,進行準備操作需要的時間可能不同。供水公司總是希望嘟嘟能找到這樣一條送水路徑,路徑上的所有管道全都準備就緒所需要的時間儘量短。嘟嘟希望你能幫助他完成這樣的一個選擇路徑的系統,以滿足供水公司的要求。另外,由於MY市的水管年代久遠,一些水管會不時出現故障導致不能使用,你的程式必須考慮到這一點。
不妨將MY市的水管網路看作一幅簡單無向圖(即沒有自環或重邊):水管是圖中的邊,水管的連線處為圖中的結點。

Input

輸入檔案第一行為3個整數:N, M, Q分別表示管道連線處(結點)的數目、目前水管(無向邊)的數目,以及你的程式需要處理的任務數目(包括尋找一條滿足要求的路徑和接受某條水管壞掉的事實)。
以下M行,每行3個整數x, y和t,描述一條對應的水管。x和y表示水管兩端結點的編號,t表示準備送水所需要的時間。我們不妨為結點從1至N編號,這樣所有的x和y都在範圍[1, N]內。
以下Q行,每行描述一項任務。其中第一個整數為k:若k=1則後跟兩個整數A和B,表示你需要為供水公司尋找一條滿足要求的從A到B的水管路徑;若k=2,則後跟兩個整數x和y,表示直接連線x和y的水管宣佈報廢(保證合法,即在此之前直接連線x和y尚未報廢的水管一定存在)。

Output

按順序對應輸入檔案中每一項k=1的任務,你需要輸出一個數字和一個回車/換行符。該數字表示:你尋找到的水管路徑中所有管道全都完成準備工作所需要的時間(當然要求最短)。

Sample Input

4 4 3
1 2 2
2 3 3
3 4 2
1 4 2
1 1 4
2 1 4
1 1 4

Sample Output

2
3

【原題資料範圍】

N ≤ 1000
M ≤ 100000
Q ≤ 100000

測試資料中宣佈報廢的水管不超過5000條;且任何時候我們考慮的水管網路都是連通的,即從任一結點A必有至少一條水管路徑通往任一結點B。

【加強版資料範圍】

N ≤ 100000
M ≤ 1000000
Q ≤ 100000
任何時候我們考慮的水管網路都是連通的,即從任一結點A必有至少一條水管路徑通往任一結點B。

題目分析:

由於這題對LCT來說刪邊十分麻煩
所以考慮離線加邊操作

先離線讀入所有詢問
構造好最終情況下的最小生成樹
然後開始倒敘處理詢問並存入棧中

對於操作1
直接查詢u到v路徑上的最大邊權並存入棧內

對於操作2
找到u到v路徑上的最大邊權
若該邊權大於待加入邊的邊權,就替換掉

最後直接輸出棧內答案

其實本體最麻煩的還是邊的編號處理
具體再程式碼中解釋

#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

int n,m,r;
struct node{int u,v,dis,num,flip;}E[1000010];
struct node2{int k,u,v,num;}q[500010];
int fa[3000010],ch[3000010][2];
int val[3000010],maxn[3000010],lzy[3000010];
int st[3000010];
int ff[3000010];
int ans[500010],cnt;

void update(int x)
{
    maxn[x]=x;
    maxn[x]=val[maxn[ch[x][0]]]>val[maxn[x]] ?maxn[ch[x][0]]:maxn[x];
    maxn[x]=val[maxn[ch[x][1]]]>val[maxn[x]] ?maxn[ch[x][1]]:maxn[x];
}

void push(int x)
{
    if(!lzy[x]) return;
    swap(ch[x][0],ch[x][1]);
    lzy[ch[x][0]]^=1; lzy[ch[x][1]]^=1;
    lzy[x]=0;
}

int isrt(int x)
{
    return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x; 
}

void rotate(int x)
{
    int y=fa[x],z=fa[y];
    int d=(ch[y][0]==x);
    if(!isrt(y))
    {
        if(ch[z][0]==y) ch[z][0]=x;
        else ch[z][1]=x;
    }
    fa[y]=x; fa[ch[x][d]]=y; fa[x]=z;
    ch[y][d^1]=ch[x][d]; ch[x][d]=y;
    update(y); update(x);
}

void splay(int x)
{
    int top=0; st[++top]=x;
    for(int i=x;!isrt(i);i=fa[i])
    st[++top]=fa[i];
    while(top) push(st[top--]);

    while(!isrt(x))
    {
        int y=fa[x],z=fa[y];
        if(!isrt(y))
        {
            if((ch[z][0]==y)^(ch[y][0]==x)) rotate(x);
            else rotate(y);
        }
        rotate(x);
    }
}

void access(int x)
{
    int t=0;
    while(x)
    {
        splay(x);
        ch[x][1]=t;
        update(x);
        t=x; x=fa[x];
    }
}

void mkrt(int x)
{
    access(x); splay(x);
    lzy[x]^=1;
}


void match(int x,int y)
{
    mkrt(x); fa[x]=y; 
}

void cut(int x,int y)
{
    mkrt(x);
    access(y); splay(y);
    fa[x]=ch[y][0]=0;
    update(y);
}

int get(int x,int y)
{
    mkrt(x);
    access(y); splay(y);
    return maxn[y];
}

bool cmp_dis(node a,node b){return a.dis<b.dis;}
bool cmp(node a,node b){return a.num<b.num;}
bool cmp_num(node a,node b)
{
    if(a.u==b.u) return a.v<b.v;
    return a.u<b.u;
}

int find_edge(int u,int v)//二分查詢邊的編號
{
    int ll=1,rr=m;
    while(ll<=rr)
    {
        int mid=ll+rr>>1;
        if(E[mid].u==u&&E[mid].v==v) return mid;
        if(E[mid].u<u||(E[mid].u==u&&E[mid].v<v)) ll=mid+1;
        else rr=mid-1;
    }
}

int find(int x)
{
    if(x==ff[x]) return x;
    else return ff[x]=find(ff[x]);
}

void kruskal()
{
    for(int i=1;i<=n;++i) ff[i]=i;
    int tot=0;
    for(int i=1;i<=m;++i)
    {
        if(E[i].flip) continue;
        int fu=find(E[i].u),fv=find(E[i].v);
        if(fu==fv) continue;
        tot++; ff[fu]=fv;
        match(E[i].u,E[i].num+n); match(E[i].v,E[i].num+n);
        if(tot==n-1) return;
    }
}

void solve()
{
    for(int i=r;i>=1;--i)
    {
        if(q[i].k==1){ ans[++cnt]=val[get(q[i].u,q[i].v)]; continue;}
        int num=get(q[i].u,q[i].v);//找到最大的邊權
        if(val[num]>val[q[i].num+n])
        {
            cut(E[num-n].u,num); cut(E[num-n].v,num);
            match(q[i].u,q[i].num+n); match(q[i].v,q[i].num+n);
        }
    }
    while(cnt) printf("%d\n",ans[cnt--]);
}

int main()
{
    n=read();m=read();r=read();
    for(int i=1;i<=m;++i)
    {
        E[i].u=read();E[i].v=read();E[i].dis=read();
        if(E[i].u>E[i].v) swap(E[i].u,E[i].v);//編號小的在前,方便二分查詢
    }

    sort(E+1,E+1+m,cmp_dis);//先按權值排序
    for(int i=1;i<=m;++i)
    {
        E[i].num=i;//設定編號
        val[i+n]=E[i].dis;
        maxn[i+n]=i+n;
    }

    sort(E+1,E+1+m,cmp_num);//按節點編號排序,方便二分查詢邊的編號
    for(int i=1;i<=r;++i)
    {
        q[i].k=read();q[i].u=read();q[i].v=read();//操作1直接儲存
        if(q[i].k==2)
        {
            if(q[i].u>q[i].v) swap(q[i].u,q[i].v);
            int k=find_edge(q[i].u,q[i].v);//找到對應邊的編號
            q[i].num=E[k].num;//記錄對應操作編號 
            E[k].flip=1;//記錄這條邊不能連線
        }
    }
    sort(E+1,E+1+m,cmp);//按編號排序,防止sort的不穩定性(相當於回到第一次排序的狀態)
    kruskal();

    solve();
    return 0;
}