题目:
http://acm.hdu.edu.cn/showproblem.php?pid=1814
题意:
有n个党派每个党派有2人,现在要选n个议员,要使每个党派都入选一人,然后有些人不能同时入选。求一个字典序最小的解
思路:
直接暴力枚举DFS,首先将所有的点都置为为染色,然后我们从第一个点开始DFS染色,我们先尝试将i染成红色(答案中的颜色1),将~i染成蓝色(2),然后dfs i的所有后继并染色,如果对于后继j没有染色,那么将j然后为红色,~j染成蓝色。如果后继j已经被染成蓝色,则说明不能选则i,如果j已经染成红色,则说明可以。那么这些后继就可以被选择。如果选择i的时候失败了,那么必定要选择~i,如果也失败,则说明无解。否则按次序选取下一个未被染色的点(这样保证了字典序最小)。时间复杂度O(nm)。
该算法虽然在多数情况下时间复杂度到不了O(NM),但是综合性能仍然不如下面的O(M)算法。不过,该算法有一个很重要的用处:求字典序最小的解!
如果原图中的同一对点编号都是连续的(01、23、45……)则可以依次尝试第0对、第1对……点,每对点中先尝试编号小的,若失败再尝试编号大的。这样一定能求出字典序最小的解(如果有解的话),因为一个点一旦被确定,则不可更改。
如果原图中的同一对点编号不连续(比如03、25、14……)则按照该对点中编号小的点的编号递增顺序将每对点排序,然后依次扫描排序后的每对点,先尝试其编号小的点,若成功则将这个点选上,否则尝试编号大的点,若成功则选上,否则(都失败)无解。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 20010;
struct edge
{
int to, next;
}g[N*10];
int cnt, head[N];
int n, m, tot;
int mark[N], tmp[N];
void add_edge(int v, int u)
{
g[cnt].to = u, g[cnt].next = head[v], head[v] = cnt++;
}
bool dfs(int v)
{
if(mark[v] == 2) return false;
if(mark[v] == 1) return true;
tmp[tot++] = v;
mark[v] = 1, mark[v^1] = 2;
for(int i = head[v]; i != -1; i = g[i].next)
{
int u = g[i].to;
if(! dfs(u)) return false;
}
return true;
}
bool solve()
{
memset(mark, 0, sizeof mark);
for(int i = 0; i < n; i++)
{
if(mark[i]) continue;
tot = 0;
if(! dfs(i))
{
for(int i = 0; i < tot; i++) mark[tmp[i]] = mark[tmp[i]^1] = 0;
if(! dfs(i^1)) return false;
}
}
return true;
}
int main()
{
int v, u;
while(~ scanf("%d%d", &n, &m))
{
cnt = 0;
memset(head, -1, sizeof head);
n *= 2;
for(int i = 1; i <= m; i++)
{
scanf("%d%d", &v, &u);
v--, u--;
add_edge(v, u^1), add_edge(u, v^1);
}
if(solve())
{
for(int i = 0; i < n; i++)
if(mark[i] == 1) printf("%d\n", i+1);
}
else printf("NIE\n");
}
return 0;
}
lrj白书上的实现
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 20010;
struct edge
{
int to, next;
}g[N*10];
int cnt, head[N];
int n, m, tot;
int tmp[N];
bool mark[N];
void add_edge(int v, int u)
{
g[cnt].to = u, g[cnt].next = head[v], head[v] = cnt++;
}
bool dfs(int v)
{
if(mark[v^1]) return false;
if(mark[v]) return true;
tmp[tot++] = v;
mark[v] = true;
for(int i = head[v]; i != -1; i = g[i].next)
{
int u = g[i].to;
if(! dfs(u)) return false;
}
return true;
}
bool solve()
{
memset(mark, 0, sizeof mark);
for(int i = 0; i < 2*n; i += 2)
if(!mark[i] && !mark[i+1])
{
tot = 0;
if(! dfs(i))
{
for(int i = 0; i < tot; i++) mark[tmp[i]] = false;
if(! dfs(i+1)) return false;
}
}
return true;
}
int main()
{
int v, u;
while(~ scanf("%d%d", &n, &m))
{
cnt = 0;
memset(head, -1, sizeof head);
for(int i = 1; i <= m; i++)
{
scanf("%d%d", &v, &u);
v--, u--;
add_edge(v, u^1), add_edge(u, v^1);
}
if(solve())
{
for(int i = 0; i < 2*n; i += 2)
if(mark[i]) printf("%d\n", i+1);
else printf("%d\n", i+2);
}
else printf("NIE\n");
}
return 0;
}