题目描述
经过千辛万苦小 A 得到了一块切糕，切糕的形状是长方体，小 A 打算拦腰将切糕切成两半分给小 B。出于美观考虑，小 A 希望切面能尽量光滑且和谐。于是她找到你，希望你能帮她找出最好的切割方案。
出于简便考虑，我们将切糕视作一个长 P、宽 Q、高 R 的长方体点阵。我们将位于第 z层中第 x 行、第 y 列上(1≤x≤P, 1≤y≤Q, 1≤z≤R)的点称为(x,y,z)，它有一个非负的不和谐值 v(x,y,z)。一个合法的切面满足以下两个条件：
与每个纵轴(一共有 P*Q 个纵轴)有且仅有一个交点。即切面是一个函数 f(x,y)，对于所有 1≤x≤P, 1≤y≤Q,我们需指定一个切割点 f(x,y),且 1≤f(x,y)≤R。
切面需要满足一定的光滑性要求，即相邻纵轴上的切割点不能相距太远。对于所有的 1≤x,x’≤P 和 1≤y,y’≤Q，若|x-x’|+|y-y’|=1，则|f(x,y)-f(x’,y’)| ≤D，其中 D 是给定的一个非负整数。 可能有许多切面f 满足上面的条件，小A 希望找出总的切割点上的不和谐值最小的那个。
输入输出格式
输入格式：
第一行是三个正整数P,Q,R，表示切糕的长P、 宽Q、高R。第二行有一个非负整数D，表示光滑性要求。接下来是R个P行Q列的矩阵，第z个 矩阵的第x行第y列是v(x,y,z) (1<=x<=P, 1<=y<=Q, 1<=z<=R)。 100%的数据满足P,Q,R<=40，0<=D<=R，且给出的所有的不和谐值不超过1000。
输出格式：
仅包含一个整数，表示在合法基础上最小的总不和谐值。
QwQ说实话 这个题一开始我连题意都没有理解，
这里是简化版的题意
可以这么来理解，就是先假设一个平面有(p*q)个点，然后对于每个点((i,j))，必须选择一个值(k)，得到一个不开心值(v(i,j,k))，相邻两个点的(k)值之差要小雨(d)，求最小化不开心值的和QwQ
其实一开始，看这个题.....没什么思路 嗯
看了题解，听过某大佬的讲解之后，才知道是个最小割
既然这样，我们就不妨先忽略这个(d)的限制条件，考虑没有限制的情况
我们不妨设源点和汇点(s,t)，增加一个附加层(r+1)层
然后(s)向第一层的点连边，(r+1)层的边向(t)连边
对于中间的点，很显然是((i,j,k)->(i,j,k+1)) 流量为(v(i,j,k)),很容易证明，这样跑最小割，每一个纵轴，会且仅会切一道。
那么加上限制.....该怎么限制呢
让删除这些边 不影响连通 不就好了咯
对于相邻的两个点((i,j))和((x,y)),我们将((i,j,k)->(x,y,k-d))流量是(inf) 就能保证QwQ当((i,j))选了k的时候，((x,y))必须要弄距离是D以内的
QwQ相反同理
然后....一定一定记得看清楚读入的顺序
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
const int maxn = 55;
const int maxm = 2e6+1e2;
const int inf = 1e9;
int a[maxn][maxn][maxn];
int point[100010],nxt[maxm],val[maxm],to[maxm];
int h[100010];
queue<int> que;
int cnt=1;
int s,t;
int dx[5]={0,1,0,-1,0};
int dy[5]={0,0,1,0,-1};
void addedge(int x,int y,int w)
{
nxt[++cnt]=point[x];
to[cnt]=y;
val[cnt]=w;
point[x]=cnt;
}
void insert(int x,int y,int w)
{
addedge(x,y,w);
addedge(y,x,0);
}
bool bfs(int s)
{
memset(h,-1,sizeof(h));
h[s]=0;
que.push(s);
while (!que.empty())
{
int x = que.front();
que.pop();
for (int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (val[i]>0 &&h[p]==-1)
{
h[p]=h[x]+1;
que.push(p);
}
}
}
if (h[t]==-1) return false;
else return true;
}
int dfs(int x,int low)
{
if (x==t || low==0) return low;
int totflow=0;
for (int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (val[i]>0 && h[p]==h[x]+1)
{
int tmp = dfs(p,min(low,val[i]));
val[i]-=tmp;
val[i^1]+=tmp;
low-=tmp;
totflow+=tmp;
if (low==0) return totflow;
}
}
if (low>0) h[x]=-1;
return totflow;
}
int solve()
{
int ans=0;
while (bfs(s))
{
ans+=dfs(s,inf);
}
return ans;
}
int p,q,r;
int d;
int main()
{
p=read(),q=read(),r=read();
d=read();
for (int k=1;k<=r;k++)
for (int i=1;i<=p;i++)
for (int j=1;j<=q;j++)
a[i][j][k]=read();
s=100004;
t=s+1;
for (int i=1;i<=p;i++)
for (int j=1;j<=q;j++)
{
insert(s,(i-1)*q+j,inf);
}
for (int i=1;i<=p;i++)
{
for (int j=1;j<=q;j++)
{
int now = (i-1)*q+j;
for (int k=1;k<=r;k++)
{
insert(now+p*q*(k-1),now+p*q*k,a[i][j][k]);
}
}
}
for (int i=1;i<=p;i++)
for (int j=1;j<=q;j++)
{
insert((i-1)*q+j+p*q*r,t,inf);
}
for (int i=1;i<=p;i++)
for (int j=1;j<=q;j++)
{
for (int o=1;o<=4;o++)
{
int x = i+dx[o],y=j+dy[o];
if (x<1 || x>p || y<1 || y>q) continue;
for (int k=d;k<=r;k++)
{
insert(p*q*k+(i-1)*q+j,p*q*(k-d)+(x-1)*q+y,inf);
}
}
}
cout<<solve();
return 0;
}