SCOI2007[蜥蜴]
2016-02-23 22:01
363 查看
[SCOI2007]蜥蜴
Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 2702 Solved: 1339
[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 200000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........
Sample Output
1HINT
100%的数据满足:1<=r, c<=20, 1<=d<=4Source
<a href="http://www.lydsy.com/JudgeOnline/problemset.php?search=Pku%202711%20Leapin" lizards'="" style="color: blue; text-decoration: none;">Pku 2711 Leapin' Lizards[Submit][Status][Discuss]
看了一下题意就反应过来是最大流了,解题思路大致是想到了,至于处理经过单点的次数,这个需要拆点处理。
剩下的建边:
所有有蜥蜴的和S联通,所有石柱能到达的点联通,可以跳出边界的和T联通。
这样建边这样的边似乎有点多.
读入的时候用scanf读所有字符串再查ASCII值,不造为毛被卡掉了- - TLE~
主要分析好了建边,就是裸的最大流了,套ISAP
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
const int inf=0x3f3f3f3f;
int n,m,d;
int a[25][25];
struct node
{
int v,cap,next;
}e[100000];
int head[1000];//拆点
int dis[1000];
int k;
int S,T;
void add(int u,int v,int cap)
{
e[k].v=v;
e[k].cap=cap;
e[k].next=head[u];
head[u]=k++;
e[k].v=u;
e[k].cap=0;
e[k].next=head[v];
head[v]=k++;
}
bool check(int x1,int y1,int x2,int y2)
{
if((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)<=d*d)
{
if(a[x2][y2]!=0)return true;
}
return false;
}
int cnt;
bool bfs()
{
memset(dis,-1,sizeof(dis));
queue<int >q;
q.push(S);
dis[S]=0;
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=head[u];~i;i=e[i].next)
{
int v=e[i].v;
if(dis[v]==-1&&e[i].cap>0)
{
dis[v]=dis[u]+1;
if(v==T)return true;
q.push(v);
}
}
}
return false;
}
int dfs(int u,int low)
{
if(u==T||low==0)return low;
int ret=0;
for(int i=head[u];~i;i=e[i].next)
{
int v=e[i].v;
if(dis[v]==dis[u]+1)
{
int f=dfs(v,min(low,e[i].cap));
e[i].cap-=f;
e[i^1].cap+=f;
low-=f;
ret+=f;
if(low==0)return ret;
}
}
dis[u]=-1;
return ret;
}
void work()
{
int ans=0;
while(bfs())
{
int x;
if(x=dfs(S,inf))ans+=x;
}
printf("%d\n",cnt-ans);
}
int main()
{
memset(head,-1,sizeof(head));
k=0;
cnt=0;
scanf("%d%d%d",&n,&m,&d);
S=0,T=n*m*2+1;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
scanf("%1d",&a[i][j]);
}
}
char C=getchar();
for(int i=1;i<=n;i++)
{
//scanf("%s",ss[i]+1);
for(int j=1;j<=m;j++)
{
C=getchar();
if(C=='L')
{
cnt++;
int id=(i-1)*m+j;
add(S,id*2-1,1);
}
}
C=getchar();
}
//能跳到的点
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(a[i][j]==0)continue;
for(int x=1;x<=n;x++)
{
for(int y=1;y<=m;y++)
{
if(i==x&&j==y)continue;
if(check(i,j,x,y))
{
int id=(i-1)*m+j;
int nid=(x-1)*m+y;
add(id*2,nid*2-1,inf);
}
}
}
}
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(a[i][j])
{
int id=(i-1)*m+j;
add(id*2-1,id*2,a[i][j]);//拆点建边
if(i-d<1||i+d>n||j-d<1||j+d>m)//和边界相连
{
add(id*2,T,inf);//和终点相连
}
}
}
}
work();
return 0;
}
相关文章推荐
- 【网络流-最大流-Dinic-建模】POJ3281 Dining:Pascal 解法
- Ford_Fulkerson算法
- 重标记与前移算法
- 网络流之增广路算法
- [笔记] 网络流-最大流 POJ-1273\HDU-4240
- 最大流/二匹配
- 最大流之预流推进
- 最大流 : FordFulkerson 算法
- [BZOJ1834][ZJOI2010][最大流][最小费用最大流]网络扩容
- POJ3281Dining解题报告
- poj_3436
- poj_1459
- PIGS
- [USACO 4.2.1] Drainage Ditches
- [wikioi 1022]覆盖
- [USACO 4.2.1] Drainage Ditches
- uva1660 最大流
- uva820 最大流模板题
- uva11082 最大流和网络流建模
- uva753 folyd EK算法