洛谷2543AHOI2005]航线规划 (树剖+线段树+割边思路)
这个题的思路还是比较巧妙的。
首先,我们发现操作只有删除和询问两种,而删除并不好维护连通性和割边之类的信息。
所以我们不妨像WC2006水管局长那样,将询问离线,然后把操作转化成加边和询问。
然后,我们会发现,若存在一条边\(x->y\),那么原本x到y的所有割边,都会变成非割边。
那意味着什么呢?
似乎加边操作,可以直接转化成区间修改。
那我们就可以首先对不涉及删除边,建一个生成树。(题目保证一定合法)
那么对于一棵树,所有的边都是割边,所以一开始所有的边的边权都是1(这里为了修改方便,我们将边权直接转化成点权了),也就是说树上除了根以外,权值都是1.
然后依次插入那些没有被删除,但是没有在生成树里面的边。每插入一条边,就涉及到一次链修改,将一条链的点的权值变成\(0\)。
然后操作中的加边也是同理。
对于询问的话,直接询问\((x,y)\)的路径和就好。
但是有一个需要注意的地方就是。由于我们边权转点权,所以\(lca\)处的点不属于路径,修改和询问的时候都需要注意的。
有些细节直接看代码吧
#include<bits/stdc++.h>
#define mk make_pair
#define pb push_back
#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+1e2;
const int maxm = 2*maxn;
struct Node{
int opt,x,y;
};
Node a[maxm];
int f[4*maxn];
int add[4*maxn];
int n,m,fa[maxn];
int dfn[maxn],size[maxn],newnum[maxn],deep[maxn];
int top[maxn],son[maxn];
int point[maxn],nxt[maxm],to[maxm];
int in[maxn],tag[maxn];
int newval[maxn];
map<pair<int,int>,int> mp;
int cnt;
void addedge(int x,int y)
{
nxt[++cnt]=point[x];
to[cnt]=y;
point[x]=cnt;
}
void dfs(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;
dfs(p,x,dep+1);
size[x]+=size[p];
if (size[p]>maxson)
{
maxson=size[p];
son[x]=p;
}
}
}
int tot;
void dfs1(int x,int chain)
{
newnum[x]=++tot;
//cout<<x<<" "<<tot<<"****"<<endl;
newval[tot]=1;
top[x]=chain;
if (!son[x]) return;
dfs1(son[x],chain);
for (int i=point[x];i;i=nxt[i])
{
int p = to[i];
if(!newnum[p])
{
dfs1(p,p);
}
}
}
void up(int root)
{
f[root]=f[2*root]+f[2*root+1];
}
void pushdown(int root,int l,int r)
{
if (add[root]!=-1)
{
add[2*root]=add[root];
add[2*root+1]=add[root];
int mid = l+r >> 1;
f[2*root]=(mid-l+1)*add[root];
f[2*root+1]=(r-mid)*add[root];
add[root]=-1;
}
}
void build(int root,int l,int r)
{
add[root]=-1;
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 y,int p)
{
if (x<=l && r<=y)
{
add[root]=p;
f[root]=(r-l+1)*p;
return;
}
pushdown(root,l,r);
int mid = l+r >> 1;
if (x<=mid) update(2*root,l,mid,x,y,p);
if (y>mid) update(2*root+1,mid+1,r,x,y,p);
up(root);
}
int query(int root,int l,int r,int x,int y)
{
if (x<=l && r<=y)
{
return f[root];
}
pushdown(root,l,r);
int mid = l+r >> 1;
int ans = 0;
if (x<=mid) ans=ans+query(2*root,l,mid,x,y);
if (y>mid) ans=ans+query(2*root+1,mid+1,r,x,y);
return ans;
}
void treeadd(int x,int y,int z)
{
while (top[x]!=top[y])
{
if (deep[top[x]]<deep[top[y]]) swap(x,y);
update(1,1,n,newnum[top[x]],newnum[x],z);
x=fa[top[x]];
}
if (deep[x]>deep[y]) swap(x,y);
int pre = query(1,1,n,newnum[x],newnum[x]);
update(1,1,n,newnum[x],newnum[y],z);
update(1,1,n,newnum[x],newnum[x],pre);
}
int treesum(int x,int y)
{
int ans=0;
while (top[x]!=top[y])
{
if (deep[top[x]]<deep[top[y]]) swap(x,y);
ans=ans+query(1,1,n,newnum[top[x]],newnum[x]);
x=fa[top[x]];
}
if (deep[x]>deep[y]) swap(x,y);
int pre = query(1,1,n,newnum[x],newnum[x]);
ans=ans+query(1,1,n,newnum[x],newnum[y]);
ans-=pre;
return ans;
}
int x[maxm],y[maxm];
int ffa[maxn];
int ans[maxm];
int find(int x)
{
if (ffa[x]!=x) ffa[x]=find(ffa[x]);
return ffa[x];
}
int main()
{
n=read(),m=read();
for (int i=1;i<=n;i++) ffa[i]=i;
for (int i=1;i<=m;i++)
{
x[i]=read(),y[i]=read();
mp[mk(x[i],y[i])]=mp[mk(y[i],x[i])]=i;
}
int tmp=0;
while (1)
{
int opt=read();
if(opt==-1) break;
a[++tmp].opt=opt;
a[tmp].x=read();
a[tmp].y=read();
}
for (int i=1;i<=tmp;i++) if (a[i].opt==0) tag[mp[mk(a[i].x,a[i].y)]]=1;
for (int i=1;i<=m;i++)
{
if (tag[i]) continue;
int f1=find(x[i]);
int f2=find(y[i]);
if (f1==f2) continue;
ffa[f1]=f2;
addedge(x[i],y[i]);
addedge(y[i],x[i]);
in[i]=1;
}
dfs(1,0,1);
dfs1(1,1);
newval[1]=0;
build(1,1,n);
int tmp1=0;
for (int i=1;i<=m;i++)
{
if (tag[i] || in[i]) continue;
treeadd(x[i],y[i],0);
}
for (int i=tmp;i>=1;i--)
{
if (a[i].opt==0)
{
treeadd(a[i].x,a[i].y,0);
}
else
{
ans[++tmp1]=treesum(a[i].x,a[i].y);
}
}
for (int i=tmp1;i>=1;i--)
{
cout<<ans[i]<<"\n";
}
return 0;
}
洛谷2543AHOI2005]航线规划 (树剖+线段树+割边思路)的更多相关文章
- 洛谷P4315 月下“毛景树”(树剖+线段树)
传送门 woc这该死的码农题…… 把每一条边转化为它连接的两点中深度较深的那一个,然后就可以用树剖+线段树对路径进行修改了 然后顺便注意在上面这种转化之后,树剖的时候不能搞$LCA$ 然后是几个注意点 ...
- BZOJ_2238_Mst_树剖+线段树
BZOJ_2238_Mst_树剖+线段树 Description 给出一个N个点M条边的无向带权图,以及Q个询问,每次询问在图中删掉一条边后图的最小生成树.(各询问间独立,每次询问不对之后的询问产生影 ...
- BZOJ_4551_[Tjoi2016&Heoi2016]树_树剖+线段树
BZOJ_4551_[Tjoi2016&Heoi2016]树_树剖+线段树 Description 在2016年,佳媛姐姐刚刚学习了树,非常开心.现在他想解决这样一个问题:给定一颗有根树(根为 ...
- 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 ...
- [CF1007D]Ants[2-SAT+树剖+线段树优化建图]
题意 我们用路径 \((u, v)\) 表示一棵树上从结点 \(u\) 到结点 \(v\) 的最短路径. 给定一棵由 \(n\) 个结点构成的树.你需要用 \(m\) 种不同的颜色为这棵树的树边染色, ...
- LOJ#3088. 「GXOI / GZOI2019」旧词(树剖+线段树)
题面 传送门 题解 先考虑\(k=1\)的情况,我们可以离线处理,从小到大对于每一个\(i\),令\(1\)到\(i\)的路径上每个节点权值增加\(1\),然后对于所有\(x=i\)的询问查一下\(y ...
- BZOJ3531-[Sdoi2014]旅行(树剖+线段树动态开点)
传送门 完了今天才知道原来线段树的动态开点和主席树是不一样的啊 我们先考虑没有宗教信仰的限制,那么就是一个很明显的树剖+线段树,路径查询最大值以及路径和 然后有了宗教信仰的限制该怎么做呢? 先考虑暴力 ...
随机推荐
- idea控制台中文乱码解决办法
也可以通过idea右下角的设置,但是properties文件是不能设置的,这个只能在file->setting->file encodings 设置
- HashSet的存储原理
HashSet的底层用哈希散列表来存储对象(默认长度为16的数组),假如: Set set=new HashSet(); set.add(obj); 内部存储过程为:定义h=obj.hashCode, ...
- Linux下SSH以及SSH秘钥
一.基于秘钥方式实现远程连接 第一步:创建密钥对(在管理端服务器上操作) 中间的输入项可以直接回车 ssh-keygen -t dsa 第二步:分发公钥(在管理端服务器执行) 这个步骤需要输入一个ye ...
- Docker(23)- 注册 docker hub 的账号
如果你还想从头学起 Docker,可以看看这个系列的文章哦! https://www.cnblogs.com/poloyy/category/1870863.html 前言 Docker Hub 是 ...
- Linu常用日志分析实战
日志结构分析 分析日志状态码所在位置为第九个 遍历取出第一行日志的每个字段 //取出第一行日志 awk 'NR==1{for(i=1;i<=NF;i++)print i"= " ...
- vue-router路由钩子
路由跳转前后,需要做某些操作,这时就可以使用路由钩子来监听路由的变化. 接收三个参数: to: Route: 即将要进入的目标路由对象 from: Route: 当前导航正要离开的路由 next: F ...
- STM32L0系列EEPROM中结构体的读取
在STM32L0中操作EEPROM本来参考了上篇操作FLASH的方法,多多少少都有些问题.我觉得可能是结构体在转换成其他变量的时候出了问题. 比如下面这段代码,在Windows上可以正常运行(使用g+ ...
- PHP中的MySQLi扩展学习(五)MySQLI_STMT对象操作
就像 PDO 中的 PDO_Statment 对象一样,MySQLI_STMT 对象也是一个预处理语句所形成的对象,专门用来操作 MySQLi 所生成的预处理语句的.其实操作方式之类也都比较相似,不外 ...
- PDO操作大数据对象
一般在数据库中,我们保存的都只是 int . varchar 类型的数据,一是因为现代的关系型数据库对于这些内容会有很多的优化,二是大部分的索引也无法施加在内容过多的字段上,比如说 text 类型的字 ...
- django使用celery搭配redis配置定时任务
已经安装环境: Python3.6 django==2.1.8(用2.2.2需要升级sqlite3) 项目名称:ceshiproject APP名称:ceshi 第一步:centos7下首先安装r ...