Codeforces 870E Points, Lines and Ready-made Titles 并查集

题意

给定平面直角坐标系上的 n n n个点。对于每个点,你可以画过这个点的一条横线,或画过这个点的一条竖线,也可以什么都不做。问有多少种不同的结果。答案模 1 0 9 + 7 10^9+7 109+7
n ≤ 100000 , ∣ x i ∣ , ∣ y i ∣ ≤ 1 0 9 n\le 100000,|xi|,|yi|\le 10^9 n100000,xi,yi109

分析

我们把每一行和每一列分别看做一个点,对于平面上的一个点(x,y),在行点x和列点y之间连一条边。
对于每个联通块,若为树,则对答案的贡献为 2 n − 1 2^n-1 2n1否则为 2 n 2^n 2n,n为连通块大小。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int N=100005;
const int MOD=1000000007;

int n,x[N],y[N],w[N],f[N*2],size[N*2],val[N*2];

int read()
{
	int x=0,f=1;char ch=getchar();
	while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}

int ksm(int x,int y)
{
	int ans=1;
	while (y)
	{
		if (y&1) ans=(LL)ans*x%MOD;
		x=(LL)x*x%MOD;y>>=1;
	}
	return ans;
}

int find(int x)
{
	if (f[x]==x) return x;
	else return f[x]=find(f[x]);
}

int main()
{
	n=read();
	for (int i=1;i<=n;i++) x[i]=read(),y[i]=read(),w[i]=x[i];
	sort(w+1,w+n+1);
	int w1=unique(w+1,w+n+1)-w-1;
	for (int i=1;i<=n;i++) x[i]=lower_bound(w+1,w+w1+1,x[i])-w;
	for (int i=1;i<=n;i++) w[i]=y[i];
	sort(w+1,w+n+1);
	w1=unique(w+1,w+n+1)-w-1;
	for (int i=1;i<=n;i++) y[i]=lower_bound(w+1,w+w1+1,y[i])-w;
	for (int i=1;i<=n*2;i++) f[i]=i,size[i]=1;
	for (int i=1;i<=n;i++)
	{
		int p=find(x[i]),q=find(y[i]+n);
		if (p!=q) f[p]=q,val[q]+=val[p]+1,size[q]+=size[p];
		else val[p]++;
	}
	int ans=1;
	for (int i=1;i<=n*2;i++)
		if (f[i]==i) ans=(LL)ans*(ksm(2,size[i])-(val[i]<size[i]))%MOD;
	printf("%d",ans);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值