总算做了一道2011以后的省选题了……
原题：
图片题面好评！
P,Q,R≤40，0≤D≤R，给出的所有的不和谐值不超过1000。
文本样例好评！
恩这个是听妹主席讲过后会写的，首先把每个点拆成链，那么割掉这个链上的某条边就表示这个点选了某个权值，边的流量就是这个点设成这个值的花费，表示要花费掉这些代价来把这条边割掉
相邻两个点的值之差<=d怎么搞呐，假设两个相邻的点x，y，那么x的第i个点往y的i-d个点连oo的边，即可，大概就像酱紫：
如果有两个被割掉的线段距离超过d，就依然存在一条路径从s到t，就要继续割
差不多是酱紫：
恩就酱
代码：
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int oo=;
int rd(){int z=,mk=; char ch=getchar();
while(ch<''||ch>''){if(ch=='-')mk=-; ch=getchar();}
while(ch>=''&&ch<=''){z=(z<<)+(z<<)+ch-''; ch=getchar();}
return z*mk;
}
struct ddd{int nxt,y,v,rvs;}e[]; int lk[],ltp=;
inline void ist(int x,int y,int z){
e[++ltp].nxt=lk[x],lk[x]=ltp,e[ltp].y=y,e[ltp].v=z,e[ltp].rvs=ltp+;
e[++ltp].nxt=lk[y],lk[y]=ltp,e[ltp].y=x,e[ltp].v=,e[ltp].rvs=ltp-;
}
const int fx[]={,-,,},fy[]={,,,-};
int n,m,h,d; int N; int s,t;
int lvl[];
int q[],hd=;
bool gtlvl(){
memset(lvl,,sizeof(lvl));
q[hd=]=s,lvl[s]=;
for(int k=;k<=hd;++k)
for(int i=lk[q[k]];i;i=e[i].nxt)if(e[i].v && !lvl[e[i].y])
lvl[e[i].y]=lvl[q[k]]+,q[++hd]=e[i].y;
return lvl[t];
}
int mxflw(int x,int y){
if(x==t) return y;
int bwl=,flw=;
for(int i=lk[x];i && bwl<y;i=e[i].nxt)if(e[i].v && lvl[e[i].y]==lvl[x]+)
if((flw=mxflw(e[i].y,min(y-bwl,e[i].v)))){
bwl+=flw;
e[i].v-=flw,e[e[i].rvs].v+=flw;
}
if(!bwl) lvl[x]=;
return bwl;
}
int dnc(){
int bwl=,flw=;
while(gtlvl())while((flw=mxflw(s,oo))) bwl+=flw;
return bwl;
}
bool chck(int x,int y,int z){
return x+fx[z]>= && x+fx[z]<=n && y+fy[z]>= && y+fy[z]<=m;
}
inline int gtid(int x,int y,int z){ return z*N+(x-)*m+y;}
int main(){//freopen("ddd.in","r",stdin);
cin>>n>>m>>h>>d; N=n*m; s=,t=N*(h+)+;
for(int k=;k<=h;++k)for(int i=;i<=n;++i)for(int j=;j<=m;++j)
ist(gtid(i,j,k-),gtid(i,j,k),rd());
for(int i=;i<=n;++i)for(int j=;j<=m;++j)for(int k=;k<;++k)if(chck(i,j,k))
for(int p=;p+d<=h;++p){
ist(gtid(i,j,p+d),gtid(i+fx[k],j+fy[k],p),oo);
ist(gtid(i+fx[k],j+fy[k],p+d),gtid(i,j,p),oo);
}
for(int i=;i<=n;++i)for(int j=;j<=m;++j) ist(s,gtid(i,j,),oo),ist(gtid(i,j,h),t,oo);
cout<<dnc()<<endl;
return ;
}