BZOJ 1066 [SCOI2007]蜥蜴

最大流。

要求逃出去的数量最多，各种状态又十分复杂，石柱还有次数限制，直接就会想到网络流的最大流。

怎么体现次数限制（石柱高度）？拆点。对于每一个点，拆成入点和出点，连边，容量为高度。建源点和汇点。源点连所有蜥蜴，容量INF。汇点连所有可以出去的地方，容量也是INF。各个能跳到的点之间也要连。跑最大流。还有，注意这题的距离不是曼哈顿距离，是平面距离。

连边的时候手残了，强行WA一次- - 大概没救了吧

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<queue>
#define N 25
#define OPP 403
using namespace std;
const int S = 401, T = 402, INF = 0x3f3f3f3f;
struct edge{int next,to,v;}e[N*N*1000];
int cur[405*2],last[405*2],h[N][N], s[N][N], tot, cnt=1, level[405*2];
void add(int a, int b, int c)
{
    e[++cnt]=(edge){last[a],b,c};
    last[a]=cnt;
}
void insert(int a, int b, int c)
{
    add(a,b,c);
    add(b,a,0);
}
int getc()
{
    char c=getchar();
    while((c<'0'||c>'9')&&c!='L'&&c!='.')c=getchar();
    if(c>='0'&&c<='9')return c-'0';
    else if(c=='L')return -1;
    else return -2;
}
int dis(int i, int j, int x, int y)
{
    return (i-x)*(i-x)+(j-y)*(j-y);
}
int dfs(int x, int flow)
{
    if(x==T)return flow;
    int use=0;
    for(int &i = cur[x]; i; i=e[i].next)
    {
        int y=e[i].to;
        if(!e[i].v || level[y]!=level[x]+1)continue;
        int w=dfs(y,min(flow-use,e[i].v));
        e[i].v-=w;
        e[i^1].v+=w;
        use+=w;
        if(use==flow)break;
    }
    return use;
}
bool bfs()
{
    queue<int> q;
    memset(level,0,sizeof(level));
    level[S]=1;
    q.push(S);
    while(!q.empty())
    {
        int x=q.front();
        q.pop();
        for(int i = last[x]; i; i=e[i].next)
        {
            int y=e[i].to;
            if(level[y] || !e[i].v)continue;
            level[y]=level[x]+1;
            if(y==T)break; 
            q.push(y);
        }
    }
    return level[T]!=0;
}
int dinic()
{
    int ret=0;
    while(bfs())
    {
        memcpy(cur,last,sizeof(cur));
        ret+=dfs(S,INF);
    }
    return ret;
}
int main()
{
    int n, m, d, sum=0;
    scanf("%d%d%d",&n,&m,&d);
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
        {
            h[i][j]=getc();
            s[i][j]=++tot;
            insert(s[i][j],s[i][j]+OPP,h[i][j]);
        }
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
        {
            if(getc()==-1)
                insert(S,s[i][j],1),sum++;//手贱连边连成INF了，吃枣药丸orz 
            if(abs(i-1)+1<=d || abs(i-n)+1<=d || abs(j-1)+1<=d || abs(j-m)+1<=d)
                insert(s[i][j]+OPP,T,INF);
            for(int x = 1; x <= n; x++)
                for(int y = 1; y <= m; y++)
                    if(dis(i,j,x,y)<=d*d)
                        insert(s[i][j]+OPP,s[x][y],INF);
        }
    int ans = dinic();
    printf("%d\n",sum-ans);
}