【POJ 3683】Priest John's Busiest Day

最新推荐文章于 2017-08-01 18:34:36 发布
最新推荐文章于 2017-08-01 18:34:36 发布
阅读量1.2k 收藏
点赞数 1
分类专栏： ~2-SAT~ 文章标签： POJ 2-SAT OI ACM
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Regina8023/article/details/42015987
版权
本文通过POJ 3683题目介绍了如何使用2-SAT算法解决一个实际问题：在一个特殊日子里，神父需要参加多对新人的婚礼，每个婚礼都有特定的开始和结束时间，且需要在婚礼开始或结束时进行特定仪式。文章详细解释了2-SAT的解决方案，包括缩点、反向图、拓扑排序和染色标记的过程，并分享了作者在理解和实现2-SAT过程中遇到的错误和感悟。
摘要由CSDN通过智能技术生成

Priest John's Busiest Day
Time Limit: 2000MS Memory Limit: 65536K
Total Submissions: 8300 Accepted: 2834 Special Judge

Description

John is the only priest in his town. September 1st is the John's busiest day in a year because there is an old legend in the town that the couple who get married on that day will be forever blessed by the God of Love. This year N couples plan to get married on the blessed day. The i-th couple plan to hold their wedding from time Si to timeTi. According to the traditions in the town, there must be a special ceremony on which the couple stand before the priest and accept blessings. The i-th couple need Diminutes to finish this ceremony. Moreover, this ceremony must be either at the beginning or the ending of the wedding (i.e. it must be either from Si to Si + Di, or from Ti- Di to Ti). Could you tell John how to arrange his schedule so that he can present at every special ceremonies of the weddings.

Note that John can not be present at two weddings simultaneously.

Input

The first line contains a integer N ( 1 ≤ N ≤ 1000). 
The next N lines contain the Si, Ti and Di. Si and Ti are in the format of hh:mm.

Output

The first line of output contains "YES" or "NO" indicating whether John can be present at every special ceremony. If it is "YES", output another N lines describing the staring time and finishing time of all the ceremonies.

Sample Input

2
08:00 09:00 30
08:15 09:00 20

Sample Output

YES
08:00 08:30
08:40 09:00

Source

POJ Founder Monthly Contest – 2008.08.31, Dagger and Facer

2-SAT可行性判断+方案输出模板题。


输出方案的方法：

缩点后将图反向，建成新图。


然后拓扑排序。


最后进行染色标记，只传递不选择的标记，不传递选择标记。



