算法原理详见 http://www.csie.ntnu.edu.tw/~u91029/Matching.html

orz 带花树很神奇,挖坑最大权匹配

#include <iostream>
#include <cstring>
#include <cstdio>
#include <deque>
using namespace std;
const int maxn = ;
deque<int> p[maxn]; //树根到x的交错路径
bool adj[maxn][maxn];
int m[maxn]; //记录匹配情况
int d[maxn]; //记录奇点和偶点
int q[maxn], *qf, *qb; //记录可延伸的偶点
int n; void label_one_side(int x, int y, int bi){
for(int i = bi+; i < p[x].size(); i++){
int z = p[x][i];
if(d[z] == ){
p[z] = p[y];
p[z].insert(p[z].end(), p[x].rbegin(), p[x].rend()-i);
d[z] = ;
*qb++ = z;
}
}
} bool BFS(int r){
for(int i = ; i < n; i++) p[i].clear();
p[r].push_back(r);
for(int i = ; i < n; i++) d[i] = -;
d[r] = ;
qf = qb = q;
*qb++ = r;
while(qf < qb){
for(int x = *qf++, y = ; y < n; y++){
if(adj[x][y] && m[y] != y){
if(d[y] == -){
if(m[y] == -){
for(int i = ; i+ < p[x].size(); i += ){
m[p[x][i]] = p[x][i+];
m[p[x][i+]] = p[x][i];
}
m[x] = y; m[y] = x;
return true;
} else {
int z = m[y];
p[z] = p[x];
p[z].push_back(y);
p[z].push_back(z);
d[y] = ; d[z] = ;
*qb++ = z;
}
} else if(d[y] == ) {
int bi = ;
while(bi < p[x].size() && bi < p[y].size() && p[x][bi] == p[y][bi]) bi++;
bi--;
label_one_side(x, y, bi);
label_one_side(y, x, bi);
}
}
}
}
return false;
} int match(){
for(int i = ; i < n; i++) m[i] = -;
int c = ;
for(int i = ; i < n; i++){
if(m[i] == -)
if(BFS(i)) c++;
else m[i] = i;
}
return c;
} int main(){
cin>>n;
int x, y;
while(cin>>x>>y) adj[x][y] = adj[y][x] = true;
match();
}

这个是缩点版本的

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N = ;
int n;
bool adj[N][N];
int p[N];
int m[N];
int d[N];
int c1[N], c2[N];
int q[N], *qf, *qb; int pp[N];
int f(int x) { return x == pp[x] ? x : pp[x] = f(pp[x]); }
void u(int x, int y) { pp[x] = y; }
int v[N]; void path(int r, int x){
if(r == x) return;
if(d[x] == ){
path(r, p[p[x]]);
int i = p[x], j = p[p[x]];
m[i] = j; m[j] = i;
} else if(d[x] == ){
path(m[x], c1[x]);
path(r, c2[x]);
int i = c1[x], j = c2[x];
m[i] = j; m[j] = i;
}
} int lca(int x, int y, int r){
int i = f(x), j = f(y);
while(i != j && v[i] != && v[j] != ){
v[i] = ; v[j] = ;
if(i != r) i = f(p[i]);
if(j != r) j = f(p[j]);
}
int b = i, z = j; if(v[j] == ) swap(b, z);
for(int i = b; i != z; i = f(p[i])) v[i] = -;
v[z] = -;
return b;
} void contract_one_side(int x, int y, int b){
for(int i = f(x); i != b; i = f(p[i])){
u(i, b);
if(d[i] == ) c1[i] = x, c2[i] = y, *qb++ = i;
}
} bool BFS(int r){
for(int i = ; i < n; i++) pp[i] = i;
for(int i = ; i < n; i++) v[i] = d[i] = -;
d[r] = ;
qf = qb = q;
*qb++ = r;
while(qf < qb){
for(int x = *qf++, y = ; y < n; y++){
if(!adj[x][y] || m[y] == y || f(x) == f(y)) continue;
if(d[y] == -){
if(m[y] == -){
path(r, x);
m[x] = y; m[y] = x;
return true;
} else {
p[y] = x; p[m[y]] = y;
d[y] = ; d[m[y]] = ;
*qb++ = m[y];
}
} else if(d[f(y)] == ){
int b = lca(x, y, r);
contract_one_side(x, y, b);
contract_one_side(y, x, b);
}
}
}
return false;
} int main(){ }

