POJ 3683 Priest John's Busiest Day

最新推荐文章于 2020-04-21 09:51:10 发布
最新推荐文章于 2020-04-21 09:51:10 发布
阅读量493 收藏
点赞数
分类专栏： ######## 图论 ####### +++++ 2-SAT +++++ POJ 文章标签： 2-SAT SCC Kosaraju
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u013569656/article/details/39028559
版权

题目大意：

题目链接

注释代码：

/*
 * Problem ID : POJ 3683 Priest John's Busiest Day
 * Author     : Lirx.t.Una
 * Language   : C++
 * Run Time   : 266 ms
 * Run Memory : 23832 KB
*/

#include <iostream>
#include <cstring>
#include <cstdio>
#include <stack>

//思路：
//如果x表示第x号婚礼在开始时举行而-x表示在结束时举行
//比如x∧y就表示婚礼x和y都在开始时举行
//但是如果从题给的事实中得出x∧y不成立（即两场婚礼会有冲突）
//则不冲突的情况只可能发生在-(x∧y)中了，如果该表达式也无法成立的话则牧师就不能完成任务了

//由此可见算法的步骤为：
//1：找出所有使i∧j不成立的(i, j)对子，其中i、j为布尔文字，i、j可以为x也可以为-x
//2：然后将i∧j化为-i∨-j形式（为了构造主合取范式），最后便是求布尔方程(a1∨b1)∧(a2∨b2)...(an^bn) = true
//3：对于每个布尔子句i∨j，可以化为蕴含表达式-i->j ∧ -j->i，这样主合取范式就全部化为蕴含式的∧
//4：对i和-i进行编号，然后对于每一对i->j做一条有向边构图进行强连通分量分解
//5：如果存在i和-i在同一个强连通分量中则布尔方程无解，否则就有解
//6：对于有解的情况，如果要求某个布尔变量x的取值
    //如果x所在的强连通分量在缩点图中的拓扑排序在-x之后的强连通分量之后就去true否则就取false

//婚礼的最大数量
#define	MAXN	1000
#define	MAXNN	2000

using namespace std;

short	s[MAXN + 1];//婚礼开始时间
short	t[MAXN + 1];//截止时间
short	d[MAXN + 1];//祈祷用时

//正向图和反向图
//其中某个点i表示i号婚礼在开始时进行，则i + n表示在i在截止时进行两种布尔状态
//其中n为婚礼数量
int		head[MAXNN + 1];
int		to[MAXNN * MAXNN + 1];
int		nxt[MAXNN * MAXNN + 1];
int		_head[MAXNN + 1];
int		_to[MAXNN * MAXNN + 1];
int		_nxt[MAXNN * MAXNN + 1];
int		e;

short	blg[MAXNN + 1];
bool	vis[MAXNN + 1];
int		nb;

stack<short>	stk;

inline int
min( int a, int b ) {

	return a < b ? a : b;
}

inline int
max( int a, int b ) {

	return a > b ? a : b;
}

inline void
addarc( int u, int v ) {

	to[e]   = v;
	nxt[e]  = head[u];
	head[u] = e;

	_to[e]   = u;
	_nxt[e]  = _head[v];
	_head[v] = e++;
}

void
dfs(int u) {

	int		v, i;

	vis[u] = true;
	for ( i = head[u]; i; i = nxt[i] )
		if ( !vis[ v = to[i] ] )
			dfs(v);
	stk.push(u);
}

void
_dfs(int u) {

	int		v, i;

	vis[u] = true;
	blg[u] = nb;

	for ( i = _head[u]; i; i = _nxt[i] )
		if ( !vis[ v = _to[i] ] )
			_dfs(v);
}

void
kosa(int nn) {

	int		i;

	for ( i = 1; i <= nn; i++ )
		if ( !vis[i] )
			dfs(i);

	memset(vis, 0, sizeof(vis));
	for ( nb = 0; !stk.empty(); stk.pop() )
		if ( !vis[ i = stk.top() ] )
			++nb, _dfs(i);
}

