UVa-12549——二分图最大匹配

题目链接：https://vjudge.net/contest/213916#problem/I

紫书图论习题11-9，和白书的二分图经典的第一个例题，UVa-11419，是一样的建模方式，可以看看白书上的解释。不过这题多了一个障碍物，也就是说不能像上面一样选择一行就将这一行（或者一列）的目标点清除。所以，我们的目的就是将图形建成没有障碍物的新图，就可以套用上面的模板了。刘汝佳的解法是，每碰见一个障碍物，就将后面的所有东西移动到下一行，这样就不会有障碍物了。但这种方法实现起来很麻烦，我没写出来，但是我看了一篇博客之后，发现了另一种建图思路：https://www.cnblogs.com/zyb993963526/p/6628791.html

现在的思路就是：反正最后我们建图只将最后有目标点的X,Y行建边，那么重新开始统计行数和列数（因为现在的行和列中如果有障碍物，那就要新建一个行或者列），对于行：遍历所有的行中的内容，如果发现有目标点，row_num++，如果之后的行中碰到了障碍物，则再换一个新行（前提是有目标点）。如对于下面这一行：

                                                                                ...*#*...

首先发现一个目标点，row_num++,flag=false；然后发现一个障碍物，如果后面有目标点，就得换新行，flag=true；最后又发现一个目标点，由于flag=true，所以row_num++。

列也是同理。最后再用一个pair存储行，列以及对应的新的行号或列号，加边的时候就利用这些数据就行了。

再详细的我也不想讲了，看代码吧，HK算法写的，快一点：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <queue>
#include <vector>
#include <map>
#include <iomanip>
#define Pair pair<int,int>
#define MP(x,y) make_pair(x,y)
using namespace std;
const int maxn=1000+50;
const int INF=0x3f3f3f3f;
Pair new_mp[maxn][maxn];
vector<int>G[maxn*2];
int uN;//左端点个数
int Mx[maxn],My[maxn];//判断是否已匹配，并记录匹配点
int dx[maxn],dy[maxn];//记录匹配点到增广路根节点的距离
int dis;
bool used[maxn];
void init()
{
    for(int i=0;i<maxn;i++)G[i].clear();
}
bool searchP()
{
    queue<int>Q;
    dis=INF;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=0;i<uN;i++)
    {
        if(Mx[i]==-1)
        {
            Q.push(i);
            dx[i]=0;
        }
    }
    while(!Q.empty())
    {
        int u=Q.front();Q.pop();
        if(dx[u]>dis)break;
        int sz=G[u].size();
        for(int i=0;i<sz;i++)
        {
            int v=G[u][i];
            if(dy[v]==-1)
            {
                dy[v]=dx[u]+1;
                if(My[v]==-1)
                    dis=dy[v];
                else
                {
                    dx[My[v]]=dy[v]+1;//v匹配的点的距离=v的距离+1
                    Q.push(My[v]);
                }

            }
        }
    }
    return dis!=INF;
}
bool DFS(int u)
{
    int sz=G[u].size();
    for(int i=0;i<sz;i++)
    {
        int v=G[u][i];
        if(!used[v]&&dy[v]==dx[u]+1)//沿着增广路dfs
        {
            used[v]=true;
            if(My[v]!=-1&&dy[v]==dis)continue;
            if(My[v]==-1||DFS(My[v]))//找到了对应的增广路
            {
                My[v]=u;
                Mx[u]=v;
                return true;
            }
        }
    }
    return false;
}
int Maxmatch()
{
    int res=0;
    memset(Mx,-1,sizeof(Mx));
    memset(My,-1,sizeof(My));
    while(searchP())//找到增广路
    {
        memset(used,false,sizeof(used));
        for(int i=0;i<uN;i++)
        {
            if(Mx[i]==-1&&DFS(i))
                res++;
        }
    }
    return res;
}

int mp[110][110];

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        for(int i=0;i<maxn;i++)G[i].clear();
        int n,m;
        scanf("%d%d",&n,&m);
        int x,y;
        memset(mp,0,sizeof(mp));
        scanf("%d",&x);
        for(int i=1;i<=x;i++)
        {
            int r,c;
            scanf("%d%d",&r,&c);
            mp[r][c]=1;
        }
        scanf("%d",&y);
        for(int i=1;i<=y;i++)
        {
            int r,c;
            scanf("%d%d",&r,&c);
            mp[r][c]=-1;
        }
        int row_num=0;//行数
        for(int i=1;i<=n;i++)
        {
            bool flag=1;//表示是否需要换新行
            for(int j=1;j<=m;j++)
            {
                if(mp[i][j]==1)
                {
                    if(flag==1)row_num++;
                    new_mp[i][j].first=row_num;
                    flag=0;
                }
                if(mp[i][j]==-1)
                    flag=1;//表示还需要再新加一行
            }
        }
        int cow_num=0;//列数
        for(int j=1;j<=m;j++)
        {
            bool flag=true;
            for(int i=1;i<=n;i++)
            {
                if(mp[i][j]==1)
                {
                    if(flag==1)cow_num++;
                    new_mp[i][j].second=cow_num;
                    flag=false;
                }
                if(mp[i][j]==-1)
                    flag=true;
            }
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {
                if(mp[i][j]==1)
                {
                    G[new_mp[i][j].first-1].push_back(new_mp[i][j].second+row_num-1);
                    G[new_mp[i][j].second+row_num-1].push_back(new_mp[i][j].first-1);
                }
            }
        }
        uN=row_num;
        printf("%d\n",Maxmatch());
    }
    return 0;
}