题解

二分答案+Dinic最大流

二分答案\(mid\)

把门拆成\(mid\)个时间点的门

相邻时间的门连一条\(inf\)的边

预处理出每个门到每个人的最短时间

为\(dis[k][i][j]\) 在\((i,j)\)的人到第\(k\)个门最短时间

然后一个人连向每个第\(dis[k][i][j]\)那个时刻的门,容量为\(1\)

然后,源点连向每个人一条容量为\(1\)的边

所有门都连向汇点一条容量为\(1\)的边(其实只要每个最后一个时刻的门连一条容量为\(mid\)的边即可)

Code

#include<bits/stdc++.h>

#define LL long long
#define RG register using namespace std;
template<class T> inline void read(T &x) {
x = 0; RG char c = getchar(); bool f = 0;
while (c != '-' && (c < '0' || c > '9')) c = getchar(); if (c == '-') c = getchar(), f = 1;
while (c >= '0' && c <= '9') x = x*10+c-48, c = getchar();
x = f ? -x : x;
return ;
}
template<class T> inline void write(T x) {
if (!x) {putchar(48);return ;}
if (x < 0) x = -x, putchar('-');
int len = -1, z[20]; while (x > 0) z[++len] = x%10, x /= 10;
for (RG int i = len; i >= 0; i--) putchar(z[i]+48);return ;
} const int N = 40, inf = 1e9; int n, m;
int fx[4] = {1, 0, -1, 0};
int fy[4] = {0, 1, 0, -1}; char ch[N][N]; int dis[N<<2][N][N], door, people; struct BB {
int x, y;
};
queue<BB> Q; void BFS(int k, int xx, int yy) {
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++)
dis[k][i][j] = inf;
dis[k][xx][yy] = 0; Q.push((BB) {xx, yy});
while (!Q.empty()) {
int x = Q.front().x, y = Q.front().y; Q.pop();
for (int w = 0; w < 4; w++) {
int i = x + fx[w], j = y + fy[w];
if (i < 1 || i > n || j < 1 || j > m || ch[i][j] != '.') continue;
if (dis[k][i][j] > dis[k][x][y]+1)
dis[k][i][j] = dis[k][x][y]+1, Q.push((BB) {i, j});
}
}
} struct node {
int to, nxt, w;
}g[N*N*N*N];
int last[N*N*N], cur[N*N*N], dep[N*N*N], s, t, gl = 1, p[N][N]; inline void add(int a, int b, int c) {
g[++gl] = (node) {b, last[a], c};
last[a] = gl;
g[++gl] = (node) {a, last[b], 0};
last[b] = gl;
}
queue<int> q;
bool bfs() {
memset(dep, 0, sizeof(dep));
q.push(s); dep[s] = 1;
while (!q.empty()) {
int u = q.front(); q.pop();
for (int i = last[u]; i; i = g[i].nxt) {
int v = g[i].to;
if (!dep[v] && g[i].w) {
dep[v] = dep[u]+1;
q.push(v);
}
}
}
return dep[t] == 0 ? 0 : 1;
} int dfs(int u, int d) {
if (u == t) return d;
for (int &i = cur[u]; i; i = g[i].nxt) {
int v = g[i].to;
if (g[i].w && dep[v] == dep[u]+1) {
int di = dfs(v, min(d, g[i].w));
if (di) {
g[i].w -= di;
g[i^1].w += di;
return di;
}
}
}
return 0;
} int Dinic() {
int ans = 0;
while (bfs()) {
for (int i = 1; i <= t; i++)
cur[i] = last[i];
while (int d = dfs(s, inf)) ans += d;
}
return ans;
} int check(int mid) {
s = door*mid+people+1, t = s+1;
memset(last, 0, sizeof(last)); gl = 1;
for (int k = 1; k <= door; k++)
for (int z = 1; z <= mid; z++) {
if (z < mid)
add((k-1)*mid+z, (k-1)*mid+z+1, inf);
else add(k*mid, t, mid);
}
for (int i = 1; i <= people; i++)
add(s, door*mid+i, 1);
for (int k = 1; k <= door; k++)
for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++)
if (ch[i][j] == '.' && dis[k][i][j] <= mid)
add(door*mid+p[i][j], (k-1)*mid+dis[k][i][j], 1);
return Dinic();
} int main() {
read(n), read(m);
for (int i = 1; i <= n; i++)
scanf("%s", ch[i]+1); for (int i = 1; i <= n; i++)
for (int j = 1; j <= m; j++)
if (ch[i][j] == 'D')
BFS(++door, i, j);
else if (ch[i][j] == '.') p[i][j] = ++people;
int l = 0, r = people, ans = 666666;
while (l <= r) {
int mid = (l + r) >> 1;
if (check(mid) == people) r = mid-1, ans = mid;
else l = mid+1;
}
if (ans == 666666) puts("impossible");
else printf("%d\n", ans);
return 0;
}