int
main() {

	int		n, nn;//婚礼数和布尔状态数
	int		hh, mm, dd;//临时接受小时、分钟、用时

	int		i, j;
	int		tmp;

	scanf("%d", &n);
	nn = n << 1;
	e  = 1;

	for ( i = 1; i <= n; i++ ) {
	
		scanf("%d:%d", &hh, &mm);
		s[i]  = hh * 60 + mm;//将时间转化为以分钟为单位

		scanf("%d:%d", &hh, &mm);
		t[i]  = hh * 60 + mm;

		scanf("%d", &dd);
		d[i] = dd;
	}

	for ( i = 1; i <= n; i++ )
		for ( j = i + 1; j <= n; j++ ) {//找出所有^表达式不成立的对子
		
			//i∧j ≡ -i∨-j ≡ i->-j∧j->-i
			if ( min( s[i] + d[i], s[j] + d[j] ) > max( s[i], s[j] ) ) {//都在开始时冲突

				addarc( i, j + n );
				addarc( j, i + n );
			}

			//i∧-j ≡ -i∨j ≡ i->j∧-j->-i
			if ( min( s[i] + d[i], t[j] ) > max( s[i], t[j] - d[j] ) ) {//开始和截止冲突
			
				addarc( i, j );
				addarc( j + n, i + n );
			}

			//-i∧j ≡ i∨-j ≡ -i->-j∧j->i
			if ( min( t[i], s[j] + d[j] ) > max( t[i] - d[i], s[j] ) ) {//截止和开始冲突
			
				addarc( i + n, j + n );
				addarc( j, i );
			}

			//-i∧-j ≡ i∨j ≡ -i->j∧-j->i
			if ( min( t[i], t[j] ) > max( t[i] - d[i], t[j] - d[j] ) ) {//都在截止时冲突
			
				addarc( i + n, j );
				addarc( j + n, i );
			}
		}

	kosa(nn);

	for ( i = 1; i <= n; i++ )
		if ( blg[i] == blg[i + n] ) {//自己和自己的非在同一个SCC中则布尔方程误解
		
			puts("NO");
			return 0;
		}

	puts("YES");
	for ( i = 1; i <= n; i++ )
		if ( blg[i] > blg[i + n] ) {//自己比自己的非的强连通分量的拓扑排序序号要大则取true（就表示在开始时举办）

			tmp = s[i] + d[i];
			printf("%02d:%02d %02d:%02d\n",
				   s[i] / 60, s[i] % 60,
				   tmp / 60, tmp %  60);
		}
		else {//否则就false的情况（在截止时举办）

			tmp = t[i] - d[i];
			printf("%02d:%02d %02d:%02d\n",
				   tmp / 60, tmp % 60,
				   t[i] / 60, t[i] %  60);
		}

	return 0;
}
无注释代码： 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <stack>

#define	MAXN	1000
#define	MAXNN	2000

using namespace std;

short	s[MAXN + 1];
short	t[MAXN + 1];
short	d[MAXN + 1];

int		head[MAXNN + 1];
int		to[MAXNN * MAXNN + 1];
int		nxt[MAXNN * MAXNN + 1];
int		_head[MAXNN + 1];
int		_to[MAXNN * MAXNN + 1];
int		_nxt[MAXNN * MAXNN + 1];
int		e;

short	blg[MAXNN + 1];
bool	vis[MAXNN + 1];
int		nb;

stack<short>	stk;

inline int
min( int a, int b ) {

	return a < b ? a : b;
}

inline int
max( int a, int b ) {

	return a > b ? a : b;
}

inline void
addarc( int u, int v ) {

	to[e]   = v;
	nxt[e]  = head[u];
	head[u] = e;

	_to[e]   = u;
	_nxt[e]  = _head[v];
	_head[v] = e++;
}

