题目大意

给定\(n\)个点\(m\)条边的无向连通图,无重边

每个点有点权

两个操作:

1.单点点权修改

2.询问从x到y的简单路径中,路径经过点的最小值的最小值时多少

(简单路径指经过每一个点至多一次)

分析

使用点双连通分量转化为树

(点双连通分量表示该连通分量中,去掉任意一个点都不回使其不连通)

这样子

当我们经过一个点双连通分量的至少两个点时,可以拿到这个点双内的最小权值

证明:

1.当只经过一个点时,经过的是一个割点,只经过一个点的情况下,不能进入该点双

2.当经过至少两个点时

①该点双只有两个点,显然成立

②该点双有\(\ge 3\)个点时,设最开始经过a,最后经过的点是b,该点双的最小值为c

当c为a或者b时,成立

否则,一定存在一条从a-c-b的路径

证明:

因为是点双,所以a-c,a-b,c-b的路径一定都是存在的

假设不存在从a-c-b的路径

说明任选一个a-c路径,所有b-c路径都与它有交

若交点有两个,则可以从a出发,走最近的交点到c,然后走另一个交点去b

否则只有一个交点,该点为b-c的必经点,删除后导致点双不连通,与定义矛盾

做法

类似圆方树的方法

点双中间加个特殊点

这篇博客挺好的

修改时不能全改(儿子多)

那就只改父亲

特判lca为特殊的情况,此时经过了该点双至少两个点

而点双中特殊点的父亲那个点并没有统计到