洛谷 P3191 [HNOI2007]紧急疏散EVACUATE(网络最大流)的更多相关文章

  1. P3191 [HNOI2007]紧急疏散EVACUATE(费用流)

    P3191 [HNOI2007]紧急疏散EVACUATE 费用流+卡常优化 我们只关心一个人通过门时的时间,在空地的行走时间可以分层维护 于是根据时间分层,到门的时候再计算代价,即代价$=$层数 每经 ...

  2. BZOJ1189: [HNOI2007]紧急疏散evacuate 二分+最大流

    1189: [HNOI2007]紧急疏散evacuate Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1132  Solved: 412[Submi ...

  3. 洛谷P3376【模板】网络最大流 ISAP

    这篇博客写得非常好呀. 传送门 于是我是DCOI这一届第一个网络流写ISAP的人了,之后不用再被YKK她们嘲笑我用Dinic了!就是这样! 感觉ISAP是会比Dinic快,只分一次层,然后不能增广了再 ...

  4. 洛谷P3376 【模板】网络最大流

    题目描述 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入输出格式 输入格式: 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行 ...

  5. 洛谷 P3376 【模板】网络最大流

    题目描述 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入输出格式 输入格式: 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行 ...

  6. 洛谷——P3376 【模板】网络最大流

    题目描述 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入输出格式 输入格式: 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行 ...

  7. 洛谷 P3376【模板】网络最大流

    题目描述 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. 输入输出格式 输入格式: 第一行包含四个正整数N.M.S.T,分别表示点的个数.有向边的个数.源点序号.汇点序号. 接下来M行每行 ...

  8. 『题解』洛谷P3376 【模板】网络最大流

    Problem Portal Portal1:Luogu Description 如题,给出一个网络图,以及其源点和汇点,求出其网络最大流. Input 第一行包含四个正整数\(N,M,S,T\),分 ...

  9. 洛谷 P3376 【模板】网络最大流 题解

    今天学了网络最大流,EK 和 Dinic 主要就是运用搜索求增广路,Dinic 相当于 EK 的优化,先用bfs求每个点的层数,再用dfs寻找并更新那条路径上的值. EK 算法 #include< ...

随机推荐

  1. 由于挂载的nfs存储目录掉下线,导致创建VM时,无法创建

    具体错误,如下截图 重新挂载存储后,在创建VM,将成功

  2. SqlServer——神奇代码1之Update

    说明:一个带有update的循环的代码.很简单,但是在QQ群里问了,应该说是很少有人注意这个问题,也就是很少有人真的理解SQL中的Update. 代码如下: if object_id('tempdb. ...

  3. Hyperledger Fabric Chaincode解析

    首先看下Blockchain结构,除了header指向下一个block的hash value外,block是由一组transaction构成, Transactions --> Blocks - ...

  4. Ubuntu 14.04 安装配置强大的星际译王(stardict)词典

    转载http://blog.csdn.net/huyisu/article/details/53437931

  5. logback.xml 实例

    <?xml version="1.0" encoding="UTF-8"?> <configuration debug="false ...

  6. HTML5之:link与title的区别

    [link]标签:外联导入样式 例1:<link rel="stylesheet" type="text/css" href="theme.cs ...

  7. (转)C# HTML解析示例---星星引发的血案

    原文地址:http://www.cnblogs.com/wurang/archive/2013/06/14/3119023.html [前言] 从CSDN转投cnBlog也有一段时间了,发现cnBlo ...

  8. html manifest 离线配置

    HTML5 引入了应用程序缓存,这意味着 web 应用可进行缓存,并可在没有因特网连接时进行访问. 应用程序缓存为应用带来三个优势: 离线浏览 - 用户可在应用离线时使用它们 速度 - 已缓存资源加载 ...

  9. 使用 Windows Phone 8 文件和 URI 关联的自动启动应用

    更详细,猛撸这里:http://msdn.microsoft.com/zh-cn/library/windowsphone/develop/jj206987(v=vs.105).aspx 在WMApp ...

  10. SpringMVC+Hibernate 项目开发之三 (创建SpringMVC项目)

    引用(很全面了):http://blog.csdn.net/dhx20022889/article/details/38041039 我只想说默认创建的项目使用的Spring版本可能不是你想要的,可以 ...