1. 程式人生 > >洛谷2472 蜥蜴(網路流)

洛谷2472 蜥蜴(網路流)

傳送門

【題目分析】

令人智熄的字串讀入操作。。。。。在BZOJ上過了然後氵谷全T?emmm。。。。。。

網路流的題難點就在於建圖,這道題還是比較明顯,首先每個點還是要拆點限制流量,上限設為石柱高度表示最多可以跳過這麼多蜥蜴。

建立一個起點s和終點t,每個有蜥蜴的石柱對應的入點就和起點s連線,s與入點間的弧容量為1(含義是該石柱只能有一個蜥蜴,如果容量大於1,就有可能在跑最大流時出現一個石柱有多個蜥蜴的情況),每個與邊界距離足夠近的石柱對應的出點和終點t連線,t與出點間的弧容量為INF,這樣t與出點間的弧容量不會影響到結果,此題的圖就建立完成了,再跑一遍最大流,最大流等於能逃脫的蜥蜴數量,輸出時用總蜥蜴數減去能逃脫的蜥蜴數即可。

【程式碼~】

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e4+10;
const int MAXM=5e5+10;
const int INF=0x3f3f3f3f;

int n,m,cnt,s,t,mind,ans;
int head[MAXN],cur[MAXN],depth[MAXN];
int nxt[MAXM],to[MAXM],w[MAXM];
int mx[MAXN],my[MAXN],tot;
char x;

int Read()
{
	int i=0,f=1;
	char c;
	for(c=getchar();(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')
	  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())
	  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void Add(int x,int y,int z)
{
	nxt[cnt]=head[x];
	head[x]=cnt;
	to[cnt]=y;
	w[cnt]=z;
	cnt++;
}

void add(int x,int y,int z)
{
	Add(x,y,z);
	Add(y,x,0);
}

bool bfs()
{
	queue<int> q;
	memset(depth,0,sizeof(depth));
	depth[s]=1;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=head[u];i!=-1;i=nxt[i])
		{
			int v=to[i];
			if(!depth[v]&&w[i])
			{
				depth[v]=depth[u]+1;
				q.push(v);
			}
		}
	}
	return depth[t]!=0;
}

int dfs(int u,int dist)
{
	if(u==t)
	  return dist;
	for(int &i=cur[u];i!=-1;i=nxt[i])
	{
		int v=to[i];
		if(depth[v]==depth[u]+1&&w[i])
		{
			int di=dfs(v,min(w[i],dist));
			if(di>0)
			{
				w[i]-=di;
				w[i^1]+=di;
				return di;
			}
		}
	}
	return 0;
}

int dinic()
{
	int ans=0;
	while(bfs())
	{
		for(int i=s;i<=t;++i)
		  cur[i]=head[i];
		while(int d=dfs(s,INF))
		  ans+=d;
	}
	return ans;
}

int dist(int x1,int y1,int x2,int y2)
{
	return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
}

int zh(int x,int y)
{
	return (x-1)*m+y;
}

int main()
{
	memset(head,-1,sizeof(head));
	scanf("%d%d%d\n",&n,&m,&mind);
	s=0,t=n*m*2+1;
	for(int i=1;i<=n;++i)
	{
		for(int j=1;j<=m;++j)
		{
		  	cin>>x;
		  	if(x!='0')
		  	{
		  	  	add(zh(i,j),zh(i,j)+n*m,x-'0');
		  	  	if(i<=mind||j<=mind||i+mind>n||j+mind>m)
		  	  	  add(zh(i,j)+n*m,t,INF);
		  	  	tot++;
		  	  	mx[tot]=i,my[tot]=j;
		  	}
		}
	}
	for(int i=1;i<=n;++i)
	{
	    for(int j=1;j<=m;++j)
	    {
	  	    cin>>x;
	  	    if(x=='L')
	  	      ans++,add(s,zh(i,j),1);
	    }
	}
	for(int i=1;i<=tot;++i)
	  for(int j=1;j<=tot;++j)
	    if(i!=j)
	    {
	    	if(dist(mx[i],my[i],mx[j],my[j])<=mind*mind)
	    	  add(zh(mx[i],my[i])+n*m,zh(mx[j],my[j]),INF);
	    }
	cout<<ans-dinic();
	return 0;
}