BZOJ 2199: [Usaco2011 Jan]奶牛议会 2-sat

2199: [Usaco2011 Jan]奶牛议会

Time Limit: 10 Sec  Memory Limit: 259 MB
Submit: 387  Solved: 250
[Submit][Status][Discuss]

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?

好久以前写的,忘了是啥了。。。

就知道是2-sat

好,那就这样吧。。。


#include<cmath>
#include<ctime>
#include<cstdio> 
#include<climits>
#include<cstring> 
#include<cstdlib> 
#include<iostream> 
#include<algorithm> 
#include<iomanip> 
#include<vector> 
#include<string>
#include<queue>  
#include<map> 
#include<set>
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<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();} 
    return f*x; 
} 
inline int getopt()
{
    char ch=getchar();
    while(ch!='N'&&ch!='Y')ch=getchar();
    if(ch=='N')return 0;
    else return 1;
}
const int N=1010;
int n,m,ecnt,last[N<<1];
bool book[N<<1];
struct EDGE{int to,nt;}e[N<<5];
inline void add(int u,int v)
{e[++ecnt]=(EDGE){v,last[u]};last[u]=ecnt;}
inline void dfs(int u)
{
    book[u]=1;
    for(int i=last[u];i;i=e[i].nt)
    if(!book[e[i].to])dfs(e[i].to);
}
inline int check(int k)
{
    memset(book,0,sizeof(book));
    dfs(k);
    for(int i=1;i<=n;i++)if(book[i<<1]&&book[i<<1|1])return 0;
    return 1;
}
int main()
{
    n=read();m=read();int x,y,opt1,opt2;char ans[N];
    for(int i=1;i<=m;i++)
    {
        x=read();opt1=getopt();
        y=read();opt2=getopt();
        add((y<<1)+(opt2^1),(x<<1)+opt1);add((x<<1)+(opt1^1),(y<<1)+opt2);
    }
    dfs(2);dfs(3);
    for(int i=1;i<=n;i++)
    {
        opt1=check(i<<1);
        opt2=check(i<<1|1);
        if(!opt1&&!opt2){puts("IMPOSSIBLE");return 0;}
        else if(opt1&&!opt2)ans[i]='N';
        else if(!opt1&&opt2)ans[i]='Y';
        else ans[i]='?';
    }
    for(int i=1;i<=n;i++)printf("%c",ans[i]);puts("");
    return 0;
} 


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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值