简要题意:
给平面上三个点集,其中第一个点集 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) E1⊆V1×(V1∪V2), E 2 ⊆ V 3 × ( V 2 ∪ V 3 ) E_2\subseteq V_3\times (V_2\cup V_3) E2⊆V3×(V2∪V3)," × \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;
}