洛谷P2472 [SCOI2007]蜥蜴——题解

题目传送门
题目大意：
有一张n*m的网格地图（习惯用n*m了），有些地方有柱子，柱子有高度，有的柱子上有蜥蜴，蜥蜴能跳到与它的欧氏距离小于等于d的柱子上（与边界的距离小于等于0可以直接跳出地图），并且起跳的柱子的高度会减1，高度减为0的柱子不能再跳到上面去。问怎样安排能使尽量多的蜥蜴跳出地图，输出最少的不能跳出去的蜥蜴个数。

---

思考过程：
题目中有一个干扰条件：“任何时刻不能有两只蜥蜴在同一个石柱上”，我们冷静思考一下就会发现，如果我们安排好了所有蜥蜴的跳跃路径，必然能保证不会有两只蜥蜴在同一时刻出现在同一个柱子上。然后这就是一个拆点网络流的裸题了。

---

具体做法：
1.将所有的点拆成入点和出点，入点向出点连一条流量为柱子高度的边
2.s（源点）向所有一开始有蜥蜴在上面的柱子的入点连一条流量为1的边
3.所有能直接跳出地图的柱子的出点向t（汇点）连一条流量为inf的边
4.相互能跳到的柱子的出点都向对方的入点连一条流量为inf的边（这个建图的时间复杂度是$O(n^4)$的，不会超）
5.跑从s到t的最大流
6.ans等于蜥蜴总数减去最大流的结果

---

代码：

#include <bits/stdc++.h>
using namespace std;

const int maxn=40,inf=1000000;
struct stu
{
    int flow,to,next;
}road[inf]; int first[maxn*maxn],cnt=1;
int n,m,d,ans=0,cnt1=0,tt;
int s,t;
char s1[25];
int dep[maxn*maxn],num[maxn][maxn];
queue <int> q;

int get(int i,int j)
{
    return (i-1)*m+j;   
}

double get_dis(int i,int j,int t1,int t2)
{
    return sqrt((i-t1)*(i-t1)+(j-t2)*(j-t2));   
}

void addedge(int x,int y,int flow)
{
    cnt++;
    road[cnt].to=y;
    road[cnt].flow=flow;
    road[cnt].next=first[x];
    first[x]=cnt;

    cnt++;
    road[cnt].to=x;
    road[cnt].flow=0;
    road[cnt].next=first[y];
    first[y]=cnt;
}

bool bfs()
{
    memset(dep,0,sizeof(dep));
    while(!q.empty()) q.pop();
    q.push(s);
    dep[s]=1;
    while(!q.empty())
    {
        int now=q.front();q.pop();
        for(int i=first[now];i;i=road[i].next)
        {
            int to=road[i].to;
            if(road[i].flow>0)
            {
                if(!dep[to])
                {
                    dep[to]=dep[now]+1;
                    q.push(to);
                }
            }
        }
    }
    if(dep[t]>0) return 1;
    return 0;
}

int dfs(int now,int flow)
{
    if(now==t) return flow;
    for(int i=first[now];i;i=road[i].next)
    {
        int to=road[i].to;
        if(road[i].flow>0&&dep[to]==dep[now]+1)
        {
            int t=dfs(to,min(flow,road[i].flow));
            if(t>0)
            {
                road[i].flow-=t;
                road[i^1].flow+=t;
                return t;   
            }
        }
    }
    return 0;
}

void max_flow()
{
    while(bfs())
    {
        while(int k=dfs(s,inf))
        {
            ans+=k;
        }
    }
}

int main()
{
    scanf("%d%d%d",&n,&m,&d);
    s=0,t=2*n*m+1,tt=n*m;
    for(int i=1;i<=n;i++)
    {
        scanf("%s",s1+1);
        for(int j=1;j<=m;j++)
        {
            num[i][j]=s1[j]-48;
            if(num[i][j]==0) continue;
            addedge(get(i,j),get(i,j)+tt,num[i][j]);    
        }
    }
    for(int i=1;i<=n;i++)
    {
        scanf("%s",s1+1);
        for(int j=1;j<=m;j++)
        {
            if(s1[j]=='L')
            {
                addedge(s,get(i,j),1);
                cnt1++; 
            }
            if(i<=d||i>=n-d+1||j<=d||j>=m-d+1)
                addedge(get(i,j)+tt,t,inf);
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            if(num[i][j]==0) continue;
            if(i<=d||i>=n-d+1||j<=d||j>=m-d+1) continue;
            for(int t1=1;t1<=n;t1++)
            {
                for(int t2=1;t2<=m;t2++)
                {
                    if(t1==i&&t2==j) continue;
                    double dis=get_dis(i,j,t1,t2);
                    if(dis<=(double)d) addedge(get(i,j)+tt,get(t1,t2),inf);
                }
            }
        }
    }
    max_flow();
    printf("%d\n",cnt1-ans);
    return 0;   
}