test/yosupo/bipartitematching.general_matching.test.cpp

• View this file on GitHub
• Last update: 2024-09-01 12:11:12+09:00
• Problem: https://judge.yosupo.jp/problem/bipartitematching
• Link: View error logs on GitHub Actions

Depends on

• 一般グラフの最大マッチング (src/Graph/general_matching.hpp)
• グラフ (src/Graph/Graph.hpp)
• CSR 表現を用いた二次元配列 他 (src/Internal/ListRange.hpp)

Code

// competitive-verifier: PROBLEM https://judge.yosupo.jp/problem/bipartitematching
// competitive-verifier: TLE 1
// competitive-verifier: MLE 64
#include <iostream>
#include "src/Graph/Graph.hpp"
#include "src/Graph/general_matching.hpp"
using namespace std;
signed main() {
 cin.tie(0);
 ios::sync_with_stdio(0);
 int L, R, M;
 cin >> L >> R >> M;
 Graph g(L + R, M);
 for (int i= 0; i < M; ++i) cin >> g[i], g[i].second+= L;
 auto [ans, _]= general_matching(g);
 cout << ans.size() << '\n';
 for (int i: ans) {
  auto [l, r]= g[i];
  cout << l << " " << r - L << '\n';
 }
 return 0;
}

#line 1 "test/yosupo/bipartitematching.general_matching.test.cpp"
// competitive-verifier: PROBLEM https://judge.yosupo.jp/problem/bipartitematching
// competitive-verifier: TLE 1
// competitive-verifier: MLE 64
#include <iostream>
#line 2 "src/Internal/ListRange.hpp"
#include <vector>
#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 2 "src/Graph/general_matching.hpp"
#include <cassert>
#line 4 "src/Graph/general_matching.hpp"
// {matching edge ids, partner(-1 if unmatched)}
std::pair<std::vector<int>, std::vector<int>> general_matching(const Graph &g, std::vector<int> partner= {}) {
 auto adj= g.adjacency_vertex(0);
 const int n= adj.size();
 std::vector<int> q, z(n), p(n);
 if (partner.empty()) partner.assign(n, -1);
 assert((int)partner.size() == n);
 std::vector<Edge> fs(n);
 auto rematch= [&](auto &&rc, int u, int v) -> void {
  int w= partner[u];
  if (partner[u]= v; w == -1 || partner[w] != u) return;
  if (auto [x, y]= fs[u]; y == -1) rc(rc, partner[w]= x, w);
  else rc(rc, x, y), rc(rc, y, x);
 };
 int t= 1;
 auto f= [&](auto &&rc, int x) -> int { return z[x] != t || p[x] == -1 ? x : (p[x]= rc(rc, p[x])); };
 auto check= [&](int r) {
  q.clear(), q.push_back(r), fs[r]= {-1, -1}, z[r]= t, p[r]= -1;
  for (int i= 0; i < (int)q.size(); ++i) {
   int x= q[i];
   for (int y: adj[x]) {
    if (y == r) continue;
    if (partner[y] == -1) return rematch(rematch, partner[y]= x, y), true;
    if (z[y] == t) {
     int u= f(f, x), v= f(f, y), w= r;
     if (u == v) continue;
     for (; u != r || v != r; fs[u]= {x, y}, u= f(f, fs[partner[u]].first)) {
      if (v != r) std::swap(u, v);
      if (fs[u].first == x && fs[u].second == y) {
       w= u;
       break;
      }
     }
     for (int a: {f(f, x), f(f, y)})
      for (; a != w; a= f(f, fs[partner[a]].first)) z[a]= t, p[a]= w, q.push_back(a);
    } else if (z[partner[y]] != t) fs[y]= {-1, -1}, fs[partner[y]]= {x, -1}, z[partner[y]]= t, p[partner[y]]= y, q.push_back(partner[y]);
   }
  }
  return false;
 };
 for (int r= n; r--;)
  if (partner[r] == -1) t+= check(r);
 q.clear();
 for (int i= 0, e= g.edge_size(); i < e; ++i)
  if (auto [u, v]= g[i]; partner[u] == v && z[u] >= 0) q.push_back(i), z[u]= z[v]= -1;
 return {q, partner};
}
#line 7 "test/yosupo/bipartitematching.general_matching.test.cpp"
using namespace std;
signed main() {
 cin.tie(0);
 ios::sync_with_stdio(0);
 int L, R, M;
 cin >> L >> R >> M;
 Graph g(L + R, M);
 for (int i= 0; i < M; ++i) cin >> g[i], g[i].second+= L;
 auto [ans, _]= general_matching(g);
 cout << ans.size() << '\n';
 for (int i: ans) {
  auto [l, r]= g[i];
  cout << l << " " << r - L << '\n';
 }
 return 0;
}

Test cases