#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <stack>
#define m 5005
using namespace std;
int cf[m],l[m],t=0,sccno[m],pre[m],n,h[m],h2[m],a[m],b[m],done[m],color[m],scc_cnt,tot=0,tot2=0,lowlink[m],id[m];
int dfs_clock=0;
stack<int> S;
struct edge
{
	int y,ne;
}e[m*200],e2[m*200];
char s[10];
int Change(char *s)
{
	int ans=0;
	ans=(s[0]-'0')*600+(s[1]-'0')*60+(s[3]-'0')*10+s[4]-'0';
	return ans;
}
bool Conflict(int a,int b,int x,int y)
{
	if (b<=x||y<=a) return false;
	return true;
}
void Add(int x,int y)
{
	tot++;
	e[tot].ne=h[x];
	e[tot].y=y;
	h[x]=tot;
}
void Add2(int x,int y)
{
	tot2++;
	e2[tot2].ne=h2[x];
	e2[tot2].y=y;
	h2[x]=tot2;
}
void dfs(int u)
{
	pre[u]=lowlink[u]=++dfs_clock;
	S.push(u);
	for (int i=h[u];i;i=e[i].ne)
	{
		int v=e[i].y;
		if (!pre[v])
		{
			dfs(v);
			lowlink[u]=min(lowlink[u],lowlink[v]);
		}
		else if (!sccno[v])
			lowlink[u]=min(lowlink[u],pre[v]);
	}
	if (pre[u]==lowlink[u])
	{
		scc_cnt++;
		for (;;)
		{
			int x=S.top();
			S.pop();
			sccno[x]=scc_cnt;
			if (x==u) break;
		}
	}
}
void Findscc()
{
	dfs_clock=scc_cnt=0;
	memset(pre,0,sizeof(pre));
	memset(sccno,0,sizeof(sccno));
	for (int i=1;i<=n*2;i++)
		if (!pre[i]) dfs(i);
}
void Topsort(int x)  //拓扑排序
{
	id[x]=-1;
	done[++t]=x;
	for (int i=h2[x];i;i=e2[i].ne)
	{
		int y=e2[i].y;
		id[y]--;
		if (!id[y]) Topsort(y);
	}
} 
void Paint(int x)  //染色，即传递不选择标记
{
	color[x]=2;
	for (int i=h2[x];i;i=e2[i].ne)
	{
		int y=e2[i].y;
		if (!color[y]) Paint(y);
	}
}
void Print(int x)
{
	if (x/60<10) printf("0");
	printf("%d:",x/60);
	if (x%60<10) printf("0");
	printf("%d",x%60);
}
int main()
{
        scanf("%d",&n);
	for (int i=1;i<=n;i++)
	{
		scanf("%s",s);
		a[i]=Change(s);
		scanf("%s",s);
		b[i]=Change(s);
		scanf("%d",&l[i]);
	}
	for (int i=1;i<=n;i++)
		for (int j=i+1;j<=n;j++)
		{
			if (Conflict(a[i],a[i]+l[i],a[j],a[j]+l[j]))
	            Add(i,j+n),Add(j,i+n);
			if (Conflict(a[i],a[i]+l[i],b[j]-l[j],b[j]))
	            Add(i,j),Add(j+n,i+n);
			if (Conflict(b[i]-l[i],b[i],a[j],a[j]+l[j]))
				Add(i+n,j+n),Add(j,i);
            if (Conflict(b[i]-l[i],b[i],b[j]-l[j],b[j]))
				Add(i+n,j),Add(j+n,i);
		}
	Findscc();
	for (int i=1;i<=n;i++)   //判断可行性
		if (sccno[i]==sccno[i+n])
		{
			printf("NO\n");
			return 0;
		}
	printf("YES\n");
	for (int i=1;i<=2*n;i++)   //建立反向后的新图
		for (int j=h[i];j;j=e[j].ne)
		{
			int y=e[j].y;
			if (sccno[y]!=sccno[i])
			{
				Add2(sccno[y],sccno[i]);
				id[sccno[i]]++;
			}
		}
	for (int i=1;i<=scc_cnt;i++)
		if (!id[i]) Topsort(i);
	for (int i=1;i<=n;i++)   //找到每一个块的对立块（由对称性可知，对立块是唯一的）
		cf[sccno[i]]=sccno[i+n],cf[sccno[i+n]]=sccno[i];
	for (int i=1;i<=scc_cnt;i++)
	{
		int x=done[i];
		if (color[x]!=0) continue;
		color[x]=1;
		Paint(cf[x]); 
	}
	for (int i=1;i<=n;i++)
		if (color[sccno[i]]==1)
			Print(a[i]),printf(" "),Print(a[i]+l[i]),printf("\n");
	    else Print(b[i]-l[i]),printf(" "),Print(b[i]),printf("\n");
	return 0;
}



感悟：

1.WA是"YES"写成了"Yes"；在拓扑排序中用到的是缩点之后(e2,h2)的图，而我还用的缩点前(e,h)的。

RE是数组开小：edge最多可能有n^2条边


2.对于2-SAT的正确性至今还不能完全的理解。。目前的理解就是：

若选择B'的后代结点A，那么A'将被标记为不选择，然后不选择标记下传给A'的后代结点B；

因为整个图是对称的，所以选A的前代结点是B'，A'的后代结点一定是B。




图的对称性：



12.19UPD:


