【CF1215F】 Radio Stations

题目

比较精妙的\(\text{2-sat}\)建图了

还是按照套路把每个电台拆成\((0/1,i)\)表示不选/选

前两种连边是板子就不解释了

考虑如何限制选择一个唯一的\(f\)，并且还能限制不选\(f\notin [l_i,r_i]\)的电台

考虑前缀优化建图，我们建\((0/1,i)\)表示在\([0,i]\)中不选/选某一个频率

对于这\(2\times m\)个点有一些显然的连边

\[(1,i)->(1,i+1) \]

即\([0,i]\)内启用一个频率则\([0,i+1]\)中必启用一个频率

\[(0,i+1)->(0,i) \]

即\([0,i+1]\)没有启用任何一个频率则\([0,i]\)必不会启用任何一个频率

注意到我们必须启用一个频率，所以连\((0,m)->(1,m)\)，这样能强制使得\(m\)的拓扑序更靠后

之后对于一个电台\(i\)，我们用\([0,l_i-1]\)和\([0,r_i]\)这两个前缀的关系来讨论一下

当\([0,l_i-1]\)内选择了一个频率时，电台\(i\)就没有办法被启用了，于是连\((1,l_i-1)->(0,i)\)，以及对称边\((1,i)->(0,l_i-1)\)

当\([0,r_i]\)中没有任何一个频率被启用时，电台\(i\)也没有办法被启用，于是连\((0,r_i)->(0,i)\)，以及对称边\((1,i)->(1,r_i)\)

根据我们连得两条对称边，我们发现当启用电台\(i\)的时候能推出\([0,l_i-1]\)没有频率被启用，\([0,r_i]\)有一个频率被启用，即被启用的频率在\([l_i,r_i]\)内，这符合题目要求

代码

#include<bits/stdc++.h>
#define re register
#define LL long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
inline int read() {
	char c=getchar();int x=0;while(c<'0'||c>'9') c=getchar();
	while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
const int maxn=2e6+1;
struct E{int v,nxt;}e[maxn*10];
int dfn[maxn],low[maxn],st[maxn],f[maxn],col[maxn],head[maxn],id[2][maxn>>1],l[maxn>>2],r[maxn>>2];
int n,m,A,B,num,__,cnt,top,col_num,ans,Ans[500005];
inline void add(int x,int y) {e[++num].v=y;e[num].nxt=head[x];head[x]=num;}
void tarjan(int x) {
	dfn[x]=low[x]=++__,st[++top]=x,f[x]=1;
	for(re int i=head[x];i;i=e[i].nxt) 
	if(!dfn[e[i].v]) tarjan(e[i].v),low[x]=min(low[x],low[e[i].v]);
	else if(f[e[i].v]) low[x]=min(low[x],dfn[e[i].v]);
	if(dfn[x]==low[x]) {
		col_num++;int now;
		do{now=st[top--];f[now]=0;col[now]=col_num;}while(x!=now);
	}
}
int main() {
	A=read(),n=read(),m=read(),B=read();
	for(re int i=1;i<=n;i++) id[0][i]=++cnt,id[1][i]=++cnt;
	for(re int i=0;i<=m;i++) id[0][i+1+n]=++cnt,id[1][i+1+n]=++cnt;
	for(re int x,y,i=1;i<=A;i++) {
		x=read(),y=read();
		add(id[0][x],id[1][y]);
		add(id[0][y],id[1][x]);
	}
	for(re int i=1;i<=n;i++) l[i]=read(),r[i]=read();
	for(re int x,y,i=1;i<=B;i++) {
		x=read(),y=read();
		add(id[1][x],id[0][y]);
		add(id[1][y],id[0][x]);
	}
	add(id[1][n+1],id[0][n+1]);
	for(re int i=1;i<=m;i++) {
		int x=i+1+n;
		add(id[0][x],id[0][x-1]);
		add(id[1][x-1],id[1][x]);
	}
	add(id[0][n+m+1],id[1][n+m+1]);
	for(re int i=1;i<=n;i++) {
		int li=l[i]+n,ri=r[i]+n+1;
		add(id[1][li],id[0][i]);
		add(id[1][i],id[0][li]);
		add(id[0][ri],id[0][i]);
		add(id[1][i],id[1][ri]);
	}
	for(re int i=1;i<=cnt;++i) if(!dfn[i]) tarjan(i);
	for(re int i=1;i<=n+1+m;i++) if(col[id[0][i]]==col[id[1][i]]) {puts("-1");return 0;}
	for(re int i=1;i<=n;i++) 
		if(col[id[1][i]]<col[id[0][i]]) Ans[++ans]=i;
	for(re int i=1;i<=m;i++) if(col[id[1][i+n+1]]<col[id[0][i+n+1]]) {
		printf("%d %d\n",ans,i);
		for(re int j=1;j<=ans;j++) printf("%d ",Ans[j]);
		return 0;
	}
}