图论--一般图带花树匹配--模板

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
using namespace std;
const int N = 250;
// 并查集维护
int belong[N];
int findb(int x) {return belong[x] == x ? x : belong[x] = findb(belong[x]);
}
void unit(int a, int b) {a = findb(a);b = findb(b);if (a != b) belong[a] = b;
}int n, match[N];
vector<int> e[N];
int Q[N], rear;
int _next[N], mark[N], vis[N];
// 朴素算法求某阶段中搜索树上两点x, y的最近公共祖先r
int LCA(int x, int y) {static int t = 0; t++;while (true) {if (x != -1) {x = findb(x); // 点要对应到对应的花上去  if (vis[x] == t)return x;vis[x] = t;if (match[x] != -1)x = _next[match[x]];else x = -1;}swap(x, y);}
}void group(int a, int p) {while (a != p) {int b = match[a], c = _next[b];// _next数组是用来标记花朵中的路径的，综合match数组来用，实际上形成了  // 双向链表，如(x, y)是匹配的，_next[x]和_next[y]就可以指两个方向了。  if (findb(c) != p) _next[c] = b;// 奇环中的点都有机会向环外找到匹配，所以都要标记成S型点加到队列中去，  // 因环内的匹配数已饱和，因此这些点最多只允许匹配成功一个点，在aug中  // 每次匹配到一个点就break终止了当前阶段的搜索，并且下阶段的标记是重  // 新来过的，这样做就是为了保证这一点。  if (mark[b] == 2) mark[Q[rear++] = b] = 1;if (mark[c] == 2) mark[Q[rear++] = c] = 1;unit(a, b); unit(b, c);a = c;}
}// 增广
void aug(int s) {for (int i = 0; i < n; i++) // 每个阶段都要重新标记  _next[i] = -1, belong[i] = i, mark[i] = 0, vis[i] = -1;mark[s] = 1;Q[0] = s; rear = 1;for (int front = 0; match[s] == -1 && front < rear; front++) {int x = Q[front]; // 队列Q中的点都是S型的  for (int i = 0; i < (int)e[x].size(); i++) {int y = e[x][i];if (match[x] == y) continue; // x与y已匹配，忽略  if (findb(x) == findb(y)) continue; // x与y同在一朵花，忽略  if (mark[y] == 2) continue; // y是T型点，忽略  if (mark[y] == 1) { // y是S型点，奇环缩点  int r = LCA(x, y); // r为从i和j到s的路径上的第一个公共节点  if (findb(x) != r) _next[x] = y; // r和x不在同一个花朵，_next标记花朵内路径  if (findb(y) != r) _next[y] = x; // r和y不在同一个花朵，_next标记花朵内路径  // 将整个r -- x - y --- r的奇环缩成点，r作为这个环的标记节点，相当于论文中的超级节点  group(x, r); // 缩路径r --- x为点  group(y, r); // 缩路径r --- y为点  }else if (match[y] == -1) { // y自由，可以增广，R12规则处理  _next[y] = x;for (int u = y; u != -1; ) { // 交叉链取反  int v = _next[u];int mv = match[v];match[v] = u, match[u] = v;u = mv;}break; // 搜索成功，退出循环将进入下一阶段  }else { // 当前搜索的交叉链+y+match[y]形成新的交叉链，将match[y]加入队列作为待搜节点  _next[y] = x;mark[Q[rear++] = match[y]] = 1; // match[y]也是S型的  mark[y] = 2; // y标记成T型  }}}
}bool g[N][N];
int main() {scanf("%d", &n);for (int i = 0; i < n; i++)for (int j = 0; j < n; j++) g[i][j] = false;// 建图，双向边  int x, y; while (scanf("%d%d", &x, &y) != EOF) {x--, y--;if (x != y && !g[x][y])e[x].push_back(y), e[y].push_back(x);g[x][y] = g[y][x] = true;}// 增广匹配  for (int i = 0; i < n; i++) match[i] = -1;for (int i = 0; i < n; i++) if (match[i] == -1) aug(i);// 输出答案  int tot = 0;for (int i = 0; i < n; i++) if (match[i] != -1) tot++;printf("%d\n", tot);for (int i = 0; i < n; i++) if (match[i] > i)printf("%d %d\n", i + 1, match[i] + 1);return 0;
}

图论--一般图带花树匹配--模板相关推荐

图论--一般带花树匹配
带花树就是说一个非二分图,图中带有奇环的图,我们不能在奇环中找增广路,因为会陷入死循环,我们可以将带花树的花(奇环)部分缩成点处理,剩下的图就是一个无奇环的图.我们再找增广路,而奇环中的的点我们可以随 ...
二分图带权匹配、最佳匹配与KM算法
---------------------以上转自ByVoid神牛博客,并有所省略. [二分图带权匹配与最佳匹配] 什么是二分图的带权匹配?二分图的带权匹配就是求出一个匹配集合,使得集合中边的权值 ...
定位匹配模板匹配地图_什么是地图匹配？
定位匹配模板匹配地图 By Marie Douriez, James Murphy, Kerrick Staley 玛丽·杜里兹(Marie Douriez),詹姆斯·墨菲(James Murph ...
“玲珑杯”ACM比赛 Round #18 C -- 图论你先敲完模板【Dp】
C -- 图论你先敲完模板 Time Limit:5s Memory Limit:256MByte Submissions:660Solved:160 DESCRIPTION 今天HHHH在操场上跑步 ...
图论--二分图最佳完美匹配（KM模板）
#include <iostream> #include <cstring> #include <cstdio>using namespace std; const ...
55.函数模板指针匹配(模板自动匹配*多的)
1 #include <iostream> 2 using namespace std; 3 4 //模板自动匹配*多的 5 template <class T> 6 void ...
OpenCV精进之路（十）：直方图匹配——模板匹配
模板匹配的作用在图像识别领域作用可大了.那什么是模板匹配? 模板匹配,就是在一幅图像中寻找另一幅模板图像最匹配(也就是最相似)的部分的技术. 说的有点抽象,下面给个例子说明就很明白了. 在上面这幅全明 ...
UOJ#80 二分图最大权匹配 [模板题]
从前一个和谐的班级,有 nlnl 个是男生,有 nrnr 个是女生.编号分别为 1,-,nl1,-,nl 和 1,-,nr1,-,nr. 有若干个这样的条件:第 vv 个男生和第 uu 个女生愿意结为 ...
HDU 2255 二分图最佳匹配模板题
题目大意: 给定每一个人能支付的房子价值,每个人最多且必须拥有一套房子,问最后分配房子可得到的最大收益抄了个别人的KM模板,就这样了... 1 #include <cstdio> 2 # ...

图论--一般图带花树匹配--模板

图论--一般图带花树匹配--模板相关推荐

最新文章

热门文章