test/hackerrank/drawing-rectangles.test.cpp

• View this file on GitHub
• Last update: 2025-08-11 23:31:06+09:00
• Problem: https://www.hackerrank.com/contests/university-codesprint-4/challenges/drawing-rectangles

Depends on

• 区間をsetで管理するやつ (src/DataStructure/RangeSet.hpp)
• 二部グラフ (src/Graph/BipartiteGraph.hpp)
• Dulmage-Mendelsohn 分解 (src/Graph/DulmageMendelsohn.hpp)
• グラフ (src/Graph/Graph.hpp)
• CSR 表現を用いた二次元配列 他 (src/Internal/ListRange.hpp)

Code

// competitive-verifier: PROBLEM https://www.hackerrank.com/contests/university-codesprint-4/challenges/drawing-rectangles
// competitive-verifier: TLE 0.5
// competitive-verifier: MLE 128
#include <iostream>
#include <vector>
#include "src/Graph/BipartiteGraph.hpp"
#include "src/Graph/DulmageMendelsohn.hpp"
#include "src/DataStructure/RangeSet.hpp"
using namespace std;
signed main() {
 cin.tie(0);
 ios::sync_with_stdio(0);
 static constexpr int N= 300'001;
 vector<RangeSet<int>> ys(N);
 int n;
 cin >> n;
 for (int i= 0; i < n; ++i) {
  int x1, y1, x2, y2;
  cin >> x1 >> y1 >> x2 >> y2;
  for (int x= x1; x <= x2; ++x) ys[x].insert(y1, y2);
 }
 BipartiteGraph bm(N, N);
 for (int x= 0; x < N; ++x)
  for (auto [y1, y2]: ys[x])
   for (int y= y1; y <= y2; ++y) bm.add_edge(x, y + N);
 DulmageMendelsohn dm(bm);
 auto vertex_cover= dm.min_vertex_cover();
 vector<int> l, r;
 for (int i: vertex_cover)
  if (i < N) l.push_back(i);
  else r.push_back(i - N);
 int L= l.size(), R= r.size();
 cout << L + R << '\n';
 cout << L << '\n';
 for (int i= 0; i < L; ++i) cout << (i ? " " : "") << l[i];
 cout << '\n';
 cout << R << '\n';
 for (int i= 0; i < R; ++i) cout << (i ? " " : "") << r[i];
 cout << '\n';
 return 0;
}

#line 1 "test/hackerrank/drawing-rectangles.test.cpp"
// competitive-verifier: PROBLEM https://www.hackerrank.com/contests/university-codesprint-4/challenges/drawing-rectangles
// competitive-verifier: TLE 0.5
// competitive-verifier: MLE 128
#include <iostream>
#include <vector>
#line 2 "src/Graph/BipartiteGraph.hpp"
#include <cassert>
#include <tuple>
#include <algorithm>
#line 4 "src/Internal/ListRange.hpp"
#include <iterator>
#include <type_traits>
#define _LR(name, IT, CT) \
 template <class T> struct name { \
  using Iterator= typename std::vector<T>::IT; \
  Iterator bg, ed; \
  Iterator begin() const { return bg; } \
  Iterator end() const { return ed; } \
  size_t size() const { return std::distance(bg, ed); } \
  CT &operator[](int i) const { return bg[i]; } \
 }
