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]={},num[N][N]={},stonenum=,xiyinum=;
int p[N*N*],next[N*N*N*N*],v[N*N*N*N*],f[N*N*N*N*],bnum=-;

bool kexing(int x,int y)
{
    if (x<||x>n||y<||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]=;
}

int dis[N*N*];

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

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

int dinic()
{
    int ans=,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*+; T=n*m*+;
    for (i=;i<=T;i++) p[i]=-;
    for (i=;i<=n;i++)
    {
        scanf("%s",s);
        for (j=;j<m;j++)
            if (s[j]!='') a[i][j+]=s[j]-'';
    }
    for (i=;i<=n;i++)
    {
        scanf("%s",s);
        for (j=;j<m;j++)
            if (s[j]=='L')
            {
                xiyinum++; a[i][j+]--;
                xiyi[xiyinum].x=i; xiyi[xiyinum].y=j+;
            }
    }
    for (i=;i<=n;i++)
        for (j=;j<=m;j++)
            if (a[i][j]>)
            {
                stonenum++;
                num[i][j]=stonenum;
            }
    for (i=;i<=n;i++)
        for (j=;j<=m;j++)
            if (a[i][j]>)
            {
                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==&&J==)&& 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=;i<=xiyinum;i++)
    {
        int x1=xiyi[i].x,y1=xiyi[i].y;
        addbian(S,i+n*m*,);
        if (bianjie(x1,y1))
            addbian(i+n*m*,T,inf);
        for (int I=-d;I<=d;I++)
            for (int J=-d;J<=d;J++)
                if (!(I==&&J==)&& I*I+J*J<=d*d)
                {
                    x=x1+I; y=y1+J;
                    if (kexing(x,y)&&a[x][y])
                        addbian(i+n*m*,num[x][y],inf);
                }
    }
    printf("%d\n",xiyinum-dinic());
}