BZOJ2199: [Usaco2011 Jan]奶牛议会（2-SAT+DP）

最新推荐文章于 2018-06-16 19:46:58 发布

原创最新推荐文章于 2018-06-16 19:46:58 发布 · 423 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#dp

2-SAT 专栏收录该内容

7 篇文章

订阅专栏

本文介绍了一种解决2-SAT问题的方法，通过构建图并利用Tarjan算法进行缩点，最后使用位运算加速拓扑排序过程。该方法能够高效判断每个变量是否能同时选择两种状态或只能选择一种状态。

传送门

题意：
2-SAT判定每个点是两种方案都可选还是只能选一种方案。

题解：
明显，两种方案都可选只有两种方案之间没有可以到达的连边。可以每个点dfs一遍判断可行性，也可以缩点后在拓扑图上dp。如果dp时用bitset会更快。

#include<bits/stdc++.h>
using namespace std;
const int Maxn=2e3+50;
int n,m,lim;char ch[3];
inline int Is(int x){return (x-1)<<1;}
inline int Isnot(int x){return (x<<1)-1;}
bitset<Maxn>t[Maxn];
vector<int>edge[Maxn];
vector<int>nowedge[Maxn];
int vis[Maxn],dfn[Maxn],id[Maxn],cnt,low[Maxn],sta[Maxn],top,ins[Maxn],ind;
inline void tarjan(int now)
{
    dfn[now]=low[now]=++ind;ins[sta[++top]=now]=1;
    for(int e=edge[now].size()-1;e>=0;e--)
    {
        int v=edge[now][e];
        if(!dfn[v])tarjan(v),low[now]=min(low[now],low[v]);
        else if(ins[v])low[now]=min(low[now],dfn[v]);
    }
    if(low[now]==dfn[now])
    {
        ++cnt;
        do
        {
            ins[sta[top]]=0;
            id[sta[top]]=cnt;
            t[cnt][sta[top]]=1;
        }while(sta[top--]!=now);
    }
}
inline void dfs(int now)
{
    vis[now]=1;
    for(int e=nowedge[now].size()-1;e>=0;e--)
    {
        int v=nowedge[now][e];
        if(!vis[v])dfs(v);
        t[now]|=t[v];
    }
}

int main()
{
    scanf("%d%d",&n,&m);lim=Isnot(n);
    for(int i=1;i<=m;i++)
    {
        int x;scanf("%d",&x);scanf("%s",ch+1);x=((ch[1]=='Y')?(Is(x)):(Isnot(x)));
        int y;scanf("%d",&y);scanf("%s",ch+1);y=((ch[1]=='Y')?(Is(y)):(Isnot(y)));
        edge[x^1].push_back(y);edge[y^1].push_back(x);
    }
    for(int i=0;i<=lim;i++)
        if(!dfn[i])tarjan(i);
    for(int i=0;i<=lim;i++)
        for(int e=edge[i].size()-1;e>=0;e--)
        if(id[edge[i][e]]!=id[i])nowedge[id[i]].push_back(id[edge[i][e]]);
    for(int i=1;i<=cnt;i++)if(!vis[i])dfs(i);
    for(int i=1;i<=n;i++)
    {
        if(t[id[Is(i)]][Isnot(i)]==1&&t[id[Isnot(i)]][Is(i)]==1){puts("IMPOSSIBLE");return 0;}
    }
    for(int i=1;i<=n;i++)
    {
        if(t[id[Is(i)]][Isnot(i)]==1)putchar('N');
        else if(t[id[Isnot(i)]][Is(i)]==1)putchar('Y');
        else putchar('?');
    }
}