对于2—SAT的新理解：

Q:为什么要把图反向呢？

A:当选择了A点，那么A'被标记为不选择，此时要把不选择的标记向前面传递（A'不选择，A'之前的一定不选择，A'之后的不一定）。


如果此时图没有反向，那么要找到A'之前的点很麻烦（现在连得边都是向后走的）；


而反向之后直接沿着边走，就相当于沿着反向边走了！


Q:为什么要拓扑排序呢？


A:我们从入度为零的点A开始染色，A标记为选择，那么A'标记为不选择；


因为图有对称性，A'的出度一定为0，只要把A’标记为不可选择就行了，此时没有其他任何点受到影响。。


然后把A和A’当做删掉了。


接下来继续按照拓扑序找到下一个点B，此时B入度为0，B'出度为0，只要把B’标记为不可选择就行了，没有必要写一个Paint()函数，


因为他后面的点之前就已经置好标记了！


附上新的代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cstdlib>
#include <cmath>
#include <stack>
#define m 5005
using namespace std;
int cf[m],l[m],t=0,sccno[m],pre[m],n,h[m],h2[m],a[m],b[m],done[m],color[m],scc_cnt,tot=0,tot2=0,lowlink[m],id[m];
int dfs_clock=0;
stack<int> S;
struct edge
{
	int y,ne;
}e[m*200],e2[m*200];
char s[10];
int Change(char *s)
{
	int ans=0;
	ans=(s[0]-'0')*600+(s[1]-'0')*60+(s[3]-'0')*10+s[4]-'0';
	return ans;
}
bool Conflict(int a,int b,int x,int y)
{
	if (b<=x||y<=a) return false;
	return true;
}
void Add(int x,int y)
{
	tot++;
	e[tot].ne=h[x];
	e[tot].y=y;
	h[x]=tot;
}
void Add2(int x,int y)
{
	tot2++;
	e2[tot2].ne=h2[x];
	e2[tot2].y=y;
	h2[x]=tot2;
}
void dfs(int u)
{
	pre[u]=lowlink[u]=++dfs_clock;
	S.push(u);
	for (int i=h[u];i;i=e[i].ne)
	{
		int v=e[i].y;
		if (!pre[v])
		{
			dfs(v);
			lowlink[u]=min(lowlink[u],lowlink[v]);
		}
		else if (!sccno[v])
			lowlink[u]=min(lowlink[u],pre[v]);
	}
	if (pre[u]==lowlink[u])
	{
		scc_cnt++;
		for (;;)
		{
			int x=S.top();
			S.pop();
			sccno[x]=scc_cnt;
			if (x==u) break;
		}
	}
}
void Findscc()
{
	dfs_clock=scc_cnt=0;
	memset(pre,0,sizeof(pre));
	memset(sccno,0,sizeof(sccno));
	for (int i=1;i<=n*2;i++)
		if (!pre[i]) dfs(i);
}
void Topsort(int x)
{
	id[x]=-1;
	done[++t]=x;
	for (int i=h2[x];i;i=e2[i].ne)
	{
		int y=e2[i].y;
		id[y]--;
		if (!id[y]) Topsort(y);
	}
}
void Print(int x)
{
	if (x/60<10) printf("0");
	printf("%d:",x/60);
	if (x%60<10) printf("0");
	printf("%d",x%60);
}
int main()
{
    scanf("%d",&n);
	for (int i=1;i<=n;i++)
	{
		scanf("%s",s);
		a[i]=Change(s);
		scanf("%s",s);
		b[i]=Change(s);
		scanf("%d",&l[i]);
	}
	for (int i=1;i<=n;i++)
		for (int j=i+1;j<=n;j++)
		{
			if (Conflict(a[i],a[i]+l[i],a[j],a[j]+l[j]))
	            Add(i,j+n),Add(j,i+n);
			if (Conflict(a[i],a[i]+l[i],b[j]-l[j],b[j]))
	            Add(i,j),Add(j+n,i+n);
			if (Conflict(b[i]-l[i],b[i],a[j],a[j]+l[j]))
				Add(i+n,j+n),Add(j,i);
            if (Conflict(b[i]-l[i],b[i],b[j]-l[j],b[j]))
				Add(i+n,j),Add(j+n,i);
		}
	Findscc();
	for (int i=1;i<=n;i++)
		if (sccno[i]==sccno[i+n])
		{
			printf("NO\n");
			return 0;
		}
	printf("YES\n");
	for (int i=1;i<=2*n;i++)
		for (int j=h[i];j;j=e[j].ne)
		{
			int y=e[j].y;
			if (sccno[y]!=sccno[i])
			{
				Add2(sccno[y],sccno[i]);
				id[sccno[i]]++;
			}
		}
	for (int i=1;i<=scc_cnt;i++)
		if (!id[i]) Topsort(i);
	for (int i=1;i<=n;i++)
		cf[sccno[i]]=sccno[i+n],cf[sccno[i+n]]=sccno[i];
	for (int i=1;i<=scc_cnt;i++)
	{
		int x=done[i];
		if (color[x]!=0) continue;
		color[x]=1;
		color[cf[x]]=2;   //原来这里是paint()函数，直接删掉
	}
	for (int i=1;i<=n;i++)
		if (color[sccno[i]]==1)
			Print(a[i]),printf(" "),Print(a[i]+l[i]),printf("\n");
	    else Print(b[i]-l[i]),printf(" "),Print(b[i]),printf("\n");
	return 0;
}


