洛谷2472 [SCOI2007]蜥蜴

Luogu

BZOJ

分析

网络流。把每根柱子拆成入点和出点。

  • 源点向每只蜥蜴连容量为$1$的边。
  • 每根柱子的入点向出点连容量为高度的边。
  • 每根柱子的出点向能到达的所有柱子的入点连容量为$\infty$的边。
  • 所有能跳出边界的柱子的出点向汇点连容量为$\infty$的边。

然后用总数减掉最大流,就做完了。

代码

//It is made by M_sea
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>
#define re register
using namespace std;

inline int read() {
    int X=0,w=1; char c=getchar();
    while (c<'0'||c>'9') { if (c=='-') w=-1; c=getchar(); }
    while (c>='0'&&c<='9') X=X*10+c-'0',c=getchar();
    return X*w;
}

const int N=20+10;
const int INF=2147483647;

int n,m,d,tot=0,s,t;
int h[N][N];
char c[N][N];

struct Edge { int v,w,nxt; };
Edge e[2000000];
int head[100000],cnt=2;

inline void addEdge(int u,int v,int w) {
    e[cnt].v=v,e[cnt].w=w,e[cnt].nxt=head[u],head[u]=cnt++;
    e[cnt].v=u,e[cnt].w=0,e[cnt].nxt=head[v],head[v]=cnt++;
}

int level[100000];

inline int bfs() {
    memset(level,0,sizeof(level));
    queue<int> Q; Q.push(s); level[s]=1;
    while (!Q.empty()) {
        int u=Q.front(); Q.pop();
        for (re int i=head[u];i;i=e[i].nxt) {
            int v=e[i].v,w=e[i].w;
            if (w&&!level[v])
                level[v]=level[u]+1,Q.push(v);
        }
    }
    return level[t]!=0;
}

inline int dfs(int u,int cpflow) {
    if (u==t) return cpflow;
    int addflow=0;
    for (re int i=head[u];i;i=e[i].nxt) {
        int v=e[i].v;
        if (e[i].w&&level[v]==level[u]+1) {
            int tmpadd=dfs(v,min(e[i].w,cpflow-addflow));
            e[i].w-=tmpadd,e[i^1].w+=tmpadd;
            addflow+=tmpadd;
        }
    }
    if (!addflow) level[u]=0;
    return addflow;
}

inline int dinic() {
    int maxflow=0;
    while (bfs()) maxflow+=dfs(s,INF);
    return maxflow;
}

inline int g(int x,int y) { return (x-1)*m+y; }

int main() {
    n=read(),m=read(),d=read(),s=0,t=2*n*m+1; int sum=0;
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j) {
            scanf("%1d",&h[i][j]);
            if (!h[i][j]) continue;
            addEdge(g(i,j),g(i,j)+n*m,h[i][j]);
            if (i-d<1||i+d>n||j-d<1||j+d>m)
                addEdge(g(i,j)+n*m,t,INF);
        }
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j) {
            c[i][j]=getchar();
            while (c[i][j]!='.'&&c[i][j]!='L') c[i][j]=getchar();
            for (re int x=1;x<=n;++x)
                for (re int y=1;y<=m;++y)
                    if (h[x][y]&&(x!=i||y!=j)&&(x-i)*(x-i)+(y-j)*(y-j)<=d*d)
                        addEdge(g(i,j)+n*m,g(x,y),INF);
            if (c[i][j]=='L') ++sum,addEdge(s,g(i,j),1);
        }
    printf("%d\n",sum-dinic());
    return 0;
}
最后修改:2019 年 09 月 24 日 10 : 13 PM

1 条评论

  1. smy

    %%%天天秒题

发表评论