UOJ #356. 【JOI2017春季合宿】Port Facility 二分图 优化建边

考虑两个货物 i,j(不妨令Ai<Aj) 当且仅当 Bi<Bj 时其二者不能同船

其他情况都不会互相影响

所以这就转化成了一个二分图模型

我们在这样的 i,j 之间连边 若无奇环 则最后的贡献为 2^(联通块数)

但是边数是 n^2 级别的 考虑优化建边

用主席树那种没想明白 还是委屈就全用了set

具体操作如下

引入0边 表示边两侧的点颜色相同 1边表示不同

这样的话 若要连接 i->a,i->b...i->z

就可以变成

1边：i->a

0边：a->b..y->z

用set维护 每次直接取出 [Aj,Bj] 之间的 Bi

利用边的性质 我们可以删去两侧都有连边的点

均摊 O(nlogn)

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<set>
#include<map>
using namespace std;

inline 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;
}
inline void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=2001000,inf=0X3f3f3f3f,mod=int(1e9)+7;

int b[N];

int last[N],ecnt;
struct EDGE{int to,nt,val;}e[N<<3];
inline void readd(int u,int v,int val)
{e[++ecnt]=(EDGE){v,last[u],val};last[u]=ecnt;}
inline void add(int u,int v,int val)
{u=b[u],v=b[v];readd(u,v,val);readd(v,u,val);}

set<int> S,T;

int a[N];
int L[N],R[N],st[N],top;
int clr[N];

void dfs(int u,int col)
{
	if(~clr[u])
	{
		if(clr[u]!=col){puts("0");exit(0);}
		return ;
	}
	clr[u]=col;
	for(int i=last[u];i;i=e[i].nt)
		dfs(e[i].to,col^e[i].val);
}

int main()
{
	int n=read();
	register int i,k,x,y;
	for(i=1;i<=n;++i)
		x=read(),y=read(),a[x]=y,b[y]=i;
	
	register set<int>::iterator it;
	S.insert(inf);T.insert(inf);
	S.insert(-inf);T.insert(-inf);
	for(i=1;i<=(n<<1);++i)
		if(y=a[i])
		{
			it=S.upper_bound(y);
			R[y]=*it;--it;L[y]=*it;
			if(*it>i) add(y,*it,1);
			it=T.upper_bound(y);--it;
			for (;*it>i;--it)
			{
				k=*it;
				if(L[k]>i) add(k,L[k],0),L[k]=-inf;
				if(R[k]<y) add(k,R[k],0),R[k]=inf;
				if(L[k]==-inf && R[k]==inf) st[++top]=k;
			}
			while(top) T.erase(st[top]),top--;
			S.insert(y);T.insert(y);
		}
	int ans(1);
	memset(clr,-1,sizeof(clr));
	for(i=1;i<=n;++i)
		if(clr[i]==-1)
			dfs(i,0),(ans<<=1)%=mod;
	cout<<ans<<endl;
	return 0;
}