jzoj5909 【NOIP2018模拟10.16】跑商（圆方树性质+树剖）

Description

一个无向图，求从i到j简单路径（无重点）上的最小点权。
n,m<=1e5

分析

考虑走到某个点双上时，假如不是同点进同点出，就可以获得这个点双上的所有贡献。
于是建出圆方树，在方点开个multiset维护点权最小值。
每次询问就是链上最值。
但是考虑修改，因为一个点可能在多个点双内，要摊一下修改和查询的复杂度。
我们使得每个方点只管儿子，不管父亲。 这样假如lca是方点，只需要额外考虑lca的父亲就行。
tarjan缩点双和缩边双经常搞混，正在探索一致的打法

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <set>
using namespace std;
const int N = 2e5 + 10;
int final[N],nex[N*2],to[N*2],tot,dep[N],top[N],size[N],hvy[N],fa[N];
int n,m,low[N],dfn[N],S[N],stm,a[N];
int gfinal[N],gnex[N*2],gto[N*2],gtot,dcnt;
multiset<int> s[N];

void link(int x,int y) {
//	printf("%d %d\n",x,y);
	to[++tot]=y,nex[tot]=final[x],final[x]=tot;
}
void linkg(int x,int y) {
	gto[++gtot]=y,gnex[gtot]=gfinal[x],gfinal[x]=gtot;
}

void tarjan(int x,int from) {
	dfn[x]=low[x]=++stm;
	S[++S[0]]=x;
	for (int i = gfinal[x]; i; i=gnex[i]){
		if (i==(from^1)) continue;
		int y = gto[i];
		if (dfn[y]==0) {
			tarjan(y,i);
			low[x]=min(low[x],low[y]);
			if (low[y] > dfn[x]) {
				link(x,y),link(y,x);
				S[0]--;
			} else if (low[y]==dfn[x]) {
				dcnt++;
				while (S[S[0]+1] != y) {
					link(S[S[0]],dcnt);
					link(dcnt,S[S[0]]);
					S[0]--;
				}
				link(x,dcnt),link(dcnt,x);
			}
		} else low[x]=min(low[x],dfn[y]);
	}
}

void dfs(int x) {
	size[x]=1;
	for (int i = final[x]; i; i=nex[i]){
		int y = to[i]; if (y != fa[x]) {
			fa[y] = x;
			dep[y] = dep[x] + 1;
			dfs(y);
			size[x]+=size[y];
			if (size[y] > size[hvy[x]])
				hvy[x]=y;
		}
	}
}

void build(int x) {
	dfn[x]=++stm;
	if (!top[x]) top[x]=x;
	if (hvy[x]) top[hvy[x]]=top[x],build(hvy[x]);
	for (int i = final[x]; i; i=nex[i]) {
		int y = to[i];
		if (y!=fa[x] && y!=hvy[x]) {
			build(y);
		}
		if (x > n) {
			s[x].insert(a[y]);
		}
	}
}

int t[N*4];
int lca(int x,int y) {
	while (top[x] != top[y]) {
		if (dep[top[x]] < dep[top[y]]) swap(x,y);
		x = fa[top[x]];
	}
	return dep[x] > dep[y] ? y : x;
}

void query(int x,int l,int r,int tl,int tr,int &ret) {
	if (l>tr || r<tl) return;
	if (tl<=l && r<=tr) {
		ret=min(ret,t[x]);
		return;
	}
	query(x*2,l,l+r>>1,tl,tr,ret);
	query(x*2+1,(l+r>>1)+1,r,tl,tr,ret);
}

char ch[10];
int query(int x,int lca) {
	int ret = 1e9;
	while (top[x] != top[lca]) {
		int mi=1e9;
		query(1,1,dcnt,dfn[top[x]],dfn[x],mi);
		ret=min(ret, mi);
		x=fa[top[x]];
	}
	int mi=1e9;
	query(1,1,dcnt,dfn[lca],dfn[x],mi);
	return min(mi,ret);
}

void change(int x,int l,int r,int tg,int v) {
	if (l==r){
		t[x]=v;
		return;
	}
	if (tg<=(l+r>>1))
		change(x*2,l,l+r>>1,tg,v);
	else
		change(x*2+1,(l+r>>1)+1,r,tg,v);
	t[x]=min(t[x*2],t[x*2+1]);
}

int Q,WW;
int main() {
	freopen("paoshang.in","r",stdin);
	freopen("paoshang.out","w",stdout);
	cin>>n>>m;
	for (int i = 1; i <= n; i++) scanf("%d",&a[i]);
	memset(t,127,sizeof t);
	gtot = 1, dcnt = n;
	for (int i = 1; i <= m; i++) {
		int u,v; scanf("%d %d",&u,&v);
		linkg(u,v),linkg(v,u);
	}
	tarjan(1,0);
	dep[1]=1; dfs(1);
	stm=0;
	build(1);

	for (int i = 1; i <= n; i++)
		change(1,1,dcnt,dfn[i],a[i]);
	for (int i = n+1; i <= dcnt; i++)  {
		change(1,1,dcnt,dfn[i],*s[i].begin());
	}

	cin>>Q;
	for (int i = 1; i <= Q; i++) {
		scanf("%s",ch);
		if (ch[0]=='Q') {
			WW++;
			int u,v; scanf("%d %d",&u,&v);
			int l = lca(u,v);
			int ans = (l > n) ? a[fa[l]] : 1e9;
			int cx = min(query(u,l),query(v,l));
			printf("%d\n",max(0,a[u]-min(ans,cx)));
		} else {
			int x,v; scanf("%d %d",&x,&v);
			change(1,1,dcnt,dfn[x],v);
			if (fa[x] > n) {
				int bl = *s[fa[x]].begin();
				s[fa[x]].erase(s[fa[x]].find(a[x]));
				s[fa[x]].insert(v);
				if (bl != *s[fa[x]].begin())
					change(1,1,dcnt,dfn[fa[x]],*s[fa[x]].begin());
			}
			a[x]=v;
		}
	}
}