题目链接:紧急疏散
这薄脊题我代码不知不觉就写长了……
这道题二分答案显然,然后用最大流\(check\)即可。设当前二分的答案为\(x\),那么把每扇门拆成\(x\)个点,第\(i\)个代表在第\(i\)个时刻从这个门走出去。然后把每个空地往可以到达的们的相应时间连边就可以了。判一下这张图是否满流即可。
然后我们就需要先求出每个空地到门的距离……注意途中不能经过另外的门,否则会被BZOJ上加强的数据给卡掉……
下面贴代码:
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout) #define rep(i,l,r) for(int i=l;i<=r;i++) #define INF 1000000000 #define maxn 200010 #define maxm 1000010 #define N 24 using namespace std; typedef long long llg; int zx[4]={1,-1,0,0},zy[4]={0,0,-1,1}; int n,m,nu[N][N],f[N*N][N*N],dor,tol,dep[maxn],d[maxn]; int head[maxn],next[maxm],to[maxm],c[maxm],tt,S,T; char s[N][N]; bool w[N*N]; void link(int x,int y,int z){ to[++tt]=y;next[tt]=head[x];head[x]=tt; to[++tt]=x;next[tt]=head[y];head[y]=tt; c[tt-1]=z; c[tt]=0; } bool bfs(){ for(int i=1;i<=T;i++) dep[i]=-1; int ld=0,rd=0; dep[d[rd++]=S]=1; while(ld!=rd){ int u=d[ld++]; for(int i=head[u],v;v=to[i],i;i=next[i]) if(c[i] && dep[v]==-1) dep[v]=dep[u]+1,d[rd++]=v; } return dep[T]!=-1; } int dfs(int u,int now){ if(!now) return 0; if(u==T) return now; int low=0,res; for(int i=head[u],v;v=to[i],i;i=next[i]) if(c[i] && dep[v]==dep[u]+1){ res=dfs(v,min(now,c[i])); low+=res; c[i]-=res; c[i^1]+=res; now-=res; } if(!low) dep[u]=-1; return low; } bool check(int x){ tt=1; S=tol*x+1; T=S+1; for(int i=1;i<=tol;i++) if(w[i]) for(int j=1;j<=x;j++){ link(tol*(j-1)+i,T,1); if(j!=x) link(tol*(j-1)+i,tol*j+i,INF); } else{ link(S,i,1); for(int j=1;j<=tol;j++) if(w[j] && f[i][j]<=x) link(i,tol*(f[i][j]-1)+j,1); } int now=0; while(bfs()) now+=dfs(S,INF); for(int i=1;i<=T;i++) head[i]=0; return now==tol-dor; } int dx[N*N],dy[N*N]; bool vis[N][N]; void getdis(int x,int y){ int u=nu[x][y]; if(w[u]) return; rep(i,1,n) rep(j,1,m) vis[i][j]=0; int ld=0,rd=0; dx[rd]=x,dy[rd++]=y; vis[x][y]=1; while(ld!=rd){ int o=dx[ld],p=dy[ld++]; for(int k=0,i,j,v;k<4;k++){ i=o+zx[k],j=p+zy[k]; v=nu[i][j]; if(i>=1 && i<=n && j>=1 && j<=m){ if(vis[i][j]) continue; if(s[i][j]!=‘X‘) f[u][v]=f[u][nu[o][p]]+1; if(s[i][j]==‘.‘) dx[rd]=i,dy[rd++]=j,vis[i][j]=1; } } } } int main(){ File("a"); scanf("%d %d",&n,&m); for(int i=1;i<=n;i++){ scanf("%s",s[i]+1); for(int j=1;j<=m;j++) if(s[i][j]!=‘X‘){ nu[i][j]=++tol; dor+=(w[tol]=(s[i][j]==‘D‘)); } } rep(i,1,tol) rep(j,1,tol) if(i!=j) f[i][j]=INF; rep(i,1,n) rep(j,1,m) if(nu[i][j]) getdis(i,j); int l=0,r=tol-dor+1,mid; while(l!=r){ mid=(l+r)>>1; if(check(mid)) r=mid; else l=mid+1; } if(l>tol-dor) printf("impossible"); else printf("%d",l); return 0; }
时间: 2024-10-23 15:56:14