P3191 [HNOI2007]紧急疏散EVACUATE(费用流)
费用流+卡常优化
我们只关心一个人通过门时的时间,在空地的行走时间可以分层维护
于是根据时间分层,到门的时候再计算代价,即代价$=$层数
每经过$1$单位时间就向下走一层
然后就是优化:
1. 删去多余点(层):只要开和人数相同的层,因为一个人的等待时间一定小等于人数
2. 能预处理的尽量预处理
3. register,快读,inline等
4. 能不用stl尽量不用(使用stl:queue=TLE),建议手写队列
复杂度$O($卡常能过$)$,详情可看code
#include<iostream>
#include<cstring>
#include<cstdio>
#define re register
using namespace std;
#define N 200005
#define M 3000005
const int d1[]={,,,-};
const int d2[]={,,-,};
char a[][];
int n,m,nm,mx,ans,tot,S,T,v[N],p[N],id[][];
bool vis[N]; int h[M],L,R;
int Cnt=,hd[N],nxt[M],ed[N],poi[M],con[M],val[M];
void adde(int x,int y,int v1,int v2){
nxt[ed[x]]=++Cnt; hd[x]=hd[x]?hd[x]:Cnt;
ed[x]=Cnt; poi[Cnt]=y; con[Cnt]=v1; val[Cnt]=v2;
}
inline void link(int x,int y,int v1,int v2){adde(x,y,v1,v2),adde(y,x,,-v2);}
#define to poi[i]
bool bfs(){//普通费用流
memset(v,,sizeof(v)); int inf=v[];
h[L=]=S; R=; vis[S]=; v[S]=;
while(L!=R){
re int x=h[L++]; vis[x]=;
if(L>=M) L=;
for(re int i=hd[x];i;i=nxt[i])
if(con[i]&&v[to]>v[x]+val[i]){
v[to]=v[x]+val[i]; p[to]=i;
if(!vis[to]){
vis[to]=,h[R++]=to;
if(R>=M) R=;
}
}
}if(v[T]==inf) return ;
tot-=; ans=max(ans,v[T]);
for(re int u=T;u!=S;u=poi[p[u]^])
con[p[u]]-=,con[p[u]^]+=;
return ;
}//因为每次流量均为1,可以省去流量数组
void draw(int x,int y){
re int p=id[x][y];
if(a[x][y]=='X') return ;
if(a[x][y]=='D') for(re int i=;i<=mx;++i) link(p+i*nm,T,,i);//分层与终点连边,代价=层数
if(a[x][y]=='.'){
link(S,p,,); ++tot;
for(re int i=;i<mx;++i){
link(p+i*nm,p+(i+)*nm,1e9,);
for(re int j=;j<;++j){
int rx=x+d1[j],ry=y+d2[j];
if(a[rx][ry]=='X') continue;
link(p+i*nm,id[rx][ry]+(i+)*nm,1e9,);//与四周连边
}
}
}
}
int main(){
scanf("%d%d",&n,&m); S=N-; T=N-; nm=n*m;
for(re int i=;i<=n;++i) scanf("%s",a[i]+);
for(re int i=;i<=n;++i)
for(re int j=;j<=m;++j)
mx+=(a[i][j]=='.'),id[i][j]=(i-)*m+j;//mx:人数,建图需要的层数
for(re int i=;i<=n;++i) for(re int j=;j<=m;++j) draw(i,j);
while(bfs()) ;
if(tot) puts("impossible");
else printf("%d",ans);
return ;
}
P3191 [HNOI2007]紧急疏散EVACUATE(费用流)的更多相关文章
- P3191 [HNOI2007]紧急疏散EVACUATE
传送门 这一题很容易想到网络流 一开始傻逼地模拟整个图每一个时间的情况,显然会爆炸 发现我们只要考虑起点到门之间的距离,不用每一步只走一格 所以直接 $BFS$ 预处理距离然后二分答案,网络流判断即可 ...
- 洛谷 P3191 [HNOI2007]紧急疏散EVACUATE(网络最大流)
题解 二分答案+Dinic最大流 二分答案\(mid\) 把门拆成\(mid\)个时间点的门 相邻时间的门连一条\(inf\)的边 预处理出每个门到每个人的最短时间 为\(dis[k][i][j]\) ...
- BZOJ1189:[HNOI2007]紧急疏散EVACUATE(最大流,枚举)
Description 发生了火警,所有人员需要紧急疏散!假设每个房间是一个N M的矩形区域.每个格子如果是'.',那么表示这是一 块空地:如果是'X',那么表示这是一面墙,如果是'D',那么表示这是 ...
- 题解 P3191 [HNOI2007]紧急疏散EVACUATE
题解 本篇题解做法为BFS+二分+最大流 二分需要撤离的时间,也就是答案(这算是一个比较套路的了) 重点在于建模(设时间为 \(tim\)): 我们将每个门拆点,拆成 \(tim\) 个,每个点向汇点 ...
- luogu P3191 [HNOI2007]紧急疏散EVACUATE
传送门 qwq这题好大力 首先可以预处理出每个人到每个门前面那个格子的最早时间,然后答案如果比最小答案大的话也是合法的,所以可以二分最终答案.检查\(mid\)是否合法就考虑每个人要去哪个门才会合法, ...
- BZOJ1189: [HNOI2007]紧急疏散evacuate 二分+最大流
1189: [HNOI2007]紧急疏散evacuate Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1132 Solved: 412[Submi ...
- 【二分答案】【最大流】[HNOI2007]紧急疏散EVACUATE
[HNOI2007]紧急疏散EVACUATE 题目描述 发生了火警,所有人员需要紧急疏散!假设每个房间是一个N M的矩形区域.每个格子如果是'.',那么表示这是一块空地:如果是'X',那么表示这是一面 ...
- Bzoj1189 [HNOI2007]紧急疏散evacuate
1189: [HNOI2007]紧急疏散evacuate Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2293 Solved: 715 Descr ...
- [HNOI2007]紧急疏散EVACUATE (湖南2007年省选)
[HNOI2007]紧急疏散EVACUATE 题目描述 发生了火警,所有人员需要紧急疏散!假设每个房间是一个N M的矩形区域.每个格子如果是'.',那么表示这是一块空地:如果是'X',那么表示这是一面 ...
随机推荐
- 零点.Net Core 接触
一.Program.cs类与Startup类 1.一切从Main开始,Main方法包含了是整个应用程序的入口 ASP.NET Core应用程序可以配置和启动主机(Host). 主机负责应用程序启动和生 ...
- 一、VS支持Vue语法
一.VS支持Vue语法
- LinuxC语言实现ATM取款机实验Socket
链接:https://pan.baidu.com/s/1sZt4qhYc6CDJVpoJHbtw-Q 提取码:53ot 复制这段内容后打开百度网盘手机App,操作更方便哦 本实验用的是Centos7t ...
- [转]Tomcat9.0安装教程 Tomcat9.0环境变量配置教程
[转]Tomcat9.0安装教程 Tomcat9.0环境变量配置教程 [转]超详细MySQL安装及基本使用教程
- Linux 查找指定内容在哪个文件中
在实际的工作中,忘记配置项放在哪个文件中时,可借助命令来查询. eg: 1.grep -r "查询内容" 文件目录 #这样查询出来的包括文件名+内容 grep -r -l ...
- mvn 打包排除test
mvn clean package compile -Dmaven.test.skip=true
- join优化
1.left outer join先执行连接操作,再将结果通过WHERE语句进行过滤 select s.ymd,s.symbol,s.price_close,d.dividend from stock ...
- Linux下统计当前文件夹下的文件个数
Linux下统计当前文件夹下的文件个数 √ ls -l |grep "^-"|wc -l
- 【HDOJ6665】Calabash and Landlord(dfs)
题意:二维平面上有两个框,问平面被分成了几个部分 x,y<=1e9 思路:分类讨论可以 但数据范围实在太小了,离散化以后随便dfs一下 #include<bits/stdc++.h> ...
- HashMap 重新学习
HashMap 重新学习 先使用 HashCode() 方法,该方法决定位置. 然后使用 equals() 方法,决定在相同位置的时候,是否覆盖. 当程序试图将一个键值对放入 HashMap 的时候, ...