BZOJ4554(Tjoi2016&Heoi2016)[游戏]--二分图最大匹配

【链接】
bzoj4554

【题目大意】
给你一个n*m的网格地图，在这张地图上最多能放上多少个炸弹能使得任意两个炸弹之间不会互相炸到。炸弹能炸到的范围是该炸弹所在的一行和一列，炸弹的威力可以穿透软石头，但是不能穿透硬石头。网格地图中’*’表示空地（可以放炸弹），’#’表示硬石头（不能放炸弹），’x’表示软石头（不能放炸弹）。

【解题报告】
初看此题大概就知道可以刷二分图最大匹配，但是却有有些无从下手。但是通过仔细思考我们可以发现硬石头将空地隔开后就这一行在硬石头两边的炸弹互相就影响不到了，列也是。所以就可以想到给行和列表上序号，对于某一行或某一列如果炸弹放在这个空地上能影响到的地方序号相同。
如图：

这样就可以轻松刷二分图最大匹配了。

Ps:因为要做到上图中的标号写起来太复杂，所以代码中真正标号并不完全相同，但是绝对正确。

#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=55,maxm=2505;
int n,m,tot,Hid,Lid,ans,who[maxm],H[maxn][maxn],L[maxn][maxn],lnk[maxm],son[maxm],nxt[maxm];
char a[maxn][maxn];
bool vis[maxm],pd[maxm];
inline int Read()
{
    int res=0;
    char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') res=res*10+ch-48,ch=getchar();
    return res;
}
inline char Read_ch()
{
    char ch=getchar();
    while (ch!='*'&&ch!='#'&&ch!='x') ch=getchar();
    return ch;
}
void Add(int x,int y)
{
    son[++tot]=y; nxt[tot]=lnk[x]; lnk[x]=tot;
}   
bool Find(int x)
{
    if (vis[x]) return 0;
    vis[x]=1;
    for (int j=lnk[x]; j; j=nxt[j])
     if (!who[son[j]]||Find(who[son[j]])) {who[son[j]]=x; return 1;}
    return 0;
}
int main()
{
    freopen("4554.in","r",stdin);
    freopen("4554.out","w",stdout);
    n=Read(); m=Read(); Hid=Lid=tot=ans=0;
    for (int i=1; i<=n; i++)
     for (int j=1; j<=m; j++) a[i][j]=Read_ch();
    memset(H,0,sizeof(H));
    memset(L,0,sizeof(L));
    for (int i=1; i<=n; i++)
    {
        Hid++;
        for (int j=1; j<=m; j++)
        {
            if (a[i][j]=='*') H[i][j]=Hid;
            if (a[i][j]=='#') Hid++;
        }
    }
    for (int j=1; j<=m; j++)
    {
        Lid++;
        for (int i=1; i<=n; i++)
        {
            if (a[i][j]=='*') L[i][j]=Lid;
            if (a[i][j]=='#') Lid++;
        }
    }
    for (int i=1; i<=n; i++)
     for (int j=1; j<=m; j++)
      if (a[i][j]=='*') Add(H[i][j],L[i][j]);
    memset(pd,0,sizeof(pd));
    for (int i=1; i<=n; i++)
     for (int j=1; j<=m; j++)
      if (a[i][j]=='*')
       {
        memset(vis,0,sizeof(vis));
        if (!pd[H[i][j]]) ans+=Find(H[i][j]),pd[H[i][j]]=1;
       }
    printf("%d",ans);
    return 0;
}