POJ 3683 Priest John's Busiest Day

最新推荐文章于 2020-04-21 09:51:10 发布
最新推荐文章于 2020-04-21 09:51:10 发布
阅读量522 收藏
点赞数
分类专栏： ######## 图论 ####### +++++ 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__
关注
专栏目录
poj3683 Priest John's Busiest Day
浅笑~如夏
01-27 216
题目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 302
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 1039
每个婚礼只能选择从开始时间开始或者终点时间结束，所以可以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 432
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 83
2-SAT 可以把每一次仪式看成变量，0/1的取值分别为开头举行和结尾举行。 转换为2-SAT接受的命题，就是看某一次仪式中有没有重合的时间段，有的话，就按照不冲突的形式连有向边。 然后跑tarjan就行啦，我们把时间全部转成分钟方便处理。。 #include <iostream> #include <cstdio> #include <stack> #incl...
【POJ3683】Priest John's Busiest Day
小蒟蒻yyb的博客
07-01 2237
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 219
传送门：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 139
题意： 题目大意：某一天结婚的人特别多但是主持婚礼的神父只有一个。婚礼时间从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...
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_3131.zip_POJ 八数码_poj
09-20
标题中的“POJ_3131.zip_POJ 八数码_poj”指的是一个与编程竞赛网站POJ（Problem Set Algorithm）相关的项目，具体是针对3131号问题的解决方案，该问题涉及到了八数码游戏。八数码游戏，又称滑动拼图，是一个经典的...
POJ2262-Goldbach's Conjecture
07-29
【标题】"POJ2262-Goldbach's Conjecture"是一个编程竞赛题目，源自北京大学的在线判题系统POJ（Problem Online Judge）。这个题目是基于著名的数学猜想——哥德巴赫猜想（Goldbach's Conjecture）设计的。 【描述...
POJ 1797 Heavy Transportation
源代码
03-02 864
题目大意： 注释代码： /* * Problem ID : POJ 1797 Heavy Transportation * Author : Lirx.t.Una * Language : GCC * Run Time : 279 ms * Run Memory
POJ 2485 Highways
源代码
03-07 683
题目大意： 注释代码： 无注释代码： 单词解释：
POJ 1860 Currency Exchange
源代码
08-20 654
题目大意： 题目链接 注释代码：
POJ 2502 Subway
源代码
08-22 653
题目大意： 题目链接 注释代码： 无注释代码： #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 3925 Minimal Ratio Tree
源代码
05-07 627
题目大意： 题目链接 注释代码： 无注释代码： 单词解释：
POJ编程题解资源集合
"这些题目是来自POJ (Programming Online Judge) 的一系列算法练习题，涵盖了基础到中等难度的编程问题。题目涉及了各种算法和数据结构，包括但不限于数学计算、字符串处理、逻辑推理、递归、动态规划、位运算等。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值