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 time Ti. 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 Di minutes 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.

题目大意：给n个婚礼，每个婚礼要举办一次祝福，这个祝福只能在婚礼的开始或结束的时候举办（大大的2SAT标志），问能否举办所有祝福，并输出祝福的时间段（任意解）。

思路：可以参考国家集训队2003年伍昱的论文，判冲连边就成。

PS：这题各种不严谨，没说Ti - Si ≥ Di，没保证不会超过24小时，反正都忽视掉是可以AC的，不忽视能不能AC我就不知道了……

  1.  #include <cstdio>
  2.  #include <cstring>
  3.  
  4.  const int MAXN = *;
  5.  const int MAXM = MAXN * MAXN * ;
  6.  
  7.  struct Topological{
  8.      int St[MAXN], c;
  9.      int n, ecnt, cnt;
  10.      int head[MAXN], order[MAXN], indeg[MAXN];
  11.      int next[MAXM], to[MAXM];
  12.  
  13.      void addEdge(int x, int y){
  14.          to[ecnt] = y; next[ecnt] = head[x]; head[x] = ecnt++;
  15.          ++indeg[y];
  16.          //printf("%d->%d\n",x,y);
  17.      }
  18.  
  19.      void init(int nn){
  20.          n = nn; ecnt = ;
  21.          memset(head, , sizeof(head));
  22.          memset(indeg,,sizeof(indeg));
  23.      }
  24.  
  25.      void build(){
  26.          c = cnt = ;
  27.          for(int i = ; i <= n; ++i)
  28.              if(indeg[i] == ) St[++c] = i;
  29.          while(c > ){
  30.              int u = St[c--]; order[cnt++] = u;
  31.              for(int p = head[u]; p; p = next[p]){
  32.                  int &v = to[p];
  33.                  --indeg[v];
  34.                  if(indeg[v] == ) St[++c] = v;
  35.              }
  36.          }
  37.      }
  38.  } T;
  39.  
  40.  struct TwoSAT{
  41.      int St[MAXN], c;
  42.      int n, ecnt, dfs_clock, scc_cnt;
  43.      int head[MAXN], sccno[MAXN], pre[MAXN], lowlink[MAXN];
  44.      int next[MAXM], to[MAXM];
  45.      int select[MAXN], sccnox[MAXN];
  46.  
  47.      void dfs(int u){
  48.          lowlink[u] = pre[u] = ++dfs_clock;
  49.          St[++c] = u;
  50.          for(int p = head[u]; p; p = next[p]){
  51.              int &v = to[p];
  52.              if(!pre[v]){
  53.                  dfs(v);
  54.                  if(lowlink[u] > lowlink[v]) lowlink[u] = lowlink[v];
  55.              }else if(!sccno[v]){
  56.                  if(lowlink[u] > pre[v]) lowlink[u] = pre[v];
  57.              }
  58.          }
  59.          if(lowlink[u] == pre[u]){
  60.              sccnox[++scc_cnt] = u;
  61.              while(true){
  62.                  int x = St[c--];
  63.                  sccno[x] = scc_cnt;
  64.                  if(x == u) break;
  65.              }
  66.          }
  67.      }
  68.  
  69.      void init(int nn){
  70.          n = nn;
  71.          ecnt = ; dfs_clock = scc_cnt = ;
  72.          memset(head,,sizeof(head));
  73.          memset(pre,,sizeof(pre));
  74.          memset(sccno,,sizeof(sccno));
  75.      }
  76.  
  77.      void addEdge(int x, int y){//x, y clash
  78.          to[ecnt] = y^; next[ecnt] = head[x]; head[x] = ecnt++;
  79.          to[ecnt] = x^; next[ecnt] = head[y]; head[y] = ecnt++;
  80.          //printf("%d<>%d\n",x,y);
  81.      }
  82.  
  83.      bool solve(){
  84.          for(int i = ; i < n; ++i)
  85.              if(!pre[i]) dfs(i);
  86.          for(int i = ; i < n; i += )
  87.              if(sccno[i] == sccno[i^]) return false;
  88.          return true;
  89.      }
  90.  
  91.      void bulid_select(){
  92.          T.init(scc_cnt);
  93.          for(int u = ; u < n; ++u){
  94.              for(int p = head[u]; p; p = next[p]){
  95.                  int &v = to[p];
  96.                  if(sccno[u] == sccno[v]) continue;
  97.                  T.addEdge(sccno[u], sccno[v]);
  98.              }
  99.          }
  100.          T.build();
  101.          memset(select,,sizeof(select));
  102.          for(int i = T.n - ; i > ; --i) {
  103.              int &x = T.order[i];
  104.              if(select[x] == -){
  105.                  select[x] = ;
  106.                  select[sccno[sccnox[x]^]] = ;
  107.              }
  108.          }
  109.      }
  110.  } G;
  111.  
  112.  const int MAXNN = ;
  113.  
  114.  int a1[MAXNN], b1[MAXNN], a2[MAXNN], b2[MAXNN], a3[MAXNN], b3[MAXNN], c[MAXNN];
  115.  
  116.  inline bool clash(int beg1, int end1, int beg2, int end2){
  117.      if(end1 <= beg2 || end2 <= beg1) return false;
  118.      return true;
  119.  }
  120.  
  121.  int main(){
  122.      int n;
  123.      while(scanf("%d", &n)!=EOF){
  124.          G.init(n*);
  125.          for(int i = ; i < n; ++i) {
  126.              scanf("%d:%d %d:%d %d", &a1[i], &a2[i], &b1[i], &b2[i], &c[i]);
  127.              a3[i] = a1[i] *  + a2[i];
  128.              b3[i] = b1[i] *  + b2[i];
  129.          }
  130.          for(int i = ; i < n; ++i) for(int j = i+; j < n; ++j) if(i != j){
  131.              if(clash(a3[i], a3[i] + c[i], a3[j], a3[j] + c[j])) G.addEdge(i*, j*);
  132.              if(clash(a3[i], a3[i] + c[i], b3[j] - c[j], b3[j])) G.addEdge(i*, j*+);
  133.              if(clash(b3[i] - c[i], b3[i], a3[j], a3[j] + c[j])) G.addEdge(i*+, j*);
  134.              if(clash(b3[i] - c[i], b3[i], b3[j] - c[j], b3[j])) G.addEdge(i*+, j*+);
  135.          }
  136.          if(G.solve()) printf("YES\n");
  137.          else {printf("NO\n"); continue;}
  138.          G.bulid_select();
  139.          for(int i = ; i < n; ++i){
  140.              //printf("%d %d\n",i*2,G.sccno[i*2]);
  141.              if(G.select[G.sccno[i*]]){
  142.                  b1[i] = a1[i]; b2[i] = a2[i] + c[i];
  143.                  while(b2[i] >= ) ++b1[i], b2[i] -= ;
  144.              } else {
  145.                  a1[i] = b1[i]; a2[i] = b2[i] - c[i];
  146.                  while(a2[i] < ) --a1[i], a2[i] += ;
  147.              }
  148.              printf("%02d:%02d %02d:%02d\n", a1[i], a2[i], b1[i], b2[i]);
  149.          }
  150.      }
  151.      return ;
  152.  }

