bzoj1066 [SCOI2007]蜥蜴（建图最大流）

把每个点拆成两个点，一个表示进入这个点，记作a，一个表示离开这个点，记作b。a->b，边权为高度，表示最多只有这么多的蜥蜴可以通过这个点。然后任意两个可以互达的石柱x,y，建边bx->ay,边权为inf，建立超级源点，向所有一开始有蜥蜴的石柱建边，边权为inf，建立超级汇点，所有可以跳出地图的点向超级汇点连边，边权为inf。然后跑dinic就好啦。

#include <bits/stdc++.h>using namespace std;#define ll long long#define N 1000#define inf 0x3f3f3f3finline int read(){    int x=0,f=1;char ch=getchar();    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();    return x*f;}int n,m,d,mp[21][21],h[N],num=1,tot=0,lev[N];struct edge{    int to,next,val;}data[170010];inline void add(int x,int y,int val){    data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val;    data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].val=0;}inline bool bfs(){    queue<int>q;memset(lev,0,sizeof(lev));    q.push(0);lev[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].to;if(lev[y]||!data[i].val) continue;            lev[y]=lev[x]+1;q.push(y);        }    }return lev[900];}int dinic(int x,int low){    if(x==900) return low;int tmp=low;    for(int i=h[x];i;i=data[i].next){        int y=data[i].to;if(lev[y]!=lev[x]+1||!data[i].val) continue;        int res=dinic(y,min(low,data[i].val));        if(!res) lev[y]=0;tmp-=res;data[i].val-=res;data[i^1].val+=res;        if(!tmp) return low;    }return low-tmp;}inline int maxflow(){    int res=0;    while(bfs()) res+=dinic(0,inf);return res;}int main(){//  freopen("a.in","r",stdin);    n=read();m=read();d=read();    for(int i=1;i<=n;++i){        char s[30];scanf("%s",s+1);        for(int j=1;j<=m;++j) mp[i][j]=s[j]-'0';    }for(int x=1;x<=n;++x)        for(int y=1;y<=m;++y){            if(!mp[x][y]) continue;            if(x-d<1||y-d<1||x+d>n||y+d>m) add((x-1)*m+y,900,inf);        }    for(int x1=1;x1<=n;++x1)        for(int y1=1;y1<=m;++y1){            if(!mp[x1][y1]) continue;add((x1-1)*m+y1+400,(x1-1)*m+y1,mp[x1][y1]);            for(int x2=1;x2<=n;++x2)                for(int y2=1;y2<=m;++y2){                    if(!mp[x2][y2]||(x1==x2&&y1==y2)||abs(x2-x1)+abs(y2-y1)>d) continue;                    add((x1-1)*m+y1,(x2-1)*m+y2+400,inf);                }        }    for(int i=1;i<=n;++i){        char s[30];scanf("%s",s+1);        for(int j=1;j<=m;++j) if(s[j]=='L') add(0,(i-1)*m+j+400,1),++tot;    }printf("%d\n",tot-maxflow());    return 0;}