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  阅读(167)  评论(0编辑  收藏  举报