CF487E Tourists + 圆方树学习笔记(圆方树+树剖+线段树+multiset)
QWQ果然我已经什么都学不会的人了。
这个题目要求的是图上所有路径的点权和!QWQ(我只会树上啊!)
这个如果是好啊
这时候就需要
圆方树!
首先在介绍圆方树之前,我们先来一点简单的前置知识
首先,我们需要知道什么是
点双联通分量
若一个无向图中的去掉任意一个节点都不会改变此图的连通性,即不存在割点,则称作点双连通图。那么一个极大的点双联通子图,就是一个双联通分量了
那么求这个方法,和普通求割点的\(tarjan\)类似
用一个栈维护所有的点
对于搜索到一个割点,然后把他的栈内部的点依次弹栈,直到这条边的\(to\)被弹出来为止(这里不能直接把当前弹出去的原因是一个点可能属于好几个点双联通,所以如果此时把这个点弹出去,是错误的
代码就直接搬了圆方树里面的,请手动忽略\(addedge\)
void tarjan(int x,int faa)
{
dfn[x]=low[x]=++tot;
st[++ttop]=x;
for (int i=point1[x];i;i=nxt1[i])
{
int p = to1[i];
if (p==faa) continue;
if (!dfn[p])
{
tarjan(p,x);
low[x]=min(low[x],low[p]);
if (low[p]>=dfn[x]) //如果这个点是割点,就把当前点双的所有点出栈,建成圆方树
{
++num;
addedge(num,x);
addedge(x,num);
do
{
addedge(st[ttop],num);
addedge(num,st[ttop]);
ttop--;
}while (st[ttop+1]!=p);
}
}
else
low[x]=min(low[x],dfn[p]);
}
}
QWQ那么圆方树是怎么一个东西呢。
实际上圆方树就是把图中所有的点双连通分量找出来,然后建一个新的方点
和所有的点双里面的点连边,然后把原图对应点双里面的边都删除(相当于建了一个新图)
然后新建的点是方点,原图的点是原点
这里直接引用了猫琨的图qwq
那么经过这么一番处理,我们会发现(详细证明请看WC的ppt)
这个由圆点和方点组成的新图是一颗树
叫做
圆方树
那么圆方树有什么性质呢?
性质
1任意两个圆点(或方点)不会相邻。
2.如果两个方点树上距离是那么一定有一个圆点连接他们,而且这个圆点是两个方点对应点双的公共点。
3,对于两个圆点 \(a,b\)之间的树上路径,
(1)路径上的所有圆点都是割点!
(2)\(a,b\)所有简单路径的并 是路径上所有方点对应的点双(包括\(a=b\)的情况)
QWQ然后我们就能够将大部分图上问题直接转化成树上问题了
那么对于这个题。
首先,询问路径的上的简单路径的点权\(min\),我们可以将方点的权值设为所连圆点的权值\(min\)(这里要开一个\(set)\),然后圆点的点权弄成原来的权值,对于每次询问,我们直接回答链\(min\)就行,这个可以用线段树+树剖来实现。
那么修改呢
QWQ
由于我们考虑到,一个圆点的点权修改,最多会影响到\(O(n)\)个级别的方点(一个点可以属于好几个点双)。
QWQ
那么该如何是好呢?
我们考虑对于每个方点,维护除去他的根(建出来树先dfs一下)的所有圆点的权值的\(min\)(由于涉及到删除什么的,所以要用\(multiset\))
那么这样每次修改,我们只需要修改对应圆点的\(fa\)的\(set\)就ok了
有一个要注意的地方就是!
如果我们对于一次询问的\(lca\),假设他的方点的话,我们还要考虑他的\(fa\)的\(val\)大小,因为我们只询问路径上的方点,不会计算这个点的。
总的来说,一般圆方树的题都要很好的运用圆方树的性质的
(虽然我根本想不到圆方树)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
#include<set>
#define mk makr_pair
#define ll long long
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
const int maxn = 2e5+5;
const int maxm = 2*maxn;
int point[maxn],nxt[maxm],to[maxm];
int point1[maxn],nxt1[maxm],to1[maxm];
int cnt1,cnt;
int n,m,k,q;
int deep[maxn];
int top[maxn],size[maxn],son[maxn];
int fa[maxn];
int f[4*maxn];
int add[4*maxn];
multiset<int> s[maxn];
int val[maxn];
int tot;
int low[maxn],dfn[maxn];
int newnum[maxn],newval[maxn];
int st[maxn],ttop;
int num,tmp;
void addedge(int x,int y)
{
// cout<<x<<" "<<y<<" *"<<endl;
nxt[++cnt]=point[x];
to[cnt]=y;
point[x]=cnt;
}
void addedge1(int x,int y)
{
nxt1[++cnt1]=point1[x];
to1[cnt1]=y;
point1[x]=cnt1;
}
void tarjan(int x,int faa)
{
dfn[x]=low[x]=++tot;
st[++ttop]=x;
for (int i=point1[x];i;i=nxt1[i])
{
int p = to1[i];
if (p==faa) continue;
if (!dfn[p])
{
tarjan(p,x);
low[x]=min(low[x],low[p]);
if (low[p]>=dfn[x]) //如果这个点是割点,就把当前点双的所有点出栈,建成圆方树
{
++num;
addedge(num,x);
addedge(x,num);
do
{
addedge(st[ttop],num);
addedge(num,st[ttop]);
ttop--;
}while (st[ttop+1]!=p);
}
}
else
low[x]=min(low[x],dfn[p]);
}
}
void dfs1(int x,int faa,int dep)
{
deep[x]=dep;
size[x]=1;
int maxson=-1;
for (int i=point[x];i;i=nxt[i])
{
int p =to[i];
if (p==faa) continue;
fa[p]=x;
dfs1(p,x,dep+1);
if (size[p]>maxson)
{
maxson=size[p];
son[x]=p;
}
size[x]+=size[p];
}
}
void dfs2(int x,int chain)
{
newnum[x]=++tmp;
newval[tmp]=val[x];
top[x]=chain;
if (!son[x]) return;
dfs2(son[x],chain);
for (int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (!newnum[p])
{
dfs2(p,p);
}
}
}
void up(int root)
{
f[root]=min(f[2*root],f[2*root+1]);
}
void build(int root,int l,int r)
{
if (l==r)
{
f[root]=newval[l];
return;
}
int mid =l + r >> 1;
build(2*root,l,mid);
build(2*root+1,mid+1,r);
up(root);
}
void update(int root,int l,int r,int x,int p)
{
if (l==r)
{
f[root]=p;
return;
}
int mid = l+r >> 1;
if (x<=mid) update(2*root,l,mid,x,p);
if (x>mid) update(2*root+1,mid+1,r,x,p);
up(root);
}
int query(int root,int l,int r,int x,int y)
{
if (x<=l && r<=y) return f[root];
int mid = l+r >> 1;
int ans=2e9;
if (x<=mid) ans=min(ans,query(2*root,l,mid,x,y));
if (y>mid) ans=min(ans,query(2*root+1,mid+1,r,x,y));
return ans;
}
int treesum(int x,int y)
{
int ans=2e9;
int xx=x,yy=y;
while (top[x]!=top[y])
{
if (deep[top[x]]<deep[top[y]]) swap(x,y);
//cout<<1<<endl;
ans=min(ans,query(1,1,num,newnum[top[x]],newnum[x]));
x=fa[top[x]];
//cout<<"***"<<endl;
}
// cout<<1<<endl;
//cout<<x<<" "<<y<<endl;
if (deep[x]>deep[y]) swap(x,y);
ans=min(ans,query(1,1,num,newnum[x],newnum[y]));
int l = x;
if (l>n) ans=min(ans,val[fa[l]]);
return ans;
}
int main()
{
n=read(),m=read(),q=read();
tmp=0;
num=n;
for (int i=1;i<=n;i++) val[i]=read();
for (int i=1;i<=m;i++)
{
int x=read(),y=read();
addedge1(x,y);
addedge1(y,x);
}
for (int i=1;i<=n;i++)
{
if(!dfn[i]) tarjan(i,0);
}
dfs1(1,0,1);
for (int i=n+1;i<=num;i++)
{
for (int j=point[i];j;j=nxt[j])
{
int p = to[j];
if (p==fa[i]) continue;
s[i].insert(val[p]);
}
if (s[i].size()==0) val[i]=2e9;
else val[i]=*(s[i].begin());
}
dfs2(1,1);
build(1,1,num);
for (int i=1;i<=q;i++)
{
char ss[10];
scanf("%s",ss+1);
if (ss[1]=='C')
{
int x=read(),y=read();
update(1,1,num,newnum[x],y);
if (fa[x]==0)
{
val[x]=y;
continue;
}
int now = *(s[fa[x]].begin());
s[fa[x]].erase(s[fa[x]].find(val[x]));
s[fa[x]].insert(y);
int tmp1 = (*(s[fa[x]].begin()));
if (tmp1!=now)
{
update(1,1,num,newnum[fa[x]],tmp1);
val[fa[x]]=tmp1;
}
val[x]=y;
}
else
{
int x=read(),y=read();
cout<<treesum(x,y)<<"\n";
}
}
return 0;
}
CF487E Tourists + 圆方树学习笔记(圆方树+树剖+线段树+multiset)的更多相关文章
- [学习笔记]Segment Tree Beats!九老师线段树
对于这样一类问题: 区间取min,区间求和. N<=100000 要求O(nlogn)级别的算法 直观体会一下,区间取min,还要维护区间和 增加的长度很不好求.... 然鹅, 从前有一个来自杭 ...
- 仙人掌&圆方树学习笔记
仙人掌&圆方树学习笔记 1.仙人掌 圆方树用来干啥? --处理仙人掌的问题. 仙人掌是啥? (图片来自于\(BZOJ1023\)) --也就是任意一条边只会出现在一个环里面. 当然,如果你的图 ...
- zkw线段树学习笔记
zkw线段树学习笔记 今天模拟赛线段树被卡常了,由于我自带常数 \(buff\),所以学了下zkw线段树. 平常的线段树无论是修改还是查询,都是从根开始递归找到区间的,而zkw线段树直接从叶子结点开始 ...
- 线段树学习笔记(基础&进阶)(一) | P3372 【模板】线段树 1 题解
什么是线段树 线段树是一棵二叉树,每个结点存储需维护的信息,一般用于处理区间最值.区间和等问题. 线段树的用处 对编号连续的一些点进行修改或者统计操作,修改和统计的复杂度都是 O(log n). 基础 ...
- 主席树学习笔记(静态区间第k大)
题目背景 这是个非常经典的主席树入门题——静态区间第K小 数据已经过加强,请使用主席树.同时请注意常数优化 题目描述 如题,给定N个整数构成的序列,将对于指定的闭区间查询其区间内的第K小值. 输入输出 ...
- BZOJ_4551_[Tjoi2016&Heoi2016]树_树剖+线段树
BZOJ_4551_[Tjoi2016&Heoi2016]树_树剖+线段树 Description 在2016年,佳媛姐姐刚刚学习了树,非常开心.现在他想解决这样一个问题:给定一颗有根树(根为 ...
- BZOJ_2238_Mst_树剖+线段树
BZOJ_2238_Mst_树剖+线段树 Description 给出一个N个点M条边的无向带权图,以及Q个询问,每次询问在图中删掉一条边后图的最小生成树.(各询问间独立,每次询问不对之后的询问产生影 ...
- BZOJ_2157_旅游_树剖+线段树
BZOJ_2157_旅游_树剖+线段树 Description Ray 乐忠于旅游,这次他来到了T 城.T 城是一个水上城市,一共有 N 个景点,有些景点之间会用一座桥连接.为了方便游客到达每个景点但 ...
- 【BZOJ5210】最大连通子块和 树剖线段树+动态DP
[BZOJ5210]最大连通子块和 Description 给出一棵n个点.以1为根的有根树,点有点权.要求支持如下两种操作: M x y:将点x的点权改为y: Q x:求以x为根的子树的最大连通子块 ...
- [LNOI2014]LCA(树剖+线段树)
\(\%\%\% Fading\) 此题是他第一道黑题(我的第一道黑题是蒲公英) 一直不敢开,后来发现是差分一下,将询问离线,树剖+线段树维护即可 \(Code\ Below:\) #include ...
随机推荐
- js日期处理多少/月/天以前
var date = new Date(); var year = date.getFullYear(); var month = (date.getMonth() + 1) > 10 ? da ...
- ecshop文件架构
/*ECShop 2.5.1 的结构图及各文件相应功能介绍ECShop2.5.1_Beta upload 的目录┣ activity.php 活动列表┣ affiche.php 广告处理文件┣ aff ...
- MySQL-存储引擎-Myisam
mysql> create table myisam_char(name char(10)) engine=myisam; Query OK, 0 rows affected (0.01 sec ...
- JavaScript之创建对象的模式
使用Object的构造函数可以创建对象或者使用对象字面量来创建单个对象,但是这些方法有一个明显的缺点:使用相同的一个接口创建很多对象,会产生大量的重复代码. (一)工厂模式 这种模式抽象了创建具体对象 ...
- SpringMVC-源码-图解
- Eclipse开发Java的简单配置
目录 Eclipse配置Java开发环境 1.配置JDK 2.切换项目的JDK版本 3.配置maven 4.配置lombok Eclipse配置Java开发环境 时隔N年重新使用Eclipse,对一些 ...
- Merchant
\(get\)二分新用法. 每道题都有答案范围提示,以前只是以为是用来提示用什么类型输出的. 从来没想过直接用它来二分. 这道题真的刷新了我的认知啊...... 整道题的复杂度是\ ...
- Python中正则表达式简介
目录 一.什么是正则表达式 二.正则表达式的基础知识 1. 原子 1)普通字符作为原子 2)非打印字符作为原子 3) 通用字符作为原子 4) 原子表 2. 元字符 1)任意匹配元字符 2)边界限制元字 ...
- 在 Docker 的 CentOS7 镜像 中安装 mysql
在 Docker 的 CentOS7 镜像 中安装 mysql 本来以为是个很简单的过程居然折腾了这么久,之前部署云服务器时也没有好好地记录,因此记录下. 特别提醒:本文的操作环境是在 Docker ...
- mybatis零碎
< < 小于号 > > 大于号 & & 和 ' ' 单 ...