【Codeforces 1215F】Radio Stations（2-SAT）

最新推荐文章于 2022-07-16 11:39:36 发布

xzyxzy?_?

最新推荐文章于 2022-07-16 11:39:36 发布

阅读量591

点赞数 3

分类专栏：算法康复训练文章标签： Codeforces 2-SAT

本文链接：https://blog.csdn.net/weixin_44231553/article/details/100927625

版权

算法康复训练专栏收录该内容

17 篇文章 1 订阅

订阅专栏

题目链接

复习一下2-SAT姿势（退役那年day2T1就是2-SAT，好像我挂了10分…

2-SAT问题：每个变量取值0/1，求一组变量取值满足一些约束关系。做法是每个变量 $x$ 开两个点 $x 1$ 和 $x 0$ ，表示变量 $x$ 选0/1。根据约束关系连边，有向边 $(a, b)$ 表示满足了 $a$ 就必须满足 $b$ 。跑一遍Tarjan，如果存在 $x$ ，使得 $x 0$ 和 $x 1$ 属于同一个强联通分量，则无解。否则对每个变量 $x$ ， $x 0$ 所属scc编号小则取0，否则取1，就是一组合法解。

关系举例：
$x$ 和 $y$ 至少选一个： $(x 0, y 1), (y 0, x 1)$
$x$ 和 $y$ 不能同时选： $(x 1, y 0), (y 1, x 0)$
选了 $x$ 就必须选 $y$ ： $(x 1, y 1), (y 0, x 0)$
总之对于 $x, y$ 的二元关系，对 $x 0, y 0, x 1, y 1$ 四个变量逐一考虑，看选了一个之后是不是有另一个必须选，有就连边。单变量的话，比如 $x$ 必须选，那就是 $(x 0, x 1)$ 。

【P.S：有向图tarjan里做low[x]=min(low[x],dfn[y])需要判断y在栈内…无向图不用…我一开始忘记了qwq】

回到这道题，若没有 $f$ 的限制就会是一个裸的2-SAT。考虑 $f$ 的限制，我们添加 $M$ 个变量，第 $i$ 个变量表示 $[f > = i]$ 。那么对于一个基站 $x$ 及其频率 $[L, R]$ ，选择 $x$ 就必须 $[f > = L]$ 取1且 $[f > = (R + 1)]$ 取0；同时对于这 $M$ 个变量之间也存在限制，若 $[f > = L + 1]$ 取1，则 $[f > = L]$ 必须取1。这个限制关系会传递下去，所以对于 $[f > = L + 1]$ 不需要考虑所有 $[f > = w] ， w < = L$ ，只需要考虑与 $[f > = L]$ 之间的限制关系即可。总复杂度仍然保持线性。

#define no(x) ((x)+n)
#define val(v) (sn+(v))
const int N=2000005;
int A,sn,sp,B,n;
vector<int>S[N];
void E(int x,int y) {S[x].push_back(y);}
int num[N],low[N],scc[N],st[N],cnt,top,inc;bool in[N];
void DFS(int x)
{
  num[x]=low[x]=++cnt;
  st[++top]=x;in[x]=1;
  for(auto &y:S[x])
    {
      if(!num[y]) DFS(y),low[x]=min(low[x],low[y]);
      else if(in[y]) low[x]=min(low[x],num[y]);
    }
  if(num[x]==low[x])
    {
      scc[x]=++inc;
      while(st[top]!=x) scc[st[top]]=inc,in[st[top]]=0,top--;
      top--;in[x]=0;
    }
}
int ans[N],u,f;
void _SAT( )
{
  int i,x,y,tot=n*2;
  for(i=1;i<=tot;i++)
    if(!num[i]) DFS(i);
  for(i=1;i<=sn;i++)
    {
      x=scc[i];y=scc[no(i)];
      if(x==y) {puts("-1");return;}
      if(x<y) ans[++u]=i;
    }
  for(i=1;i<=sp;i++)
    {
      x=scc[val(i)];y=scc[no(val(i))];
      if(x==y) {puts("-1");return;}
      if(x<y) f=i;
    }
  printf("%d %d\n",u,f);
  for(i=1;i<=u;i++) printf("%d ",ans[i]);
}
int main( )
{
  read(A);read(sn);read(sp);read(B);n=sn+sp;
  int i,j,x,y,l,r;
  for(i=1;i<=A;i++)
    {
      read(x),read(y);
      E(no(x),y),E(no(y),x); //x和y必须选至少一个
    }
  for(x=1;x<=sn;x++)
    {
      read(l),read(r),r++;
      E(x,val(l)),E(no(val(l)),no(x));  //选x必须选f>=l
      if(r<=sp) E(x,no(val(r))),E(val(r),no(x)); //选x则不能选f>=(r+1)
    }
  for(i=1;i<=B;i++)
    {
      read(x);read(y);
      E(x,no(y));E(y,no(x));  //x和y不能同时选
    }
  for(i=1;i<sp;i++)
    {
      l=val(i);r=val(i+1);
      E(r,l);E(no(l),no(r));  //选f>=l+1必须选f>=l
    }
  _SAT( );
  return 0;
}