[HEOI2016/TJOI2016] 游戏

题目描述

在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

    一类网络流套路。
    首先把每行每列不含'#’的那些极长段处理出来，显然每一段只能最多放一个，最后再建个二分图匹配一下就好啦。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=3005;
#define pb push_back
vector<int> g[maxn];
struct lines{
	int to,flow,cap;
}l[maxn*maxn];
int S,t=-1,T,d[maxn],cur[maxn];
bool v[maxn];

inline void add(int from,int to,int cap){
	l[++t]=(lines){to,0,cap},g[from].pb(t);
	l[++t]=(lines){from,0,0},g[to].pb(t);
}

inline bool BFS(){
	memset(v,0,sizeof(v));
	queue<int> q; q.push(S);
	v[S]=1,d[S]=0;
	int x; lines e;
	
	while(!q.empty()){
		x=q.front(),q.pop();
		for(int i=g[x].size()-1;i>=0;i--){
			e=l[g[x][i]];
			if(e.flow<e.cap&&!v[e.to]) v[e.to]=1,d[e.to]=d[x]+1,q.push(e.to);
		}
	}
	
	return v[T];
}

int dfs(int x,int A){
	if(x==T||!A) return A;
	int flow=0,f,sz=g[x].size();
	for(int &i=cur[x];i<sz;i++){
		lines &e=l[g[x][i]];
		if(d[x]+1==d[e.to]&&(f=dfs(e.to,min(A,e.cap-e.flow)))){
			A-=f,flow+=f;
			e.flow+=f,l[g[x][i]^1].flow-=f;
			if(!A) break;
		}
	}
	
	return flow;
}

inline int max_flow(){
	int an=0;
	while(BFS()){
		memset(cur,0,sizeof(cur));
		an+=dfs(S,1<<30);
	}
	return an;
}

int n,m,L,cnt,R[55][55],C[55][55];
char s[233][233];

inline void build(){
	for(int i=1;i<=n;i++){
		R[i][1]=++cnt;
		for(int j=2;j<=m;j++){
			if(s[i][j]=='#') R[i][j]=++cnt;
			else R[i][j]=cnt;
		}
	}
	
	L=cnt;
	
	for(int i=1;i<=m;i++){
		C[1][i]=++cnt;
		for(int j=2;j<=n;j++){
			if(s[j][i]=='#') C[j][i]=++cnt;
			else C[j][i]=cnt;
		}
	}	
	
	T=cnt+1;
	
	for(int i=1;i<=L;i++) add(S,i,1);
	for(int i=L+1;i<T;i++) add(i,T,1);
	
	for(int i=1;i<=n;i++)
	    for(int j=1;j<=m;j++) if(s[i][j]=='*') add(R[i][j],C[i][j],1);
}

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%s",s[i]+1);
	
	build();	
	
	printf("%d
",max_flow());
	
	return 0;
}