[题解] BZOJ 3144 [HNOI2013] 切糕

Chlience

于 2018-03-22 21:36:54 发布

阅读量277

点赞数

分类专栏：【图论】网络流文章标签：省选网络流

本文链接：https://blog.csdn.net/Force_CHL/article/details/79660367

版权

【图论】网络流专栏收录该内容

12 篇文章 0 订阅

订阅专栏

题目描述 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希望找出总的切割点上的不和谐值最小的那个,即 $\sum 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
样例解释
第一组样例中最佳切面的f为f(1,1)=f(2,1)=2,f(1,2)=f(2,2)=1.
第二组样例中最佳切面的f为f(1,1)=f(2,1)=f(1,2)=f(2,2)=1.
数据范围及提示 Data Size & Hint
100%的数据满足 $P,Q,R<=40,0<=D<=R$ ,且给出的所有的不和谐值不超过1000

Solution

将一个切糕切成两半?
切的时候还要满足全局最优?
这不是最小割么(逃

emm,像各位一样的dalao一定都看出来了吧

显然,切糕上的每个切割点将其变为网络流上的一条边,流量即为不和谐值,在切糕的顶上放一个s,在切糕的底下放一个t,跑最小割(最大流)即可

但是我们又要满足一定的光滑性要求,那就让切了两个不满足要求之间的点还是能够联通就好了(无法割开),所以说对与上下左右不满足平滑度的点连上一条容量为INF的边就好了

将点转化为边的方式要仔细看哦!

代码如下

#include <bits/stdc++.h>
using namespace std;
#define INF 2100000000
struct Edge {
    int to;
    int next;
    int flow;
}e[2000010];
int head[200010];
int cur[200010];
int f[200010],d[200010];
int p,q,r,D,s,t,tot=1;
int g[50][50][50],cnt;
int move[4][2]={{1,0},{-1,0},{0,1},{0,-1}};
int read() {
    int ans=0,flag=1;
    char ch=getchar();
    while((ch>'9' || ch<'0') && ch!='-') ch=getchar();
    if(ch=='-') flag=-1,ch=getchar();
    while(ch>='0' && ch<='9') ans=ans*10+ch-'0',ch=getchar();
    return ans*flag;
}
void addedge(int u,int v,int w) {
    e[++tot].to=v;
    e[tot].flow=w;
    e[tot].next=head[u];
    head[u]=tot;

    e[++tot].to=u;
    e[tot].flow=0;
    e[tot].next=head[v];
    head[v]=tot;
}
int dfs(int now,int flow) {
    if(now==t)
        return flow;
    int use=0;
    for(int i=cur[now];i;i=e[i].next) {
        cur[now]=i;
        if(e[i].flow>0 && d[e[i].to]+1==d[now]) {
            int temp=dfs(e[i].to,min(flow-use,e[i].flow));
            use+=temp;
            e[i].flow-=temp;
            e[i^1].flow+=temp;
            if(flow==use)
                return use;
        }
    }
    cur[now]=head[now];

    if(!(--f[d[now]]))
        d[s]=cnt+2;
    ++f[++d[now]];
    return use;
}
int main() {
    p=read(),q=read(),r=read();
    D=read();
    for(int i=1;i<=r+1;i++)
        for(int j=1;j<=p;j++)
            for(int k=1;k<=q;k++)
                g[i][j][k]=++cnt;
    s=0;t=cnt+1;
    for(int i=1;i<=p;i++)
        for(int j=1;j<=q;j++) {
            addedge(s,g[1][i][j],INF);
            addedge(g[r+1][i][j],t,INF);
        }
    for(int i=1;i<=r;i++)
        for(int j=1;j<=p;j++)
            for(int k=1;k<=q;k++) {
                int x=read();
                addedge(g[i][j][k],g[i+1][j][k],x);
            }
    for(int i=D+1;i<=r;i++)
        for(int j=1;j<=p;j++)
            for(int k=1;k<=q;k++)
                for(int l=0;l<4;l++)
                    if(g[i-D][j+move[l][0]][k+move[l][1]])
                        addedge(g[i][j][k],g[i-D][j+move[l][0]][k+move[l][1]],INF);
    int ans=0;
    f[0]=cnt+2;
    while(d[s]<cnt+2)
        ans+=dfs(s,1<<30);
    printf("%d\n",ans);
    return 0;
}