BZOJ3144【HNOI2013】切糕

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=3144

【分析】

    所有人都说这是个经典的最小割模型......蒟蒻泪流满面......

    建一个r+1层的，每层都是p*q的图，从源出发向第一层连inf的边，从第r+1层上每个点出发向汇连inf的边，对i,j,k与i,j,k+1之间连v(i,j,k)的边。

    如果没有高度限制，那么这么做就是在每一条纵轴上选一个点。

    因为有高度限制，所以对于每个层数大于d的点，要向x,y坐标相邻的，层数小d的点连一条inf的边。

    以p=1,q=2,r=5,d=2为例。（先忽略从右边较高的点连向左边较低的点的那些边）

    

    下面是源，上面是汇。如图所示，如果我把图中的红边割掉，那么图中的蓝边也就没有用了。如果此时把绿边割掉，那么还是存在一条从源到汇的路径；如果割绿色的边上面的边，那么源到汇就没有路径可达了。

    仔细理解一下，也就是在左边选了一条边以后，右边的高度小于当前边的边就不能选了。如果把右上向左下连的边也连起来，那么高度就完全限制住了。

    然后跑一下最大流就好了。

【代码】

#include <cstdio>
#include <algorithm>
using namespace std;
const int maxl=45,maxn=75000,maxm=maxn<<4,inf=~0U>>1,mov[4][2]={{0,1},{1,0},{0,-1},{-1,0}};

int head[maxn],next[maxm],E[maxm],F[maxm],Ecnt;
int id[maxl][maxl][maxl];
int p,q,r,d;
int s,e;

inline void Add_Edge(int x,int y,int f) {next[++Ecnt]=head[x];head[x]=Ecnt;E[Ecnt]=y;F[Ecnt]=0;next[++Ecnt]=head[y];head[y]=Ecnt;E[Ecnt]=x;F[Ecnt]=f;}

void Init()
{
	scanf("%d%d%d%d",&p,&q,&r,&d);
	int cnt=0,x;
	Ecnt=1;
	s=++cnt;
	for (int i=1;i<=r+1;i++)
	for (int j=1;j<=p  ;j++)
	for (int k=1;k<=q  ;k++)
	id[i][j][k]=++cnt;
	e=++cnt;
	
	for (int i=1;i<=r;i++)
	for (int j=1;j<=p;j++)
	for (int k=1;k<=q;k++)
	scanf("%d",&x),Add_Edge(id[i][j][k],id[i+1][j][k],x);
	
	for (int i=1;i<=p;i++)
	for (int j=1;j<=q;j++)
	Add_Edge(s,id[1][i][j],inf),Add_Edge(id[r+1][i][j],e,inf);
	
	#define nx (j+mov[mv][0])
	#define ny (k+mov[mv][1])
	for (int i=d+1;i<=r+1;i++)
	for (int j=1;j<=p;j++)
	for (int k=1;k<=q;k++)
	for (int mv=0;mv<4;mv++)
	if (id[i][nx][ny])
	Add_Edge(id[i][j][k],id[i-d][nx][ny],inf);
	#undef nx
	#undef ny
}

int Q[maxn],Vis[maxn],h[maxn];
bool BFS(int time)
{
	int top=1,tail=1;
	Q[1]=s;Vis[1]=time;h[1]=0;
	while (top<=tail)
	{
		int x=Q[top++],hx=h[x];
		for (int i=head[x];i;i=next[i]) if (Vis[E[i]]!=time && F[i^1])
			h[E[i]]=h[x]+1,Vis[E[i]]=time,Q[++tail]=E[i];
		if (Vis[e]==time) break;
	}
	return Vis[e]==time;
}

int DFS(int x,int f)
{
	if (x==e) return f;
	int res=f;
	for (int i=head[x];i && res;i=next[i]) if (h[E[i]]==h[x]+1 && F[i^1])
	{
		int d=DFS(E[i],min(F[i^1],res));
		res-=d;F[i]+=d;F[i^1]-=d;
	}
	if (f==res) h[x]=-1;
	return f-res;
}

void Dinic()
{
	int res=0,time=0;
	while (BFS(++time)) res+=DFS(s,inf);
	printf("%d\n",res);
}

int main()
{
	Init();
	Dinic();
	return 0;
}