bzoj 1066: [SCOI2007] 蜥蜴

    这道题还是挺好想的，但我一开始还是想错了……

    把每个石柱拆成两个点，一个入度，一个出度，两个点连一条容量为高度的边，这样就可以限制从此石柱上经过的蜥蜴的数量。关于蜥蜴是否单独成点，我是单独当成了一个点，貌似做麻烦了，可以直接源点连石柱，但那样我想会不会造成一些问题，貌似也没有。

    虽然很水，但还是调了很久。主要问题出在建图上，我把一个点拆成了高度个点，这样无法达到上面说的限制蜥蜴经过的数量这个功能，所以ＷＡ了很久，看了题解，才突然明白，这么搞不行……

    代码如下：

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <queue>
#define N 25
#define inf 1<<30
using namespace std;

struct sss
{
    int x,y;
}xiyi[N*N];
int n,m,S,T,d;
int a[N][N]={0},num[N][N]={0},stonenum=0,xiyinum=0;
int p[N*N*5],next[N*N*N*N*6],v[N*N*N*N*6],f[N*N*N*N*6],bnum=-1;

bool kexing(int x,int y)
{
    if (x<1||x>n||y<1||y>m) return false;
    else return true;
}

void addbian(int x,int y,int flow)
{
    bnum++; next[bnum]=p[x];
    p[x]=bnum; v[bnum]=y; f[bnum]=flow;
    bnum++; next[bnum]=p[y];
    p[y]=bnum; v[bnum]=x; f[bnum]=0;
}

int dis[N*N*5];

bool BFS()
{
    int i,j,k;
    queue<int> q;
    for (i=1;i<=T;i++) dis[i]=-1;
    dis[S]=0; q.push(S);
    while (!q.empty())
    {
        j=q.front(); q.pop(); k=p[j];
        while (k!=-1)
        {
            if (f[k]&&dis[v[k]]==-1)
            {
                dis[v[k]]=dis[j]+1;
                q.push(v[k]);
            }
            k=next[k];
        }
    }
    if (dis[T]==-1) return false;
    else return true;
}

int DFS(int now,int change)
{
    int i,j,k,flow=0;
    if (now==T||change==0) return change;
    k=p[now];
    while (k!=-1)
    {
        if (f[k]&&dis[v[k]]==dis[now]+1&&(i=DFS(v[k],min(change,f[k])))>0)
        {
            f[k]-=i; f[k^1]+=i;
            flow+=i; change-=i;
            if (change==0) break;
        }
        k=next[k];
    }
    dis[now]=-1;
    return flow;
}

int dinic()
{
    int ans=0,flow;
    while (BFS())
        while (flow=DFS(S,inf))
            ans+=flow;
    return ans;
}

bool bianjie(int x,int y)
{
    if (x<=d||y<=d) return true;
    else if (n-x<d||m-y<d) return true;
    else return false;
}

int main()
{
    int i,j,k,x,y;
    char s[N];
    scanf("%d%d%d",&n,&m,&d);
    S=n*m*3+1; T=n*m*3+2;
    for (i=1;i<=T;i++) p[i]=-1;
    for (i=1;i<=n;i++)
    {
        scanf("%s",s);
        for (j=0;j<m;j++)
            if (s[j]!='0') a[i][j+1]=s[j]-'0';
    }
    for (i=1;i<=n;i++)
    {
        scanf("%s",s);
        for (j=0;j<m;j++)
            if (s[j]=='L')
            {
                xiyinum++; a[i][j+1]--;
                xiyi[xiyinum].x=i; xiyi[xiyinum].y=j+1;
            }
    }
    for (i=1;i<=n;i++)
        for (j=1;j<=m;j++)
            if (a[i][j]>0)
            {
                stonenum++;
                num[i][j]=stonenum;
            }
    for (i=1;i<=n;i++)
        for (j=1;j<=m;j++)
            if (a[i][j]>0)
            {
                addbian(num[i][j],num[i][j]+n*m,a[i][j]);
                if (bianjie(i,j))
                    addbian(num[i][j]+n*m,T,inf);
                for (int I=-d;I<=d;I++)
                    for (int J=-d;J<=d;J++)
                    if (!(I==0&&J==0)&& I*I+J*J<=d*d)
                    {
                        x=i+I; y=j+J;
                        if (kexing(x,y)&&a[x][y])
                            addbian(num[i][j]+n*m,num[x][y],inf);
                    }
            }
    for (i=1;i<=xiyinum;i++)
    {
        int x1=xiyi[i].x,y1=xiyi[i].y;
        addbian(S,i+n*m*2,1);
        if (bianjie(x1,y1))
            addbian(i+n*m*2,T,inf);
        for (int I=-d;I<=d;I++)
            for (int J=-d;J<=d;J++)
                if (!(I==0&&J==0)&& I*I+J*J<=d*d)
                {
                    x=x1+I; y=y1+J;
                    if (kexing(x,y)&&a[x][y])
                        addbian(i+n*m*2,num[x][y],inf);
                }
    }
    printf("%d\n",xiyinum-dinic());
}