_LR(ListRange, iterator, T);
_LR(ConstListRange, const_iterator, const T);
#undef _LR
template <class T> struct CSRArray {
 std::vector<T> dat;
 std::vector<int> p;
 size_t size() const { return p.size() - 1; }
 ListRange<T> operator[](int i) { return {dat.begin() + p[i], dat.begin() + p[i + 1]}; }
 ConstListRange<T> operator[](int i) const { return {dat.cbegin() + p[i], dat.cbegin() + p[i + 1]}; }
};
template <template <class> class F, class T> std::enable_if_t<std::disjunction_v<std::is_same<F<T>, ListRange<T>>, std::is_same<F<T>, ConstListRange<T>>, std::is_same<F<T>, CSRArray<T>>>, std::ostream &> operator<<(std::ostream &os, const F<T> &r) {
 os << '[';
 for (int _= 0, __= r.size(); _ < __; ++_) os << (_ ? ", " : "") << r[_];
 return os << ']';
}
#line 3 "src/Graph/Graph.hpp"
struct Edge: std::pair<int, int> {
 using std::pair<int, int>::pair;
 Edge &operator--() { return --first, --second, *this; }
 int to(int v) const { return first ^ second ^ v; }
 friend std::istream &operator>>(std::istream &is, Edge &e) { return is >> e.first >> e.second, is; }
};
struct Graph: std::vector<Edge> {
 size_t n;
 Graph(size_t n= 0, size_t m= 0): vector(m), n(n) {}
 size_t vertex_size() const { return n; }
 size_t edge_size() const { return size(); }
 size_t add_vertex() { return n++; }
 size_t add_edge(int s, int d) { return emplace_back(s, d), size() - 1; }
 size_t add_edge(Edge e) { return emplace_back(e), size() - 1; }
#define _ADJ_FOR(a, b) \
 for (auto [u, v]: *this) a; \
 for (size_t i= 0; i < n; ++i) p[i + 1]+= p[i]; \
 for (int i= size(); i--;) { \
  auto [u, v]= (*this)[i]; \
  b; \
 }
#define _ADJ(a, b) \
 vector<int> p(n + 1), c(size() << !dir); \
 if (!dir) { \
  _ADJ_FOR((++p[u], ++p[v]), (c[--p[u]]= a, c[--p[v]]= b)) \
 } else if (dir > 0) { \
  _ADJ_FOR(++p[u], c[--p[u]]= a) \
 } else { \
  _ADJ_FOR(++p[v], c[--p[v]]= b) \
 } \
 return {c, p}
 CSRArray<int> adjacency_vertex(int dir) const { _ADJ(v, u); }
 CSRArray<int> adjacency_edge(int dir) const { _ADJ(i, i); }
#undef _ADJ
#undef _ADJ_FOR
};
#line 6 "src/Graph/BipartiteGraph.hpp"
// [0, L) is left, [L, n) is right
struct BipartiteGraph: Graph {
 size_t L;
 BipartiteGraph() {}
 BipartiteGraph(size_t L, size_t R, size_t m= 0): Graph(L + R, m), L(L) {}
 size_t left_size() const { return L; }
 size_t right_size() const { return this->n - L; }
};
std::vector<int> paint_two_colors(const CSRArray<int> &adj) {
 const int n= adj.size();
 std::vector<int> col(n, -1);
 for (int s= n; s--;)
  if (col[s] == -1) {
   std::vector<int> q= {s};
   for (int i= col[s]= 0, v; i < (int)q.size(); ++i)
    for (int u: adj[v= q[i]])
     if (int c= col[v]; col[u] == c) return {};
     else if (col[u] == -1) col[u]= c ^ 1, q.push_back(u);
  }
 return col;
}
std::vector<int> paint_two_colors(const Graph &g) { return paint_two_colors(g.adjacency_vertex(0)); }
// { BipartiteGraph , original to new, new to original }
// {{},{},{}} if not bipartite
std::tuple<BipartiteGraph, std::vector<int>, std::vector<int>> graph_to_bipartite(const Graph &g, std::vector<int> color= {}) {
 if (color.empty()) color= paint_two_colors(g);
 if (color.empty()) return {};
 const int n= g.vertex_size(), m= g.edge_size();
 std::vector<int> a(n), b(n);
 int l= 0, r= n;
 for (int i= n; i--;) b[a[i]= color[i] ? --r : l++]= i;
 BipartiteGraph bg(l, n - l, m);
 for (int i= m; i--;) {
  auto [u, v]= g[i];
  bg[i]= std::minmax(a[u], a[v]);
 }
 return {bg, a, b};
}
namespace _bg_internal {
std::vector<int> _bm(int L, const CSRArray<int> &adj, std::vector<int> &m) {
 std::vector<int> a, p, q(L);
 for (bool u= true; u;) {
  u= false, a.assign(L, -1), p.assign(L, -1);
  int t= 0;
  for (int l= L; l--;)
   if (m[l] == -1) q[t++]= a[l]= p[l]= l;
  for (int i= 0; i < t; ++i)
   if (int l= q[i], x; m[a[l]] == -1)
    for (int r: adj[l]) {
     if (x= m[r]; x == -1) {
      for (u= true; r != -1; l= p[l]) m[r]= l, std::swap(m[l], r);
      break;
     }
     if (p[x] == -1) a[q[t++]= x]= a[p[x]= l];
    }
 }
 return a;
}
}
template <bool lexical= false> std::pair<std::vector<int>, std::vector<int>> bipartite_matching(const BipartiteGraph &bg, std::vector<int> partner= {}) {
 const int L= bg.left_size(), M= bg.edge_size();
 if (partner.empty()) partner.assign(bg.vertex_size(), -1);
 assert(partner.size() == bg.vertex_size());
 {
  CSRArray<int> adj{std::vector<int>(M), std::vector<int>(L + 1)};
  for (auto [l, r]: bg) ++adj.p[l];
  for (int i= 0; i < L; ++i) adj.p[i + 1]+= adj.p[i];
  for (auto [l, r]: bg) adj.dat[--adj.p[l]]= r;
  if constexpr (lexical) {
   for (int l= L; l--;) std::sort(adj[l].begin(), adj[l].end());
   _bg_internal::_bm(L, adj, partner);
   std::vector<char> a(L, 1);
   for (int l= 0; l < L; ++l)
    if (int r= partner[l], v= l; r != -1) {
     std::vector<int> p(L, partner[v]= partner[r]= -1), c(adj.p.begin(), adj.p.begin() + L);
     for (p[v]= -2;;) {
      if (c[v] == adj.p[v + 1]) v= p[v];
      else if (int u= partner[r= adj.dat[c[v]++]]; u == -1) {
       for (; r != -1; v= p[v]) partner[r]= v, std::swap(partner[v], r);
       break;
      } else if (a[u] && p[u] == -1) p[u]= v, v= u;
     }
     a[l]= 0;
    }
  } else _bg_internal::_bm(L, adj, partner);
 }
 std::vector<int> c;
 std::vector<char> p(L);
 for (int i= 0; i < M; ++i)
  if (auto [l, r]= bg[i]; partner[l] == r && !p[l]) c.push_back(i), p[l]= 1;
 return {c, partner};
}
#line 3 "src/Graph/DulmageMendelsohn.hpp"
#include <numeric>
#line 5 "src/Graph/DulmageMendelsohn.hpp"
class DulmageMendelsohn {
 int L;
 std::vector<int> b, m, a;
 CSRArray<int> dag[2];
public:
 DulmageMendelsohn(const BipartiteGraph &bg): L(bg.left_size()) {
  auto adj= bg.adjacency_vertex(0);
  const int n= adj.size();
  m.assign(n, -1), b.assign(n, -3), a= _bg_internal::_bm(L, adj, m);
  std::vector<int> q(n - L);
  int t= 0, k= 0;
  for (int l= L; l--;)
   if (a[l] != -1)
    if (b[l]= -1; m[l] != -1) b[m[l]]= -1;
  for (int r= n; r-- > L;)
   if (m[r] == -1) b[q[t++]= r]= 0;
  for (int i= 0, r, w; i < t; ++i)
   for (int l: adj[r= q[i]])
    if (b[l] == -3)
     if (b[l]= 0, w= m[l]; w != -1 && b[w] == -3) b[q[t++]= w]= 0;
  t= 0;
  {
   std::vector<int> c(adj.p.begin(), adj.p.begin() + L);
   for (int l= L; l--;)
    if (int v= l; b[v] == -3)
     for (b[v]= -2; v >= 0;) {
      if (c[v] == adj.p[v + 1]) a[t++]= v, v= b[v];
      else if (int w= m[adj.dat[c[v]++]]; b[w] == -3) b[w]= v, v= w;
     }
  }
  for (int i= 0, e= 0, r; t--;)
   if (int s= a[t], p= m[s]; b[p] == -3)
    for (b[q[e++]= p]= b[s]= ++k; i < e; ++i)
     for (int l: adj[r= q[i]])
      if (b[m[l]] == -3) b[q[e++]= m[l]]= b[l]= k;
  ++k;
  for (int l= L; l--;)
   if (b[l] == -1)
    if (b[l]= k; m[l] != -1) b[m[l]]= k;
  a.assign(k + 2, 0);
  for (int i= n; i--;) ++a[b[i]];
  for (int i= 0; i <= k; ++i) a[i + 1]+= a[i];
  for (int i= n; i--;) m[--a[b[i]]]= i;
  Graph h(k + 1);
  for (auto [l, r]: bg)
   if (b[l] != b[r]) h.add_edge(b[l], b[r]);
  std::sort(h.begin(), h.end()), h.erase(std::unique(h.begin(), h.end()), h.end()), dag[0]= h.adjacency_vertex(1), dag[1]= h.adjacency_vertex(-1);
 }
 size_t size() const { return a.size() - 1; }
 ConstListRange<int> block(int k) const { return {m.cbegin() + a[k], m.cbegin() + a[k + 1]}; }
 int operator()(int i) const { return b[i]; }
 std::vector<int> min_vertex_cover(std::vector<int> ord= {}) const {
  if (ord.empty()) ord.resize(b.size()), std::iota(ord.begin(), ord.end(), 0);
  std::vector<char> z(size(), -1);
  std::vector<int> q(size()), vc;
  z[0]= 1, z.back()= 0;
  for (int v: ord) {
   int c= (v >= L), k= b[v], s= z[k];
   if (s == -1) {
    auto &adj= dag[z[q[0]= k]= s= !c];
    for (int i= 0, t= 1; i < t; ++i)
     for (int u: adj[q[i]])
      if (z[u] == -1) z[u]= s, q[t++]= u;
   }
   if (c ^ s) vc.push_back(v);
  }
  return vc;
 }
};
#line 3 "src/DataStructure/RangeSet.hpp"
#include <set>
#line 5 "src/DataStructure/RangeSet.hpp"
#include <limits>
#line 7 "src/DataStructure/RangeSet.hpp"
template <class Int, bool merge= true> class RangeSet {
 struct ClosedSection {
  Int l, r;
  Int length() const { return r - l + 1; }
  bool operator<(const ClosedSection &cs) const { return l < cs.l || (l == cs.l && r > cs.r); }
  operator bool() const { return l <= r; }
  friend std::ostream &operator<<(std::ostream &os, const ClosedSection &cs) { return cs ? os << "[" << cs.l << "," << cs.r << "]" : os << "∅"; }
 };
 std::set<ClosedSection> mp;
public:
 RangeSet() {
  constexpr Int INF= std::numeric_limits<Int>::max() / 2;
  mp.insert({INF, INF}), mp.insert({-INF, -INF});
 }
 ClosedSection covered_by(Int l, Int r) const {
  assert(l <= r);
  if (auto it= std::prev(mp.upper_bound(ClosedSection{l, l})); it->l <= l && r <= it->r) return *it;
  return {1, 0};
 }
 ClosedSection covered_by(Int x) const { return covered_by(x, x); }
 ClosedSection covered_by(const ClosedSection &cs) const { return covered_by(cs.l, cs.r); }
 size_t size() const { return mp.size() - 2; }
 auto begin() const { return std::next(mp.begin()); }
 auto end() const { return std::prev(mp.end()); }
 Int insert(Int l, Int r) {
  assert(l <= r);
  auto it= std::prev(mp.upper_bound(ClosedSection{l, l}));
  Int sum= 0, x= it->l, y= it->r;
  if (x <= l && r <= y) return sum;
  if (x <= l && l <= y + merge) sum+= y - (l= x) + 1, it= mp.erase(it);
  else std::advance(it, 1);
  for (; it->r < r; it= mp.erase(it)) sum+= it->r - it->l + 1;
  if (x= it->l, y= it->r; x - merge <= r && r <= y) sum+= (r= y) - x + 1, mp.erase(it);
  return mp.insert({l, r}), r - l + 1 - sum;
 }
 Int insert(Int x) { return insert(x, x); }
 Int insert(const ClosedSection &cs) { return insert(cs.l, cs.r); }
 Int erase(Int l, Int r) {
  assert(l <= r);
  auto it= std::prev(mp.upper_bound(ClosedSection{l, l}));
  Int sum= 0, x= it->l, y= it->r;
  if (x <= l && r <= y) {
   if (mp.erase(it); x < l) mp.insert({x, l - 1});
   if (r < y) mp.insert({r + 1, y});
   return r - l + 1;
  }
  if (x <= l && l <= y) {
   if (x < l) mp.insert({x, l - 1});
   sum+= y - l + 1, it= mp.erase(it);
  } else std::advance(it, 1);
  for (; it->r <= r; it= mp.erase(it)) sum+= it->r - it->l + 1;
  if (x= it->l, y= it->r; x <= r && r <= y)
   if (sum+= r - x + 1, mp.erase(it); r < y) mp.insert({r + 1, y});
  return sum;
 }
 Int erase(Int x) { return erase(x, x); }
 Int erase(const ClosedSection &cs) { return erase(cs.l, cs.r); }
 Int mex(Int x) const {
  auto cs= covered_by(x);
  return cs ? cs.r + 1 : x;
 }
 friend std::ostream &operator<<(std::ostream &os, const RangeSet &rs) {
  os << "[";
  for (auto it= rs.begin(); it != rs.end(); ++it) os << (it == rs.begin() ? "" : ",") << *it;
  return os << "]";
 }
};
#line 9 "test/hackerrank/drawing-rectangles.test.cpp"
using namespace std;
signed main() {
 cin.tie(0);
 ios::sync_with_stdio(0);
 static constexpr int N= 300'001;
 vector<RangeSet<int>> ys(N);
 int n;
 cin >> n;
 for (int i= 0; i < n; ++i) {
  int x1, y1, x2, y2;
  cin >> x1 >> y1 >> x2 >> y2;
  for (int x= x1; x <= x2; ++x) ys[x].insert(y1, y2);
 }
 BipartiteGraph bm(N, N);
 for (int x= 0; x < N; ++x)
  for (auto [y1, y2]: ys[x])
   for (int y= y1; y <= y2; ++y) bm.add_edge(x, y + N);
 DulmageMendelsohn dm(bm);
 auto vertex_cover= dm.min_vertex_cover();
 vector<int> l, r;
 for (int i: vertex_cover)
  if (i < N) l.push_back(i);
  else r.push_back(i - N);
 int L= l.size(), R= r.size();
 cout << L + R << '\n';
 cout << L << '\n';
 for (int i= 0; i < L; ++i) cout << (i ? " " : "") << l[i];
 cout << '\n';
 cout << R << '\n';
 for (int i= 0; i < R; ++i) cout << (i ? " " : "") << r[i];
 cout << '\n';
 return 0;
}