solution

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cctype>
#include <functional>
#include <algorithm>
#include <queue>
using namespace std;
const int M=2e5+7;
const int INF=1e9+7; inline int ri(){
int x=0;bool f=1;char c=getchar();
for(;!isdigit(c);c=getchar()) if(c=='-') f=0;
for(;isdigit(c);c=getchar()) x=x*10+c-48;
return f?x:-x;
} int n,m,q;
int val[M]; struct vec{
int g[M],te;
struct edge{
int y,nxt;
edge(int _y=0,int _nxt=0){y=_y,nxt=_nxt;}
}e[M<<1];
vec(){memset(g,0,sizeof g);te=0;}
inline void push(int x,int y){e[++te]=edge(y,g[x]);g[x]=te;}
inline void push2(int x,int y){push(x,y);push(y,x);}
inline int& operator () (int x){return g[x];}
inline edge& operator [] (int x){return e[x];}
}e; struct Heap{
priority_queue<int,vector<int>,greater<int> >q,ers;
inline void ins(int d){q.push(d);}
inline void del(int d){ers.push(d);}
int top(){
while(!ers.empty()&&q.top()==ers.top()) q.pop(),ers.pop();
return q.top();
}
}; namespace Seg{
int a[524299];
int n; inline void pushup(int x){a[x]=min(a[x<<1],a[x<<1|1]);} void ins(int x,int l,int r,int to,int d){
if(l==r){
a[x]=d;
return;
}
int mid=l+r>>1;
if(to<=mid) ins(x<<1,l,mid,to,d);
else ins(x<<1|1,mid+1,r,to,d);
pushup(x);
} int get(int x,int l,int r,int tl,int tr){
if(tl<=l&&r<=tr) return a[x];
int mid=l+r>>1;
if(tr<=mid) return get(x<<1,l,mid,tl,tr);
if(mid<tl) return get(x<<1|1,mid+1,r,tl,tr);
return min(get(x<<1,l,mid,tl,tr),get(x<<1|1,mid+1,r,tl,tr));
} void ins(int x,int d){ins(1,1,n,x,d);}
int get(int x,int y){return get(1,1,n,x,y);}
int get(int x){return get(1,1,n,x,x);}
} namespace Tr{
vec e;
Heap h[M]; int pre[M],dep[M];
int sz[M],son[M];
int top[M],tdfn;
int tid[M],pid[M]; void dfs1(int x){
int p,y;
sz[x]=1; son[x]=0;
for(p=e(x);p;p=e[p].nxt)
if((y=e[p].y)!=pre[x]){
dep[y]=dep[x]+1;
pre[y]=x;
dfs1(y);
sz[x]+=sz[y];
if(sz[y]>sz[son[x]]) son[x]=y;
}
} void dfs2(int x){
int p,y;
pid[tid[x]=++tdfn]=x;
if(y=son[x]){
top[y]=top[x];
dfs2(y);
}
for(p=e(x);p;p=e[p].nxt)
if((y=e[p].y)!=pre[x]&&y!=son[x]){
top[y]=y;
dfs2(y);
}
} void split(){
pre[1]=0;dep[1]=1;
dfs1(1);
top[1]=1;tdfn=0;
dfs2(1);
} int LCA(int x,int y){
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
x=pre[top[x]];
}
return dep[x]<dep[y]?x:y;
} void mdf(int x,int pr,int nw){
Seg::ins(tid[x],nw);
if(pre[x]==0) return;
if(pr>0) h[pre[x]].del(pr);
h[pre[x]].ins(nw);
Seg::ins(tid[pre[x]],h[pre[x]].top());
} int get(int x,int y){
int lca=LCA(x,y),res=INF;
for(;dep[top[x]]>dep[lca];x=pre[top[x]])
res=min(res,Seg::get(tid[top[x]],tid[x]));
if(dep[x]>dep[lca])
res=min(res,Seg::get(tid[lca]+1,tid[x]));
for(;dep[top[y]]>=dep[lca];y=pre[top[y]])
res=min(res,Seg::get(tid[top[y]],tid[y]));
if(dep[y]>=dep[lca])
res=min(res,Seg::get(tid[lca],tid[y]));
if(lca>n&&pre[lca]) res=min(res,Seg::get(tid[pre[lca]]));
return res;
} void push(int x,int y){e.push2(x,y);}
} namespace G{
vec e; int dfn[M],low[M],tdfn;
int stack[M],Top;
int cnt; void BCC(int x,int fr){
int p,y;
dfn[x]=low[x]=++tdfn;
stack[++Top]=x;
for(p=e(x);p;p=e[p].nxt)
if((y=e[p].y)!=fr){
if(dfn[y]) low[x]=min(low[x],dfn[y]);
else{
BCC(y,x);
low[x]=min(low[x],low[y]);
if(low[y]>=dfn[x]){
cnt++;
Tr::push(x,cnt);
do{Tr::push(cnt,stack[Top]);}while(y!=stack[Top--]);
}
}
}
} void push(int x,int y){e.push2(x,y);}
} int main(){ int i,x,y; char s[9];
n=ri(),m=ri(),q=ri();
for(i=1;i<=n;i++) val[i]=ri(); for(i=1;i<=m;i++) G::push(ri(),ri()); G::cnt=n;
G::BCC(1,0);
Seg::n=G::cnt; Tr::split(); for(i=1;i<=n;i++) Tr::mdf(i,0,val[i]); for(i=1;i<=q;i++){
scanf("%s",s);
if(s[0]=='C'){
x=ri(), y=ri();
Tr::mdf(x,val[x],y);
val[x]=y;
}
else if(s[0]=='A'){
x=ri(),y=ri();
printf("%d\n",Tr::get(x,y));
}
} return 0;
}