void
dfs(int u) {

	int		v, i;

	vis[u] = true;
	for ( i = head[u]; i; i = nxt[i] )
		if ( !vis[ v = to[i] ] )
			dfs(v);
	stk.push(u);
}

void
_dfs(int u) {

	int		v, i;

	vis[u] = true;
	blg[u] = nb;

	for ( i = _head[u]; i; i = _nxt[i] )
		if ( !vis[ v = _to[i] ] )
			_dfs(v);
}

void
kosa(int nn) {

	int		i;

	for ( i = 1; i <= nn; i++ )
		if ( !vis[i] )
			dfs(i);

	memset(vis, 0, sizeof(vis));
	for ( nb = 0; !stk.empty(); stk.pop() )
		if ( !vis[ i = stk.top() ] )
			++nb, _dfs(i);
}

int
main() {

	int		n, nn;
	int		hh, mm, dd;

	int		i, j;
	int		tmp;

	scanf("%d", &n);
	nn = n << 1;
	e  = 1;

	for ( i = 1; i <= n; i++ ) {
	
		scanf("%d:%d", &hh, &mm);
		s[i]  = hh * 60 + mm;

		scanf("%d:%d", &hh, &mm);
		t[i]  = hh * 60 + mm;

		scanf("%d", &dd);
		d[i] = dd;
	}

	for ( i = 1; i <= n; i++ )
		for ( j = i + 1; j <= n; j++ ) {
		
			if ( min( s[i] + d[i], s[j] + d[j] ) > max( s[i], s[j] ) ) {

				addarc( i, j + n );
				addarc( j, i + n );
			}

			if ( min( s[i] + d[i], t[j] ) > max( s[i], t[j] - d[j] ) ) {
			
				addarc( i, j );
				addarc( j + n, i + n );
			}

			if ( min( t[i], s[j] + d[j] ) > max( t[i] - d[i], s[j] ) ) {
			
				addarc( i + n, j + n );
				addarc( j, i );
			}

			if ( min( t[i], t[j] ) > max( t[i] - d[i], t[j] - d[j] ) ) {
			
				addarc( i + n, j );
				addarc( j + n, i );
			}
		}

	kosa(nn);

	for ( i = 1; i <= n; i++ )
		if ( blg[i] == blg[i + n] ) {
		
			puts("NO");
			return 0;
		}

	puts("YES");
	for ( i = 1; i <= n; i++ )
		if ( blg[i] > blg[i + n] ) {

			tmp = s[i] + d[i];
			printf("%02d:%02d %02d:%02d\n",
				   s[i] / 60, s[i] % 60,
				   tmp / 60, tmp %  60);
		}
		else {

			tmp = t[i] - d[i];
			printf("%02d:%02d %02d:%02d\n",
				   tmp / 60, tmp % 60,
				   t[i] / 60, t[i] %  60);
		}

	return 0;
}


