題目鏈接：點(diǎn)擊查看

題目大意：給出一個(gè)n*n的實(shí)驗(yàn)室，其中有一些安全的核反應(yīng)堆，我們視為墻，還有一個(gè)出了故障的核反應(yīng)堆，會(huì)在t秒之后爆炸，出故障的反應(yīng)堆會(huì)擴(kuò)散毒氣，毒氣在每個(gè)時(shí)間單位都會(huì)向相鄰的方格擴(kuò)散一個(gè)單位，現(xiàn)在有數(shù)個(gè)科學(xué)家以及數(shù)個(gè)安全膠囊，如果科學(xué)家可以在毒氣擴(kuò)散以及核反應(yīng)堆爆炸之前到達(dá)安全氣囊，則可以獲救，其余情況都會(huì)死亡，現(xiàn)在規(guī)定科學(xué)家的動(dòng)作如下：

每個(gè)時(shí)間單位都可以向相鄰的格子移動(dòng)一個(gè)單位

不可以進(jìn)入毒氣的區(qū)域

每個(gè)格子同時(shí)可以容納無(wú)數(shù)個(gè)科學(xué)家

每個(gè)安全膠囊只能容納一個(gè)科學(xué)家

在滿足以上所有條件的情況下，求在核反應(yīng)堆爆炸之前，最多能有多少個(gè)科學(xué)家獲救

題目分析：這個(gè)題目一方面是要求一個(gè)最大值，換種說(shuō)法就是要求一種最優(yōu)方案，但是這個(gè)方案中的每一個(gè)科學(xué)家的決策都會(huì)對(duì)其他科學(xué)家產(chǎn)生影響，所以不能單純的考慮最短路的問(wèn)題，既然是要在全局情況下考慮最優(yōu)方案，我們不妨用網(wǎng)絡(luò)流來(lái)解決，在讀入數(shù)據(jù)后可以預(yù)處理出dp[time][x][y][xx][yy]數(shù)組，表示在time時(shí)間時(shí)可以從點(diǎn)(x,y)到達(dá)點(diǎn)(xx,yy)，同時(shí)可以更簡(jiǎn)化一下，也就是預(yù)處理出在給定時(shí)間t之前，能否從點(diǎn)(x,y)到達(dá)點(diǎn)(xx,yy)，這樣就可以直接建圖了：

源點(diǎn)->每個(gè)科學(xué)家，流量為當(dāng)前格子科學(xué)家的數(shù)量

每個(gè)科學(xué)家->可到達(dá)的安全膠囊，流量為無(wú)窮大

每個(gè)安全膠囊->匯點(diǎn)，流量為當(dāng)前格子安全膠囊的數(shù)量

最大流的建圖比較簡(jiǎn)單，這道題目麻煩就麻煩在該如何預(yù)處理，可以用dfs和bfs，以及一點(diǎn)簡(jiǎn)單的動(dòng)態(tài)規(guī)劃

有一個(gè)細(xì)節(jié)需要注意一下，就是當(dāng)某個(gè)格子(x,y)在第k秒時(shí)才被毒氣擴(kuò)散到，那么這個(gè)格子是可以進(jìn)入的，但如果這個(gè)格子在k秒之前就被毒氣擴(kuò)散到了，那么這個(gè)格子才是不可行的，這個(gè)細(xì)節(jié)對(duì)應(yīng)著dfs中何時(shí)該跳出

代碼：