141MS

后记：

上面那个输出解的方法略显麻烦了，看到了一个新的方法，就是若sccno[2i]＜sccno[2i+1]就选i，否则选i的对立面，不用求拓扑排序了。省了一大堆代码量。

证明：

对于任意点a、b（设~a、~b分别为他们的对立面），rank为他们的拓扑序号
那么我们选择a当且仅当rank[a]＞rank[~a]，选择b当且仅当rank[b]＞rank[~b]
假设a和b有矛盾，那么有a到~b有路径（由对称性b到~a有路径），则rank[~b]＞rank[a]（rank[~a]＞rank[b]）
联合上式，有rank[a]＞rank[~a]＞rank[b]＞rank[~b]＞rank[a]矛盾
假设不成立，a和b没有矛盾
所以对于任意点a，选择a当且仅当rank[a]＞rank[~a]是合法的

当我们用tarjan求强联通分量的时候，实际上求出来的scnno[]的值便是一个逆序的拓扑排序值。
仔细想想这没什么问题，因为再tarjan求强联通分量的时候，对于任意点x，它的后继顶点一定会比x先编入强联通分量，那么x的后继一定要比x的scnno值要小。

代码（125MS）：

  1.  #include <cstdio>
  2.  #include <cstring>
  3.  
  4.  const int MAXN = *;
  5.  const int MAXM = MAXN * MAXN * ;
  6.  
  7.  struct TwoSAT{
  8.      int St[MAXN], c;
  9.      int n, ecnt, dfs_clock, scc_cnt;
  10.      int head[MAXN], sccno[MAXN], pre[MAXN], lowlink[MAXN];
  11.      int next[MAXM], to[MAXM];
  12.  
  13.      void dfs(int u){
  14.          lowlink[u] = pre[u] = ++dfs_clock;
  15.          St[++c] = u;
  16.          for(int p = head[u]; p; p = next[p]){
  17.              int &v = to[p];
  18.              if(!pre[v]){
  19.                  dfs(v);
  20.                  if(lowlink[u] > lowlink[v]) lowlink[u] = lowlink[v];
  21.              }else if(!sccno[v]){
  22.                  if(lowlink[u] > pre[v]) lowlink[u] = pre[v];
  23.              }
  24.          }
  25.          if(lowlink[u] == pre[u]){
  26.              ++scc_cnt;
  27.              while(true){
  28.                  int x = St[c--];
  29.                  sccno[x] = scc_cnt;
  30.                  if(x == u) break;
  31.              }
  32.          }
  33.      }
  34.  
  35.      void init(int nn){
  36.          n = nn;
  37.          ecnt = ; dfs_clock = scc_cnt = ;
  38.          memset(head,,sizeof(head));
  39.          memset(pre,,sizeof(pre));
  40.          memset(sccno,,sizeof(sccno));
  41.      }
  42.  
  43.      void addEdge(int x, int y){//x, y clash
  44.          to[ecnt] = y^; next[ecnt] = head[x]; head[x] = ecnt++;
  45.          to[ecnt] = x^; next[ecnt] = head[y]; head[y] = ecnt++;
  46.          //printf("%d<>%d\n",x,y);
  47.      }
  48.  
  49.      bool solve(){
  50.          for(int i = ; i < n; ++i)
  51.              if(!pre[i]) dfs(i);
  52.          for(int i = ; i < n; i += )
  53.              if(sccno[i] == sccno[i^]) return false;
  54.          return true;
  55.      }
  56.  } G;
  57.  
  58.  const int MAXNN = ;
  59.  
  60.  int a1[MAXNN], b1[MAXNN], a2[MAXNN], b2[MAXNN], a3[MAXNN], b3[MAXNN], c[MAXNN];
  61.  
  62.  inline bool clash(int beg1, int end1, int beg2, int end2){
  63.      if(end1 <= beg2 || end2 <= beg1) return false;
  64.      return true;
  65.  }
  66.  
  67.  int main(){
  68.      int n;
  69.      while(scanf("%d", &n)!=EOF){
  70.          G.init(n*);
  71.          for(int i = ; i < n; ++i) {
  72.              scanf("%d:%d %d:%d %d", &a1[i], &a2[i], &b1[i], &b2[i], &c[i]);
  73.              a3[i] = a1[i] *  + a2[i];
  74.              b3[i] = b1[i] *  + b2[i];
  75.          }
  76.          for(int i = ; i < n; ++i) for(int j = i+; j < n; ++j) if(i != j){
  77.              if(clash(a3[i], a3[i] + c[i], a3[j], a3[j] + c[j])) G.addEdge(i*, j*);
  78.              if(clash(a3[i], a3[i] + c[i], b3[j] - c[j], b3[j])) G.addEdge(i*, j*+);
  79.              if(clash(b3[i] - c[i], b3[i], a3[j], a3[j] + c[j])) G.addEdge(i*+, j*);
  80.              if(clash(b3[i] - c[i], b3[i], b3[j] - c[j], b3[j])) G.addEdge(i*+, j*+);
  81.          }
  82.          if(G.solve()) printf("YES\n");
  83.          else {printf("NO\n"); continue;}
  84.          for(int i = ; i < n; ++i){
  85.              if(G.sccno[i * ] < G.sccno[i *  + ]) {
  86.                  b1[i] = a1[i]; b2[i] = a2[i] + c[i];
  87.                  while(b2[i] >= ) ++b1[i], b2[i] -= ;
  88.              } else {
  89.                  a1[i] = b1[i]; a2[i] = b2[i] - c[i];
  90.                  while(a2[i] < ) --a1[i], a2[i] += ;
  91.              }
  92.              printf("%02d:%02d %02d:%02d\n", a1[i], a2[i], b1[i], b2[i]);
  93.          }
  94.      }
  95.      return ;
  96.  }

