Loading [Contrib]/a11y/accessibility-menu.js
$\newcommand{\O}{\mathrm{O}}$ My Algorithm : kopricky アルゴリズムライブラリ

kopricky アルゴリズムライブラリ

Bipartite Matching Algorithm

コードについての説明(個人的メモ)

二部グラフの最大マッチングを求めるアルゴリズム. 例えば $n$ 人の人と $m$ 個の仕事があり, それぞれの人について行うことのできる仕事が分かっていて $1$ 人が高々 $1$ つの仕事を行うというような状況を考えます. このとき最大いくつの仕事を行うことができるか? というような問題をこのアルゴリズムを用いて高速に解くことができます.
人と仕事に対応する $n + m$ 個の頂点と仕事ができるできないに対応する 人 → 仕事 の有向辺からなる二部グラフを考える. ここで超頂点 $s,t$ を source と sink として $s$ → 人, 仕事 → $t$ の辺を先ほどの二部グラフに加えた有向グラフ上で $s$ → $t$ の最大流問題を解く. このときの答えが問題の答えと一致している.
最大流の計算に Dinic 法 を用いると時間計算量は普通に考えると $\O (n^2 m)$ であるがこの問題に対しては $\O(m \sqrt{n})$ と高速に動作することが分かっている.
(注)
以下の実装で $t$ に向かう辺のみは level の下がる向きに進んでおり厳密には最短とは限らない増加道を通って更新をしているがそのような辺が saturated した後にできる逆辺はそれ以降の $s$ -> $t$ 増加道に関与しないため "増加道がフェーズごとに狭義単調増加すること" および "同じフェーズで saturated されてできた逆辺を通る最短の増加道は存在しないこと" が保たれ正当である.

(コンストラクタ)
BM$(u, v)$ : $u$ が人の数, $v$ が仕事の数
(関数)
add_edge$(u,v)$ : 人 $u$ と仕事 $v$ の間に辺を追加する
solve(): 最大マッチングを求める(メンバ変数 asg に結果が格納される)
minimum_vertex_cover(): (solve() を呼び出した後に呼び出すことで)最小頂点被覆を求める

時間計算量: $\O (m \sqrt{n})$

コード

  1. class BM {
  2. private:
  3. const int U, V;
  4. vector<vector<int> > G;
  5. vector<int> level, que, prv, rasg;
  6. int bfs(){
  7. int last = -1;
  8. fill(level.begin(), level.end(), -1);
  9. int qh = 0, qt = 0;
  10. for(int i = 0; i < U; ++i){
  11. if(asg[i] < 0) level[i] = 0, que[qt++] = i, prv[i] = -1;
  12. }
  13. while(qh < qt){
  14. const int u = que[qh++];
  15. if(u >= U){
  16. const int v = rasg[u - U];
  17. if(v >= 0) level[v] = level[u] + 1, que[qt++] = v, prv[v] = u;
  18. else last = u;
  19. }else{
  20. for(const int v : G[u]){
  21. if(level[v] < 0) level[v] = level[u] + 1, que[qt++] = v, prv[v] = u;
  22. }
  23. }
  24. }
  25. return last;
  26. }
  27. bool dfs(const int u){
  28. const int tmp = level[u];
  29. level[u] = -1;
  30. if(u >= U){
  31. if(rasg[u - U] < 0) return true;
  32. else return dfs(rasg[u - U]);
  33. }else{
  34. for(const int v : G[u]){
  35. if(tmp < level[v]){
  36. if(dfs(v)){
  37. asg[u] = v - U, rasg[v - U] = u;
  38. return true;
  39. }
  40. }
  41. }
  42. }
  43. return false;
  44. }
  45. void hint_search(int cur, int& flow){
  46. ++flow;
  47. while(cur >= 0){
  48. level[cur] = -1;
  49. if(cur >= U) asg[prv[cur]] = cur - U, rasg[cur - U] = prv[cur];
  50. cur = prv[cur];
  51. }
  52. }
  53. public:
  54. BM(const int u, const int v)
  55. : U(u), V(v), G(U + V), level(U + V), que(U + V), prv(U + V), rasg(V, -1), asg(U, -1){}
  56. void add_edge(const int from, const int to){
  57. G[from].push_back(U + to);
  58. }
  59. // asg に左側頂点がどの右側頂点とマッチングされるかが格納される
  60. vector<int> asg;
  61. int solve(){
  62. int flow = 0;
  63. for(;;){
  64. const int cur = bfs();
  65. if(cur < 0) break;
  66. hint_search(cur, flow);
  67. for(int i = 0; i < U; ++i){
  68. if(asg[i] < 0) flow += dfs(i);
  69. }
  70. }
  71. return flow;
  72. }
  73. // solve() を呼び出した後に呼び出す(左側頂点 i の添字は U(右側頂点数) + i とする)
  74. vector<int> minimum_vertex_cover(){
  75. vector<int> mvc;
  76. for(int i = 0; i < U; ++i){
  77. if(level[i] < 0) mvc.push_back(i);
  78. }
  79. for(int i = U; i < U + V; ++i){
  80. if(level[i] >= 0) mvc.push_back(i);
  81. }
  82. return mvc;
  83. }
  84. };

verify 用の問題

yosupo さんの Library-Checker : Matching on Bipartite Graph 提出コード