【校内模拟】坦克（LCT）

最新推荐文章于 2023-02-14 19:37:30 发布
最新推荐文章于 2023-02-14 19:37:30 发布
阅读量200 收藏
点赞数
分类专栏： LCT 校内模拟
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zxyoi_dreamer/article/details/102922113
版权

简要题意见此：here

题解：

其实不是平面图也可以做，化边为点，边的权值设为删除时间，然后双指针就行了。

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	template<typename T>T get(){
		char c;T num;
		while(!isdigit(c=gc()));num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}inline int gi(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

#define y1 y_1
#define y2 y_2

cs int N=1e5+7;

int n,m1,m2,tot,ct1,ct2,ec;
struct edge{int u,v;}e[N<<2|3];

int vl[N*7],tim[N*7];
namespace LCT{
	cs int N=::N*7;
#define lc(u) son[u][0]
#define rc(u) son[u][1]
	int son[N][2],fa[N],sum[N],sub[N];
	bool rev[N];int mx[N];
	inline void rever(int u){rev[u]^=1;std::swap(lc(u),rc(u));}
	inline void pushdn(int u){
		if(rev[u]){
			if(lc(u))rever(lc(u));
			if(rc(u))rever(rc(u));
			rev[u]=0;
		}
	}
	inline void pushup(int u){
		sum[u]=sum[lc(u)]+sum[rc(u)]+vl[u]+sub[u];mx[u]=u;
		if(tim[mx[lc(u)]]>tim[mx[u]])mx[u]=mx[lc(u)];
		if(tim[mx[rc(u)]]>tim[mx[u]])mx[u]=mx[rc(u)];
	}
	inline bool get(int u){return rc(fa[u])==u;}
	inline bool isrt(int u){return lc(fa[u])!=u&&rc(fa[u])!=u;}
	inline void Rotate(int u){
		int p=fa[u],pp=fa[p],d=get(u);
		if(!isrt(p))son[pp][get(p)]=u;
		son[p][d]=son[u][!d],son[u][!d]=p;
		fa[p]=u,fa[u]=pp,fa[son[p][d]]=p;
		pushup(p);pushup(u);
	}
	void Splay(int u){
		static int q[N],qn;q[qn=1]=u;
		for(int re p=u;!isrt(p);p=fa[p])q[++qn]=fa[p];
		for(int re i=qn;i;--i)pushdn(q[i]);
		for(int re p=fa[u];!isrt(u);Rotate(u),p=fa[u])
		if(!isrt(p))Rotate(get(u)==get(p)?p:u);pushup(u); 
	}
	void access(int u){
		for(int re ch=0;u;u=fa[ch=u]){
			Splay(u);if(rc(u))sub[u]+=sum[rc(u)];
			if(ch)sub[u]-=sum[ch];rc(u)=ch;
			pushup(u);
		}
	}
	void makert(int u){access(u);Splay(u),rever(u);}
	int find_rt(int u){
		access(u);Splay(u);
		while(pushdn(u),lc(u))u=lc(u);
		Splay(u);return u;
	}
	bool iscon(int u,int v){makert(u);return find_rt(v)==u;}
	void link(int u,int v){
		makert(u);access(v);Splay(v);
		fa[u]=v;sub[v]+=sum[u];pushup(v);
	}
	void cut(int u,int v){
		makert(u);access(v);Splay(v);
		fa[lc(v)]=0,lc(v)=0,pushup(v),pushup(u);
	}
	int qsz(int u){access(u);Splay(u);return sum[u];}
	int find_e(int u,int v){
		makert(u),access(v),Splay(v);
		return mx[v];
	}
#undef rc
#undef lc
}

int y1[N],y2[N];
int p1[N],p2[N],po[N];
std::vector<int> e1[N],e2[N];
std::vector<int> nd1[N],nd2[N];
int b[N];
signed main(){
#ifdef zxyoi
	freopen("mt.in","r",stdin);
#endif
	n=gi(),m1=gi(),m2=gi();
	tot=n+m1+n+m2;
	for(int re i=1;i<=m1;++i)p1[i]=++tot;
	for(int re i=1;i<=m2;++i)p2[i]=++tot;
	for(int re i=1;i<=n;++i)po[i]=++tot,vl[tot]=1;
	for(int re i=1;i<=tot;++i)LCT::pushup(tot);
	for(int re i=1;i<=m1;++i)gi(),b[i]=y1[i]=gi();
	std::sort(b+1,b+m1+1);ct1=std::unique(b+1,b+m1+1)-b-1;
	for(int re i=1;i<=m1;++i){
		y1[i]=std::lower_bound(b+1,b+ct1+1,y1[i])-b;
		nd1[y1[i]].push_back(i);
	}
	for(int re i=1;i<n+m1;++i){
		int t=gi(),u=gi(),v=gi();
		if(t==1){
			e[++ec]=(edge){po[u],p1[v]};
			e1[v].push_back(ec);tim[ec]=y1[v];
			LCT::link(ec,po[u]);LCT::link(ec,p1[v]);
		}else {
			e[++ec]=(edge){p1[u],p1[v]};
			if(y1[u]>y1[v])e1[u].push_back(ec),tim[ec]=y1[u];
			else e1[v].push_back(ec),tim[ec]=y1[v];
			LCT::link(ec,p1[u]);LCT::link(ec,p1[v]);
		}
	}
	for(int re i=1;i<=m2;++i)gi(),b[i]=y2[i]=gi();
	std::sort(b+1,b+m2+1);ct2=std::unique(b+1,b+m2+1)-b-1;
	for(int re i=1;i<=m2;++i){
		y2[i]=std::lower_bound(b+1,b+ct2+1,y2[i])-b;
		nd2[y2[i]].push_back(i);
	}
	for(int re i=1;i<n+m2;++i){
		int t=gi(),u=gi(),v=gi();
		if(t==1){
			e[++ec]=(edge){po[u],p2[v]};
			e2[v].push_back(ec);
		}else {
			e[++ec]=(edge){p2[u],p2[v]};
			if(y2[u]>y2[v])e2[u].push_back(ec);
			else e2[v].push_back(ec);
		}
	}
	ll ans=0,all=(ll)m1*m2;
	for(int re i1=ct1,i2=1,now=0;i1;--i1){
		for(int re u:nd1[i1])
		for(int re j:e1[u]){
			if(LCT::iscon(j,e[j].u))LCT::cut(j,e[j].u);
			if(LCT::iscon(j,e[j].v))LCT::cut(j,e[j].v);
		}
		while(i2<=ct2&&LCT::qsz(po[1])<n){
			for(int re u:nd2[i2])
			for(int re j:e2[u]){
				if(LCT::iscon(e[j].u,e[j].v)){
					int p=LCT::find_e(e[j].u,e[j].v);
					LCT::cut(p,e[p].u);
					LCT::cut(p,e[p].v);
				}
				LCT::link(j,e[j].u);
				LCT::link(j,e[j].v);
			}
			now+=nd2[i2++].size();
		}
		ans+=(ll)nd1[i1].size()*(m2-now);
	}
	if(ans==0)cout<<"0\n";
	else if(ans==all)cout<<"1\n";
	else {
		ll g=std::__gcd(ans,all);
		cout<<ans/g<<"/"<<all/g<<"\n";
	}
	return 0;
}

再来一份没有化边为点的，跑得居然更慢：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	template<typename T>T get(){
		char c;T num;
		while(!isdigit(c=gc()));num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num; 
	}inline int gi(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

#define y1 y_1
#define y2 y_2

cs int N=1e5+7;

int n,m1,m2,tot,ct1,ct2,ec;
struct edge{int u,v;}e[N<<2|3];

int tim[N<<2|3],vl[N*3];
namespace LCT{
	cs int N=::N*3;
#define lc(u) son[u][0]
#define rc(u) son[u][1]
	int son[N][2],fa[N],sum[N],sub[N];
	bool rev[N];int mx[N],pr[N],sf[N];
	inline void rever(int u){rev[u]^=1;
		std::swap(lc(u),rc(u));
		std::swap(pr[u],sf[u]);
	}
	void pushdn(int u){
		if(rev[u]){
			if(lc(u))rever(lc(u));
			if(rc(u))rever(rc(u));
			rev[u]=0;
		}
	}
	void pushup(int u){
		sum[u]=sum[lc(u)]+sum[rc(u)]+vl[u]+sub[u];mx[u]=0;
		(lc(u)&&tim[pr[u]]>tim[mx[u]])&&(mx[u]=pr[u]);
		(rc(u)&&tim[sf[u]]>tim[mx[u]])&&(mx[u]=sf[u]);
		(tim[mx[lc(u)]]>tim[mx[u]])&&(mx[u]=mx[lc(u)]);
		(tim[mx[rc(u)]]>tim[mx[u]])&&(mx[u]=mx[rc(u)]);
	}
	inline bool get(int u){return rc(fa[u])==u;}
	inline bool isrt(int u){return lc(fa[u])!=u&&rc(fa[u])!=u;}
	inline void Rotate(int u){
		int p=fa[u],pp=fa[p],d=get(u);
		if(!isrt(p))son[pp][get(p)]=u;
		son[p][d]=son[u][!d],son[u][!d]=p;
		fa[p]=u,fa[u]=pp,fa[son[p][d]]=p;
		pushup(p);pushup(u);
	}
	void Splay(int u){
		static int q[N],qn;q[qn=1]=u;
		for(int re p=u;!isrt(p);p=fa[p])q[++qn]=fa[p];
		for(int re i=qn;i;--i)pushdn(q[i]);
		for(int re p=fa[u];!isrt(u);Rotate(u),p=fa[u])
		if(!isrt(p))Rotate(get(p)==get(u)?p:u);
	}
	void access(int u){
		for(int re ch=0;u;u=fa[ch=u]){
			Splay(u);if(rc(u))sub[u]+=sum[rc(u)],rc(u)=0;
			if(ch){
				sub[u]-=sum[ch];
				while(pushdn(ch),lc(ch))ch=lc(ch);
				Splay(ch);sf[u]=pr[ch];rc(u)=ch;
			}pushup(u);
		}
	}
	void makert(int u){access(u);Splay(u);rever(u);}
	int find_rt(int u){
		access(u);Splay(u);
		while(pushdn(u),lc(u))u=lc(u);
		Splay(u);return u;
	}
	bool iscon(int u,int v){makert(u);return find_rt(v)==u;}
	void link(int u,int v,int e){
		makert(u),access(v),Splay(v);
		fa[u]=v,sub[v]+=sum[u];pr[u]=sf[v]=e;
		pushup(u),pushup(v);
	}
	void cut(int u,int v){
		makert(u);access(v);Splay(v);
		fa[lc(v)]=0,lc(v)=0,pushup(v),pushup(u);
	}
	int qsz(int u){access(u);Splay(u);return sum[u];}
	int find_e(int u,int v){
		makert(u),access(v),Splay(v);
		return mx[v];
	}
#undef lc
#undef rc
}

bool del[N<<2|3];
int y1[N],y2[N];
int p1[N],p2[N],po[N];
std::vector<int> e1[N],e2[N];
std::vector<int> nd1[N],nd2[N]; 
int b[N];

signed main(){
#ifdef zxyoi
	freopen("mt.in","r",stdin);
#endif
	n=gi(),m1=gi(),m2=gi();tot=0;
	for(int re i=1;i<=m1;++i)p1[i]=++tot;
	for(int re i=1;i<=m2;++i)p2[i]=++tot;
	for(int re i=1;i<=n;++i)po[i]=++tot,vl[tot]=1;
	for(int re i=1;i<=tot;++i)LCT::pushup(i);
	for(int re i=1;i<=m1;++i)gi(),b[i]=y1[i]=gi();
	std::sort(b+1,b+m1+1);ct1=std::unique(b+1,b+m1+1)-b-1;
	for(int re i=1;i<=m1;++i){
		y1[i]=std::lower_bound(b+1,b+ct1+1,y1[i])-b;
		nd1[y1[i]].push_back(i); 
	}
	for(int re i=1;i<n+m1;++i){
		int t=gi(),u=gi(),v=gi();
		if(t==1){
			e[++ec]=(edge){po[u],p1[v]};
			e1[v].push_back(ec);tim[ec]=y1[v];
			LCT::link(po[u],p1[v],ec);
		}else {
			e[++ec]=(edge){p1[u],p1[v]};
			if(y1[u]>y1[v])e1[u].push_back(ec),tim[ec]=y1[u];
			else e1[v].push_back(ec),tim[ec]=y1[v];
			LCT::link(p1[u],p1[v],ec);
		}
	}
	for(int re i=1;i<=m2;++i)gi(),b[i]=y2[i]=gi();
	std::sort(b+1,b+m2+1);ct2=std::unique(b+1,b+m2+1)-b-1;
	for(int re i=1;i<=m2;++i){
		y2[i]=std::lower_bound(b+1,b+ct2+1,y2[i])-b;
		nd2[y2[i]].push_back(i);
	}
	for(int re i=1;i<n+m2;++i){
		int t=gi(),u=gi(),v=gi();
		if(t==1){
			e[++ec]=(edge){po[u],p2[v]};
			e2[v].push_back(ec);
		}else {
			e[++ec]=(edge){p2[u],p2[v]};
			if(y2[u]>y2[v])e2[u].push_back(ec);
			else e2[v].push_back(ec);
		}
	}
	ll ans=0,all=(ll)m1*m2;
	for(int re i1=ct1,i2=1,now=0;i1;--i1){
		for(int re u:nd1[i1])
		for(int re j:e1[u]){
			if(!del[j])LCT::cut(e[j].u,e[j].v);
		}
		while(i2<=ct2&&LCT::qsz(po[1])<n){
			for(int re u:nd2[i2])
			for(int re j:e2[u]){
				if(LCT::iscon(e[j].u,e[j].v)){
					int p=LCT::find_e(e[j].u,e[j].v);
					LCT::cut(e[p].u,e[p].v);
					del[p]=true;
				}
				LCT::link(e[j].v,e[j].u,j);
			}
			now+=nd2[i2++].size();
		}
		ans+=(ll)nd1[i1].size()*(m2-now);
	}
	if(ans==0)cout<<"0\n";
	else if(ans==all)cout<<"1\n";
	else {
		ll g=std::__gcd(ans,all);
		cout<<ans/g<<"/"<<all/g<<"\n";
	}
	return 0;
}
zxyoi_dreamer
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
动态树LCT BySemiWaker
03-05
动态树LCT
LCT Unlock MTK V1.8 CN
04-08
《LCT Unlock MTK V1.8 CN：MTK芯片解锁工具详解》 在移动通信技术领域，MTK（MediaTek）是一家知名的芯片制造商，其生产的芯片广泛应用于各种智能手机和平板电脑。LCT Unlock MTK V1.8 CN是专为MTK芯片设计的一款...
[校内模拟题1]
weixin_30338497的博客
07-25 114
竞赛时间：????年??月??日??:??-??:?? 题目名称 算 游 运 名称 sum tour lucky 输入 sum.in tour.in lucky.in 输出 sum.out tour.out lucky.out 每个测试点时限 1 秒 1 秒 1 秒 内存限制 512MB 512MB 512MB 测试点数目 10 10 10 每个测试点分值 1...
NOIP模拟题 斐波那契数列
weixin_30820077的博客
11-03 147
题目大意 给定长度为$n$序列$A$，将它划分成尽可能少的若干部分，使得任意部分内两两之和均不为斐波那契数列中的某一项。 题解 不难发现$2\times 10^9$之内的斐波那契数不超过$50$个 先求出第$i$个数之前最后一个能和第$i$个数相加为斐波那契数的位置$last_i$。 考虑每一部分$[l,r]$只需满足$\max\{last_i\}<l(i\in [l,r...
[校内模拟]喜欢最最痛
时光真疯狂, 我一路执迷于匆忙.
06-27 489
Description 有一棵n个点的树，边有边权。 有m条额外边，第i条的边权为ai，对于每个i∈[0,m]，问加入前i条额外边后，从1出发经过所有树边至少一次最后回到1的最短路径的长度。 n,m<=10^5 Solution 容易发现这题分为两部分，先求出an[i]表示选择树上i条不相交的链的长度最大值，然后对于每个前缀i，求一个j，使得an[j]-前i个a的前j小最大 an的话可以考虑...
LCT
andyc_03的博客
04-25 111
【算法简介】 又是一个大数据结构！ 简单的说，LCT是维护了一个森林，每个树都是一棵splay，由虚边连接着各个树 支持操作： 查询、修改链上的信息（最值，总和等） 随意指定原树的根（即换根） 动态连边、删边 动态维护连通性 性质： 每一个Splay维护的是一条从上到下按在原树中深度严格递增的路径，且中序遍历Splay得到的每个点的深度序列严格递增。 边分为实边和虚边，实边包含在Splay中，而虚边总是由一棵Splay指向另一个节点（指向该Splay中中序遍历最靠前的点在原树中的父亲）。
NOIP模拟赛8.16 C、长跑（LCT）
cqbzcsq的博客
08-16 187
长跑 内存限制：128 MiB 时间限制：1000 ms 【题目描述】 某校开展了同学们喜闻乐见的阳光长跑活动。为了能“为祖国健康工作五十 年”，同学们纷纷离开寝室，离开教室，离开实验室，到操场参加 3000 米长跑运 动。一时间操场上熙熙攘攘，摩肩接踵，盛况空前。 为了让同学们更好地监督自己，学校推行了刷卡机制。 学校中有 n 个地点，用 1 到 n 的整数表示，每个地点设有若干个刷卡机...
计算机的术语中lcT,技术|模拟 Altair 8800 计算机
weixin_42497828的博客
07-08 679
Altair 8800 是 1975 年发布的自建家用电脑套件。Altair 基本上是第一台个人电脑(PC)，虽然 PC 这个名词好几年前就出现了。对 Dell、HP 或者 Macbook 而言它是亚当(或者夏娃)。有些人认为为 Z80(与 Altair 的 Intel 8080 密切相关的处理器)编写仿真器真是太棒了，并认为它需要一个模拟 Altair 的控制面板。所以如果你想知道 1975 年...
【校内模拟】子串（SAM）（LCT）（主席树）
zxyoi_dreamer的博客（不定期诈尸）
09-30 263
题解： 最后实现注意一下字符串边界，LCT稍微考虑一下可以不用写cut。 代码： #include<bits/stdc++.h> #define ll long long #define re register #define cs const namespace IO{ cs int Rlen=1<<22|1; char buf[Rlen],*p1,*p2; i...
校内模拟赛 旅行(by NiroBC)
weixin_30369041的博客
04-02 93
题意： 　　n个点的无向图，Q次操作，每次操作可以连接增加一条边，询问两个点之间有多少条边是必经之路。如果不连通，输出-1。 分析： 　　首先并查集维护连通性，每次加入一条边后，如果这条边将会连接两个联通块，那么lct连接两个点，边权化为点权，新增一个点，点权为1。否则，构成了环，环上的边都变为0，lct维护覆盖标记。询问就是对一条链进行询问。 　　离线+树剖的做法：从前往后建出树，如果出...
动态树（LCT）
最新发布
DY-hzx的博客
02-14 615
关于LCT的介绍
lct.zip_LCT
09-19
【标题】"lct.zip_LCT" 涉及到的是一个名为 "lct.zip" 的压缩文件，其中包含了与"LCT"相关的资料。LCT通常代表“逻辑控制技术”（Logical Control Technology）或者“线性电流变换器”（Linear Current Transformer...
lct_LCT_
09-29
标题中的“lct_LCT_”很可能代表“Link-Cut Tree”，这是一项在计算机科学领域，特别是算法设计中广泛使用的数据结构。Link-Cut Tree（链接切割树）是一种自顶向下的动态树数据结构，主要用于解决与树相关的操作，如...
lct.zip_lct文件
09-24
标题中的“lct.zip_lct文件”指的是一个压缩包，其中包含了名为“lct”的文件。这个“lct”文件可能是一个库文件，用于在特定的编程环境或项目中提供预定义的功能，使得开发者能够更方便地进行软件开发，而无需从零...
【SCOI2018】Tree（LCT）
zxyoi_dreamer的博客（不定期诈尸）
04-05 517
描述 在大小为 N 的树上，点从 1 到 N 标号，第 i 个点有权值 Ai，现在需要支持两 种操作： 第一种操作格式为“1 U” ，表示询问从 U 出发的简单路径，经过的点权值 之和的最大值； 第二种操作格式为“2 U V” ，表示将 U 的权值修改为 V。 输入 第一行两个整数 N 和 M，表示树的大小和操作数； 第二行 N-1 个整数，第 i 个整数 Pi(1<=Pi<=i)表示...
2018.12.12【BZOJ5192】【洛谷P4271】New Barns（LCT维护直径）
zxyoi_dreamer的博客（不定期诈尸）
12-12 415
BZOJ传送门 洛谷传送门 解析： 许久没有写LCT，手都有点生了。 思路： 这道题是显然的需要维护直径的，因为有结论，树上距离一个点的最远点一定是某个直径的端点。证明很显然，就不证了。 那么怎么维护动态的直径呢？ 我们可以用LCT维护一个森林，对于每棵树记录直径的两个端点，加入一个新点的时候更新，我们有如下结论：新的直径与原来直径至少有一个端点重合。 那么LCT乱搞一波就做完了。 代码： #...
【模板】LCT求最小生成树
zxyoi_dreamer的博客（不定期诈尸）
10-07 371
参考题目：洛谷P3366 解析： 联赛完后统一更模板题解析。 代码： #include&lt;bits/stdc++.h&gt; using namespace std; #define ll long long #define re register #define gc getchar #define pc putchar #define cs const inline int geti...
【校内模拟】小B的图（贪心）（MST）（LCT）
zxyoi_dreamer的博客（不定期诈尸）
02-11 351
简要题意： 给一张图，有两种边，第一种的权值为 v+xv+xv+x，其中 xxx 是变量，第二种的权值为 v−xv-xv−x。 现在给出若干 xxx 的值，请你求出MST的大小。 分别只保留两种边的图仍然连通。 题解： 由于本身两种图是连通的，所以在新MST上的边只可能是原来两种边的MST中的边。 将 v+xv+xv+x 拿来做MST，建立LCT。考虑用 v−xv-xv−x 来替换。 注意替换顺...
动态树问题：LCT与Splay应用
其中，Link-Cut-Tree (LCT) 是解决这类问题的关键技术，由传奇算法学家Tarjan提出，用于处理树的动态更新，而不同于一般的动态树概念。 LCT 建立在Splay树的基础上，但并非简单的套用Splay树的模板。Splay树本身是...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值