Hdu 5352 MZL's City (多重匹配)

题目链接：

　　Hdu 5352 MZL's City

题目描述：

　　有n各节点，m个操作。刚开始的时候节点都是相互独立的，一共有三种操作：

　　1：把所有和x在一个连通块内的未重建过的点全部重建。

　　2：建立一条双向路(x,y)

　　3：又发生了地震，p条路被毁。

　　问最后最多有多少个节点被重建，输出重建节点的最小字典序。

解题思路：

　　这几天正好在学匹配，但是昨天下午还是没有看出来这个是匹配题目。看了题解扪心自问了自己三次，是不是傻。就是把每个需要重建的节点x拆成k个点，然后对每个拆分后的点和与拆点在同一连通块里面的点建边，然后按照倒序进行匹配，保证字典序最大。

　　但是匹配好像要比网络流慢好多，改天还是去学一下新姿势的好。

#include <vector>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int maxn = 205;
const int N = 205*500;
vector <int> G[N];
int ans[maxn*2], p[maxn], used[maxn], vis[maxn];
int n, m, k, num, nu, maps[maxn][maxn];
void dfs (int u)
{//求连通块内的点
    vis[u] = 1;
    p[num ++] = u;
    for (int i=1; i<=n; i++)
        if (!vis[i] && maps[u][i])
            dfs (i);
}
bool Find (int u)
{
    for (int i=0; i<G[u].size(); i++)
    {
        int v = G[u][i];
        if (!vis[v])
        {
            vis[v] = 1;
            if (!used[v] || Find(used[v]))
            {
                used[v] = u;
                return true;
            }
        }
    }
    return false;
}
int hungry ()
{
    int res = 0;
    memset (used, 0, sizeof(used));
    for (int i=nu-1; i>=0; i--)//倒序求最大匹配
        for (int j=i*k; j<(i+1)*k; j++)
        {
            memset (vis, 0, sizeof(vis));
            if (Find(j))
            {
                res ++;
                ans[i] ++;
            }
        }
    return res;
}
int main ()
{
    int t;
    scanf ("%d", &t);
    while (t --)
    {
        scanf ("%d %d %d", &n, &m, &k);
        memset (maps, 0, sizeof(maps));
        for (int i=0; i<N; i++)
            G[i].clear();
        nu = 0;
        while (m --)
        {
            int type, x, y, z;
            scanf ("%d", &type);
            if (type == 1)
            {
                scanf ("%d", &x);
                num = 0;
                memset (vis, 0, sizeof(vis));
                dfs (x);
                for (int i=0; i<num; i++)//建边
                    for (int j=nu*k; j<(nu+1)*k; j++)//拆点
                        G[j].push_back(p[i]);
                nu ++;
            }
            else if (type == 2)
            {
                scanf ("%d %d", &x, &y);
                maps[x][y] = maps[y][x] = 1;
            }
            else
            {
                scanf ("%d", &z);
                while (z --)
                {
                    scanf ("%d %d", &x, &y);
                    maps[x][y] = maps[y][x] = 0;
                }
            }
        }
        memset (ans, 0, sizeof(ans));
        int res = hungry ();
        printf ("%d
", res);
        for (int i=0; i<nu; i++)
            printf ("%d%c", ans[i], i==nu-1?'
':' ');
    }
    return 0;
}