題目：https://www.lydsy.com/JudgeOnline/problem.php?id=2007

一個(gè)點(diǎn)的高度一定不是0就是1。答案一定形如一個(gè)左上角的連通塊全是0的點(diǎn)、一個(gè)右下角的連通塊全是1的點(diǎn)。

注意從東到西還有從南到北的邊也有用！因?yàn)椴灰欢ㄊ且粋€(gè)階梯形的，還可以拐來(lái)拐去，只是一定是兩個(gè)連通塊罷了。

所以最小割一下那個(gè)分界線就行了。但會(huì)TLE。

#include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int N=500*501+5,M=4000*501+5,INF=2e6+5; int n,t,bh[505][505],hd[N],xnt=1,cur[N],to[M],nxt[M],cap[M]; int dfn[N],q[N],he,tl; int Mn(int a,int b){return a<b?a:b;} int rdn() {int ret=0;bool fx=1;char ch=getchar();while(ch>'9'||ch<'0'){if(ch=='-')fx=0;ch=getchar();}while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();return fx?ret:-ret; } void add(int x,int y,int z) {to[++xnt]=y;nxt[xnt]=hd[x];hd[x]=xnt;cap[xnt]=z;to[++xnt]=x;nxt[xnt]=hd[y];hd[y]=xnt;cap[xnt]=0; } bool bfs() {memset(dfn,0,sizeof dfn);dfn[0]=1;q[he=tl=1]=0;while(he<=tl){int k=q[he++];for(int i=hd[k],v;i;i=nxt[i])if(cap[i]&&!dfn[v=to[i]])dfn[v]=dfn[k]+1,q[++tl]=v;}return dfn[t]; } int dinic(int cr,int flow) {if(cr==t)return flow;int use=0;for(int& i=cur[cr],v;i;i=nxt[i])if(cap[i]&&dfn[v=to[i]]==dfn[cr]+1){int tmp=dinic(v,Mn(flow-use,cap[i]));if(!tmp)dfn[v]=0;use+=tmp;cap[i]-=tmp;cap[i^1]+=tmp;if(use==flow)return use;}return use; } int main() {n=rdn();for(int i=0;i<=n;i++)for(int j=0;j<=n;j++)bh[i][j]=t++;t--; int d;for(int i=0;i<=n;i++)for(int j=1;j<=n;j++)d=rdn(),add(bh[i][j-1],bh[i][j],d);for(int i=1;i<=n;i++)for(int j=0;j<=n;j++)d=rdn(),add(bh[i-1][j],bh[i][j],d);for(int i=0;i<=n;i++)for(int j=1;j<=n;j++) d=rdn(),add(bh[i][j],bh[i][j-1],d);for(int i=1;i<=n;i++)for(int j=0;j<=n;j++)d=rdn(),add(bh[i][j],bh[i-1][j],d);int ans=0;while(bfs())memcpy(cur,hd,sizeof hd),ans+=dinic(0,INF);printf("%d\n",ans);return 0; } View Code

可以轉(zhuǎn)成最短路。注意邊的方向。

學(xué)習(xí)了學(xué)長(zhǎng)的不顯式建圖的方法。大概 dis[ ][ ] 記錄的就是從起點(diǎn)走到格子的距離，再記4個(gè) dis[ ][ ] 表示它的周圍4條邊的容量，之類的。

#include<cstdio> #include<cstring> #include<algorithm> #include<queue> using namespace std; const int N=505,INF=2e6+5; int n,dis[5][N][N],ans=INF;bool vis[N][N]; struct Node{int x,y,dis;Node(int a=0,int b=0,int d=0):x(a),y(b),dis(d) {}bool operator< (const Node &b)const{return dis>b.dis;} }; priority_queue<Node> q; int Mn(int a,int b){return a<b?a:b;} int rdn() {int ret=0;bool fx=1;char ch=getchar();while(ch>'9'||ch<'0'){if(ch=='-')fx=0;ch=getchar();}while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();return fx?ret:-ret; } void add(int x,int y,int d) {if(d<dis[4][x][y])dis[4][x][y]=d,q.push(Node(x,y,d)); } void dj() {for(int i=1;i<=n;i++)add(1,i,dis[0][1][i]);for(int i=1;i<=n;i++)add(i,n,dis[1][i][n+1]);while(q.size()){int x=q.top().x,y=q.top().y,d=q.top().dis; q.pop();if(vis[x][y])continue; vis[x][y]=1;if(x<n)add(x+1,y,d+dis[0][x+1][y]);//x+1(up)if(y>1)add(x,y-1,d+dis[1][x][y]);if(x>1)add(x-1,y,d+dis[2][x-1][y]);//x-1(dn)if(y<n)add(x,y+1,d+dis[3][x][y]);}for(int i=1;i<=n;i++)ans=Mn(ans,dis[4][i][1]+dis[1][i][1]);for(int i=1;i<=n;i++)ans=Mn(ans,dis[4][n][i]+dis[0][n+1][i]);//n+1 } int main() {n=rdn();int d=n+1;for(int i=1;i<=d;i++)for(int j=1;j<=n;j++)dis[0][i][j]=rdn();//upfor(int i=1;i<=n;i++)for(int j=1;j<=d;j++)dis[1][i][j]=rdn();//leftfor(int i=0;i<=n;i++)//0~n & 1~n !!!for(int j=1;j<=n;j++)dis[2][i][j]=rdn();//dnfor(int i=1;i<=n;i++)for(int j=0;j<=n;j++)dis[3][i][j]=rdn();//rightmemset(dis[4],0x3f,sizeof dis[4]);dj();printf("%d\n",ans);return 0; } View Code

轉(zhuǎn)載于:https://www.cnblogs.com/Narh/p/10165510.html

總結(jié)

以上是生活随笔為你收集整理的bzoj 2007 [Noi2010]海拔——最小割转最短路的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

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