[bzoj3144][HNOI2013] 切糕 最小割

3144:[HNOI2013]切糕

题目描述 Description
经过千辛万苦小 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)。一个合法的切面满足以下两个条件：
1. 与每个纵轴(一共有 P*Q 个纵轴)有且仅有一个交点。即切面是一个函数 f(x,y)，对于所有 1≤x≤P, 1≤y≤Q,我们需指定一个切割点 f(x,y),且 1≤f(x,y)≤R。
2. 切面需要满足一定的光滑性要求，即相邻纵轴上的切割点不能相距太远。对于所有的 1≤x,x’≤P 和 1≤y,y’ ≤Q，若|x-x’|+|y-y’|=1，则|f(x,y)-f(x’,y’)| ≤D，其中 D 是给定的一个非负整数。
可能有许多切面f 满足上面的条件，小A 希望找出总的切割点上的不和谐值最小的那个，即 ∑v(x,y, f(x,y))最小。

输入描述 Input Description
输入文件第一行是三个正整数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)。

输出描述 Output Description
输出仅包含一个整数，表示在合法基础上最小的总不和谐值。
样例输入 Sample Input
input1
2 2 2
1
6 1
6 1
2 6
2 6
input2
2 2 2
0
5 1
5 1
2 5
2 5
样例输出 Sample Output
output1
6
output2
12
数据范围及提示 Data Size & Hint
100%的数据满足P,Q,R≤40，0≤D≤R，且给出的所有的不和谐值不超过1000。

这波很亏！！！！
反正错了很多道，都是邻接表的错！！！！！
最小割模型，多看几遍题就知道了

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
const int INF = 1000000000;
const int MX = 100000 + 5;
const int MAXN = 1000000 + 5;
using namespace std;
int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int P,Q,R,T,D,cnt=1; 
int dx[4] = {0,0,1,-1},dy[4] = {1,-1,0,0};
int head[MX], cur[MX], h[MX], q[MX];
struct edge {
    int to, next, v;
}e[MAXN];
int f[45][45][45];
void adde( int u, int v, int w ) {
    e[++cnt].to=v;e[cnt].next=head[u];head[u]=cnt;e[cnt].v=w;
    e[++cnt].to=u;e[cnt].next=head[v];head[v]=cnt;e[cnt].v=0;
}
int p(int x,int y,int z){
    if(z==0)return 0;
    return (z-1)*P*Q+(x-1)*Q+y;
}
void build() {
    for( int i = 1; i <= P; i++ ) 
        for( int j = 1; j <= Q; j++ ) {
            for( int k = 1; k <= R; k++ ) {
                adde( p(i,j,k-1), p(i,j,k), f[i][j][k]);
                if( k > D )
                for( int jk = 0; jk < 4; jk++ ) {
                    int xx = i + dx[jk]; int yy = j + dy[jk];
                    if( xx < 1 || yy < 1 || xx > P || yy > Q ) continue;
                    adde( p( i, j, k), p( xx, yy, k-D), INF );
                }
            }
            adde( p( i, j, R ), T, INF);
        }
}
bool bfs() {
    int hd = 0,tl = 1;
    memset(h,-1,sizeof(h));
    q[0] = 0;h[0] = 0;
    while( hd != tl) {
        int now = q[hd];hd++;
        for( int i = head[now]; i; i = e[i].next)
            if( h[e[i].to] == -1 && e[i].v ) {
                h[e[i].to] = h[now]+1;
                q[tl++] = e[i].to;
            }
    }
    return h[T] != -1;
}
int dfs( int x, int f ) {
    if( x == T ) return f;
    int w,used = 0;
    for( int i = cur[x]; i; i = e[i].next )
        if( h[e[i].to] == h[x]+1 ){
            w = dfs( e[i].to, min( f - used, e[i].v));
            used += w; e[i].v -= w; e[i^1].v += w;
            if( e[i].v ) cur[x] = i;
            if( used == f ) return f;
        }
    if( !used ) h[x] = -1;
    return used;
}
int dinic() {
    int tmp = 0;
    while( bfs() ) {
        for( int i = 0; i <= T; i++ ) cur[i] = head[i];
        tmp += dfs( 0, INF );
    }
    return tmp;
}
int main() {
    P = read(); Q = read(); R = read();
    T = P*Q*R+1;
    D = read();
    for( int i = 1; i <= R; i++ )
        for( int j = 1;j <= P; j++ )
            for( int k = 1; k <= Q; k++ )
                f[j][k][i] = read();
    build();
    printf("%d\n",dinic());
    return 0;
}