#include<iostream> #include<cstdlib> #include<string> #include<cstring> #include<cstdio> #include<algorithm> #include<climits> #include<cmath> #include<cctype> #include<stack> #include<queue> #include<list> #include<vector> #include<set> #include<map> #include<sstream> using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=210;struct Edge {int to,w,next; }edge[N*N];//邊數(shù)const int b[4][2]={0,1,0,-1,1,0,-1,0};int n,t,st=N-1,ed=st-1,cap[15][15],sci[15][15],maze[65][15][15],head[N],cnt;//maze[time][x][y]:在time時(shí)間下點(diǎn)(x,y)的狀態(tài):0->安全 1->墻 2->毒氣bool vis[15][15],dp[65][15][15][15][15];//dp[time][x][y][xx][yy]:在time時(shí)間時(shí)能否從點(diǎn)(x,y)到達(dá)點(diǎn)(xx,yy)void addedge(int u,int v,int w) {edge[cnt].to=v;edge[cnt].w=w;edge[cnt].next=head[u];head[u]=cnt++;edge[cnt].to=u;edge[cnt].w=0;//反向邊邊權(quán)設(shè)置為0edge[cnt].next=head[v];head[v]=cnt++; }int d[N],now[N];//深度 當(dāng)前弧優(yōu)化bool bfs(int s,int t)//尋找增廣路 {memset(d,0,sizeof(d));queue<int>q;q.push(s);now[s]=head[s];d[s]=1;while(!q.empty()){int u=q.front();q.pop();for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;if(d[v])continue;if(!w)continue;d[v]=d[u]+1;now[v]=head[v];q.push(v);if(v==t)return true;}}return false; }int dinic(int x,int t,int flow)//更新答案 {if(x==t)return flow;int rest=flow,i;for(i=now[x];i!=-1&&rest;i=edge[i].next){int v=edge[i].to;int w=edge[i].w;if(w&&d[v]==d[x]+1){int k=dinic(v,t,min(rest,w));if(!k)d[v]=0;edge[i].w-=k;edge[i^1].w+=k;rest-=k;}}now[x]=i;return flow-rest; }int solve(int st,int ed) {int ans=0,flow;while(bfs(st,ed))while(flow=dinic(st,ed,inf))ans+=flow;return ans; }int get_id(int x,int y,int k)//science:0 capsule:1 {return (x-1)*n+y+k*n*n; }void dfs(int x,int y,int xx,int yy,int step) {vis[xx][yy]=true;dp[step][x][y][xx][yy]=true;if(step>=t||maze[step][xx][yy])return;for(int k=0;k<4;k++){int nx=xx+b[k][0];int ny=yy+b[k][1];if(nx<=0||ny<=0||nx>n||ny>n)continue;if(vis[nx][ny])continue;if(maze[step][nx][ny])continue; dfs(x,y,nx,ny,step+1);} }void init() {memset(now,0,sizeof(now));memset(head,-1,sizeof(head));cnt=0;for(int time=1;time<=t;time++)//首先預(yù)處理出時(shí)間與(x,y)的關(guān)系{for(int x=1;x<=n;x++){for(int y=1;y<=n;y++){if(maze[time-1][x][y])maze[time][x][y]=maze[time-1][x][y];else{for(int k=0;k<4;k++){int xx=x+b[k][0];int yy=y+b[k][1];if(xx<=0||yy<=0||xx>n||yy>n)continue;if(maze[time-1][xx][yy]==2)maze[time][x][y]=2;}}}}}for(int i=1;i<=n;i++)//預(yù)處理出dp數(shù)組{for(int j=1;j<=n;j++){if(sci[i][j]){memset(vis,false,sizeof(vis));dfs(i,j,i,j,0);}}}for(int i=1;i<=t;i++)for(int x=1;x<=n;x++)for(int y=1;y<=n;y++)for(int xx=1;xx<=n;xx++)for(int yy=1;yy<=n;yy++)//將dp數(shù)組的狀態(tài)統(tǒng)一到時(shí)間t下，方便操作dp[i][x][y][xx][yy]|=dp[i-1][x][y][xx][yy]; }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);scanf("%d%d",&n,&t);for(int i=1;i<=n;i++){char s[15];scanf("%s",s+1);for(int j=1;j<=n;j++){if(isdigit(s[j]))sci[i][j]=s[j]-'0';else if(s[j]=='Y')maze[0][i][j]=1;else if(s[j]=='Z')maze[0][i][j]=2;}}for(int i=1;i<=n;i++){char s[15];scanf("%s",s+1);for(int j=1;j<=n;j++){if(isdigit(s[j]))cap[i][j]=s[j]-'0';}}init();for(int i=1;i<=n;i++){for(int j=1;j<=n;j++){if(sci[i][j])addedge(st,get_id(i,j,0),sci[i][j]);if(cap[i][j])addedge(get_id(i,j,1),ed,cap[i][j]);}}for(int x=1;x<=n;x++)for(int y=1;y<=n;y++)for(int xx=1;xx<=n;xx++)for(int yy=1;yy<=n;yy++)if(sci[x][y]&&cap[xx][yy]&&dp[t][x][y][xx][yy])addedge(get_id(x,y,0),get_id(xx,yy,1),inf);printf("%d\n",solve(st,ed));return 0; }

?

總結(jié)

以上是生活随笔為你收集整理的CodeForces - 78E Evacuation(最大流)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。