Env	Name	Status	Elapsed	Memory
g++-13	augmented_cycle_00	AC	896 ms	13 MB
g++-13	augmented_cycle_01	AC	510 ms	13 MB
g++-13	augmented_cycle_02	TLE	1005 ms	0 MB
g++-13	cycle_00	AC	66 ms	13 MB
g++-13	cycle_01	AC	76 ms	13 MB
g++-13	example_00	AC	6 ms	3 MB
g++-13	issue1068_00	AC	6 ms	3 MB
g++-13	issue1068_large_00	AC	38 ms	11 MB
g++-13	issue1068_large_01	AC	38 ms	11 MB
g++-13	issue1068_large_02	AC	41 ms	11 MB
g++-13	issue1068_large_03	AC	39 ms	11 MB
g++-13	issue1068_large_04	AC	38 ms	11 MB
g++-13	issue1068_large_05	AC	41 ms	11 MB
g++-13	issue1124_00	AC	45 ms	9 MB
g++-13	issue1124_01	AC	84 ms	9 MB
g++-13	issue1124_02	AC	31 ms	9 MB
g++-13	kuhn_killer_00	AC	42 ms	13 MB
g++-13	line_00	AC	44 ms	13 MB
g++-13	line_01	AC	46 ms	13 MB
g++-13	line_random_00	AC	69 ms	13 MB
g++-13	line_random_01	AC	79 ms	13 MB
g++-13	many_paths_00	AC	302 ms	12 MB
g++-13	many_paths_01	AC	440 ms	12 MB
g++-13	many_paths_02	AC	194 ms	12 MB
g++-13	many_smalls_00	AC	37 ms	10 MB
g++-13	many_smalls_01	AC	36 ms	10 MB
g++-13	max_random_00	AC	50 ms	13 MB
g++-13	max_random_01	AC	50 ms	13 MB
g++-13	max_random_02	AC	51 ms	13 MB
g++-13	random_00	AC	10 ms	6 MB
g++-13	random_01	AC	28 ms	8 MB
g++-13	random_02	AC	26 ms	8 MB
g++-13	random_03	AC	27 ms	7 MB
g++-13	random_04	AC	13 ms	7 MB
g++-13	random_05	AC	15 ms	6 MB
g++-13	random_06	AC	43 ms	10 MB
g++-13	random_07	AC	15 ms	7 MB
g++-13	random_08	AC	57 ms	8 MB
g++-13	random_09	AC	30 ms	8 MB
g++-13	unique_matching_00	AC	890 ms	13 MB
g++-13	unique_matching_01	AC	953 ms	12 MB
g++-13	unique_matching_02	AC	974 ms	12 MB
g++-13	unique_matching_03	AC	659 ms	12 MB
g++-13	unique_matching_04	AC	506 ms	12 MB
clang++-18	augmented_cycle_00	AC	930 ms	13 MB
clang++-18	augmented_cycle_01	AC	522 ms	12 MB
clang++-18	augmented_cycle_02	TLE	1005 ms	0 MB
clang++-18	cycle_00	AC	70 ms	13 MB
clang++-18	cycle_01	AC	80 ms	13 MB
clang++-18	example_00	AC	6 ms	4 MB
clang++-18	issue1068_00	AC	5 ms	4 MB
clang++-18	issue1068_large_00	AC	37 ms	11 MB
clang++-18	issue1068_large_01	AC	37 ms	11 MB
clang++-18	issue1068_large_02	AC	41 ms	11 MB
clang++-18	issue1068_large_03	AC	39 ms	11 MB
clang++-18	issue1068_large_04	AC	40 ms	11 MB
clang++-18	issue1068_large_05	AC	44 ms	11 MB
clang++-18	issue1124_00	AC	48 ms	9 MB
clang++-18	issue1124_01	AC	90 ms	9 MB
clang++-18	issue1124_02	AC	31 ms	9 MB
clang++-18	kuhn_killer_00	AC	42 ms	13 MB
clang++-18	line_00	AC	48 ms	13 MB
clang++-18	line_01	AC	45 ms	13 MB
clang++-18	line_random_00	AC	68 ms	13 MB
clang++-18	line_random_01	AC	81 ms	13 MB
clang++-18	many_paths_00	AC	320 ms	12 MB
clang++-18	many_paths_01	AC	467 ms	12 MB
clang++-18	many_paths_02	AC	216 ms	12 MB
clang++-18	many_smalls_00	AC	36 ms	10 MB
clang++-18	many_smalls_01	AC	37 ms	10 MB
clang++-18	max_random_00	AC	50 ms	13 MB
clang++-18	max_random_01	AC	50 ms	13 MB
clang++-18	max_random_02	AC	50 ms	13 MB
clang++-18	random_00	AC	10 ms	6 MB
clang++-18	random_01	AC	28 ms	8 MB
clang++-18	random_02	AC	27 ms	8 MB
clang++-18	random_03	AC	27 ms	7 MB
clang++-18	random_04	AC	13 ms	7 MB
clang++-18	random_05	AC	15 ms	6 MB
clang++-18	random_06	AC	44 ms	10 MB
clang++-18	random_07	AC	13 ms	7 MB
clang++-18	random_08	AC	57 ms	8 MB
clang++-18	random_09	AC	30 ms	8 MB
clang++-18	unique_matching_00	AC	971 ms	13 MB
clang++-18	unique_matching_01	AC	999 ms	12 MB
clang++-18	unique_matching_02	TLE	1005 ms	0 MB
clang++-18	unique_matching_03	AC	697 ms	12 MB
clang++-18	unique_matching_04	AC	555 ms	12 MB

Back to top page