[hit 1917] Peaceful Commission(2-sat)

标标准准的2-sat问题。
题意是有n组人,每组两人。编号从1开始,所以第i组的人的编号分别是2i-1和2i。
给出一些关系,表明a和b不和。不和的人不能同时被选中。问有没有办法在符合要求的情况下,从每组中选出一个人来。
这道题不仅要求是否有解,还要输出一组可行解,任意可行解均能AC。
模板题,总结一下解题过程:
    1.建图
    2.求有向图的强连通分量
    3.判断是否存在解
    4.将每个强连通分量缩成一个点,此时原图成为向无环图
    5.对该有向无环图进行拓扑排序
    6.逆拓扑序列进行筛选、染色
    7.输出结果

还是略繁琐的。

#include<iostream>
#include<vector>
#include<cstdio>
#include<algorithm>
using namespace std;

#define maxn 20000

vector<int> tt[maxn]; // 原有向图

int n, m;

int vis[maxn], tag[maxn];
int low[maxn], dfn[maxn];
bool instack[maxn];
int stap[maxn];
int ind, stop, cnt;
int belong[maxn];

vector<int> mp[maxn]; // 存联通分量内的结点
vector<int> ord; // 存拓扑序列

int color[maxn]; // 非着色结点即为答案

int min(int a, int b)
{
return a < b ? a : b;
}

void tarjan(int u) // 求强连通分量
{
dfn[u] = low[u] = ind++;
instack[u] = true;
stap[++stop] = u;
for(int i = 0; i < tt[u].size(); i++)
{
int v = tt[u][i];
if(!dfn[v])
{
tarjan(v);
low[u] = min(low[u], low[v]);
}
else if(instack[v])
low[u] = min(low[u], dfn[v]);
}
if(low[u] == dfn[u])
{
cnt++;
int v;
do
{
v = stap[stop--];
instack[v] = false;
belong[v] = cnt;
}while(v != u);
}
}

void dfs(int u) //拓扑排序
{
vis[u] = 1;
for(int i = 0; i < mp[u].size(); i++)
{
if(!vis[mp[u][i]])
dfs(mp[u][i]);
}
ord.push_back(u);
}

int main()
{
int a, b;
while(~scanf("%d%d", &n, &m))
{
for(int i = 0; i <= 2*n; i++)
tt[i].clear();
// 建边
for(int i = 1; i <= m; i++)
{
scanf("%d%d", &a, &b);
tt[a].push_back(((b-1)^1)+1);
tt[b].push_back(((a-1)^1)+1);
}

memset(low, 0, sizeof(low));
memset(dfn, 0, sizeof(dfn));
memset(instack, false, sizeof(instack));
memset(belong, 0, sizeof(belong));
ind = stop = cnt = 0;
// 求强连通分量
for(int i = 1; i <= 2*n; i++)
{
if(!dfn[i])
tarjan(i);
}
// 判断是否有解
int flag = 0;
for(int i = 1; i <= n; i++)
{
if(belong[2*i-1] == belong[2*i])
{
flag = 1;
break;
}
}
if(flag)
{
printf("NIE\n");
continue;
}

for(int i = 0; i <= 2*n; i++)
mp[i].clear();
// 缩点
for(int i = 1; i <= 2*n; i++)
{
for(int j = 0; j < tt[i].size(); j++)
{
int v = tt[i][j];
// 这个结点的邻边不在这个节点的强连通分量中
// 那么就要建立一条有向边,即删除同一分量的边
if(belong[v] != belong[i])
mp[belong[i]].push_back(belong[v]);

}
}

memset(vis, 0, sizeof(vis));
ord.clear();

// 对缩点后的有向图进行拓扑排序
for(int i = 1; i <= cnt; i++)
if(!vis[i]) dfs(i);

for(int i = 0; i <= cnt; i++)
mp[i].clear();

// 把同个连通分量内的结点放在一起
for(int i = 1; i <= 2*n; i++)
mp[belong[i]].push_back(i);

memset(vis, 0, sizeof(vis));
memset(color, 0, sizeof(color));

//逆向拓扑序列进行染色
for(int i = 0; i < ord.size(); i++)
{
int u = ord[i];
// 未访问过,这个连通分量全选
if(!vis[u])
{
for(int j = 0; j < mp[u].size(); j++)
{
// 找到此联通分量中每个结点的另一半
// 然后染色,染色的结点必定在某个分量中
// 所以那个分量全部不要
int cur = ((mp[u][j]-1)^1)+1;
color[cur] = 1;
vis[belong[cur]] = 1;
}
}
}
for(int i = 1; i <= 2*n; i++)
if(!color[i])
printf("%d\n", i);
}
return 0;
}



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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值