[题解] BZOJ 3144 [HNOI2013] 切糕

题目描述 Description
经过千辛万苦小A得到了一块切糕,切糕的形状是长方体,小A打算拦腰将切糕切成两半分给小B.出于美观考虑,小A希望切面能尽量光滑且和谐.于是她找到你,希望你能帮她找出最好的切割方案.

出于简便考虑,我们将切糕视作一个长 P P 、宽Q、高 R R 的长方体点阵.我们将位于第z层中第 x x 行、第y列上 (1xP,1yQ,1zR) ( 1 ≤ x ≤ P , 1 ≤ y ≤ Q , 1 ≤ z ≤ R ) 的点称为 (x,y,z) ( x , y , z ) ,它有一个非负的不和谐值 v(x,y,z) v ( x , y , z ) .一个合法的切面满足以下两个条件:

1.与每个纵轴(一共有 PQ P ∗ Q 个纵轴)有且仅有一个交点.即切面是一个函数 f(x,y) f ( x , y ) ,对于所有 1xP,1yQ 1 ≤ x ≤ P , 1 ≤ y ≤ Q ,我们需指定一个切割点 f(x,y) f ( x , y ) ,且 1f(x,y)R 1 ≤ f ( x , y ) ≤ R .

2.切面需要满足一定的光滑性要求,即相邻纵轴上的切割点不能相距太远.对于所有的 1x,xP 1 ≤ x , x ′ ≤ P 1y,yQ 1 ≤ y , y ′ ≤ Q ,若 |xx|+|yy|=1 | x − x ′ | + | y − y ′ | = 1 ,则 |f(x,y)f(x,y)|D | f ( x , y ) − f ( x ′ , y ′ ) | ≤ D ,其中 D D 是给定的一个非负整数.

可能有许多切面f满足上面的条件,小A希望找出总的切割点上的不和谐值最小的那个,即 v(x,y,f(x,y)) ∑ v ( x , y , f ( x , y ) ) 最小

输入描述 Input Description
第一行是三个正整数 P,Q, P , Q , R,表示切糕的长 P P 、 宽Q、高 R R .第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q列的矩阵,第 z z 个 矩阵的第x行第 y 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 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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值