[BZOJ1066] [SCOI2007] 蜥蜴 - 网络流

1066: [SCOI2007]蜥蜴

Time Limit: 1 Sec   Memory Limit: 162 MB
Submit: 2937   Solved: 1467
[ 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

<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] 
HOME   Back

---

한국어   中文   فارسی   English   ไทย 

版权所有 ©2008-2012 大视野在线测评 |  湘ICP备13009380号 |  站长统计

Based on opensource project  hustoj.

        一眼网络流啊QAQ。根据题目的含义，石柱的高度就是在这个石柱上能够通过的蜥蜴数量。

        因此我们要考虑如何限制通过的蜥蜴数量。非常显然的，我们将一个石柱拆成两个点，相当于即将落地和即将起跳的过程，中间连一条流量为石柱高度的边，就能做到限制通过的蜥蜴数。随后从源到有蜥蜴的点连一条流量为1的边，代表这个点有一条蜥蜴。

        对于两两间曼哈顿距离小于d的石柱a和b，我们在这两个石柱拆出来的点的每一对即将起跳和即将落地的点上连一条双向的边，流量为INF，代表从石柱a跳到石柱b，或从石柱b跳到石柱a。如果石柱高度为，那么即使蜥蜴跳过去了也并没有什么用 = = 所以这样连边是成立的。

        最后对于到边界距离小于d的石柱，连一条边到汇，代表跳出地图，这样最大流就等于最多跳出地图的蜥蜴数。

        最后用蜥蜴总数相减即可。

#include "stdio.h"
#include "iostream"
#define inf 0x7fffffff
#define rep(f,a,b) for(f=a;f<=b;f++)
#define travel(x,f) for(f=head[x];f;f=e[f].next)
#define bk(x) x%2?x+1:x-1
using namespace std;
const int N=805,M=80005;
int read(){
	int v=0; char ch=getchar();
	while (ch<'0'||ch>'9') ch=getchar();
	while (ch<='9'&&ch>='0') { v=v*10+ch-'0'; ch=getchar();}
	return v;
}
char ma[25][25],td[25][25];
int S,T,head[N+M],n,m,tmp,le,cnt;
struct edge {
    int next,to,flow;
    void ins(int a,int b,int &c){
        to=a,flow=b;
        next=c; c=tmp;
    }
} e[2*M];
inline void add(int u,int v,int f){
    ++tmp; e[tmp].ins(v,f,head[u]);
    ++tmp; e[tmp].ins(u,0,head[v]);
}
int maxflow,q[N+M],l,r,d[N+M],inq[N+M];
bool bfs(){
    l=1,r=0; q[++r]=S; int i,u,v;
    rep(i,S,T) d[i]=inf,inq[i]=0;
    d[S]=0; inq[S]=1;
    while(l<=r){ u=q[l++]; inq[u]=0;
        travel(u,i){ v=e[i].to;
            if (e[i].flow&&d[v]>d[u]+1){
                d[v]=d[u]+1;
                if(!inq[v]){
                    q[++r]=v;
                    inq[v]=1;
                }
            }
        }
    }
    return d[T]!=inf;
}
int dfs(int v,int mf){
    int i,cost=0,p,to;
    if(v==T) return mf;
    travel(v,i){
        if(mf==0) return cost; to=e[i].to;
        if(e[i].flow&&d[to]==d[v]+1){
            p=dfs(e[i].to,min(mf,e[i].flow));
            e[i].flow-=p; e[bk(i)].flow+=p;
            mf-=p; cost+=p;
        }
    }
    if(cost==0) d[v]=-1;
    return cost;
}
int dinic(){
    while (bfs()) 
        maxflow+=dfs(S,inf);
    return maxflow;
}
int main(){
    n=read(),m=read(); le=read(); 
    S=0; T=2*n*m+1; int i,j,k,t;
    rep(i,1,n) scanf("%s",ma[i]+1);
    rep(i,1,n) scanf("%s",td[i]+1);
    rep(i,1,n) rep(j,1,m){ int q=(i-1)*m+j;
        if(td[i][j]=='L') add(S,q,1),++cnt;
        if(ma[i][j]>'0'){ 
            add(q,q+n*m,ma[i][j]-'0');
            if(i<=le||j<=le||i+le>n||j+le>m)
                add(q+n*m,T,inf);
            rep(k,max(i-le,1),min(i+le,n))
            rep(t,max(j-le+max(i-k,k-i),1),min(j+le-max(i-k,k-i),m))
                if(ma[k][t]>'0')
                add(q+n*m,(k-1)*m+t,inf);
        }
    }
    printf("%d\n",cnt-dinic());
    return 0;
}