【题解】P2472 [SCOI2007]蜥蜴

tips:我习惯把边的容量叫成流量……见谅吧

这道题刚学完网络流就开始肝了……一直到上上周才肝完（

果然我是全机房最蒟蒻之人（

---

这是最经典的拆点题目之一了……吧？

看看题面：

石柱高度，每跳一次容量-1 -> 网络流中边的流量

可这是节点啊！！！

那怎么办？不做了

拆点！

对于每一个点u，我们拆成点u和点u’，从u到u’连一条流量为石柱高度的边

有蜥蜴的石柱->从原点s到对应石柱的u连一条流量为1的边

能相互到达的石柱：

对于一对石柱u,v，从u’（出点）到v（入点）连一条流量为inf的边

能出去的石柱->从u’到汇点t连一条流量为inf的边

---

坑点：

1.拆点时注意了，不要让你的原点和某个点共用一个编号（

2.确实有一个能让你的代码量简化一点的连汇点的方式，但是……容易出锅……

---

代码qwq！

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>
using namespace std;

struct note
{
    int t;
    int next;
    int val;
};

const int inf=1e9;
note e[1000010];
int head[1000010],sta[1000010],cur[1000010];
int n,m,s,t,cnt=1;

inline void add(int x,int y,int val)
{
    e[++cnt].t=y;
    e[cnt].val=val;
    e[cnt].next=head[x];
    head[x]=cnt;
}

inline void add_edge(int x,int y,int val)
{
    add(x,y,val);
    add(y,x,0);
}

bool bfs(int s)
{
    memset(sta,0,sizeof(sta));
    memcpy(cur,head,sizeof(head));
    sta[s]=1;
    queue<int> q;
    q.push(s);
    while(!q.empty())
    {
        int p=q.front();
        q.pop();
        for(int i=head[p];i+1;i=e[i].next)
        {
            int u=e[i].t;
            if(!sta[u]&&e[i].val!=0)
            {
                sta[u]=sta[p]+1;
                q.push(u);
            }
        }
    }
    return sta[t];
}

int dfs(int x,int min_n)
{
    if(x==t)
    {
        return min_n;
    }
    for(int i=cur[x];i+1;i=e[i].next)
    {
        cur[x]=i;
        int u=e[i].t;
        if(e[i].val>0&&sta[u]==sta[x]+1)
        {
            int temp=dfs(u,min(min_n,e[i].val));
            if(temp>0)
            {
                e[i].val-=temp;
                e[i^1].val+=temp;
                return temp;
            }
        }
    }
    return 0;
}

int dinic()
{
    int ret=0;
    while(bfs(s))
    {
        int temp=dfs(s,inf);
        while(temp)
        {
            ret+=temp;
            temp=dfs(s,inf);
        }
    }
    return ret;
}

int stone[410][3];
int xi[410][2];
int r,c,d,sum;
int con_s,con_x;
bool g[30][30];

bool check(int x,int y,int x1,int y1)
{
    double dis=sqrt(abs(x-x1)*abs(x-x1)+abs(y-y1)*abs(y-y1));
    if(dis<=(double)d)
    {
        return 1;
    }
    return 0;
}

int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d%d\n",&r,&c,&d);
    for(int i=1;i<=r;i++)
    {
        for(int j=1;j<=c;j++)
        {
            char ch;
            cin>>ch;
            if(ch!='0')
            {
                stone[++sum][0]=i;
                stone[sum][1]=j;
                stone[sum][2]=ch-'0';
            }
        }
    }
    con_s=sum;
    sum=0;
    for(int i=1;i<=r;i++)
    {
        for(int j=1;j<=c;j++)
        {
            char ch;
            cin>>ch;
            if(ch=='L')
            {
                xi[++sum][0]=i;
                xi[sum][1]=j;
                g[i][j]=1;
            }
        }
    }
    con_x=sum;
    t=con_s+con_s+1;
    for(int i=1;i<=con_s;i++)
    {
        int x=stone[i][0];
        int y=stone[i][1];
        add_edge(i,con_s+i,stone[i][2]);
        if(x<=d||r-x<d||y<=d||c-y<d)
        {
            add_edge(con_s+i,t,inf);
        }
        for(int j=1;j<=con_s;j++)
        {
            if(j==i)
            {
                continue;
            }
            int x1=stone[j][0];
            int y1=stone[j][1];
            if(check(x,y,x1,y1))
            {
                add_edge(con_s+i,j,inf);
            }
        }
        if(g[x][y])
        {
            add_edge(s,i,1);
        }
    }
    printf("%d\n",con_x-dinic());
    return 0;
}
/*
5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........
*/