bzoj2594 [Wc2006]水管局长数据加强版——LCT
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2594
时间倒序一下,就是 魔法森林 那道题;
有个不解的地方,是 access 里面关于 pushup 的地方,两种写法都可以,但不明白没有注释掉的写法为什么也可以。
代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int const maxn=,maxm=1e6+;
int n,m,Q,val[maxn],mx[maxn],fa[maxn],pre[maxn],c[maxn][],sta[maxn],rev[maxn];
struct N{int hd,to,w,id; bool b;}edge[maxm];
struct T{int u,v,ans,f,id;}q[];
bool cmp(N x,N y){return x.w<y.w;}
bool cmp2(N x,N y){return x.hd<y.hd || (x.hd==y.hd && x.to<y.to);}
bool cmpp(N x,N y){return x.id<y.id;}
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
int rd()
{
int ret=,f=; char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=-; ch=getchar();}
while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
return ret*f;
}
int get(int x,int y)
{
int l=,r=m;
while(l<=r)
{
int mid=((l+r)>>);
if(edge[mid].hd==x && edge[mid].to==y)return mid;
if(edge[mid].hd<x || (edge[mid].hd==x && edge[mid].to<y))l=mid+;
else r=mid-;
}
}
bool isroot(int x)
{
return c[pre[x]][]!=x && c[pre[x]][]!=x;
}
void pushup(int x)
{
int ls=c[x][],rs=c[x][];
mx[x]=x;
if(val[mx[ls]]>val[mx[x]])mx[x]=mx[ls];
if(val[mx[rs]]>val[mx[x]])mx[x]=mx[rs];
}
void reverse(int x)
{
int ls=c[x][],rs=c[x][];
if(rev[x])
{
rev[ls]^=; rev[rs]^=; rev[x]=;
swap(c[x][],c[x][]);
}
}
void rotate(int x)
{
int y=pre[x],z=pre[y],d=(c[y][]==x);
if(!isroot(y))c[z][c[z][]==y]=x;
pre[x]=z; pre[y]=x; pre[c[x][!d]]=y;
c[y][d]=c[x][!d]; c[x][!d]=y;
pushup(y); pushup(x);
}
void splay(int x)
{
int top=; sta[++top]=x;
for(int i=x;!isroot(i);i=pre[i])sta[++top]=pre[i];
for(int i=top;i;i--)reverse(sta[i]);
for(;!isroot(x);rotate(x))
{
int y=pre[x],z=pre[y];
if(isroot(y))continue;
((c[y][]==x)^(c[z][]==y))?rotate(x):rotate(y);
}
pushup(x);//和下面配套
}
void access(int x)
{
for(int t=;x;c[x][]=t,t=x,x=pre[x])splay(x);//这样为什么也可以?
}
//void access(int x)
//{
// int t=0;
// while(x)
// {
// splay(x);c[x][1]=t;pushup(x);t=x;x=pre[x];
// }
//}
void makeroot(int x)
{
access(x); splay(x); rev[x]^=;
}
void link(int x,int y)
{
makeroot(x); pre[x]=y;//
}
void cut(int x,int y)
{
makeroot(x); access(y); splay(y);
pre[x]=; c[y][]=;
}
int query(int x,int y)
{
makeroot(x); access(y); splay(y); return mx[y];
}
int main()
{
n=rd(); m=rd(); Q=rd();
for(int i=;i<=n;i++)fa[i]=i;
for(int i=;i<=m;i++)
{
edge[i].hd=rd(); edge[i].to=rd(); edge[i].w=rd();
if(edge[i].hd>edge[i].to)swap(edge[i].hd,edge[i].to);
}
sort(edge+,edge+m+,cmp);//w
for(int i=;i<=m;i++)
{
edge[i].id=i;
val[i+n]=edge[i].w; mx[i+n]=i+n;//
}
sort(edge+,edge+m+,cmp2);//u,v
for(int i=;i<=Q;i++)
{
q[i].f=rd(); q[i].u=rd(); q[i].v=rd();
if(q[i].f==)
{
if(q[i].u>q[i].v)swap(q[i].u,q[i].v);
int t=get(q[i].u,q[i].v);
// edge[t].b=1; q[i].id=edge[i].id;
edge[t].b=; q[i].id=edge[t].id;
}
}
sort(edge+,edge+m+,cmpp);//w
for(int i=,tot=;i<=m;i++)
{
if(edge[i].b)continue;//kruskal 连不删的边
int u=edge[i].hd,v=edge[i].to;
if(find(u)!=find(v))
{
fa[find(u)]=find(v);
link(u,i+n); link(v,i+n);
tot++;
if(tot==n-)break;
}
}
for(int i=Q;i;i--)
{
int u=q[i].u,v=q[i].v;
if(q[i].f==)q[i].ans=val[query(u,v)];
else
{
int k=q[i].id,t=query(u,v);
if(edge[k].w<val[t])
{
cut(t,edge[t-n].hd); cut(t,edge[t-n].to);
link(u,k+n); link(v,k+n);
}
}
}
for(int i=;i<=Q;i++)
if(q[i].f==)printf("%d\n",q[i].ans);
return ;
}
bzoj2594 [Wc2006]水管局长数据加强版——LCT的更多相关文章
- [bzoj2594][Wc2006]水管局长数据加强版 (lct)
论蒟蒻的自我修养T_T.. 和noi2014魔法森林基本一样...然而数据范围大得sxbk...UPD:这题如果用lct判联通的话可能会被卡到O(mlogm)..所以最好还是用并查集吧 一开始数组开太 ...
- BZOJ2594 [Wc2006]水管局长数据加强版 LCT kruskal
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ2594 题意概括 N个点的图,M条带权边.(N<=100000,M<=1000000) ...
- [bzoj2594][Wc2006]水管局长数据加强版——lct+离线
Brief Description 您有一个无向带权图,您需要支持两种操作. 询问两个点之间的最大权最小路径. 删除一条边. Algorithm Design 我们首先提出一个猜想:最优路径一定在原图 ...
- BZOJ2594: [Wc2006]水管局长数据加强版
题解: 裸LCT+离线+二分+MST... 代码:(几乎摘抄自hzwer) #include<cstdio> #include<cstdlib> #include<cma ...
- BZOJ 2594: [Wc2006]水管局长数据加强版 [LCT kruskal]
2594: [Wc2006]水管局长数据加强版 Time Limit: 25 Sec Memory Limit: 128 MBSubmit: 2917 Solved: 918[Submit][St ...
- BZOJ 2594: [Wc2006]水管局长数据加强版( LCT )
离线然后就是维护加边的动态MST, Link cut tree秒掉..不过我写+调了好久...时间复杂度O(NlogN + MlogM) ------------------------------- ...
- BZOJ2594 [Wc2006]水管局长数据加强版 【LCT维护最小生成树】
题目 SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的 ...
- [BZOJ2594] [Wc2006]水管局长数据加强版(LCT + kruskal + 离线)
传送门 WC这个题真是丧心病狂啊,就是想学习一下怎么处理边权,给我来了这么一个破题! ORZ hzwer 临摹黄学长代码233 但还是复杂的一匹 理一下思路吧 题目大意:给定一个无向图,多次删除图中的 ...
- BZOJ 2594: [Wc2006]水管局长数据加强版 (LCT维护最小生成树)
离线做,把删边转化为加边,那么如果加边的两个点不连通,直接连就行了.如果联通就找他们之间的瓶颈边,判断一下当前边是否更优,如果更优就cut掉瓶颈边,加上当前边. 那怎么维护瓶颈边呢?把边也看做点,向两 ...
随机推荐
- //……关于HTTP与HTTPS
图1 图2 第一张访问域名http://www.tsinghua.edu.cn,谷歌浏览器提示不安全链接,第二张是https://www.12306.cn/index,浏览器显示安全,为什么会这样子呢 ...
- python之抽象 2014-4-6
#抽象 8.40am-1.懒惰即美德2.抽象和结构3.创建函数 内建的callable 函数可以判定函数是否可以调用 >>> import math >>> x=1 ...
- hihoCoder#1054 滑动解锁
原题地址 回溯搜索 对于每个待枚举的点,检查: 1. 度数检查:是否违反了出度入度限制.因为生成的路径除了首尾节点外,其他节点的出度和入度只能为2 2. 共线检查:是否违反了共线条件.即跨越了尚未枚举 ...
- Django开发:(3.1)ORM:单表操作
MVC或者MVC框架中包括一个重要的部分,就是ORM,它实现了数据模型与数据库的解耦,即数据模型的设计不需要依赖于特定的数据库,通过简单的配置就可以轻松更换数据库,这极大的减轻了开发人员的工作量,不需 ...
- node的express框架,核心第三方模块body-parser 获取我们所有post请求传过来数据
- 安装 body-parser模块- npm install body-parser -S - 调用- let bodyParser=require('body-parser'); - 设置中间件- ...
- HttpClientUtil 工具类 实现跨域请求数据
package com.xxx.common.util; import java.io.IOException; import java.net.URI; import java.util.Array ...
- JAVA分布式架构
- mkswap,swapon, swapoff命令:创建交换分区
linux支持虚拟内存,用作虚拟内存的硬盘部分被称为交互空间(swap space),虚拟内存是指使用磁盘当作内存的扩展,这样可用内存的大小就相应的增大了.内核会将暂时不用的内存块的内容写到硬盘上,从 ...
- GDAL源码编译
转自阿Fai, GDAL源码编译 在这里,我使用源码编译出C#可以使用的dll静态文件. 一.简单的编译 1.简单的认识 首先进入GDAL的源代码目录,可以看到有几个sln为后缀的文件名,比如make ...
- [转]一个完整的Installshield安装程序实例
@import url("http://files.cnblogs.com/files/go-jzg/vs.css"); --> Installshield安装程序实例—基本 ...