[HNOI2013] 切糕

题目描述：

切蛋糕！

题目分析：

这数据范围，还有求最小的值，基本上最小割没跑了.
先不考虑限制.
我们新建一层虚拟蛋糕
从下面一层向上面一层连容量为花费的边
从s向1层连容量为inf的边，代表不可割
从虚拟层向t连容量为inf的边，代表不可割
然后…就不会啦qwqqqq
怎么考虑限制哇
我们从本点向四周高度为h-d的点连容量为inf的边就可以限制住了！
为什么呢？

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

题目链接：

BZOJ 3144
COGS 2398
Luogu 3227

Ac 代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <queue>
#define inf 0x7fffffff
const int maxm=500*50*50;
int head[maxm],to[maxm<<1],cap[maxm<<1],net[maxm<<1];
int cnt=1;
inline void addedge(int u,int v,int c)
{
	cnt++;
	to[cnt]=v,cap[cnt]=c,net[cnt]=head[u],head[u]=cnt;
	cnt++;
	to[cnt]=u,cap[cnt]=0,net[cnt]=head[v],head[v]=cnt;
}
namespace Maxflow{
	int deep[maxm];
	std::queue <int> dl;
	inline bool BFS(int s,int t)
	{
		memset(deep,-1,sizeof(deep));
		deep[s]=0;
		dl.push(s);
		while(!dl.empty())
		{
			int now=dl.front();
			dl.pop();
			for(int i=head[now];i;i=net[i])
			if(cap[i]&&(deep[to[i]]==-1))
			 deep[to[i]]=deep[now]+1,dl.push(to[i]);
		}
		return deep[t]!=-1;
	}
	int DFS(int now,int flow,int t)
	{
		if(now==t) return flow;
		int used=0,w;
		for(int i=head[now];i;i=net[i])
		if(cap[i]&&deep[to[i]]==deep[now]+1)
		{
			w=DFS(to[i],std::min(cap[i],flow-used),t);
			cap[i]-=w,cap[i^1]+=w;
			used+=w;
			if(flow==used) return flow;
		}
		if(!used) deep[now]=-1;
		return used;
	}
	inline int Dinic(int s,int t)
	{
		int maxflow=0;
		while(BFS(s,t)) maxflow+=DFS(s,inf,t);
		return maxflow;
	} 
}
int P,Q,R,D;
int id[70][70][70],tot;
int dx[]={0,0,0,1,-1},dy[]={0,-1,1,0,0};
bool judge(int x,int y){return (x>=1)&&(x<=P)&&(y>=1)&&(y<=Q);}
int main()
{
	//freopen("nutcake.in","r",stdin);
	//freopen("nutcake.out","w",stdout);
	scanf("%d%d%d%d",&P,&Q,&R,&D);
	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]=++tot;
	int s=0,t=tot+10;
	for(int i=1;i<=R;i++)
	 for(int j=1;j<=P;j++)
	  for(int k=1;k<=Q;k++)
	   {
	   	 int x;
	   	 scanf("%d",&x);
	   	 addedge(id[i][j][k],id[i+1][j][k],x);
	   }
	for(int i=1;i<=P;i++)
	 for(int j=1;j<=Q;j++)
	  addedge(s,id[1][i][j],inf),addedge(id[R+1][i][j],t,inf);
	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 l=1;l<=4;l++)
        {
        	int xx=j+dx[l],yy=k+dy[l];
        	if(!judge(xx,yy)) continue;
        	addedge(id[i][j][k],id[i-D][xx][yy],inf);
        }
    printf("%d\n",Maxflow::Dinic(s,t));
    return 0;
}