题目描述:
切蛋糕!
题目分析:
这数据范围,还有求最小的值,基本上最小割没跑了.
先不考虑限制.
我们新建一层虚拟蛋糕
从下面一层向上面一层连容量为花费的边
从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;
}