cf 487E Tourist的更多相关文章

  1. cf C. Tourist Problem

    http://codeforces.com/contest/340/problem/C #include <cstdio> #include <cstring> #includ ...

  2. CF 732F Tourist Reform——v-SCC+dfs

    题目:http://codeforces.com/contest/732/problem/F 给无向图定向使得从每个点出发能去的点数最小值最大. SCC.点内部dfs定向.点间以siz最大的为起点反向 ...

  3. Solution -「CF 487E」Tourists

    \(\mathcal{Description}\)   Link.   维护一个 \(n\) 个点 \(m\) 条边的简单无向连通图,点有点权.\(q\) 次操作: 修改单点点权. 询问两点所有可能路 ...

  4. CF数据结构练习

    1. CF 438D The Child and Sequence 大意: n元素序列, m个操作: 1,询问区间和. 2,区间对m取模. 3,单点修改 维护最大值, 取模时暴力对所有>m的数取 ...

  5. 732F Tourist Reform

    // CF 732F Tourist Reform // 思路:两遍tarjan // 找强联通分量 #include <bits/stdc++.h> using namespace st ...

  6. 2018 桂林ccpc现场赛 总结

    Day 0 5个小时的火车,坐的昏昏欲睡.桂林站出来没有地铁,而是出租车排成长队依次上车,也算是某种意义上的地铁吧.到了酒店才发现学校那边又给我们换了,又拖着行李找新的酒店,途中路过一家餐馆,所有人都 ...

  7. Note -「圆方树」学习笔记

    目录 圆方树的定义 圆方树的构造 实现 细节 圆方树的运用 「BZOJ 3331」压力 「洛谷 P4320」道路相遇 「APIO 2018」「洛谷 P4630」铁人两项 「CF 487E」Touris ...

  8. C. Tourist Problem 2021.3.29 晚vj拉题 cf 1600 纯数学题

    拉题链接  https://vjudge.net/contest/430219#overview 原题链接  https://codeforces.com/problemset/problem/340 ...

  9. ORA-00494: enqueue [CF] held for too long (more than 900 seconds) by 'inst 1, osid 5166'

    凌晨收到同事电话,反馈应用程序访问Oracle数据库时报错,当时现场现象确认: 1. 应用程序访问不了数据库,使用SQL Developer测试发现访问不了数据库.报ORA-12570 TNS:pac ...

随机推荐

  1. Linux Shell 几个特殊符号命令 & 、&& 、 ||

    & 放在启动参数后面表示设置此进程为后台进程 默认情况下,进程是前台进程,这时就把Shell给占据了,我们无法进行其他操作,对于那些没有交互的进程,很多时候,我们希望将其在后台启动,可以在启动 ...

  2. SpringMVC 多视图解析器 跳转问题

    在SpringMVC的配置文件中加入以下配置: <!--  下面红色的配置必须要在--> <mvc:default-servlet-handler /> <bean id ...

  3. 6- vue django restful framework 打造生鲜超市 -完成商品列表页(下)

    Vue+Django REST framework实战 搭建一个前后端分离的生鲜超市网站 Django rtf 完成 商品列表页下 drf中的request和response drf对于django的 ...

  4. ZendFramework-2.4 源代码 - ViewManager类图

  5. HOJ 13845 Atomic Computer有向无环图的动态规划

    考虑任意一个数字,任何一个都会有奇怪的..性质,就是一个可以保证不重复的方案——直接简单粗暴的最高位加数字..于是,如同上面的那个题:+1.-1.0 但是考虑到65536KB的标准内存限制,会得出一个 ...

  6. Storm: 遇到问题总结

    1.没有ack : kafkaspout id 重复导致每次读最新没有数据. 2.由于storm提供的读取kafka的enternal工具存在bug,导致重复读取数据,致使数据不准确.storm bu ...

  7. 3 View - Request对象

    1.HttpReqeust对象 服务器接收到http协议的请求后,会根据报文创建HttpRequest对象 视图函数的第一个参数是HttpRequest对象 在django.http模块中定义了Htt ...

  8. heap&stack的区别

    1. heap (堆)是一个可动态申请的内存空间,一般所有创建的对象和数组都放在这里.stack (栈)是一个先进后出的数据结构,通常用于保存方法(函数)中的参数,局部变量.stack (栈)的空间小 ...

  9. daily algorithm 判断链表是否有环

    public static boolean isLoopLink(ListNode head) { if (head == null) { return false; } ListNode fast ...

  10. 【LoadRunner】loadrunner常见问题汇总

    LoadRunner常见问题1.LR 脚本为空的解决方法: 1.去掉ie设置中的第三方支持取消掉 2.在系统属性-高级-性能-数据执行保护中,添加loadrunner安装目录中的vugen.exe文件 ...