【最大流】BZOJ 1066 [SCOI2007]蜥蜴 题解

（传送门）

题目大意

已知n只跳跃距离为D的蜥蜴被困在石柱林中，它们可以通过在石柱上跳跃来逃脱，但每个石柱都有一个通过的限制，问最少有多少只蜥蜴最终困在石柱林中。

解题分析

ans可以转换成总蜥蜴数-最大逃脱蜥蜴数。

每个石柱明显说有通过的限制，那么网络流建模，每个石柱拆点，流量限制为通过限制，建立超级源和超级汇，等等… 总之不错的网络流建模题。

解题代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 805
#define maxe 28005
#define INF (((1<<30)-1)<<1)+1
using namespace std;
int n,m,d,k,s,t,num,tot,ans,son[maxe],nxt[maxe],cap[maxe],flow[maxe],lnk[maxn],lst[maxn],dst[maxn],que[maxn];
bool vs[maxn];
void _add(int x,int y,int z){
    son[++tot]=y; cap[tot]=z; flow[tot]=0; nxt[tot]=lnk[x]; lnk[x]=tot;
    son[++tot]=x; cap[tot]=0; flow[tot]=0; nxt[tot]=lnk[y]; lnk[y]=tot;
}
char getc(){char ch=getchar(); while (('0'>ch||ch>'3')&&ch!='L'&&ch!='.') ch=getchar(); return ch;}
int geti(int id,int x,int y){return id*n*m+(x-1)*m+y;}
void _init(){
    freopen("lizard.in","r",stdin);
    freopen("lizard.out","w",stdout);
    scanf("%d%d%d",&n,&m,&d); tot=1; t=k=2*n*m+2; s=k-1;
    memset(lnk,0,sizeof(lnk));
    memset(nxt,0,sizeof(nxt));
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++){
            int x=getc()-'0';
            if (x) _add(geti(0,i,j),geti(1,i,j),x); //石柱的通过限制
        }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            if (getc()=='L') {num++; _add(s,geti(0,i,j),1);} //超级源向出事蜥蜴位置连边
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            for (int x=max(1,i-d);x<=min(n,i+d);x++)
                for (int y=max(1,j-d);y<=min(m,j+d);y++)
                    if ((i-x)*(i-x)+(j-y)*(j-y)<=d*d&&(i!=x||j!=y))
                        _add(geti(1,i,j),geti(0,x,y),INF); //相邻石柱连边
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            if (i<=d||j<=d||n-i+1<=d||m-j+1<=d) _add(geti(1,i,j),t,INF); //超级汇
}
bool _bfs(){
    memset(vs,0,sizeof(vs));
    int hed=0,til=1; vs[t-1]=1; dst[t-1]=0; que[1]=t-1;
    while (hed!=til){
        int x=que[++hed];
        for (int j=lnk[x];j;j=nxt[j])
            if (cap[j]>flow[j]&&(!vs[son[j]])){
                dst[son[j]]=dst[x]+1;
                vs[son[j]]=1; que[++til]=son[j];
                if (son[j]==k) return true;
            }
    }
    return false;
}
int _dfs(int x,int t,int now){
    if (now==0||x==t) return now; int tem=0;
    for (int j=lst[x];j;lst[x]=j=nxt[j])
        if (dst[son[j]]==dst[x]+1){
            int ew=_dfs(son[j],t,min(now,cap[j]-flow[j]));
            if (ew){
                flow[j]+=ew; flow[j^1]-=ew; now-=ew; tem+=ew;
                if (!now) break;
            }
        }
    return tem;
}
void _solve(){
    ans=0;
    while (_bfs()){
        for (int i=1;i<=k;i++) lst[i]=lnk[i];
        ans+=_dfs(s,t,INF);
    }
    printf("%d",num-ans);
}
int main()
{
    _init();
    _solve();
    return 0;
}