poj 3648 2-sat 输出任意一组解模板
转载地址:http://blog.csdn.net/qq172108805/article/details/7603351
/*
2-sat问题,题意:有对情侣结婚,请来n-1对夫妇,算上他们自己共n对,编号为0~~n-1,他们自己编号为0
所有人坐在桌子两旁,新娘不想看到对面的人有夫妻关系或偷奸关系,若有解,输出一组解,无解输出bad luck
思路:
1.根据偷奸关系建图(1h和2h有偷奸关系,建边1h->2w 2h->1w)
2.求强连通分量
3.判断有无解(任一对夫妇不在同一强连通分量中,有解;否则无解)
4.缩点建图(建反向图)
5.拓扑排序
6.由底向上求解(由于上面建的是反向图,所以自顶(无入度)向下)
至此,选出来的是做新娘对面的,要求输出,坐在新娘一边的
*/
#include<stdio.h>
#include<string.h>
#include<queue>
#include<stack>
using namespace std;
struct edge
{
int yong;
int v[1000000];
int next[1000000];
int head[1000000];
edge()
{
clear();
}
void clear()
{
yong=1;
memset(head,0,sizeof(head));
}
void add(int u,int w)
{
v[yong]=w;
next[yong]=head[u];
head[u]=yong;
yong++;
}
}e1,e2;
int dfn[1000],low[1000];
int n,m,index,scc;
queue<int>q;
stack<int>ss;
int ins[1000],belong[1000],dui[1000],ind[1000],fang[1000];
void tarjan(int u)//tarjan求强连通分量
{
int i,v;
dfn[u]=low[u]=index++;
ins[u]=1;
ss.push(u);//之前误写成队列
for(i=e1.head[u];i;i=e1.next[i])
{
v=e1.v[i];
if(dfn[v]==0)
{
tarjan(v);
if(low[v]<low[u])
low[u]=low[v];
}else if(ins[v]&&dfn[v]<low[u])
low[u]=dfn[v];
}
if(low[u]==dfn[u])
{
do{
v=ss.top();
ss.pop();
ins[v]=0;
belong[v]=scc;//标记所属强连通分量的标号
}while(v!=u);
scc++;
}
}
void topsort()//拓扑排序 和 求解
{
int v,u;
while(!q.empty())
q.pop();
int i;
for(i=0;i<scc;i++)
{
if(ind[i]==0)//0入度的节点入队
q.push(i);
} while(!q.empty())
{
v=q.front();
q.pop();
if(!fang[v])//0入度的节点涂色
{//fang[i]表示标号为i的强连通分量被选择 或 删除
fang[v]=1;//表示选择
fang[dui[v]]=2;//表示删除
} for(i=e2.head[v];i;i=e2.next[i])
{
u=e2.v[i];
if(--ind[u]==0)//减入度,0入度的入队
{
q.push(u);
}
}
}
} int main()
{
int i,j,a,b,v;
char c1,c2;
while(scanf("%d%d",&n,&m),m+n)
{//初始化
e1.clear();
e2.clear();
//建原图
for(i=1;i<=m;i++)
{
scanf("%d%c %d%c",&a,&c1,&b,&c2);
e1.add(c1=='h'?a*2+1:a*2,c2=='h'?b*2:b*2+1);
e1.add(c2=='h'?b*2+1:b*2,c1=='h'?a*2:a*2+1);//这里处理要仔细
}
e1.add(0,1);//新郎必选
//求强连通分量
index=1;//dfn[]标记
scc=0;//强连通分量个数,从0开始计数
memset(dfn,0,sizeof(dfn));
memset(ins,0,sizeof(ins));
for(i=0;i<2*n;i++)
{
if(!dfn[i])
tarjan(i);
}
//判断有无解
for(i=0;i<n;i++)
{//belong[i]表示i所属的强连通分量的标号
if(belong[2*i]==belong[2*i+1])//若有夫妇在同一强连通分量中,无解
break;
dui[belong[2*i+1]]=belong[2*i];//记录配偶所属的强连通分量的标号,求解的时候有用
dui[belong[2*i]]=belong[2*i+1];
}
if(i!=n)
{
printf("bad luck\n");
continue;
}
//构建反向缩点图
memset(ind,0,sizeof(ind));
for(i=0;i<2*n;i++)
{
for(j=e1.head[i];j;j=e1.next[j])
{
v=e1.v[j];
if(belong[i]!=belong[v])
{
e2.add(belong[v],belong[i]);//注意,这里的顺序
ind[belong[i]]++;//统计入度
}
}
}
//拓扑排序 和 求解 同时进行了
memset(fang,0,sizeof(fang));
topsort();
//输出解
for(i=1;i<n;i++)
{
if(fang[belong[2*i]]==1)
printf("%dh ",i);
else printf("%dw ",i);
}
printf("\n");
}
return 0;
}
poj 3648 2-sat 输出任意一组解模板的更多相关文章
- HDU 1507 Uncle Tom's Inherited Land*(二分匹配,输出任意一组解)
Uncle Tom's Inherited Land* Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (J ...
- poj 3683 2-sat问题,输出任意一组可行解
/* 2sat问题 输出任意一组可行解 */ #include<stdio.h> #include<string.h> #include<stdlib.h> #in ...
- 洛谷 P4174 [NOI2006]最大获利 && 洛谷 P2762 太空飞行计划问题 (最大权闭合子图 && 最小割输出任意一组方案)
https://www.luogu.org/problemnew/show/P4174 最大权闭合子图的模板 每个通讯站建一个点,点权为-Pi:每个用户建一个点,点权为Ci,分别向Ai和Bi对应的点连 ...
- Wedding (poj 3648 2-SAT 输出随意一组解)
Language: Default Wedding Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9004 Accept ...
- 2-sat 输出任意一组可行解&拓扑排序+缩点 poj3683
Priest John's Busiest Day Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 8170 Accept ...
- POJ - 3648 Wedding (2-SAT 输出解决方案)
题意:有N-1对夫妇和1对新郎新娘要出席婚礼,这N对人要坐在走廊两侧.要求每对夫妇要坐在不同侧.有M对人有通奸关系,对于这一对人,不能同时坐在新娘对面(新娘新郎也可能和别人有通奸关系).求如何避免冲突 ...
- SAM I AM UVA - 11419(最小顶点覆盖+输出一组解)
就是棋盘问题输出一组解 https://blog.csdn.net/llx523113241/article/details/47759745 http://www.matrix67.com/blog ...
- POJ 3648 Wedding(2-SAT的模型运用+DFS | Tarjan)
Wedding Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 10427 Accepted: 3170 Specia ...
- POJ 3648 Wedding (2-SAT,经典)
题意:新郎和新娘结婚,来了n-1对夫妻,这些夫妻包括新郎之间有通奸关系(包括男女,男男,女女),我们的目地是为了满足新娘,新娘对面不能坐着一对夫妻,也不能坐着有任何通奸关系的人,另外新郎一定要坐新娘对 ...
随机推荐
- windows deintall 12c client
1.unintall: close all oracle app C:\app\client\CICadmin\product\12.1.0\client_1\deinstall deinstall ...
- D. Winter Is Coming 贪心(好题)
http://codeforces.com/contest/747/problem/D 大概的思路就是找到所有两个负数夹着的线段,优先覆盖最小的长度.使得那时候不用换鞋,是最优的. 但是这里有个坑点, ...
- python 使用 Pyscript 调试 报错
UnicodeEncodeError: 'ascii' codec can't encode characters in position 13-16: ordinal not in range(12 ...
- 【工具】sublime使用技巧
Ctrl+N 新建一个编辑区,Ctrl+Shift+C 或!加 Ctrl+E新建一个骨架完好的文件. Ctrl+Shift+P开启命令模式,sshtml 切换html语法. esc退出,Ctrl+`打 ...
- Red Hat Linux常用命令
1.查看机器型号 [root@local ~]# dmidecode | grep "Product Name" Product Name: VMware Virtual Plat ...
- CentOS 6.4 linux下编译安装MySQL5.6.14
CentOS 6.4下通过yum安装的MySQL是5.1版的,比较老,所以就想通过源代码安装高版本的5.6.14. 正文: 一:卸载旧版本 使用下面的命令检查是否安装有MySQL Server rpm ...
- spring cloud 概念
微服务构架需要使用场景: 1.可以将一个系统拆分成几个系统. 2.每个子系统可以部署多个应用,多个应用之间可以使用负载均衡. 3.需要一个服务注册中心,所有的服务都在一个注册中心注册,负载均衡也是通过 ...
- mysql中对order by的函数substring_index() , find_in_set()使用
题目是这样的:sql = "select * from table where id in(3,1,2,5)"; 怎样使得查询的结果按照 3 ,1 , 2, 5来排序: ...
- 几个net命令
A.显示当前工作组服务器列表 net view,当不带选项使用本命令时,它就会显示当前域或网络上的计算机上的列表. 比如:查看这个IP上的共享资源,就可以 C:\\>net view 192 ...
- ansible API(开发应用)
7. ansible API(开发应用) 官网链接