【题解】Uoj#30 Tourist(广义圆方树+树上全家桶)

名字听起来很霸气其实算法很简单....

仙人掌上的普通圆方树是普及题,但是广义圆方树虽然很直观但是有很多地方值得深思

说一下算法的流程:

  • 对于所有点强连通分量(强联通,意味着要找极大的那个),建立一个虚点\(u\),然后把环内所有边断开,紧接着让环内所有点向这个虚点连边。可以看出对于每一个大小为\(s\)的SCC,我们导出了一个点数为\(s+1\)边数为\(n\)的图且联通,所以圆方树是树。
  • 为了方便讨论,对于每个桥加个虚点。虚点维护SCC内所有点的信息

正确性是显然的,因为我如果要从环的某一点a走到另一点b,我一定会在这颗圆方树上经过虚点。虚点可以用不同的方法维护点权信息。

然而不好维护边权,所以我们的做法是直接把边看做一个点建图...

用Tarjan魔改一下就能找了

这一题就是模板题,直接建出来树剖即可。

然后树上的点权修改如果单次修改和度数有关是\(O(n)\)的,一个常见讨论是在父亲处打tag,此时为了方便讨论加的规则就派上了用(少讨论很多东西)

//@winlere
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<set>
#define mid ((l+r)>>1)
#define lef l,mid,pos<<1
#define rgt mid+1,r,pos<<1|1
#define getchar() (__c==__ed?(__ed=__buf+fread(__c=__buf,1,1<<21,stdin),*__c++):*__c++)

using namespace std;  typedef long long ll;   char __buf[1<<21],*__c=__buf,*__ed=__buf;
inline int qr(){
    int ret=0,f=0,c=getchar();
    while(!isdigit(c))f|=c==45,c=getchar();
    while(isdigit(c)) ret=ret*10+c-48,c=getchar();
    return f?-ret:ret;
}

const int maxn=1e5+5;
vector<int> eg[maxn],et[maxn<<1];
int w[maxn<<1];
int dfn[maxn<<1],low[maxn],stk[maxn],top;
int arc[maxn<<1],fi[maxn<<1],d[maxn<<1],r[maxn<<1];
int seg[maxn<<3],siz[maxn<<1],son[maxn<<1];
int n,m,q,cnt;

struct Pri{
    multiset<int> s;
    int Top;
    Pri(){Top=1e9+1;s.clear();}
    inline void insert(int x){s.insert(x);Top=*s.begin();}
    inline void del(int x){s.erase(s.find(x));Top=*s.begin();}
    inline int top(){return Top;}
}p[maxn];

void add(vector<int>*e,int fr,int to){
    e[fr].push_back(to);
    e[to].push_back(fr);
}

void dfs(int now,int last){
    stk[++top]=now; dfn[now]=low[now]=++*dfn;
    for(auto t:eg[now]){
        if(!dfn[t]){
            dfs(t,now);
            if(low[t]>=dfn[now]){
                ++cnt; int temp;
                add(et,now,cnt);
                do temp=stk[top--],add(et,cnt,temp),p[cnt-n].insert(w[temp]);
                while(temp!=t);
            }
            low[now]=min(low[now],low[t]);
        }else low[now]=min(low[now],dfn[t]);
    }
}

void build(int l,int r,int pos){
    if(l==r) return seg[pos]=w[arc[l]],void();
    build(lef); build(rgt);
    seg[pos]=min(seg[pos<<1],seg[pos<<1|1]);
}

void upd(int v,int p,int l,int r,int pos){
    if(p<l||p>r) return;
    if(l==r) return seg[pos]=v,void();
    upd(v,p,lef); upd(v,p,rgt);
    seg[pos]=min(seg[pos<<1],seg[pos<<1|1]);
}

int que(int L,int R,int l,int r,int pos){
    if(L>r||R<l) return 1e9;
    if(L<=l&&r<=R) return seg[pos];
    return min(que(L,R,lef),que(L,R,rgt));
}

void dfs1(int now,int last){
    r[now]=last; d[now]=d[last]+1;
    siz[now]=1;
    for(auto t:et[now])
        if(!siz[t])
            dfs1(t,now),siz[now]+=siz[t],son[now]=siz[son[now]]>siz[t]?son[now]:t;
}

void dfs2(int now,int last){
    fi[now]=last; dfn[now]=++*dfn; arc[*dfn]=now;
    if(son[now]) dfs2(son[now],last);
    for(auto t:et[now])
        if(!dfn[t]) dfs2(t,t);
}

int que(int u,int v){
    int ret=1e9+1;
    while(fi[u]^fi[v]){
        if(d[fi[u]]<d[fi[v]]) swap(u,v);
        ret=min(ret,que(dfn[fi[u]],dfn[u],1,cnt,1));
        u=r[fi[u]];
    }
    if(d[u]<d[v]) swap(u,v);
    ret=min(que(dfn[v],dfn[u],1,cnt,1),ret);
    if(v>n) ret=min(ret,w[r[v]]);
    return ret;
}

int main(){
    cnt=n=qr(); m=qr(); q=qr();
    for(int t=1;t<=n;++t) w[t]=qr();
    for(int t=1,a,b;t<=m;++t)
        a=qr(),b=qr(),add(eg,a,b);
    for(int t=1;t<=n;++t) if(!dfn[t]) dfs(t,0);
    memset(dfn,0,sizeof dfn);
    dfs1(1,0); dfs2(1,1);
    for(int t=n+1;t<=cnt;++t) w[t]=p[t-n].top();
    build(1,cnt,1);
    while(q--){
        char c=getchar();
        while(!isalpha(c)) c=getchar();
        int u=qr(),v=qr();
        if(c=='A') printf("%d\n",que(u,v));
        else{
            upd(v,dfn[u],1,cnt,1);
            if(r[u]){
                p[r[u]-n].del(w[u]);
                p[r[u]-n].insert(v);
                upd(p[r[u]-n].top(),dfn[r[u]],1,cnt,1);
            }
            w[u]=v;
        }
    }
    return 0;
}

12-21 18:40