题目链接: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;
}