BZOJ 4554: [Tjoi2016&Heoi2016]游戏

Description

在2016年，佳缘姐姐喜欢上了一款游戏，叫做泡泡堂。简单的说，这个游戏就是在一张地图上放上若干个炸弹，看是否能炸到对手，或者躲开对手的炸弹。在玩游戏的过程中，小H想到了这样一个问题：当给定一张地图，在这张地图上最多能放上多少个炸弹能使得任意两个炸弹之间不会互相炸到。炸弹能炸到的范围是该炸弹所在的一行和一列，炸弹的威力可以穿透软石头，但是不能穿透硬石头。给定一张n*m的网格地图:其中代表空地，炸弹的威力可以穿透，可以在空地上放置一枚炸弹。x代表软石头，炸弹的威力可以穿透，不能在此放置炸弹。#代表硬石头，炸弹的威力是不能穿透的，不能在此放置炸弹。例如：给出1*4的网格地图*xx，这个地图上最多只能放置一个炸弹。给出另一个1*4的网格地图x#，这个地图最多能放置两个炸弹。现在小H任意给出一张n*m的网格地图，问你最多能放置多少炸弹

Input

第一行输入两个正整数n,m,n表示地图的行数，m表示地图的列数。1≤n,m≤50。接下来输入n行m列个字符，代表网
格地图。*的个数不超过n*m个

Output

输出一个整数a，表示最多能放置炸弹的个数

Sample Input

4 4

#***

*#**

**#*

xxx#

Sample Output

5

分析

第一次写拆点网络流，调了半天，拍也拍不出错，结果发现，边集数组开小了。

代码

#include <bits/stdc++.h>

#define INF 0x7fffffff
#define N 55

struct NOTE
{
    int to,next,c;
}e[N * N * 2];

int cnt = 1;
int ans;
int next[N * N];

void add(int x,int y,int f)
{
    e[++cnt].to = y;
    e[cnt].c = f;
    e[cnt].next = next[x];
    next[x] = cnt;

    e[++cnt].to = x;
    e[cnt].c = 0;
    e[cnt].next = next[y];
    next[y] = cnt;
}

int s,t;
int dis[N * N],cur[N * N];

bool bfs()
{
    for (int i = s; i <= t; i++)
        dis[i] = 0;
    dis[s] = 1;
    std::queue<int> Q;
    Q.push(s);
    while (!Q.empty())
    {
        int u = Q.front();
        Q.pop();
        for (int i = next[u]; i; i = e[i].next)
        {
            if (e[i].c && !dis[e[i].to])
            {
                dis[e[i].to] = dis[u] + 1;
                if (e[i].to == t)
                    return 1;
                Q.push(e[i].to);
            }
        }
    }
    return 0;
}

int dfs(int x,int maxf)
{
    if (x == t || !maxf)
        return maxf;
    int ret = 0;
    for (int &i = cur[x]; i; i = e[i].next)
        if (e[i].c && dis[e[i].to] == dis[x] + 1)
        {
            int f = dfs(e[i].to,std::min(e[i].c,maxf - ret));
            ret += f;
            e[i].c -= f;
            e[i ^ 1].c += f;
            if (ret == maxf)
                break;
        }
    return ret;
}

void dinic()
{
    while (bfs())
    {
        for (int i = s; i <= t; i++)
            cur[i] = next[i];
        ans += dfs(s,INF);
    }
}

int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    char ch[N];
    int mapx[N][N],mapy[N][N],map[N][N];
    int x = 0,y = 0;
    for (int i = 1; i <= n; i++)
    {
        scanf("%s",ch + 1);
        for (int j = 1; j <= m; j++)
            if (ch[j] == '*')
                map[i][j] = 0;
                else
                if (ch[j] == 'x')
                    map[i][j] = 1;
                    else map[i][j] = 2;
    }
    for (int i = 1; i <= n; i++)
    {
        for (int j = 1; j <= m; j++)
        {
            if (!mapx[i][j] && map[i][j]!=2)
            {
                x++;
                int k=j;
                while (k<=m && !mapx[i][k] && map[i][k]!=2)
                    mapx[i][k]=x,k++;
            }
        }
    }
    y = x;
    for (int j = 1; j <= m; j++)
        for (int i = 1; i <= n; i++)
            if (!mapy[i][j] && map[i][j]!=2)
            {
                y++;
                int k = i;
                while (k <= n && !mapy[k][j] && map[k][j]!=2)
                    mapy[k][j] = y,k++;
            }
    s = 0;
    t = y + 1;
    for (int i = 1; i <= x; i++)
        add(s,i,1);
    for (int i = x + 1; i <= y; i++)
        add(i,t,1);
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= m; j++)
            if (!map[i][j])
                add(mapx[i][j],mapy[i][j],1);

    dinic();

    printf("%d\n",ans);
}