(时间短了~)


Regina8023
关注
专栏目录
poj.rar_MáS
09-21
题意：给一个无向图，求从起点s到终点e尽力n条边的路径最小值。 边的数量2<=t<=100(每条边可以重复遍历) 算法：利用dp求出 path[ l ] [ i ] [ k ]=MIN(path[ l ][ i ] [ k ] , path[ l ][ i ] [ j ] + path...
【POJ3683】Priest John's Busiest Day
小蒟蒻yyb的博客
07-01 2228
POJ 3683 题目John is the only priest in his town. September 1st is the John’s busiest day in a year because there is an old legend in the town that the couple who get married on that day will be forever blessed b
POJ 3683 2-sat问题，不过要输出路径用topsort
bingshen的专栏
02-24 1119
Priest John's Busiest DayTime Limit: 2000MS Memory Limit: 65536KTotal Submissions: 3985 Accepted: 1394 Special Judge<br />Description<br />John is the only priest in his town. September 1st is the John's busiest day in a year because there is an old legend
poj 3683 2-sat
渐近自由
08-16 368
题意： 有n对新人打算在同一天举行结婚仪式，只有一个神父。 第i对新人的结婚仪式时间为 sti 到 edi，在其仪式开始或者结束时，需要进行一个用时为di的特别仪式。 即 si 到 si + di，或者 edi - di 到 edi，这两个时间段，该时间段需要神父在场。 判断是否可以通过合理安排每个特别仪式在开始或者结束矩形，使得神父能够出席所有特殊仪式。 解析： 挑战程序设计
POJ 3683 2-SAT 输出可行方案
Hermit_Inwind的博客
08-01 424
http://poj.org/problem?id=3683 1，根据矛盾建边 2，tarjan求scc判断可行性 3，根据scc缩点建反向边 4，拓扑排序，按顺序输出选择情况 #include #include #include #include #include #include using namespace std; const int maxn = 2005
POJ3683 Priest John's Busiest Day
weixin_30872157的博客
11-23 66
POJ3683 Priest John's Busiest Day Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 5196 Accepted: 1802 Special Judge Description John is the only prie...
POJ2262-Goldbach's Conjecture
07-29
【标题】"POJ2262-Goldbach's Conjecture"是一个编程竞赛题目，源自北京大学的在线判题系统POJ（Problem Online Judge）。这个题目是基于著名的数学猜想——哥德巴赫猜想（Goldbach's Conjecture）设计的。 【描述...
POJ2388-Who's in the Middle
07-31
《POJ2388-Who's in the Middle》是一个典型的计算机编程竞赛题目，主要考察参赛者的算法设计和实现能力。这个题目源自北京大学的在线评测系统POJ（PKU Online Judge），是编程爱好者和学生提升编程技能的重要平台之...
POJ2528-Mayor's posters 测试数据
06-27
【标题】"POJ2528-Mayor's posters"是一个经典的编程竞赛题目，源自2003年Alberta Collegiate Programming Contest。该比赛每年都会吸引众多编程爱好者参加，旨在锻炼参赛者的算法设计和实现能力。题目“市长的海报...
poj 4003 Bob’s Race.md
11-17
poj 4003 Bob’s Race.md
poj 3683 Priest John's Busiest Day 【2-sat 经典建图输出一组解】【有向图tarjan + 反向拓扑 + 染色】
世界很大
07-24 844
Priest John's Busiest Day Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 8652   Accepted: 2954   Special Judge Description John is the only priest in his
poj3683
fighting
07-15 545
Priest John's Busiest Day Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 8640   Accepted: 2947   Special Judge Description John is the only priest in his
poj3683 Priest John's Busiest Day【2-sat n+m算法】
禅境花园
07-13 397
#include<stdio.h> #include<algorithm> #include<vector> #include<string.h> #include<math.h> using namespace std; const int N=1024; int mark[2*N],s[2*N],c,n,st[N][2],ed[N][2],t[N]; vector<int>g[2*N];int
POJ 3683 2-sat
Dan__ge的博客
10-28 486
点击打开链接 题意：给了n个人的结婚时间，然后要有一个D的时间完成一个特殊活动，这个活动只能在开始或者结束之前完成，牧师需要主持这个特殊活动，问牧师是否可以主持n个活动不冲突，若可以输出一种可能的情况 思路：挑战程序设计上的例题，对于活动A和B来说，若A开始B开始冲突的话，那么!A V !B为真，然后还有开始和结束，结束和开始，结束和结束这四种情况的式子，然后连边跑scc即可，然后输出的时候c
POJ 3683 Priest John's Busiest Day
Salvare219
05-16 460
题目链接 : POJ3683 Priest John's Busiest Day Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 6797   Accepted: 2318   Special Judge Description John is the o
poj 3683 2-sat问题，输出任意一组可行解
weixin_30267697的博客
11-10 97
/* 2sat问题 输出任意一组可行解 */ #include<stdio.h> #include<string.h> #include<stdlib.h> #include<algorithm> #include<queue> #include<vector> using namespace std; #define...
POJ 3683 2-SAT +输出方案
lentty1452的专栏
11-03 599
输出方案好烦啊 代码超级多~第一次写这种代码，照着别人的代码打的。。   建图： 对于婚礼i和婚礼j，i表示开始主持，i`表示结束主持，j也是一样 枚举每一对不同的i和j， 如果i和j冲突，连接ij` 如果i和j`冲突，连接ij 如果i`和j冲突，连接i`j` 如果i`和j冲突，连接i`j   输出方案： 我们可以把处在同一个强连通分量中的点和边缩
2-SAT——5.0（poj3683 Priest John's Busiest Day）
laoda扯一扯
08-08 2746
poj3683 Priest John's Busiest Day 题意：一个小镇里面只有一个牧师，现在有些新人要结婚，需要牧师分别去主持一个仪式，给出每对新人婚礼的开始时间 s 和结束时间 t ，还有他们俩的这个仪式需要的时间（每对新人需要的时间长短可能不同） d ，牧师可以
POJ编程题解资源集合
"这些题目是来自POJ (Programming Online Judge) 的一系列算法练习题，涵盖了基础到中等难度的编程问题。题目涉及了各种算法和数据结构，包括但不限于数学计算、字符串处理、逻辑推理、递归、动态规划、位运算等。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值