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;
}
}
}