洛谷 P2825 [HEOI2016] 游戏

题目描述

在2016年，佳缘姐姐喜欢上了一款游戏，叫做泡泡堂。简单的说，这个游戏就是在一张地图上放上若干个炸弹，看是否能炸到对手，或者躲开对手的炸弹。在玩游戏的过程中，小H想到了这样一个问题：当给定一张地图，在这张地图上最多能放上多少个炸弹能使得任意两个炸弹之间不会互相炸到。炸弹能炸到的范围是该炸弹所在的一行和一列，炸弹的威力可以穿透软石头，但是不能穿透硬石头。给定一张n*m的网格地图:其中代表空地，炸弹的威力可以穿透，可以在空地上放置一枚炸弹。x代表软石头，炸弹的威力可以穿透，不能在此放置炸弹。#代表硬石头，炸弹的威力是不能穿透的，不能在此放置炸弹。例如：给出1*4的网格地图*xx，这个地图上最多只能放置一个炸弹。给出另一个1*4的网格地图x#，这个地图最多能放置两个炸弹。现在小H任意给出一张n*m的网格地图，问你最多能放置多少炸弹。

输入输出格式

输入格式：
第一行输入两个正整数n,m,n表示地图的行数，m表示地图的列数。1≤n,m≤50。接下来输入n行m列个字符，代表网格地图。*的个数不超过n*m个。

输出格式：
输出一个整数a，表示最多能放置炸弹的个数

输入输出样例

输入样例#1：
4 4

∗∗∗

∗#∗∗
∗∗#∗
xxx#
输出样例#1：
5

---

【分析】
是个神奇的模型…参见poj 2226 Muddy Fields

---

【代码】

//[HEOI2016]游戏 
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#define inf 1e9+7
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
queue <int> q;
const int mxn=100005;
int n,m,s,t,nx,ny,cnt;
bool map[2555];
char c[55][55];
int head[mxn],x[55][55],y[55][55],dis[mxn];
struct edge {int to,flow,next;} f[mxn<<1];
inline void add(int u,int v,int flow)
{
    f[++cnt].to=v,f[cnt].next=head[u],f[cnt].flow=flow,head[u]=cnt;
    f[++cnt].to=u,f[cnt].next=head[v],f[cnt].flow=0,head[v]=cnt;
}
inline bool bfs()
{
    memset(dis,-1,sizeof dis);
    q.push(s);
    dis[s]=0;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u];i;i=f[i].next)
        {
            int v=f[i].to,flow=f[i].flow;
            if(dis[v]==-1 && flow>0)
              dis[v]=dis[u]+1,q.push(v);
        }
    }
    return dis[t]>0;
}
inline int find(int u,int low)
{
    int i,j,k,sum=0,a=0;
    if(u==t) return low;
    for(int i=head[u];i;i=f[i].next)
    {
        int v=f[i].to,flow=f[i].flow;
        if(dis[v]==dis[u]+1 && flow>0 && (a=find(v,min(flow,low-sum))))
        {
            sum+=a;
            f[i].flow-=a;
            if(i&1) f[i+1].flow+=a;
            else f[i-1].flow+=a;
        }
    }
    if(!sum) dis[u]=-1;
    return sum;
}
int main()
{
    int i,j,k,ans=0;
    scanf("%d%d",&n,&m);
    fo(i,1,n) scanf("%s",c[i]+1);
    fo(i,1,n) fo(j,1,m) x[i][j]=(j==1||c[i][j]=='#')?++nx:nx; ny=nx+1;
    fo(j,1,m) fo(i,1,n) y[i][j]=(i==1||c[i][j]=='#')?++ny:ny;
    fo(i,1,n) fo(j,1,m) if(c[i][j]=='*') add(x[i][j],y[i][j],1);
    s=0,t=ny+1;
//  fo(i,1,n)
//  {
//      fo(j,1,m)
//        cout<<x[i][j];
//      cout<<endl;
//  }
    memset(map,0,sizeof map);
    fo(i,1,n) fo(j,1,m) if(c[i][j]=='*' && !map[x[i][j]]) add(s,x[i][j],1),map[x[i][j]]=1;
    memset(map,0,sizeof map);
    fo(i,1,n) fo(j,1,m) if(c[i][j]=='*' && !map[y[i][j]]) add(y[i][j],t,1),map[y[i][j]]=1;
    while(bfs())
      ans+=find(s,inf);
    printf("%d\n",ans);
    return 0;
}
/*
4 4
#***
*#**
**#*
xxx#
*/