CF487E Tourists

传送门

很显然可以用圆方树转成一棵树,然后就可以愉快的上树剖了
接下来考虑方点,当路径经过方点的时候,这个方点的贡献应该是它所表示的点双内的所有点的最小值
然后我们可以考虑将每个方点用一个multiset存下当前点双内的所有点权
但是发现修改时需要枚举包含当前圆点相邻的所有方点,当该图是菊花图时,就会被卡成\(O(n^2)\)
然后只好去看题解,居然发现志同道合的大佬,然后就借鉴了一下他的处理方法
就是我们可以对于每个方点的multiset只存下它的子节点的权值
那么每次修改就只需要修改当前圆点的父节点就好了
但是这样做的弊端就是:当查询的两个点的lca是一个方点的时候,lca父节点的权值也要考虑进来(显然,lca代表的点双包括它的父节点)
代码:

#include<cstdio>
#include<iostream>
#include<set>
#include<algorithm>
using namespace std;
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
#define rg register
const int maxn=2e5+10;char p[3];
int dep[maxn],f[maxn],top[maxn],n,m,k,w[maxn],tot,size[maxn];
int dfn[maxn],low[maxn],tmp,st[maxn],tp,id[maxn],num,nid[maxn];
multiset<int>d[maxn];
struct segment_tree{int l,r,mn;}s[maxn*4];
multiset<int>::iterator it;
struct oo
{
    int cnt,pre[maxn*2],nxt[maxn*2],h[maxn];
    void add(int x,int y)
    {
        pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt;
        pre[++cnt]=x,nxt[cnt]=h[y],h[y]=cnt;
    }
}a,b;
void tarjan(int x,int fa)
{
    dfn[x]=low[x]=++tmp;st[++tp]=x;
    for(rg int i=a.h[x];i;i=a.nxt[i])
        if(a.pre[i]!=fa)
        {
            if(!dfn[a.pre[i]])
            {
                tarjan(a.pre[i],x);
                low[x]=min(low[x],low[a.pre[i]]);
                if(dfn[x]<=low[a.pre[i]])
                {
                    tot++;int z;
                    do{z=st[tp--],b.add(z,tot);}while(z!=a.pre[i]);
                    b.add(x,tot);
                }
            }
            else low[x]=min(low[x],dfn[a.pre[i]]);
        }
}
void dfs(int x,int fa)
{
    size[x]=1;
    for(rg int i=b.h[x];i;i=b.nxt[i])
        if(b.pre[i]!=fa)
        {
            f[b.pre[i]]=x,dep[b.pre[i]]=dep[x]+1;
            if(x>n)d[x].insert(w[b.pre[i]]);
            dfs(b.pre[i],x),size[x]+=size[b.pre[i]];
        }
}
void dfs2(int x,int f)
{
    top[x]=f,id[x]=++num;nid[num]=x;int k=0;
    for(rg int i=b.h[x];i;i=b.nxt[i])
        if(dep[b.pre[i]]>dep[x]&&size[b.pre[i]]>size[k])k=b.pre[i];
    if(!k)return ;dfs2(k,f);
    for(rg int i=b.h[x];i;i=b.nxt[i])
        if(dep[b.pre[i]]>dep[x]&&k!=b.pre[i])dfs2(b.pre[i],b.pre[i]);
}
void update(int x){s[x].mn=min(s[x<<1].mn,s[x<<1|1].mn);}
void build(int x,int l,int r)
{
    s[x].l=l,s[x].r=r;int mid=(l+r)>>1;
    if(l==r){s[x].mn=nid[l]<=n?w[nid[l]]:*d[nid[l]].begin();return ;}
    build(x<<1,l,mid),build(x<<1|1,mid+1,r);
    update(x);
}
void change(int x,int a,int b,bool c)
{
    if(s[x].l==s[x].r)
    {
        if(c)s[x].mn=*d[nid[a]].begin();
        else s[x].mn=b;
        return ;
    }
    int mid=(s[x].l+s[x].r)>>1;
    if(a<=mid)change(x<<1,a,b,c);
    else change(x<<1|1,a,b,c);
    update(x);
}
int get(int x,int a,int b)
{
    if(a<=s[x].l&&b>=s[x].r)return s[x].mn;
    int mid=(s[x].l+s[x].r)>>1,ans=1e9;
    if(a<=mid)ans=min(ans,get(x<<1,a,b));
    if(b>mid)ans=min(ans,get(x<<1|1,a,b));
    return ans;
}
int qsum(int x,int y)
{
    int ans=1e9;
    while(top[x]!=top[y])
    {
        if(dep[top[x]]<dep[top[y]])swap(x,y);
        ans=min(ans,get(1,id[top[x]],id[x]));
        x=f[top[x]];
    }
    if(id[x]>id[y])swap(x,y);
//  printf("%d %d\n",x,y);
    ans=min(ans,get(1,id[x],id[y]));
    if(x>n)ans=min(w[f[x]],ans);
    return ans;
}
int main()
{
    read(n),read(m),read(k),tot=n;
    for(rg int i=1;i<=n;i++)read(w[i]);
    for(rg int i=1,x,y;i<=m;i++)read(x),read(y),a.add(x,y);
    tarjan(1,0),dfs(1,0),dfs2(1,1),build(1,1,num);
    for(rg int i=1,x,y;i<=k;i++)
    {
        scanf("%s",p+1),read(x),read(y);
        if(p[1]=='C')
        {
            if(f[x])d[f[x]].erase(w[x]),d[f[x]].insert(y),change(1,id[f[x]],y,1);
            w[x]=y,change(1,id[x],y,0);
        }
        else printf("%d\n",qsum(x,y));
    }
}
posted @ 2019-03-28 19:26 蒟蒻--lichenxi 阅读(...) 评论(...) 编辑 收藏