码迷,mamicode.com
首页 > 编程语言 > 详细

算法复习——2—sat(bzoj2199)

时间:2017-09-05 10:03:58      阅读:288      评论:0      收藏:0      [点我收藏+]

标签:script   false   line   main   包含   ring   div   tar   lib   

题目:

Description

由于对Farmer John的领导感到极其不悦,奶牛们退出了农场,组建了奶牛议会。议会以“每头牛 都可以获得自己想要的”为原则,建立了下面的投票系统: M只到场的奶牛 (1 <= M <= 4000) 会给N个议案投票(1 <= N <= 1,000) 。每只 奶牛会对恰好两个议案 B_i and C_i (1 <= B_i <= N; 1 <= C_i <= N)投 出“是”或“否”(输入文件中的‘Y‘和‘N‘)。他们的投票结果分别为VB_i (VB_i in {‘Y‘, ‘N‘}) and VC_i (VC_i in {‘Y‘, ‘N‘})。 最后,议案会以如下的方式决定:每只奶牛投出的两票中至少有一票和最终结果相符合。 例如Bessie给议案1投了赞成‘Y‘,给议案2投了反对‘N‘,那么在任何合法的议案通过 方案中,必须满足议案1必须是‘Y‘或者议案2必须是‘N‘(或者同时满足)。 给出每只奶牛的投票,你的工作是确定哪些议案可以通过,哪些不能。如果不存在这样一个方案, 输出"IMPOSSIBLE"。如果至少有一个解,输出: Y 如果在每个解中,这个议案都必须通过 N 如果在每个解中,这个议案都必须驳回 ? 如果有的解这个议案可以通过,有的解中这个议案会被驳回 考虑如下的投票集合: - - - - - 议案 - - - - - 1 2 3 奶牛 1 YES NO 奶牛 2 NO NO 奶牛 3 YES YES 奶牛 4 YES YES 下面是两个可能的解: * 议案 1 通过(满足奶牛1,3,4) * 议案 2 驳回(满足奶牛2) * 议案 3 可以通过也可以驳回(这就是有两个解的原因) 事实上,上面的问题也只有两个解。所以,输出的答案如下: YN?

Input

* 第1行:两个空格隔开的整数:N和M * 第2到M+1行:第i+1行描述第i只奶牛的投票方案:B_i, VB_i, C_i, VC_i

Output

* 第1行:一个含有N个字符的串,第i个字符要么是‘Y‘(第i个议案必须通过),或者是‘N‘ (第i个议案必须驳回),或者是‘?‘。 如果无解,输出"IMPOSSIBLE"。

Sample Input


3 4
1 Y 2 N
1 N 2 N
1 Y 3 Y
1 Y 2 Y


Sample Output

YN?

HINT

 

Source

题解:

2-sat问题的模板题目,先说2-sat问题的基本解法:

一些问题可以转成布尔方程来求解····

我们的目的是将其布尔方程的每个文字拆开成两点,一点表示其本身,一点表示它的非,比如a就拆成a与┐a,并且将各种运算符号转化为只含有^(与)和->(A->B表示A为真则B为真)的形式,比如∨转化为┐a -> b ^ ┐b -> a  ,a一定为真就转换为  ┐a->a  的形式,然后将->转换成边,两边连上对应的点。

如果a与┐a在最后建成的图的同一个强连通分量里···那么布尔方程有解

如果a所在强连通分量的拓扑序在┐a所在强连通分量的拓扑序之后,那么a为真,之前a为假,如果相等则真假均可以取。这里求拓扑序直接用tarjian即可,先找到的强连通分量的拓扑序一定更大

以上就是基本知识

但这道题有点特殊·····因为包含a拓扑序与┐a相等的情况要判断····用tarjian的话有点麻烦···

但n很小···直接dfs即可····若a可以到达┐a,则说明a可能与┐a在同一强连通分量或者a所在强连通分量的拓扑序小于等于┐a的拓扑序

代码:

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cctype>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;
const int N=2005;
const int M=10005;
int first[N],next[M],go[M],tot=1,n,m;
bool visit[N];
inline void comb(int a,int b)
{
  next[++tot]=first[a],first[a]=tot,go[tot]=b;
}
inline int tran(int a)
{
  return (a%2==1)?a+1:a-1;
}
inline int R()
{
  char c;
  int f=0;
  for(c=getchar();c<0||c>9;c=getchar());
  for(;c<=9&&c>=0;c=getchar())
    f=(f<<3)+(f<<1)+c-0;
  return f;
}
inline void dfs(int u)
{
  visit[u]=true;
  for(int e=first[u];e;e=next[e])
  {  
    if(!visit[go[e]])
      dfs(go[e]);
  }
}
inline bool jud(int u)
{
  memset(visit,false,sizeof(visit));
  dfs(u);
  if(!visit[tran(u)])  return true; 
  else return false;
}
int main()
{
  //freopen("a.in","r",stdin);
  n=R(),m=R();
  char s[5],t[5];
  int a,b;
  for(int i=1;i<=m;i++)
  {
    scanf("%d%s%d%s",&a,s,&b,t);
    int t1,t2;
    if(s[0]==Y)
      t1=a*2-1;
    else
      t1=a*2;
    if(t[0]==Y)
      t2=b*2-1;
    else
      t2=b*2;
    comb(tran(t2),t1);
    comb(tran(t1),t2);
  }
  for(int i=1;i<=n;i++)
  {
    bool flag1=jud(i*2-1);
    bool flag2=jud(i*2);
    if(!flag1&&!flag2)  {cout<<"IMPOSSIBLE"<<endl;return 0;}
    else if(!flag1)  cout<<"N";
    else if(!flag2)  cout<<"Y";
    else cout<<"?";
  }
  return 0;
}

 

算法复习——2—sat(bzoj2199)

标签:script   false   line   main   包含   ring   div   tar   lib   

原文地址:http://www.cnblogs.com/AseanA/p/7476747.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!