[Luogu P3227] [BZOJ 3144] [HNOI2013]切糕

洛谷传送门

BZOJ传送门

题目描述

经过千辛万苦小 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 希望找出总的切割点上的不和谐值最小的那个。

输入输出格式

输入格式：

第一行是三个正整数$P,Q,R$，表示切糕的长$P$、 宽$Q$、高$R$。第二行有一个非负整数$D$，表示光滑性要求。接下来是$R$个$P$行$Q$列的矩阵，第$z$个 矩阵的第$x$行第$y$列是$v(x,y,z)$ ($1\leq x\leq P, 1\leq y\leq Q, 1\leq z\leq R$)。 100%的数据满足$P,Q,R\leq 40，0\leq D\leq R$，且给出的所有的不和谐值不超过$1000$。

输出格式：

仅包含一个整数，表示在合法基础上最小的总不和谐值。

输入输出样例

输入样例#1：

2  2 2
1
6  1
6  1
2  6
2  6

输出样例#1：

6

说明

最佳切面的$f$为$f(1,1)=f(2,1)=2,f(1,2)=f(2,2)=1$

解题分析

最小割经典模型，和这道题如出一辙。

如果没有高度差的限制的话， 我们直接选每个高为$R$的单位立方体上的最小值即可。换句话说， 我们可以把同一水平位置所有高度上的点依次连在一起， 边的权值为其不和谐值， 求出最小割即可。

有了高度限制， 我们就必须保证相邻两个格子超出范围的部分没有被割掉。 怎么办？ 设当前点为$(x,y,z)$， 相邻点的平面坐标为$(a,b)$，直接从$(x,y,z)$向$(a,b,z-D)$连边即可。 这样如果割在了$(a,b,z-D)$之前的话整个图仍然是联通的。 同时为了使这种边不被割掉， 其流量应该设为$INF$。

此题复杂度甚为玄学， 看似完全无法过但却跑的飞快…不过仍需加上当前弧优化才能$AC$。

代码如下：

#include <cstdio>
#include <cstring>
#include <cctype>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <queue>
#define R register
#define IN inline
#define gc getchar()
#define W while
#define MX 200050
#define S 0
#define T 65000
#define INF 100000000
#define base 45
template <class TT>
IN void in(TT &x)
{
    x = 0; R char c = gc;
    W (!isdigit(c)) c = gc;
    W (isdigit(c))
    x = (x << 1) + (x << 3) + c - 48, c = gc;
}
int n, m, h, cnt = -1, del, siz;
int head[MX], layer[MX], val[45][45][45], cpy[MX];
std::queue <int> q;
struct Edge {int to, fl, nex;} edge[MX << 1];
IN void add(R int from, R int to, R int fl)
{
    edge[++cnt] = {to, fl, head[from]}, head[from] = cnt;
    edge[++cnt] = {from, 0, head[to]}, head[to] = cnt;
}
IN bool BFS()
{
    std::memset(layer, 0, sizeof(layer));
    layer[S] = 1; q.push(S);
    R int now;
    W (!q.empty())
    {
        now = q.front(); q.pop();
        for (R int i = head[now]; ~i; i = edge[i].nex)
        {
            if(edge[i].fl && !layer[edge[i].to])
            layer[edge[i].to] = layer[now] + 1, q.push(edge[i].to);
        }
    }
    return layer[T];
}
int DFS(R int now, R int val)
{
    if(now == T) return val;
    int lef = val, buf;
    for (R int &i = head[now]; ~i; i = edge[i].nex)
    {
        if(edge[i].fl && layer[edge[i].to] == layer[now] + 1)
        {
            buf = DFS(edge[i].to, std::min(lef, edge[i].fl));
            if(!buf) continue;
            lef -= buf; edge[i].fl -= buf, edge[i ^ 1].fl += buf;
            if(lef <= 0) return val;
        }
    }
    return val - lef;
}
int Dinic()
{
    std::memcpy(cpy, head, 66000 * sizeof(int));
    int ret = 0;
    W (BFS()) ret += DFS(S, INF), std::memcpy(head, cpy, 66000 * sizeof(int));
    return ret;
}
IN int getid(R int height, R int hang, R int lie) {return (height - 1) * siz + (hang - 1) * m + lie;}
int main(void)
{
    int id, bd, a, b, c;
    std::memset(head, -1, sizeof(head));
    in(n), in(m), in(h), in(del); siz = n * m;
    for (R int i = 1; i <= h; ++i)
        for (R int j = 1; j <= n; ++j)
                for (R int k = 1; k <= m; ++k)
                in(val[i][j][k]);
    for (R int i = 1; i <= h; ++i)
    {
        for (R int j = 1; j <= n; ++j)
        for (R int k = 1; k <= m; ++k)
        {
            id = getid(i, j, k);
            if(i == 1) add(S, id, INF);
            if(i != h) add(id, getid(i + 1, j, k), val[i][j][k]);
            else add(id, T, val[i][j][k]);
            if(i > del)
            {
                if(j > 1) add(id, getid(i - del, j - 1, k), INF);
                if(j < n) add(id, getid(i - del, j + 1, k), INF);
                if(k > 1) add(id, getid(i - del, j, k - 1), INF);
                if(k < m) add(id, getid(i - del, j, k + 1), INF);
            }
        }
    }
    printf("%d", Dinic());
}