UOJ 79 带花树入门

从前一个和谐的班级，所有人都是搞OI的。有 nn 个是男生，有 00 个是女生。男生编号分别为 1,…,n1,…,n。

现在老师想把他们分成若干个两人小组写动态仙人掌，一个人负责搬砖另一个人负责吐槽。每个人至多属于一个小组。

有若干个这样的条件：第 vv 个男生和第 uu 个男生愿意组成小组。

请问这个班级里最多产生多少个小组？
输入格式

第一行两个正整数，n,mn,m。保证 n≥2n≥2。

接下来 mm 行，每行两个整数 v,uv,u 表示第 vv 个男生和第 uu 个男生愿意组成小组。保证 1≤v,u≤n1≤v,u≤n，保证 v≠uv≠u，保证同一个条件不会出现两次。
输出格式

第一行一个整数，表示最多产生多少个小组。

接下来一行 nn 个整数，描述一组最优方案。第 vv 个整数表示 vv 号男生所在小组的另一个男生的编号。如果 vv 号男生没有小组请输出 00。
样例一
input

10 20
9 2
7 6
10 8
3 9
1 10
7 1
10 9
8 6
8 2
8 1
3 1
7 5
4 7
5 9
7 8
10 4
9 1
4 8
6 3
2 5

output

5
9 5 6 10 2 3 8 7 1 4

样例二
input

5 4
1 5
4 2
2 1
4 3

output

2
2 1 4 3 0

限制与约定

1≤n≤5001≤n≤500，1≤m≤1247501≤m≤124750。

时间限制：1s1s

空间限制：256MB256MB

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>

using namespace std;

#define LL long long
#define pb push_back
#define Set(a, v) memset(a, v, sizeof(a))
#define For(i, a, b) for(int i = (a); i <= (int)(b); i++)
#define Forr(i, a, b) for(int i = (a); i >= (int)(b); i--)

#define MAXN (500+5)

struct BlossomTree{
    vector<int> G[MAXN];
    queue<int> q;
    int n;
    int pre[MAXN], link[MAXN], col[MAXN], fa[MAXN];
    bool vis[MAXN];

    void init(int rn){
        n = rn;
        Set(link, 0);
        For(i, 0, n) G[i].clear();
    }

    void Addedges(int u, int v){
        G[u].pb(v); G[v].pb(u);
    }

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

    void Combine(int u, int lca){
        while(u != lca){
            int f = link[u], g = pre[f];
            if(find(g) != lca) pre[g] = f; //put the rever edge
            if(col[f] == 1) col[f] = 2, q.push(f);
            fa[find(u)] = find(f); fa[find(f)] = find(g);

            u = pre[link[u]];
        }
    }

    void Contract(int u, int v){
        Set(vis, 0);

        int lca = 0;
        for(int i = u; i; i = pre[link[i]]) i = find(i), vis[i] = true;
        for(int i = v; i; i = pre[link[i]]){
            i = find(i);
            if(vis[i]){lca = i; break;}
        }

        if(find(u) != lca) pre[u] = v;
        if(find(v) != lca) pre[v] = u;
        Combine(u, lca); Combine(v, lca);
    }

    void Bfs(int s){
        while(!q.empty()) q.pop();
        Set(col, 0); Set(pre, 0);
        For(i, 1, n) fa[i] = i;

        q.push(s);
        while(!q.empty()){
            int now = q.front(); q.pop();

            For(i, 0, G[now].size()-1){
                int v = G[now][i];
                if(find(now)==find(v) || link[now]==v || col[v]==1) continue;

                if(col[v] == 2) Contract(now, v);  
                else if(link[v]){
                    pre[v] = now; col[v] = 1;
                    col[link[v]] = 2;
                    q.push(link[v]);
                }else{
                    pre[v] = now;

                    int g = link[now], f = now, o = v;
                    while(o){
                        link[f] = o; link[o] = f;
                        o = g; f = pre[o]; g = link[f];
                    }
                    return;
                }
            }
        }
    }

    int MaxMatch(){
        For(i, 1, n) if(!link[i]) Bfs(i);

        int ret = 0;
        For(i, 1, n) if(link[i]) ret++;
        return ret/2;
    }

    void print(){
        printf("%d", link[1]);
        For(i, 2, n) printf(" %d", link[i]);
        printf("\n");
    }
}BT;

int main(){
    int n, m;
    scanf("%d%d", &n, &m);

    BT.init(n);
    For(i, 1, m){
        int u, v;
        scanf("%d%d", &u, &v);
        BT.Addedges(u, v);
    }

    printf("%d\n", BT.MaxMatch());
    BT.print();

    return 0;
}