poj3057 Evacuation 二分图匹配

首先BFS出每个出口到每个人的距离,然后枚举时间。

因为每一秒只能出一个人,就相当于是一个人与一个<门,时间>的二元组来匹配。

枚举时间的时候逐渐地的增加右边<门,时间>的数量,直到匹配数等于人数为止。


#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cstdlib>
#include <cmath>
#include <cctype>
#include <queue>
#include <stack>
#include <set>
#include <map>
#define pb push_back
#define mp make_pair
#define N 14
using namespace std;
int n , m , ans;
char g[N][N];

int dis[48][144] , dcnt , pcnt , id[N][N];

int dx[] = {1 , 0 , -1 , 0} , dy[] = {0 , 1 , 0 , -1};
bool in(int x , int y) {return x > 0 && x <= n && y > 0 && y <= m;}
int d[N][N];
void BFS(int num , int x , int y)
{
  queue<pair<int ,int> > q;
  memset(d , -1 , sizeof(d));
  memset(dis[num] , -1 , sizeof(dis[num]));
  q.push(mp(x , y)) , d[x][y] = 0;

  while (!q.empty())
  {
    x = q.front().first , y = q.front().second , q.pop();
    for (int i = 0 ; i < 4 ; ++ i) if (in(x + dx[i] , y + dy[i]) && !~d[x + dx[i]][y + dy[i]] && g[x + dx[i]][y + dy[i]] == '.')
    {
      d[x + dx[i]][y + dy[i]] = d[x][y] + 1;
      dis[num][id[x + dx[i]][y + dy[i]]] = d[x][y] + 1;
      q.push(mp(x + dx[i] , y + dy[i]));
    }
  }
}
int re[4 * N][N * N];
bool f[4 * N][N * N];

bool find(int x)
{
  for (int i = 1 ; i <= dcnt ; ++ i) if (~dis[i][x])
    for (int j = ans ; j >= dis[i][x] ; -- j)
      if (!f[i][j])
      {
        f[i][j] = 1;
        if (!re[i][j] || find(re[i][j]))
        {
          re[i][j] = x;
          return 1;
        }
      }
  return 0;
}

void work()
{
  int i , j , k , sum = 0;
  scanf("%d%d",&n,&m);
  pcnt = dcnt = ans = 0;
  memset(g , 0 , sizeof(g));
  memset(id , 0 , sizeof(id));

  for (i = 1 ; i <= n ; ++ i)
    scanf("%s" , g[i] + 1);
  for (i = 1 ; i <= n ; ++ i)
    for (j = 1 ; j <= m ; ++ j)
      if (g[i][j] == '.')
        id[i][j] = ++ pcnt;
  for (i = 1 ; i <= n ; ++ i)
    for (j = 1 ; j <= m ; ++ j)
      if (g[i][j] == 'D')
      {
        ++ dcnt;
        BFS(dcnt , i , j);
      }

  if (!pcnt) {puts("0");return;}
  for (ans = 1 ; ans <= n * m ; ++ ans)
  {
    memset(re , 0 , sizeof(re)) , sum = 0;
    for (i = 1 ; i <= pcnt ; ++ i)
    {
      memset(f , 0 , sizeof(f));
      sum += find(i);
      if (sum == pcnt)
      {
        printf("%d\n" , ans);
        return;
      }
    }
  }
 // hehe:
  puts("impossible");
}

int main()
{
  int _; scanf("%d",&_); while (_--)
  //while (scanf("%d%d%d%d%d",&p , &n , &m , &a , &b) , p)
    work();
  return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值