Test cases

Env	Name	Status	Elapsed	Memory
g++-13	00	AC	54 ms	56 MB
g++-13	01	AC	55 ms	56 MB
g++-13	02	AC	61 ms	56 MB
g++-13	03	AC	59 ms	56 MB
g++-13	04	AC	61 ms	60 MB
g++-13	05	AC	57 ms	56 MB
g++-13	06	AC	54 ms	56 MB
g++-13	07	AC	54 ms	56 MB
g++-13	08	AC	54 ms	56 MB
g++-13	09	AC	55 ms	56 MB
g++-13	10	AC	58 ms	57 MB
g++-13	11	AC	63 ms	60 MB
g++-13	12	AC	62 ms	60 MB
g++-13	13	AC	61 ms	60 MB
g++-13	14	AC	63 ms	60 MB
g++-13	15	AC	67 ms	60 MB
g++-13	16	AC	71 ms	60 MB
g++-13	17	AC	82 ms	60 MB
g++-13	18	AC	93 ms	61 MB
g++-13	19	AC	118 ms	62 MB
g++-13	20	AC	160 ms	69 MB
g++-13	21	AC	182 ms	72 MB
g++-13	22	AC	163 ms	70 MB
g++-13	23	AC	205 ms	75 MB
g++-13	24	AC	54 ms	56 MB
g++-13	25	AC	55 ms	56 MB
g++-13	26	AC	56 ms	56 MB
g++-13	27	AC	57 ms	56 MB
g++-13	28	AC	58 ms	56 MB
g++-13	29	AC	60 ms	57 MB
g++-13	30	AC	64 ms	57 MB
g++-13	31	AC	73 ms	58 MB
g++-13	32	AC	84 ms	59 MB
g++-13	33	AC	115 ms	62 MB
g++-13	34	AC	157 ms	69 MB
g++-13	35	AC	156 ms	69 MB
g++-13	36	AC	160 ms	69 MB
g++-13	37	AC	168 ms	70 MB
g++-13	38	AC	168 ms	70 MB
g++-13	39	AC	73 ms	57 MB
g++-13	40	AC	204 ms	57 MB
g++-13	41	AC	64 ms	56 MB
g++-13	42	AC	71 ms	57 MB
g++-13	43	AC	85 ms	58 MB
g++-13	44	AC	80 ms	58 MB
g++-13	45	AC	128 ms	58 MB
g++-13	46	AC	304 ms	58 MB
g++-13	47	AC	74 ms	58 MB
g++-13	48	AC	120 ms	57 MB
g++-13	49	AC	240 ms	56 MB
g++-13	50	AC	54 ms	56 MB
clang++-18	00	AC	51 ms	56 MB
clang++-18	01	AC	52 ms	56 MB
clang++-18	02	AC	57 ms	56 MB
clang++-18	03	AC	57 ms	56 MB
clang++-18	04	AC	60 ms	60 MB
clang++-18	05	AC	54 ms	56 MB
clang++-18	06	AC	53 ms	56 MB
clang++-18	07	AC	54 ms	56 MB
clang++-18	08	AC	53 ms	56 MB
clang++-18	09	AC	54 ms	56 MB
clang++-18	10	AC	55 ms	57 MB
clang++-18	11	AC	60 ms	60 MB
clang++-18	12	AC	61 ms	60 MB
clang++-18	13	AC	61 ms	60 MB
clang++-18	14	AC	61 ms	60 MB
clang++-18	15	AC	64 ms	60 MB
clang++-18	16	AC	69 ms	60 MB
clang++-18	17	AC	79 ms	60 MB
clang++-18	18	AC	90 ms	61 MB
clang++-18	19	AC	112 ms	62 MB
clang++-18	20	AC	151 ms	69 MB
clang++-18	21	AC	181 ms	72 MB
clang++-18	22	AC	159 ms	70 MB
clang++-18	23	AC	201 ms	75 MB
clang++-18	24	AC	52 ms	56 MB
clang++-18	25	AC	53 ms	56 MB
clang++-18	26	AC	54 ms	56 MB
clang++-18	27	AC	54 ms	56 MB
clang++-18	28	AC	56 ms	56 MB
clang++-18	29	AC	57 ms	56 MB
clang++-18	30	AC	59 ms	57 MB
clang++-18	31	AC	71 ms	58 MB
clang++-18	32	AC	82 ms	59 MB
clang++-18	33	AC	111 ms	62 MB
clang++-18	34	AC	152 ms	69 MB
clang++-18	35	AC	149 ms	68 MB
clang++-18	36	AC	160 ms	69 MB
clang++-18	37	AC	165 ms	70 MB
clang++-18	38	AC	165 ms	70 MB
clang++-18	39	AC	75 ms	57 MB
clang++-18	40	AC	257 ms	57 MB
clang++-18	41	AC	62 ms	57 MB
clang++-18	42	AC	72 ms	57 MB
clang++-18	43	AC	82 ms	58 MB
clang++-18	44	AC	82 ms	58 MB
clang++-18	45	AC	144 ms	58 MB
clang++-18	46	AC	367 ms	58 MB
clang++-18	47	AC	75 ms	58 MB
clang++-18	48	AC	139 ms	57 MB
clang++-18	49	AC	221 ms	56 MB
clang++-18	50	AC	53 ms	56 MB

Back to top page