[最大流] BZOJ1066: [SCOI2007]蜥蜴

题意

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

题解

近似于裸的最大流。
注意到石柱的减小其实就是点访问次数的限制。转化为最大流模型就是点有流量限制。拆点即可。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn=10005,maxe=1000005;
struct Edge {
    int from,to,cap,flow;
    Edge(int x1=0,int x2=0,int x3=0,int x4=0):from(x1),to(x2),cap(x3),flow(x4){}
} Es[maxe];
int N,n,m,d[maxn],nxt[maxe],S,T,fir[maxn],pos[maxn],tot=1;
queue <int> que;
void add(int x,int y,int w){
    Es[++tot]=Edge(x,y,w,0);
    nxt[tot]=fir[x]; fir[x]=tot;
    Es[++tot]=Edge(y,x,0,0);
    nxt[tot]=fir[y]; fir[y]=tot;
}
bool BFS(){
    memset(d,63,sizeof(d)); int INF=d[0];
    while(!que.empty()) que.pop(); 
    que.push(S); d[S]=0; 
    while(!que.empty()){
        int x=que.front(); que.pop();
        for(int j=fir[x];j;j=nxt[j]) 
        if(d[Es[j].to]==INF&&Es[j].cap>Es[j].flow){
            d[Es[j].to]=d[x]+1; que.push(Es[j].to);
        }
    }
    return d[T]!=INF;
}
int find_DFS(int x,int flow){
    if(x==T||flow==0) return flow;
    int res=0,f;
    for(int &j=pos[x];j;j=nxt[j]){
        if(d[x]+1==d[Es[j].to]&&(f=find_DFS(Es[j].to,min(flow,Es[j].cap-Es[j].flow)))>0){
            Es[j].flow+=f; Es[j^1].flow-=f; 
            res+=f; flow-=f; if(flow==0) break;
        }
    }
    return res;
}
int GetMaxFlow() {
    int MaxFlow=0,now;
    while(BFS()){
        for(int i=1;i<=N;i++) pos[i]=fir[i]; 
        MaxFlow+=find_DFS(S,1e+9);
    }
    return MaxFlow;
}
int a[30][30],all;
double D; 
inline bool is_edge(int x,int y){ return x==0||x==n+1||y==0||y==m+1; }
inline int id(int x,int y){
    if(is_edge(x,y)) return T;
    return (x-1)*(2*m)+y*2-1;
}
inline bool check(int x,int y,int xx,int yy){
    return (x-xx)*(x-xx)+(y-yy)*(y-yy)<=D*D&&!(x==xx&&y==yy);
}
int main(){
    scanf("%d%d%lf",&n,&m,&D);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
        char ch=getchar(); while(!('0'<=ch&&ch<='9')) ch=getchar();
        add(id(i,j),id(i,j)+1,(a[i][j]=ch-'0'));
     }
    N=n*m*2; S=++N; T=++N;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
        char ch=getchar(); while(ch!='.'&&ch!='L') ch=getchar();
        if(ch=='L') 
         add(S,id(i,j),1), all++;
     }
    for(int i1=0;i1<=n+1;i1++)
     for(int j1=0;j1<=m+1;j1++) if(a[i1][j1]||is_edge(i1,j1))
      for(int i2=1;i2<=n;i2++)
       for(int j2=1;j2<=m;j2++) if(a[i2][j2])
        if(check(i1,j1,i2,j2)){
            add(id(i2,j2)+1,id(i1,j1),1e+9);
            //printf("edge: %d %d -> %d %d\n",i2,j2,is_edge(i1,j1)?0:i1,is_edge(i1,j1)?0:j1);
        }
    printf("%d\n",all-GetMaxFlow());
    return 0;
}