__Lirx_t_Una__
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
poj3683 Priest John's Busiest Day
浅笑~如夏
01-27 196
题目2-sat问题，要求要输出方案。先考虑建模，把每个婚礼都拆成两个部分。对于两个婚礼，A，B，用A，B表示开始，A’，B’表示结束，如果A与B,B’都冲突的话，A就不能选，连向A’，否则只有一个冲突的话，A连B’,A’连B。这样图就建好了，那么如何输出方案呢，一般tarjan后跑一个拓扑序就好了，但是有更好的方案，tarjan相当于就是一个拓扑了，最后只要A与A’谁的强连通分量编号小就选谁就好了。
POJ-3683 Priest John's Busiest Day
flx413的博客
05-04 277
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...
poj 3683 Priest John's Busiest Day
work_freedom的专栏
05-02 994
每个婚礼只能选择从开始时间开始或者终点时间结束，所以可以2-sat。 建图：第i个婚礼  从开始时间开始用点 2*i 标记， 在终点时间结束用2*i+1 标记， 比较婚礼 i和j ，如果 2*i 和 2*j 有冲突， 则连一条 2*i到2*j+1的边(其他同理）...... 输出解： 先将图缩点，然后建立一个逆向的图，在拓扑输出解。因为如果正向的图，选了一个点，它能到达的点都需要选，而如果逆向
POJ3683 Priest John's Busiest Day
Jelly_acm的博客
10-09 387
Priest John’s Busiest DayTime Limit: 2000MS Memory Limit: 65536K Total Submissions: 10647 Accepted: 3626 Special JudgeDescriptionJohn is the only priest in his town. September 1st is
POJ 3683 Priest John's Busiest Day (算竞进阶习题)
weixin_30556959的博客
05-10 71
2-SAT 可以把每一次仪式看成变量，0/1的取值分别为开头举行和结尾举行。 转换为2-SAT接受的命题，就是看某一次仪式中有没有重合的时间段，有的话，就按照不冲突的形式连有向边。 然后跑tarjan就行啦，我们把时间全部转成分钟方便处理。。 #include <iostream> #include <cstdio> #include <stack> #incl...
【POJ3683】Priest John's Busiest Day
小蒟蒻yyb的博客
07-01 2134
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 Priest John's Busiest Day 【经典2-SAT】
Waiting For You
08-08 159
传送门：POJ 3683 Priest John's Busiest Day Time Limit: 2000MS Memory Limit: 65536K Problem Description John is the only priest in his town. September 1st is the John’s busiest day in a year because t...
Poj 3683 Priest John's Busiest DayPriest Day
weixin_45539557的博客
04-21 107
题意： 题目大意：某一天结婚的人特别多但是主持婚礼的神父只有一个。婚礼时间从s开始到e结束，神父必须在s到s+d或者e-d到e这段时间内在。给定了n个婚礼的s，e，d，求一种方案能使得神父主持所有的婚礼。 明显是2-SAT问题，把靠近前面区间定义为0，后面定义为1，则n^2先连边建图，判断是否存在则是i与i+n是否在同一scc 难点在于构造方案： 第一个方法首先缩点，其实是要在原图上从0出节点开始...
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...
POJ2262-Goldbach's Conjecture
07-29
北大POJ2262-Goldbach's Conjecture 解题报告+AC代码
POJ2388-Who's in the Middle
07-31
北大POJ2388-Who's in the Middle 解题报告+AC代码
POJ2528-Mayor's posters 测试数据
06-27
POJ2528-Mayor's posters 测试数据。数据来源：Alberta Collegiate Programming Contest 2003.10.18 – 问题G
poj 4003 Bob’s Race.md
11-17
poj 4003 Bob’s Race.md
POJ 1797 Heavy Transportation
源代码
03-02 840
题目大意： 注释代码： /* * Problem ID : POJ 1797 Heavy Transportation * Author : Lirx.t.Una * Language : GCC * Run Time : 279 ms * Run Memory
POJ 2485 Highways
源代码
03-07 662
题目大意： 注释代码： 无注释代码： 单词解释：
POJ 1860 Currency Exchange
源代码
08-20 632
题目大意： 题目链接 注释代码：
POJ 2502 Subway
源代码
08-22 620
题目大意： 题目链接 注释代码： 无注释代码： #include #include #include #include #define INF 1E15 #define MAXN 202 #define V_W 10 #define V_S 40 using namespace std; struct Point { int x, y; friend i
POJ 1511 Invitation Cards
源代码
03-05 599
题目大意： 注释代码： 无注释代码： 单词解释： theater：n, 电影院，剧院 invitation：n, 邀请
POJ1159
最新发布
06-01
POJ1159是一个动态规划问题，也被称为“Palindrome”。题目描述如下： 给定一个字符串，你需要将它分割成若干个子串，使得每个子串都是回文串。求最少需要分割几次。 例如，对于字符串“abcbm”，最少需要分割一次...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值