Augmenting Path Algorithm : 一般图最大匹配的更多相关文章

  1. [转]带花树,Edmonds's matching algorithm,一般图最大匹配

    看了两篇博客,觉得写得不错,便收藏之.. 首先是第一篇,转自某Final牛 带花树……其实这个算法很容易理解,但是实现起来非常奇葩(至少对我而言). 除了wiki和amber的程序我找到的资料看着都不 ...

  2. Maximum Cardinality Bipartite Matching: Augmenting Path Algorithm

    http://www.csie.ntnu.edu.tw/~u91029/Matching.html int nx,ny; int mx[N],my[N]; bool vy[N]; bool g[N][ ...

  3. [转载]Maximum Flow: Augmenting Path Algorithms Comparison

    https://www.topcoder.com/community/data-science/data-science-tutorials/maximum-flow-augmenting-path- ...

  4. UOJ79 一般图最大匹配

    题目描述 从前一个和谐的班级,所有人都是搞OI的.有 nn 个是男生,有 00 个是女生.男生编号分别为 1,-,n1,-,n. 现在老师想把他们分成若干个两人小组写动态仙人掌,一个人负责搬砖另一个人 ...

  5. HDOJ 4687 Boke and Tsukkomi 一般图最大匹配带花树+暴力

    一般图最大匹配带花树+暴力: 先算最大匹配 C1 在枚举每一条边,去掉和这条边两个端点有关的边.....再跑Edmonds得到匹配C2 假设C2+2==C1则这条边再某个最大匹配中 Boke and ...

  6. HDU 4687 Boke and Tsukkomi (一般图最大匹配)【带花树】

    <题目链接> 题目大意: 给你n个点和m条边,每条边代表两点具有匹配关系,问你有多少对匹配是冗余的. 解题分析: 所谓不冗余,自然就是这对匹配关系处于最大匹配中,即该匹配关系有意义.那怎样 ...

  7. 【UOJ#79】一般图最大匹配(带花树)

    [UOJ#79]一般图最大匹配(带花树) 题面 UOJ 题解 带花树模板题 关于带花树的详细内容 #include<iostream> #include<cstdio> #in ...

  8. ZOJ 3316 Game 一般图最大匹配带花树

    一般图最大匹配带花树: 建图后,计算最大匹配数. 假设有一个联通块不是完美匹配,先手就能够走那个没被匹配到的点.后手不论怎么走,都必定走到一个被匹配的点上.先手就能够顺着这个交错路走下去,最后一定是后 ...

  9. 【Learning】带花树——一般图最大匹配

    一般图最大匹配--带花树 问题 ​ 给定一个图,求该图的最大匹配.即找到最多的边,使得每个点至多属于一条边. ​ 这个问题的退化版本就是二分图最大匹配. ​ 由于二分图中不存在奇环,偶环对最大匹配并无 ...

随机推荐

  1. Vue报错 [Vue warn]: Cannot find element

    在前端开发全面进入前端的时代 作为一个合格的前端开发工作者 框架是不可或缺的Vue React Anguar 作为前端小白,追随大佬的脚步来到来到博客园,更新现在正在学习的Vue 注 : 相信学习Vu ...

  2. tp3.2报错;syntax error, unexpected 'function' (T_FUNCTION), expecting identifier (T_STRING) or \\ (T_NS_SEPARATOR)

    出错原因:这个是php版本问题,laravel5.1的php版本要求是PHP >= 5.5.9,切换一下PHP版本就行.

  3. 微信小程序 提示框延时跳转

    wx.showToast({ title: '成功', icon: 'success', duration: 2000, success:function(){ console.log('haha') ...

  4. 软件的按契约设计(DbC---Design by Contract)

    一.DbC基本概念 DbC的思想源于商业活动中商家和用户的行为(义务和利益关系),双方都要遵守一个契约(合同),交易才能完成. 商家与用户的契约关系如下: 1. 商家必须提供某种产品(义务),并有权获 ...

  5. Kubernetes-apiserver

    Kubernetes API服务器为API对象验证和配置数据,这些对象包含Pod.Service.ReplicationController等等.API Server提供REST操作以及前端到集群的共 ...

  6. 39-Role以及Claims授权

    asp.net core多鼓励使用claims授权 1-使用role授权 在类或方法上贴上Roles,这样就知道有user的角色才可以访问 [Authorize(Roles="user&qu ...

  7. 一种简单实用的双向电平转换电路3.3V-5V

    当你使用3.3V的单片机的时候,电平转换就在所难免了,经常会遇到3.3转5V或者5V转3.3V的情况,这里介绍一个简单的电路,他可以实现两个电平的相互转换(注意是相互哦,双向的,不是单向的!).电路十 ...

  8. P1509 找啊找啊找GF

    P1509 找啊找啊找GF 题目背景 "找啊找啊找GF,找到一个好GF,吃顿饭啊拉拉手,你是我的好GF.再见." "诶,别再见啊..." 七夕...七夕...七 ...

  9. Android StateListDrawable的坑

    有问题的代码: StateListDrawable background = new StateListDrawable(); CircleDrawable pressedDrawable = new ...

  10. 获得通讯录并拨打电话 Android

    由于通讯录在手机里是以数据库贮存的 所以我们可以通过getContentResolver来获得通讯录 ,这个方法返回一个游标的数据类型,通过moveToNext()方法来获取所有的手机号码信息, 当然 ...