POJ 3683 Priest John's Busiest Day（2-SAT 并输出解）的更多相关文章

  1. POJ 3683 Priest John's Busiest Day（2-SAT+方案输出）

    Priest John's Busiest Day Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 10010   Accep ...

  2. POJ 3683 Priest John's Busiest Day / OpenJ_Bailian 3788 Priest John's Busiest Day(2-sat问题)

    POJ 3683 Priest John's Busiest Day / OpenJ_Bailian 3788 Priest John's Busiest Day(2-sat问题) Descripti ...

  3. POJ 3683 Priest John's Busiest Day （2-SAT）

    Priest John's Busiest Day Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 6900   Accept ...

  4. poj - 3683 - Priest John's Busiest Day（2-SAT）

    题意:有N场婚礼,每场婚礼的开始时间为Si,结束时间为Ti,每场婚礼有个仪式,历时Di,这个仪式要么在Si时刻开始,要么在Ti-Di时刻开始,问能否安排每场婚礼举行仪式的时间,使主持人John能参加所 ...

  5. POJ 3683 Priest John's Busiest Day (2-SAT)

    题意:有n对新人要在同一天结婚.结婚时间为Ti到Di,这里有时长为Si的一个仪式需要神父出席.神父可以在Ti-(Ti+Si)这段时间出席也可以在(Di-Si)-Si这段时间.问神父能否出席所有仪式,如 ...

  6. POJ 3683 Priest John's Busiest Day (2-SAT，常规)

    题意: 一些人要在同一天进行婚礼,但是牧师只有1个,每一对夫妻都有一个时间范围[s , e]可供牧师选择,且起码要m分钟才主持完毕,但是要么就在 s 就开始,要么就主持到刚好 e 结束.因为人数太多了 ...

  7. POJ 3683 Priest John's Busiest Day

    2-SAT简单题,判断一下两个开区间是否相交 #include<cstdio> #include<cstring> #include<cmath> #include ...

  8. POJ 3683 Priest John's Busiest Day[2-SAT 构造解]

    题意: $n$对$couple$举行仪式,有两个时间段可以选择,问是否可以不冲突举行完,并求方案 两个时间段选择对应一真一假,对于有时间段冲突冲突的两人按照$2-SAT$的规则连边(把不冲突的时间段连 ...

  9. POJ 3683 Priest John's Busiest Day 【2-Sat】

    这是一道裸的2-Sat,只要考虑矛盾条件的判断就好了. 矛盾判断: 对于婚礼现场 x 和 y,x 的第一段可以和 y 的第一段或者第二段矛盾,同理,x 的第二段可以和 y 的第一段或者第二段矛盾,条件 ...

随机推荐

  1. 20145330《Java程序设计》第四次实验报告

    20145330<Java程序设计>第四次实验报告 实验四 Android环境搭建 实验内容 1.搭建Android环境 2.运行Android 3.修改代码,能输出学号 实验步骤 搭建A ...

  2. nmap的script参数列表

    在新的nmap版本中,添加了script功能的使用.在nmap的安装目录的share/nmap/scripts中,已经有将61个写好的脚本提供. 具体的用法可以参考:http://nmap.org/b ...

  3. ios修改产品名

    在创建项目的时候,会设置一个项目名,以后生成的APP名字也就是这个了,但由于某种原因,我想修改APP名字,也就是屏幕程序图标下面显示的那个,这该怎么办呢? 下面有三种方法都可以: 修改Product ...

  4. NodeJS学习笔记之Connect中间件模块（一）

    NodeJS学习笔记之Connect中间件模块(一) http://www.jb51.net/article/60430.htm NodeJS学习笔记之Connect中间件模块(二) http://w ...

  5. WebApp之Meta标签

    <meta name="apple-touch-fullscreen" content="yes">"添加到主屏幕“后,全屏显示 < ...

  6. cookie的设置、获取和删除封装

    在我们为了去完成数据储存,有时会用到cookie,简单封装一下cookie: <!DOCTYPE html> <html lang="en"> <he ...

  7. Android-----工程文件目录介绍

  8. Jquery，YUI这个著名js库名称作用的理解

    看廖雪峰大神的教程,其中讲到变量作用域问题.在命名空间中,写到:因为全局变量绑到了window上,不同的js文件访问相同全局变量或者定义了相同名字的顶层函数,都会造成命名冲突,并且很难被发现. 减少冲 ...

  9. Vim 常见操作

    1.复制所有内容 按照此顺序敲即可:gg(光标定位到文件头)  V(选中该行) G(选中该行开始到最后一行结尾) y(复制选中内容) 2.粘贴所有内容 正常模式下,敲p即可,如果遇到粘贴内容不全,需要 ...

  10. NeuSoft(2)添加系统调用

    1.下载内核 apt-get install linux-source 在/usr/src下 2.解压内核 cd /usr/src tar -jxvf linux-source-3.2.0.tar.b ...