码迷,mamicode.com
首页 > 其他好文 > 详细

【题解】Uoj79一般图最大匹配

时间:2018-06-29 18:29:05      阅读:214      评论:0      收藏:0      [点我收藏+]

标签:its   bsp   getchar   节点   for   开始   void   浅谈   c++   

  带花树裸题,感觉带花树强强……不会的勿看此文,解释的可能不对,只是给自己看的!!!如题,带花树即为求一般图最大匹配算法(匈牙利与dinic为二分图最大匹配)。推荐论文:2015年《浅谈图的匹配算法及其应用》(长郡中学    ——陈胤伯)。论文当中有对于带花树算法的详细解析,在这里只想记录一下算法的基本流程:

——————————————————————————

  \(id\) : 记录一个点为奇点/偶点(0偶1奇)。

  \(fa\) : 并查集记录一个点属于哪一个点为根的花。

——————————————————————————

  我们对于每一个没有匹配到的点进行 bfs,默认该点为一个偶点。

  当我们遍历到一个未访问过的点(在此次bfs中):

    -若该点未匹配:

    Yes!我们找到了一条新的增广路。顺着 \( pre \) 数组的指引反向增广,修改匹配的对象。(当前点为\(x\), \(y = pre[x]\),\(z = match[y]\), 则应使 \(x\), \(y\) 成为匹配边,\(z\) 点继续增广(过程同上,一直到增广不下去了为止))。

    -若该点已匹配:

    将它的匹配点标记为偶点,推进队列。

  遍历到一个访问过的点(在此次bfs中):

    -若该点是一个奇点:

    找到的是一个偶环,无视。

    -若该点是一个偶点:

    如果这两个点本来就在同一朵花中,无视。否则进行缩花。我们找到这两点在bfs树上的 lca: \(k\);将这两点之间连上边从\(x\) 开始向上遍历所有的匹配边,找到一朵花的根节点就修改其花根为 \(k\); 若找到的点有匹配点,则将匹配点也标记为偶点并压入队列。

  总之是个玄学算法……

#include <bits/stdc++.h>
using namespace std;
#define maxn 300000
int n, m, match[maxn];
int pre[maxn], id[maxn];
int timer, fa[maxn], ans;
int vis[maxn];
queue <int> q;

struct edge
{
    int cnp = 1, head[maxn], to[maxn], last[maxn];
    void add(int u, int v)
    {
        to[cnp] = v, last[cnp] = head[u], head[u] = cnp ++;
        to[cnp] = u, last[cnp] = head[v], head[v] = cnp ++;
    }
}E1;

int read()
{
    int x = 0, k = 1;
    char c;
    c = getchar();
    while(c < 0 || c > 9) { if(c == -) k = -1; c = getchar(); }
    while(c >= 0 && c <= 9) x = x * 10 + c - 0, c = getchar();
    return x * k;
}

int find(int x) { return x == fa[x] ? x : fa[x] = find(fa[x]); }

void Clear(queue<int> &q) {
    queue <int> empty;
    swap(empty, q);
}

int LCA(int x, int y)
{
    timer ++;
    while(vis[x] != timer) 
    {
        if(x)
        {
            x = find(x);
            if(vis[x] == timer) return x;
            vis[x] = timer;
            if(match[x]) x = find(pre[match[x]]);
            else x = 0;
        }
        swap(x, y);
    }
    return x;
}

void Change(int x, int y, int k)
{
    while(find(x) != k)
    {
        pre[x] = y; int z = match[x];
        if(id[z] == 1) id[z] = 0, q.push(z);
        if(find(z) == z) fa[z] = k;
        if(find(x) == x) fa[x] = k;
        y = z; x = pre[y];
    }
}

bool bfs(int u)
{
    for(int i = 1; i <= n; i ++) id[i] = -1, fa[i] = i;
    Clear(q); q.push(u); id[u] = 0;
    while(!q.empty())
    {
        int u = q.front(); q.pop();
        for(int i = E1.head[u]; i; i = E1.last[i])
        {
            int v = E1.to[i];
            if(id[v] == -1)
            {
                pre[v] = u; id[v] = 1;
                if(match[v])
                {
                    id[match[v]] = 0; q.push(match[v]);
                    continue;
                }
                int last, t, now = v;
                while(now)
                {
                    t = pre[now]; last = match[t];
                    match[t] = now, match[now] = t;
                    now = last;
                }
                return 1;
            }
            else if(!id[v] && find(u) != find(v))
            {
                int lca = LCA(u, v);
                Change(u, v, lca), Change(v, u, lca);
            }
        }
    }
    return 0;
}

int main()
{
    n = read(), m = read();
    for(int i = 1; i <= m; i ++)
    {
        int x = read(), y = read();
        E1.add(x, y);
    }
    for(int i = 1; i <= n; i ++) if(!match[i] && bfs(i)) ans ++;
    printf("%d\n", ans);
    for(int i = 1; i <= n; i ++) printf("%d ", match[i]);
    return 0;
}

 

【题解】Uoj79一般图最大匹配

标签:its   bsp   getchar   节点   for   开始   void   浅谈   c++   

原文地址:https://www.cnblogs.com/twilight-sx/p/9244779.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!