HDU 3036 Escape 网格图多人逃生 网络流||二分匹配 建图技巧
题意：
每一个' . '有一个姑娘， E是出口，'.'是空地 , 'X‘ 是墙。
每秒钟每一个姑娘能够走一步（上下左右）
每秒钟每一个出口仅仅能出去一个人
给定n*m的地图， 时限T
问全部姑娘是否能在T秒内逃生，若能输出最小值，不能输出"impossible"
思路：
显然是二分答案+网络流判可行。
由于每一个出口每秒钟仅仅能出去一个人，那么就把每一个出口按时间拆点，则T秒钟就拆成T个点。
网络流建图
1、源点 到 每一个姑娘 建流量为1的边。
2、若某姑娘到 a出口须要时间为 t秒，则建一条流量为1的边 连向a出口拆点为t秒的点。
3、每一个出口的全部拆点向汇点连一条流量为1的边。
4、对于每一个出口u的x秒拆点，向u的x+1秒的拆点连一条流量为inf的边（表示从x秒来的人假设x秒还从u出不去，能够在u等到x+1秒出去）
二分一下 第二点中的 t 秒(即答案），推断最大流是否等于人数
二分匹配建图：
枚举时间（由于时间最大仅仅有12*12）
在原图的基础上加上以下的边：
在i时间内若能出去则那姑娘向全部能出去的 出口时间拆点连边。
再在原来的匹配上继续增广，累计最大匹配数。
当最大匹配数==人数时则是最小的时间。
数据较水能够用点非主流的建图卡过去，预计仅仅有30组数据。
网络流代码：
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
using namespace std; #define ll int #define N 20050
#define M 105000
#define inf 10737418
struct Edge{
ll from, to, cap, nex;
}edge[M*4];//注意这个一定要够大 不然会re 还有反向弧 ll head[N], edgenum;
void add(ll u, ll v, ll cap){
Edge E = { u, v, cap, head[u]};
edge[ edgenum ] = E;
head[u] = edgenum ++; Edge E2= { v, u, 0, head[v]};
edge[ edgenum ] = E2;
head[v] = edgenum ++;
}
ll sign[N];
bool BFS(ll from, ll to){
memset(sign, -1, sizeof(sign));
sign[from] = 0; queue<ll>q;
q.push(from);
while( !q.empty() ){
int u = q.front(); q.pop();
for(ll i = head[u]; i!=-1; i = edge[i].nex)
{
ll v = edge[i].to;
if(sign[v]==-1 && edge[i].cap)
{
sign[v] = sign[u] + 1, q.push(v);
if(sign[to] != -1)return true;
}
}
}
return false;
}
ll Stack[N], top, cur[N];
ll dinic(ll from, ll to){
ll ans = 0;
while( BFS(from, to) )
{
memcpy(cur, head, sizeof(head));
ll u = from; top = 0;
while(1)
{
if(u == to)
{
ll flow = inf, loc;//loc 表示 Stack 中 cap 最小的边
for(ll i = 0; i < top; i++)
if(flow > edge[ Stack[i] ].cap)
{
flow = edge[Stack[i]].cap;
loc = i;
} for(ll i = 0; i < top; i++)
{
edge[ Stack[i] ].cap -= flow;
edge[Stack[i]^1].cap += flow;
}
ans += flow;
top = loc;
u = edge[Stack[top]].from;
}
for(ll i = cur[u]; i!=-1; cur[u] = i = edge[i].nex)//cur[u] 表示u所在能增广的边的下标
if(edge[i].cap && (sign[u] + 1 == sign[ edge[i].to ]))break;
if(cur[u] != -1)
{
Stack[top++] = cur[u];
u = edge[ cur[u] ].to;
}
else
{
if( top == 0 )break;
sign[u] = -1;
u = edge[ Stack[--top] ].from;
}
}
}
return ans;
} void init(){memset(head,-1,sizeof head);edgenum = 0;} char mp[15][15];
int n, m, T;
int Hash(int x,int y){return x*m+y;} vector<int>E,P;
int dis[150][150], step[4][2]={1,0,-1,0,0,1,0,-1};
bool vis[150][150];
void bfs(int sx,int sy){
int start = Hash(sx,sy);
memset(vis, 0, sizeof vis);
vis[sx][sy] = 1;
dis[start][start] = 0;
queue<int>qx,qy; while(!qx.empty())qx.pop(); while(!qy.empty())qy.pop();
qx.push(sx), qy.push(sy);
while(!qx.empty()){
int x = qx.front(), y = qy.front();
qx.pop(); qy.pop();
for(int i = 0; i < 4; i++){
int dx = x + step[i][0], dy = y + step[i][1];
if(!(0<=dx&&dx<n&&0<=dy&&dy<m))continue;
if(vis[dx][dy] || mp[dx][dy]!='.')continue;
vis[dx][dy] = 1;
dis[Hash(dx,dy)][start] = dis[start][Hash(dx,dy)] = dis[start][Hash(x,y)]+1;
qx.push(dx); qy.push(dy);
}
}
}
bool ok(int TIME){
init();
int from = N-2, to = N-1;
for(int i = 0; i < P.size(); i++)add(from, P[i], 1);
for(int i = 0; i < P.size(); i++)
{
for(int j = 0; j < E.size(); j++)
if(dis[P[i]][E[j]]<=TIME) add(P[i],j*150+150+dis[P[i]][E[j]],1);
}
for(int i = 0; i < E.size(); i++)
for(int j = 1; j <= TIME; j++)
{
add(i*150+150+j, to, 1);
if(j!=TIME)
add(i*150+150+j,i*150+150+j+1,inf);
}
return dinic(from,to)==P.size();
}
int main(){
//freopen("date.in","r+",stdin);
//freopen("ans.out","w+",stdout); int i, j;
while(~scanf("%d %d %d",&n,&m,&T)){
E.clear(); P.clear();
memset(dis, 0, sizeof dis);
memset(mp, 0, sizeof mp);
for(i=0;i<n;i++)scanf("%s",mp[i]); for(i=0;i<n;i++)for(j=0;j<m;j++)
{
if(mp[i][j]=='E')E.push_back(Hash(i,j));
else if(mp[i][j]=='.')P.push_back(Hash(i,j));
}
if(P.size()==0){puts("0");continue;} if(E.size()==0){puts("impossible");continue;} for(i = 0; i < E.size(); i++)bfs(E[i]/m, E[i]%m);
int l = 1, r = 256, ans = inf;
while(l<=r)
{
int mid = (l+r)>>1;
if(ok(mid))ans = min(ans, mid), r = mid-1;
else l = mid+1;
}
if(T<ans || ans==inf){puts("impossible");continue;}
else printf("%d\n",ans);
}
return 0;
}
/*
6 12 100000
......E....E
E.EEE...E...
...E....EE..
.E.E.......E
.....E......
.E...E...EE.
11 3 100000
X..
...
.E.
EE.
E.E
...
...
.E.
E..
..E
E..
7 6 100000
.E..E.
...E..
......
...E.E
...X..
...EE.
....E.
11 10 100000
.E........
......E...
........E.
EE....E..E
...E...E..
......XE.E
..........
.........X
.........E
.EE.......
..EE.E.E.E
7 3 100000
...
..E
...
..X
EE.
...
.X.
3 2 100000
..
..
EE
6 7 100000
..E....
.......
.....E.
...EE..
.......
.......
6 8 100000
..EEEE..
...E....
........
........
.......E
E......E 10 8 100000
E.E..X..
........
E...E...
E..E....
........
.....E..
......EX
........
.E......
..E..EE. 5 5 100
.....
XXXXX
EEEEE
.....
XXXXX 1 1 1
E 1 1 1
P
1 1 1
X 2 2 1
E.
.X 2 2 2
E.
.X 8 1 1000
.
E
.
.
.
.
.
. 3 4 100
E...
....
...E 3 4 100
E...
.E..
...E 4 5 100
E....
.E...
...E.
..... 4 5 100
.....
.E...
.....
..... 12 12 10000000
...E.....E.E
E.E.E.....E.
.......E....
............
......E.....
............
...E.....E.E
E.E.E.....E.
.......E....
............
......E.....
............ 12 12 6
...E.....E.E
E.E.E.....E.
.......E....
............
......E.....
............
...E.....E.E
E.E.E.....E.
.......E....
............
......E.....
............ 12 12 10000000
E...........
............
............
............
............
............
............
............
............
............
............
............ 12 12 10000000
E...........
............
............
............
............
.....E......
............
............
............
............
............
...........E 3 3 10
...
.E.
... */
二分匹配代码：
#include <set>
#include <map>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std; /* '0. Macros, Preprocessers, */
#pragma comment(linker, "/STACK:102400000,102400000")
#ifdef _WIN32
typedef __int64 int64;
#define OD(_TYPE) printf("%I64d\n", _TYPE)
#define ODC(_CASE, _TYPE) printf("Case %d: %I64d\n", _CASE++, _TYPE);
#else
typedef long long int64;
#define OD(_TYPE) printf("%lld\n", _TYPE)
#define ODC(_CASE, _TYPE) printf("Case %d: %lld\n", _CASE++, _TYPE);
#endif
#define LSON l,mid,id<<1
#define RSON mid+1,r,id<<1|1
#define MM (l+r)>>1 /* '1. Constants, */
const double EPS = 1e-8;
const double PI = acos(-1.0);
const int INF = 0x3f3f3f3f; /* '2. Coding Area, */
const int N = 15;
const int M = 205;
const int dir[4][2] = { { -1, 0 }, { 0, 1 }, { 1, 0 }, { 0, -1 } };
struct STATE {
int x, y, t;
STATE() { }
STATE(int x, int y, int t) : x(x), y(y), t(t) { }
};
int uN, vN;
vector<int> g[M];
int linker[M * M]; // 每一个v相应的u匹配
bool used[M * M];
char mp[N][N];
int dis[M][M]; // 每一个女孩到出口的最短时间
bool vis[N][N]; // 訪问过或没有
int r, c, T;
int gid, eid;
map<pair<int, int>, int> mpg, mpe; // 每一个点的编号 bool dfs(int u) {
for (unsigned int i = 0; i < g[u].size(); i++) {
int v = g[u][i];
if (used[v]) continue;
used[v] = true;
if (linker[v] == -1 || dfs(linker[v])) {
linker[v] = u;
return true;
}
}
return false;
} int hungary() {
int res = 0;
memset(linker, -1, sizeof(linker));
for (int u = 0; u < uN; u++) {
memset(used, 0, sizeof(used));
if (dfs(u))
++res;
}
return res;
} void bfs(int x, int y) {
queue<STATE> q;
int eid = mpe[make_pair(x, y)];
memset(vis, 0, sizeof(vis));
vis[x][y] = 1;
q.push(STATE(x, y, 0));
while (!q.empty()) {
STATE now = q.front();
q.pop();
for (int i = 0; i < 4; i++) {
int dx = now.x + dir[i][0];
int dy = now.y + dir[i][1];
if (dx >= 0 && dx < r && dy >= 0 && dy < c) {
if (mp[dx][dy] != '.') continue;
if (vis[dx][dy]) continue;
int gid = mpg[make_pair(dx, dy)];
q.push(STATE(dx, dy, now.t + 1));
dis[gid][eid] = now.t + 1;
vis[dx][dy] = 1;
}
}
}
} bool build_and_run(int limit) {
for (int i = 0; i < M; i++) g[i].clear();
uN = gid, vN = eid * limit;
for (int i = 0; i < gid; i++) {
for (int j = 0; j < eid; j++) {
for (int k = dis[i][j]; k <= limit; k++) {
g[i].push_back(j * limit + k);
}
}
}
int ans = hungary();
if (ans >= gid) return true;
return false;
} void gao() {
gid = eid = 0;
mpg.clear(); mpe.clear();
memset(dis, 0x3f, sizeof(dis));
for (int i = 0; i < r; i++) {
scanf("%s", mp[i]);
for (int j = 0; j < c; j++) {
if (mp[i][j] == 'E') mpe[make_pair(i, j)] = eid++;
if (mp[i][j] == '.') mpg[make_pair(i, j)] = gid++;
}
} for (int i = 0; i < r; i++) {
for (int j = 0; j < c; j++) {
if (mp[i][j] == 'E') {
bfs(i, j);
}
}
} int L = 1, R = 256, ans = INF;
while (L <= R) {
int mid = (L + R) >> 1;
if (build_and_run(mid)) {
ans = mid;
R = mid - 1;
} else {
L = mid + 1;
}
}
if (ans == INF || ans > T) {
printf("impossible\n");
} else {
printf("%d\n", ans);
}
} int main() {
while (~scanf("%d%d%d", &r, &c, &T)) {
memset(g, 0, sizeof(g));
gao();
}
return 0;
}
HDU 3036 Escape 网格图多人逃生 网络流||二分匹配 建图技巧的更多相关文章
- TTTTTTTTTTTTTTTTT POJ 2226 草地覆木板 二分匹配 建图
Muddy Fields Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 9754 Accepted: 3618 Desc ...
- TTTTTTTTTTTTTTTTTT POJ 2724 奶酪消毒机 二分匹配 建图 比较难想
Purifying Machine Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 5004 Accepted: 1444 ...
- hdu 4857 逃生 拓扑排序+逆向建图
逃生 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Problem Descr ...
- Food HDU - 4292 网络流 拆点建图
http://acm.hdu.edu.cn/showproblem.php?pid=4292 给一些人想要的食物和饮料,和你拥有的数量,问最多多少人可以同时获得一份食物和一份饮料 写的时候一共用了2种 ...
- hdu 4725 The Shortest Path in Nya Graph （最短路+建图）
The Shortest Path in Nya Graph Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K ...
- 图论--网络流--最小割 HDU 2485 Destroying the bus stations（最短路+限流建图）
Problem Description Gabiluso is one of the greatest spies in his country. Now he's trying to complet ...
- poj 1149 Pigs 网络流-最大流 建图的题目（明天更新）-已更新
题目大意:是有M个猪圈,N个顾客,顾客要买猪,神奇的是顾客有一些猪圈的钥匙而主人MIRKO却没有钥匙,多么神奇?顾客可以在打开的猪圈购买任意数量的猪,只要猪圈里有足够数量的猪.而且当顾客打开猪圈后mi ...
- POJ 1274 The Perfect Stall、HDU 2063 过山车（最大流做二分匹配）
The Perfect Stall Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 24081 Accepted: 106 ...
- poj 3281 Dining 网络流-最大流-建图的题
题意很简单:JOHN是一个农场主养了一些奶牛,神奇的是这些个奶牛有不同的品味,只喜欢吃某些食物,喝某些饮料,傻傻的John做了很多食物和饮料,但她不知道可以最多喂饱多少牛,(喂饱当然是有吃有喝才会饱) ...
随机推荐
- memory引擎的索引失效一例
memory引擎的索引失效一例 memory引擎的索引失效一例
- [置顶] 简单解析linux下进程通信方法
linux下的进程通信手段基本上是从Unix平台上的进程通信手段继承而来的.而对Unix发展做出重大贡献的两大主力AT&T的贝尔实验室及BSD(加州大学伯克利分校的伯克利软件发布中心)在进程间 ...
- (step6.3.2)hdu 1068(Girls and Boys——二分图的最大独立集)
题目大意:第一行输入一个整数n,表示有n个节点.在接下来的n行中,每行的输入数据的格式是: 1: (2) 4 6 :表示编号为1的人认识2个人,他们分别是4.6: 求,最多能找到多少个人,他们互不认识 ...
- 居然还有FindFirstChangeNotification函数
http://download.csdn.net/download/sololie/5966243
- PHP+lghttpd + postgresql 环境搭建
PHP+lghttpd + postgresql 环境搭建 Linux 下PHP环境搭建 安装环境:PHP+lghttpd + postgresql 1:yum install lighttpd 红旗 ...
- oracle 之 内存—鞭辟近里（三）
oracle 之 内存—鞭辟近里(三) 今天是2013-07-08,今天晚上突然接到一个电话,我的外甥问我的qq是多少,我感觉很吃惊,他长大了.在他现在这个年龄就开始接触网络,我难免有少许担心,希望他 ...
- FZU2176（二维线段树+dfs）
传送门:easy problem 题意:给定一棵n个节点以1为根的树,初始每个节点的值为0,现在我们要在树上进行一些操作,操作有两种类型. 1 x val 表示对以x为根的子树的每个点进行加权操作(我 ...
- Wix打包系列（一）如何使用wix制作安装程序
原文:Wix打包系列(一)如何使用wix制作安装程序 最近由于项目需要,需要给客户制作安装程序,一开始使用vs2005自带的打包工程来打包,但用了一段时间发现vs打包太死板,而且使用起来问题很多.收费 ...
- poj3280（区间dp）
题目连接:http://poj.org/problem?id=3280 题意:给定一个长度为m(m<=2000)的小写字母字符串,在给定组成该字符串的n(n<=26)个字符的添加和删除费用 ...
- 组队赛第二场：字符串哈希+DP
长春现场赛 HDU 4821 思路:这题周赛的时候没做出来,有点可惜了.要是当时记起来unsigned long long自己主动取模,然后提醒一下大帝的话,后续大帝就能过了. 唉,导致让他取了好多 ...