BZOJ 1066 POJ 2711 [SCOI2007]蜥蜴

1066: [SCOI2007]蜥蜴

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 2944  Solved: 1471
[Submit][Status][Discuss]

Description

　　在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃
到边界外。 每行每列中相邻石柱的距离为1，蜥蜴的跳跃距离是d，即蜥蜴可以跳到平面距离不超过d的任何一个石
柱上。石柱都不稳定，每次当蜥蜴跳跃时，所离开的石柱高度减1（如果仍然落在地图内部，则到达的石柱高度不
变），如果该石柱原来高度为1，则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个
石柱上。

Input

　　输入第一行为三个整数r，c，d，即地图的规模与最大跳跃距离。以下r行为石竹的初始状态，0表示没有石柱
，1~3表示石柱的初始高度。以下r行为蜥蜴位置，“L”表示蜥蜴，“.”表示没有蜥蜴。

Output

　　输出仅一行，包含一个整数，即无法逃离的蜥蜴总数的最小值。

Sample Input

5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........

Sample Output

1

HINT

100%的数据满足：1<=r, c<=20, 1<=d<=4

Source

【题目分析】

   

参考了mhr的ppt才做出来orzzzzzzzz。。。感觉网络流里拆点的方法很有趣啊

【代码】

#include <cstdio>#include <cstring>#include <string>#include <iostream>#include <queue>#define INF 10000000#define M(a) memset(a,0,sizeof a)#define F(i,j,k) for (int i=j;i<=k;++i)using namespace std;int list[51][51],num[51][51];int n,m,d,en=0;int S=0,T,sum=0,all=0;int u[100001],v[100001],h[100001],ne[100001],f[100001],map[100001];inline void add(int a,int b,int r){u[en]=a;v[en]=b;ne[en]=h[a];f[en]=r;h[a]=en++;u[en]=b;v[en]=a;ne[en]=h[b];f[en]=0;h[b]=en++;}inline bool tell(){memset(map,-1,sizeof map);queue <int> q;q.push(S);map[S]=0;while (!q.empty()){int u=q.front(); q.pop();for (int i=h[u];i!=-1;i=ne[i]){if (map[v[i]]==-1&&f[i]){map[v[i]]=map[u]+1;q.push(v[i]);}}}if (map[T]!=-1) return true;else return false;}inline int zeng (int k,int now){if (k==T) return now;int r=0;for (int i=h[k];i!=-1&&now>r;i=ne[i]){if (map[k]+1==map[v[i]]&&f[i]!=0){int t=zeng(v[i],min(now-r,f[i]));f[i]-=t;f[i^1]+=t;r+=t;}}if (!r) map[k]=-1;return r;}inline int dinic(){int r=0,t;while (tell()) while (t=zeng(S,INF)) r+=t;return r; }inline bool ok(int a,int b,int c,int e){if ((c-a)*(c-a)+(e-b)*(e-b)<=d*d) return true;return false;}int main(){memset(h,-1,sizeof h);memset(ne,-1,sizeof ne);memset(u,-1,sizeof u);memset(v,-1,sizeof v);scanf("%d%d%d",&n,&m,&d);F(i,1,n)F(j,1,m){char ch;cin>>ch;list[i][j]=ch-'0';if (list[i][j]!=0){num[i][j]=++sum;}}T=sum*2+1;F(i,1,n)F(j,1,m){char ch;cin>>ch;if (ch=='L') {add(S,num[i][j]*2,1);all++;}}F(i,1,n)F(j,1,m) if ((i<=d||j<=d||i+d>n||j+d>m)&&(num[i][j])){add(num[i][j]*2+1,T,INF);}F(i,1,n)F(j,1,m) if (num[i][j]) {add(num[i][j]*2,num[i][j]*2+1,list[i][j]);}F(i,1,n)F(j,1,m){if (num[i][j]){F(x,min(i-d,1),min(i+d,n)) F(y,min(j-d,1),min(j+d,m))  if (num[x][y]&&(i!=x||j!=y)) if (ok(i,j,x,y)) {  add(num[i][j]*2+1,num[x][y]*2,INF);  }}}cout<<all-dinic()<<endl;}