luogu2472&&bzoj1066 [SCOI2007]蜥蜴

(http://www.elijahqi.win/2017/11/23/luogu2472bzoj1066-scoi2007%E8%9C%A5%E8%9C%B4/)

题目背景

07四川省选

题目描述

在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃到边界外。

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

输入输出格式

输入格式：

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

输出格式：

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

输入输出样例

输入样例#1： 复制

5 8 2

00000000

02000000

00321100

02000000

00000000

……..

……..

..LLLL..

……..

……..

输出样例#1： 复制

1

说明

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

建边 每个柱子拆点 用来限制跳跃次数

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cmath>
#define inf 0x3f3f3f3f
using namespace std;
int h[1000],num=1;
int r,c,d,ans1,mp[22][22],level[1000];char s1[22];
struct node{
int x,y,next,z;
}data[644000];
inline void insert1(int x,int y,int z){
data[++num].y=y;data[num].z=z;data[num].next=h[x];h[x]=num;data[num].x=x;
data[++num].y=x;data[num].z=0;data[num].next=h[y];h[y]=num;data[num].x=y;
}
inline int calc(int i,int j,int op){return (i-1)*c*2+(j-1)*2+op;}
inline bool bfs(){
memset(level,0,sizeof(level));queue<int>q;q.push(0);level[0]=1;
while (!q.empty()){
int x=q.front();q.pop();
for (int i=h[x];i;i=data[i].next){
int y=data[i].y,z=data[i].z;
if (level[y]||!z) continue;q.push(y);
level[y]=level[x]+1;if (y==900) return true;
}
}return false;
}
inline bool check(int x,int y,int x1,int y1){
if (!mp[x1][y1]) return false;
if ((x1-x)*(x1-x)+(y1-y)*(y1-y)>d*d) return false;else return true;
}
int dfs(int x,int s){
if (x==900) return s;int ss=s;
for (int i=h[x];i;i=data[i].next){
int y=data[i].y,z=data[i].z;
if (level[x]+1==level[y]&&z){
int xx=dfs(y,min(z,s));
s-=xx;data[i].z-=xx;data[i^1].z+=xx;
if (!s) return ss;
}
}return ss-s;
}
int main(){
freopen("2472.in","r",stdin);
scanf("%d%d%d",&r,&c,&d);
for (int i=1;i<=r;++i){
scanf("%s",s1+1);
for (int j=1;j<=c;++j) mp[i][j]=s1[j]-'0';
}
for (int i=1;i<=r;++i){
scanf("%s",s1+1);
for (int j=1;j<=c;++j) if (s1[j]=='L') insert1(0,calc(i,j,1),1),ans1++;
}
for (int i=1;i<=r;++i){
for (int j=1;j<=c;++j){
if (!mp[i][j]) continue;
for (int x=1;x<=r;++x){
for (int y=1;y<=c;++y){
if (i==x&&j==y) continue;
if (check(i,j,x,y)){
insert1(calc(i,j,2),calc(x,y,1),inf);
}
}
}
if (i<=d||abs(r-i+1)<=d||j<=d||abs(c-j+1)<=d)insert1(calc(i,j,2),900,inf);
}
}int ans=0;
for (int i=1;i<=r;++i)
for (int j=1;j<=c;++j)
if (mp[i][j])insert1(calc(i,j,1),calc(i,j,2),mp[i][j]);
//  for (int i=2;i<=num;++i) printf("%d %d %d\n",data[i].x,data[i].y,data[i].z);
while(bfs())
ans+=dfs(0,inf);
printf("%d",ans1-ans);
return 0;
}