[cf1215f]Radio Stations(2-SAT)

题意就是有p个点,可以选取或不取。有n对关系,满足两者至少取一个。又有m对关系,满足两个不能同时取。同时每个点有区间[li,ri],最终要满足存在这样一个数属于所有取到点的区间,也就是说他们的交集不是空集,并且已知li、ri∈[1,M]。

如果不看最后区间的限制条件的话,就是裸的2-SAT了。对于n对关系,将不取a向取b连边、不取b向取a连边;对于m对关系,取a向不取b连边、取b向不取a连边。

然后我们考虑如何通过连边描述每个点的区间。一个点的左端点所限制的其实是:如果最后选取的数小于li,那么是不能取这个点的;如果取了这个点,那么最终的数一定大于等于li。右端点同理。由于左右断点范围已知且不大,所以我们再建立2*(M+2)个点(从0到M+1),分别表示最终选择的数大于等于x和小于等于x。当li==x时,取i向大于等于x连边、小于等于x-1向不取i连边。对于ri同理,大家不妨自己推导一下。同时小于等于x向小于等于x+1连边、大于等于x+1向大于等于x连边。而且需要注意的是,对于所有端点的l和r应分别做排序。

其他的大概也就没什么了。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=400010,NN=1600010;
struct edge{
	int y,next;
}data[10*N];
struct node{
	int x,id;
}l1[N],r1[N];
int n,p,M,m,flag,n1,n2,dfn[NN],low[NN],sta[NN],top,f[NN],h[NN],num,num1,num2,num3,a1[N],l,r;
bool insta[NN];
inline int read(){
	int x=0;char ch=getchar();
	while(ch<'0'||ch>'9')ch=getchar();
	while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
	return x;
}
inline void addedge(int u,int v){
	data[++num].y=v,data[num].next=h[u],h[u]=num;
}
inline int min1(int i,int j){return i<j?i:j;}
inline bool cmp1(node i,node j){return i.x<j.x;}
inline bool cmp2(node i,node j){return i.x>j.x;}
inline bool cmp3(node i,node j){return i.id<j.id;}
void tarjan(int u){
	dfn[u]=low[u]=++num1,sta[++top]=u,insta[u]=true;
	for(int i=h[u],v;i!=-1;i=data[i].next){
		v=data[i].y;
		if(!dfn[v]){
			tarjan(v);low[u]=min1(low[v],low[u]);
		}else if(insta[v]&&dfn[v]<low[u])low[u]=dfn[v];
	}
	if(dfn[u]==low[u]){
		++num2;
		while(dfn[sta[top]]!=low[sta[top]]){
			f[sta[top]]=num2;insta[sta[top--]]=false;
		}
		f[u]=num2;insta[sta[top--]]=false;
	}
}
int main(){
	n=read(),p=read(),M=read(),m=read();
	memset(h,-1,sizeof h);num=0;
	for(int i=1,x,y;i<=n;++i){
		x=read(),y=read();
		addedge(x<<1|1,y<<1),addedge(y<<1|1,x<<1);
	}
	for(int i=1;i<=p;++i){
		l1[i].x=read(),r1[i].x=read();
		l1[i].id=r1[i].id=i;
	}
	n1=p<<1|1;n2=2+M;
	sort(l1+1,l1+p+1,cmp2);sort(r1+1,r1+p+1,cmp1);
	l=r=1;
	for(int i=1;i<=M;++i){
		while(r<=p&&r1[r].x==i){addedge(r1[r].id<<1,n1+n2+i+1);addedge(n1+i+2,r1[r].id<<1|1);++r;}
		addedge(n1+i+1,n1+i);addedge(n1+n2+i,n1+n2+i+1);
	}
	addedge(n1+M+2,n1+M+1);addedge(n1+n2+M+1,n1+n2+M+2);
	for(int i=M;i>=1;--i){
		while(l<=p&&l1[l].x==i){addedge(l1[l].id<<1,n1+i+1);addedge(n1+n2+i,l1[l].id<<1|1);++l;}
	}
	for(int i=1,x,y;i<=m;++i){
		x=read(),y=read();
		addedge(x<<1,y<<1|1),addedge(y<<1,x<<1|1);
	}
	sort(r1+1,r1+p+1,cmp3);
	memset(dfn,0,sizeof dfn);
	for(int i=2;i<=n1+2*n2;++i)if(!dfn[i])tarjan(i);
	flag=1;
	for(int i=1;i<=p;++i)if(f[i<<1]==f[i<<1|1]){flag=0;break;}
	if(!flag){printf("-1");return 0;}
	r=M,num3=0;
	for(int i=1;i<=p;++i)
		if(f[i<<1]<f[i<<1|1]){r=min1(r,r1[i].x);a1[++num3]=i;}
	printf("%d %d\n",num3,r);
	for(int i=1;i<=num3;++i)printf("%d ",a1[i]);
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值