BZOJ2150【最小路径覆盖】

最新推荐文章于 2018-10-17 07:55:00 发布

Lethelody

最新推荐文章于 2018-10-17 07:55:00 发布

阅读量437

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/Lethelody/article/details/44865785

版权

网络流专栏收录该内容

15 篇文章 0 订阅

订阅专栏

/* I will wait for you */

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iostream>
#include<fstream>
#include<vector>
#include<queue>
#include<deque>
#include<set>
#include<map>
#include<string>
#define make make_pair
#define fi first
#define se second

using namespace std;

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;

const int maxn=10010;
const int maxm=1000000;
const int maxs=100;
const int INF=1<<29;
const int P=1000000007;
const double error=1e-9;

struct edge{ int v,a,next; }e[maxm];

int s,t,n,m,c,r,w,cnt=0,ans=0;
int head[maxn],dis[maxn],cur[maxn],q[maxm],a[maxs][maxs];

void addedge(int u,int v,int a)
{
	e[cnt]=(edge){v,a,head[u]},head[u]=cnt++;
	e[cnt]=(edge){u,0,head[v]},head[v]=cnt++;
}

bool bfs()
{
	memset(dis,-1,sizeof(dis));
	int he=1,ta=1;dis[q[1]=s]=0;
	
	while(he<=ta)
	{
		int u=q[he++];
		for(int i=head[u],v;i!=-1;i=e[i].next)
			if(e[i].a&&dis[v=e[i].v]==-1) dis[v]=dis[u]+1,q[++ta]=v;
	}
	return dis[t]!=-1;
}

int dfs(int u,int a)
{
	if(u==t) return a;
	int f,flow=0;
	for(int i=cur[u];a&&i!=-1;i=e[i].next)
	{
		int v=e[i].v;
		if(e[i].a&&dis[v]==dis[u]+1&&(f=dfs(v,min(a,e[i].a)))>0)
		{
			e[i].a-=f,e[i^1].a+=f,a-=f,flow+=f;
			if(e[i].a) cur[u]=i;
		}
	}
	if(!flow) dis[u]=-1;
	return flow;
}

void dinic()
{
	while(bfs()) memcpy(cur,head,sizeof(head)),ans-=dfs(s,INF);
}

int main()
{
	scanf("%d%d%d%d",&n,&m,&r,&c);
	
	w=n*m;s=2*w;t=2*w+1;
	
	for(int i=0;i<n;i++)
	{
		char s[100];scanf("%s",s);
		for(int j=0;j<m;j++) a[i][j]=s[j]=='.'?1:0;
	}
		
	memset(head,-1,sizeof(head));
	
	for(int i=0;i<n;i++) for(int j=0;j<m;j++) if(a[i][j]) 
	{
		int u=i*m+j;
		if(i+r<n)
		{
			if(j+c<m&&a[i+r][j+c]) addedge(u,w+(i+r)*m+(j+c),1);
			if(j-c>=0&&a[i+r][j-c]) addedge(u,w+(i+r)*m+(j-c),1);
		}
		if((i+r>=n||r!=c)&&i+c<n)
		{
			if(j+r<m&&a[i+c][j+r]) addedge(u,w+(i+c)*m+(j+r),1);
			if(j-r>=0&&a[i+c][j-r]) addedge(u,w+(i+c)*m+(j-r),1);
		}
	}
	
	for(int i=0;i<n;i++) for(int j=0;j<m;j++) 
		if(a[i][j]) ans++,addedge(s,i*m+j,1),addedge(w+i*m+j,t,1);
	
	dinic();printf("%d\n",ans);

	return 0;
}