【校内模拟】坦克(平面图生成树)(并查集)(双指针)

简要题意:

给平面上三个点集,其中第一个点集 V 1 V_1 V1中的点全部在第一象限,第二个点集 V 2 V_2 V2中的点全部在 x x x 轴上,第三个点集 V 3 V_3 V3中的点全部在第四象限。

然后给两个无向边集,其中 E 1 ⊆ V 1 × ( V 1 ∪ V 2 ) E_1\subseteq V_1\times (V_1\cup V_2) E1V1×(V1V2) E 2 ⊆ V 3 × ( V 2 ∪ V 3 ) E_2\subseteq V_3\times (V_2\cup V_3) E2V3×(V2V3)," × \times ×"表示集合的笛卡尔积,且上述两个边集分别成树。

保证所有边和点形成的图是平面图

等概率在 V 1 V_1 V1中选择一个点,把所有纵坐标不小于这个点的点全部删除,等概率在 V 3 V_3 V3中选择一个点,把所有纵坐标不大于这个点的点全部删除。

询问此时剩余点的导出子图中,所有 V 2 V_2 V2中的点连通的概率。

题解:

其实双指针的思路还是挺显然的,问题在于如何维护连通性。

可以化边为点,边上维护删除时间然后上LCT。

由于要维护大量的边,常数十分窒息

题目给了一个性质,就是说这张图保证是一个平面图

容易想到平面图最短路转对偶图最小割

注意到极左点和极右点连通的时候,相当于一个从极上到极下的割。

好其实这些都还不是最重要的。

我们考虑按照离 x x x轴的距离加点的时候对连通性的影响。

由于这是一个平面图,这个影响可以表示为区间之间的包含关系。那么删掉一个点会造成的影响就是 i i i i + 1 i+1 i+1不连通。而加一个点的影响也可以描述为 i i i i + 1 i+1 i+1从不连通变为连通。

对于上下两个点集用带权并查集算一下中间点集区间连通的时候的 y y y坐标,然后双指针就完事了。


代码:

#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;

cs int N=1e5+7;
#define y1 y_1
#define y2 y_2
int n,m1,m2;
int b1[N],b2[N];
int y1[N],y2[N];

int fa[N<<1|1],L[N<<1|1];int id[N];
std::vector<int> G[N<<1|1];

int gf(int u){while(u!=fa[u])u=fa[u]=fa[fa[u]];return u;}
void get_area(int m,int *bel,int *y){
	for(int re i=1;i<=n+m;++i)G[i].clear(),fa[i]=i;
	for(int re i=1;i<=m;++i)gi(),y[i]=gi();
	for(int re i=1;i<n+m;++i){
		int t=gi(),u=gi(),v=gi();
		if(t==1)G[v+n].push_back(u);
		else if(y[u]>y[v]||(y[u]==y[v]&&v<u))G[u+n].push_back(v+n);
		else G[v+n].push_back(u+n);
	}
	for(int re i=1;i<=n;++i)L[i]=i;
	for(int re i=n+1;i<=n+m;++i)L[i]=0;
	for(int re i=1;i<=m1;++i)id[i]=i+n;
	std::sort(id+1,id+m1+1,[&](int i,int j){
		i-=n,j-=n;return y[i]<y[j]||(y[i]==y[j]&&i<j);
	});
	for(int re i=1;i<=m1;++i){
		int u=id[i];std::vector<int> cc;
		for(int re v:G[u]){int p=gf(v);cc.push_back(p);fa[p]=u;}
		std::sort(cc.begin(),cc.end(),[](int i,int j){
			return L[i]<L[j];
		});
		for(int re i=1;i<cc.size();++i)bel[L[cc[i]]]=y[u-n];
		L[u]=L[cc.front()];
	}
}

signed main(){
#ifdef zxyoi
	freopen("mt.in","r",stdin);
#endif
	n=gi(),m1=gi(),m2=gi();
	get_area(m1,b1,y1);
	get_area(m2,b2,y2);
	std::sort(y1+1,y1+m1+1,std::greater<int>());
	std::sort(y2+1,y2+m2+1);
	for(int re i=1;i<n;++i)id[i]=i+1;
	std::sort(id+1,id+n,[](int i,int j){return b1[i]>b1[j];});
	ll ans=0,all=(ll)m1*m2;int ban=0;
	for(int re i=1,j=1,k=0;i<=m1;++i){
		while(j<n&&b1[id[j]]>=y1[i])ban=std::max(ban,b2[id[j++]]);
		while(k<m2&&ban>=y2[k+1])++k;
		ans+=m2-k;
	}
	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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值