沉迷Link-Cut tree无法自拔之:[BZOJ2594][Wc2006]水管局长数据加强版
来自蒟蒻 \(Hero \_of \_Someone\) 的 \(LCT\) 学习笔记
$
$
这应该算是道套路题吧, 如果将图中的边转换成点, 再将边权变点权, 就可以用 \(LCT\) 来维护了
这道题的基本做法就是, 用 \(LCT\) 来动态地维护最小生成树,
如果这样做的话, 题目中要求的删边操作就不太好搞,
但是既然只有删边操作的话, 我们就可以考虑离线处理,
将不会被删除的边先加进图中跑 \(Kruskal\) , 然后化删边为添边, 倒着来处理每一组询问.
$
$
此外, 用 \(LCT\) 维护 \(MST\) , 就是在添边的时候如果遇到环且环上最长的边边权大于当前边, 就将最大边 \(cut\) , 再将当前边添入
//made by Hero_of_Someone
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#define N (100010)
#define M (1000010)
#define RG register
using namespace std;
inline int gi(){ RG int x=0; RG char ch=getchar(); while(ch<'0'||ch>'9') ch=getchar();
while('0'<=ch&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar(); return x; }
int n,m,q;
struct Edge{
int u,v,w;
bool operator<(const Edge&a)const{
return u==a.u?v<a.v:u<a.u;
}
bool operator<(const pair<int,int>& a)const{
return u==a.first?v<a.second:u<a.first;
}
}e[M];
struct Que{ int op,x,y; }Q[N];
//-----------------------------------------------------------
int Fa[N];
inline int Find(RG int x){ return Fa[x]==x?x:Fa[x]=Find(Fa[x]); }
struct EEE{
int u,v,w,id;
bool operator<(const EEE& a)const{ return w<a.w; }
void operator=(const Edge& a){
u=a.u,v=a.v,w=a.w,id=(&a)-e;
}
}E[M];
int num,head[N],nxt[N<<1],to[N<<1],w[N<<1]; //w为该边原来的编号
inline void add(int u,int v,int d){
nxt[++num]=head[u];to[num]=v;w[num]=d;head[u]=num;
nxt[++num]=head[v];to[num]=u;w[num]=d;head[v]=num;
}
bool v[M];
inline void kruskal(){
RG int top=0;
for(RG int i=1;i<=n;i++) Fa[i]=i;
for(RG int i=1;i<=m;i++){
if(v[i]) continue;
E[++top]=e[i];
}
sort(E+1,E+top+1);
for(RG int i=1;i<=top;i++){
RG int x=E[i].u,y=E[i].v;
RG int fx=Find(x),fy=Find(y);
if(fx==fy) continue;
Fa[fx]=fy; add(x,y,E[i].id);
}
}
inline void init(){
n=gi(),m=gi(),q=gi();
for(RG int i=1;i<=m;i++){
e[i].u=gi(),e[i].v=gi(),e[i].w=gi();
}
sort(e+1,e+m+1);
for(RG int i=1;i<=q;i++){
Q[i].op=gi(),Q[i].x=gi(),Q[i].y=gi();
if(Q[i].op==2){
RG int V=lower_bound(e+1,e+m+1,make_pair(Q[i].x,Q[i].y))-e;
Q[i].x=V; v[V]=1;
}
}q++;
kruskal();
}
//---------------------------------------------------------------
int val[N+M],Max[N+M];
int ch[N+M][2],fa[N+M],rev[N+M];
inline void cur(int x,int y){ val[x]=Max[x]=y; }
int vis[N];
inline void dfs(RG int x,RG int f,RG int d){
if(vis[x]) return ;
vis[x]=1; cur(x,0);
if(f){ cur(n+d,d); fa[n+d]=f; fa[x]=n+d; } //将边转换成点插入树中
for(RG int i=head[x];i;i=nxt[i])
dfs(to[i],x,w[i]);
}
inline bool cnm(int x,int y){ return e[x].w<e[y].w; }
inline void up(int x){
Max[x]=max(Max[ch[x][0]],Max[ch[x][1]],cnm);
Max[x]=max(Max[x],val[x],cnm);
}
inline void reverse(int x){
swap(ch[x][0],ch[x][1]);
rev[x]^=1;
}
inline void down(int x){
if(!rev[x]) return ;
reverse(ch[x][0]);
reverse(ch[x][1]);
rev[x]=0;
}
inline bool is_root(int x){ return ch[fa[x]][0]!=x && x!=ch[fa[x]][1]; }
inline bool lr(int x){ return x==ch[fa[x]][1]; }
inline void rotate(int x){
RG int y=fa[x],z=fa[y],k=lr(x);
if(!is_root(y)) ch[z][lr(y)]=x;
fa[x]=z; fa[ch[x][k^1]]=y; fa[y]=x;
ch[y][k]=ch[x][k^1]; ch[x][k^1]=y;
up(y); up(x);
}
int st[N+M];
inline void splay(int x){
RG int y=x,top=0;
while(1){
st[++top]=y;
if(is_root(y)) break;
y=fa[y];
}
for(RG int i=top;i;i--) down(st[i]);
while(!is_root(x)){
if(!is_root(fa[x])) rotate(lr(x)^lr(fa[x])?x:fa[x]);
rotate(x);
}
}
inline void access(int x){
RG int y=0;
while(x){ splay(x);
ch[x][1]=y; fa[y]=x;
up(x); y=x; x=fa[x];
}
}
inline void make_root(int x){
access(x); splay(x); reverse(x);
}
inline int query(int x,int y){
make_root(x); access(y); splay(y);
return Max[y];
}
inline int find(int x){
while(fa[x]) x=fa[x];
return x;
}
inline void link(int x,int y){
if(find(x)==find(y)) return ;
make_root(x); fa[x]=y;
}
inline void cut(int x,int y){
make_root(x); access(y); splay(y);
if(ch[y][0]==x) ch[y][0]=0,fa[x]=0,up(y);
}
inline void Insert(int id){
RG int x=e[id].u,y=e[id].v;
if(find(x)==find(y)){
RG int tmp=query(x,y);
if(e[tmp].w<=e[id].w) return ;
cut(n+tmp,e[tmp].u);
cut(n+tmp,e[tmp].v);
}
cur(n+id,id);
link(x,n+id);
link(y,n+id);
}
int Ans[N],top;
inline void work(){
for(RG int i=1;i<=n;i++) dfs(i,0,0);
while(--q){
if(Q[q].op==2) Insert(Q[q].x);
else Ans[++top]=e[query(Q[q].x,Q[q].y)].w;
}
for(RG int i=top;i;i--) printf("%d\n",Ans[i]);
}
int main(){ init(); work(); return 0; }
沉迷Link-Cut tree无法自拔之:[BZOJ2594][Wc2006]水管局长数据加强版的更多相关文章
- BZOJ2594: [Wc2006]水管局长数据加强版
题解: 裸LCT+离线+二分+MST... 代码:(几乎摘抄自hzwer) #include<cstdio> #include<cstdlib> #include<cma ...
- [bzoj2594][Wc2006]水管局长数据加强版 (lct)
论蒟蒻的自我修养T_T.. 和noi2014魔法森林基本一样...然而数据范围大得sxbk...UPD:这题如果用lct判联通的话可能会被卡到O(mlogm)..所以最好还是用并查集吧 一开始数组开太 ...
- BZOJ2594 [Wc2006]水管局长数据加强版 【LCT维护最小生成树】
题目 SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的 ...
- BZOJ2594 [Wc2006]水管局长数据加强版 LCT kruskal
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ2594 题意概括 N个点的图,M条带权边.(N<=100000,M<=1000000) ...
- [bzoj2594][Wc2006]水管局长数据加强版——lct+离线
Brief Description 您有一个无向带权图,您需要支持两种操作. 询问两个点之间的最大权最小路径. 删除一条边. Algorithm Design 我们首先提出一个猜想:最优路径一定在原图 ...
- bzoj2594 [Wc2006]水管局长数据加强版——LCT
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2594 时间倒序一下,就是 魔法森林 那道题: 有个不解的地方,是 access 里面关于 p ...
- [BZOJ2594] [Wc2006]水管局长数据加强版(LCT + kruskal + 离线)
传送门 WC这个题真是丧心病狂啊,就是想学习一下怎么处理边权,给我来了这么一个破题! ORZ hzwer 临摹黄学长代码233 但还是复杂的一匹 理一下思路吧 题目大意:给定一个无向图,多次删除图中的 ...
- 【bzoj2594】 Wc2006—水管局长数据加强版
http://www.lydsy.com/JudgeOnline/problem.php?id=2594 (题目链接) 题意 给出一个带边权的无向简单,要求维护两个操作,删除${u,v}$之间的连边: ...
- 【bzoj2594】[Wc2006]水管局长数据加强版
真是神题 当时调了几天没调出来 后来没管了 当时把fread去掉就TLE,加上就RE 一直在底下跟网上的程序拍,尝试各种优化常数都没用 拍出几组不一样的,发现我是对的,醉了,网上那个是怎么过的 记一下 ...
随机推荐
- HNOI2019 多边形 polygon
HNOI2019 多边形 polygon https://www.luogu.org/problemnew/show/P5288 这题镪啊... 首先堆结论: 显然终止状态一定是所有边都连向n了 根据 ...
- ABP+AdminLTE+Bootstrap Table权限管理系统第六节--abp控制器扩展及json封装以及6种处理时间格式化的方法
返回总目录:ABP+AdminLTE+Bootstrap Table权限管理系统一期 一,控制器AbpController 说完了Swagger ui 我们再来说一下abp对控制器的处理和json的封 ...
- slurm.conf系统初始配置
#slurm集群配置 ##集群名称 ClusterName=myslurm ##主控制器的主机名 ControlMachine=node11 ##主控制器的IP地址 ControlAddr=192.1 ...
- 【nodejs】让nodejs像后端mvc框架(asp.net mvc)一orm篇【如EF般丝滑】typeorm介绍(8/8)
文章目录 前情概要 在使用nodejs开发过程中,刚好碰到需要做一个小工具,需要用到数据库存储功能.而我又比较懒,一个小功能不想搞一个nodejs项目,又搞一个后端项目.不如直接在nodejs里面把对 ...
- 基于Nginx+Keepalived的LB服务监控(邮件报警)
IDC两台机器上部署了Nginx+Keepalived主从模式的LB代理负载层,现在需要对LB进行每日巡检和服务监控,利用SendEmail邮件监控. 0)SendEmail部署 参考:http:// ...
- Pupet自动化管理环境部署记录
废话不多说了,下面记录下Puppet在Centos下的部署过程: puppet是什么puppet是一种基于ruby语言开发的Lnux.Unix.windows平台的集中配置管理系统.它使用自有的pup ...
- Coolest Ski Route-不定起点和终点----在有向变的情况下---求最长路
这题最开始给你了N个点,M条边,边是单向边,问不指定起点和终点,最长路是什么??? 脑补一下,不定起点和终点的最短路,用弗洛伊德算法搞一搞,但是...那个垃圾算法的复杂度是N^3的,但是这个算法的M高 ...
- A. A Prank
题意 有数列从小到大排列,都是不同范围1~ 1000,问你最多去掉多少个数字还能复原 由于wrong很多发所以写一下 链接 [http://codeforces.com/contest/1062/pr ...
- 第三个spring冲刺第2天
今天我们有了计时功能的实现,并且在考虑如何使得计时器美观好看达到我们的要求,对此我们换了不同的背景,时钟框,效果还有待查看.
- Oracle PLSQL 客户端 连接Oracle12.2 出现权限问题的解决办法以及绿色版Oracle客户端的使用.
1. 同事反馈登录最新的oracle12.2 的数据库时登录不上报错: 2. 记得当时查过资料, Oracle 在12.2 增加了客户端连接数据库的加密级别 比较早的oracle客户端比如11.2.0 ...