【BZOJ1066】最大流

题目通俗易懂
分析：
一开始没看出来是一个最大流。
不过…
每个柱子跳一下就会变矮…这不就是流量吗…
遂最大流。
每个柱子拆成柱底和柱顶，容量就是高度。
然后距离为d的柱子连起来，
起点连青蛙，容量为1，因为一只青蛙跳完了就没了。
可跳出去的点连汇点。
最大流就行了。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1010;
const int INF = 1e9+7;
int n,m;
struct edge
{
	int to;
	int cap;
	int rev;
};
vector<edge> G[maxn];
int depth[maxn],arc[maxn];
int M[maxn][maxn];
void add_edge(int from,int to,int cost)
{
	G[from].push_back((edge){to,cost,G[to].size()});
	G[to].push_back((edge){from,0,G[from].size()-1});
}
void bfs(int s)
{
	queue<int> q;
	memset(depth,-1,sizeof(depth));
	depth[s] = 0;
	q.push(s);
	while(!q.empty())
	{
		int v = q.front();
		q.pop();
		for(int i=0;i<G[v].size();i++)
		{
			edge &e = G[v][i];
			if(e.cap>0 && depth[e.to]<0)
			{
				depth[e.to] = depth[v]+1;
				q.push(e.to);
			}
		}
	}
}
int dfs(int v,int t,int f)
{
	if(v==t) return f;
	for(int &i=arc[v];i<G[v].size();i++)
	{
		edge &e = G[v][i];
		if(e.cap>0 && depth[e.to]==depth[v]+1)
		{
			int d = dfs(e.to,t,min(f,e.cap));
			if(d>0)
			{
				e.cap -= d;
				G[e.to][e.rev].cap += d;
				return d;
			}
		}
	}
	return 0;
}
int dinic(int s,int t)
{
	int flow = 0;
	while(true)
	{
		bfs(s);
		if(depth[t]<0) return flow;
		memset(arc,0,sizeof(arc));
		int f;
		while((f=dfs(s,t,INF))>0)
		{
			flow += f;
		}
	}
}
int dis(int x1,int y1,int x2,int y2)
{
	return abs(x1-x2)+abs(y1-y2);
}
int hah(int x,int y)
{
	return (x-1)*m+y;
}
int main()
{
	int cnt = 0;
	int d;
	scanf("%d%d%d",&n,&m,&d);
	int s = 0,t = n*m+n*m+1;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			char x;
			cin>>x;
			int t = x-'0';
			M[i][j] = t;
			if(t)
			{
				add_edge(hah(i,j),hah(i,j)+n*m,t);
			}
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			char x;
			cin>>x;
			if(x=='L')
			{
				cnt++;
				add_edge(s,hah(i,j),1);
			}
		}
	}
	for(int x1=1;x1<=n;x1++)
	{
		for(int y1=1;y1<=m;y1++)
		{
			for(int x2=1;x2<=n;x2++)
			{
				for(int y2=1;y2<=m;y2++)
				{
					if(M[x1][y1])
					{
						if(x1<=d || y1<=d || x1>=n-d+1 || y1>=m-d+1) add_edge(hah(x1,y1)+n*m,t,INF);
					} 
					if(M[x1][y1] && M[x2][y2])
					{ 
						if(x1==x2 && y1==y2) continue;
						if(dis(x1,y1,x2,y2)<=d) 
						{
							add_edge(hah(x2,y2)+n*m,hah(x1,y1),INF);
							add_edge(hah(x1,y1)+n*m,hah(x2,y2),INF);
						}
					}
				}
			}
		}
	}
	printf("%d\n",cnt-dinic(s,t));
	return 0;
}