题目描述 Description
经过千辛万苦小A得到了一块切糕,切糕的形状是长方体,小A打算拦腰将切糕切成两半分给小B.出于美观考虑,小A希望切面能尽量光滑且和谐.于是她找到你,希望你能帮她找出最好的切割方案.
出于简便考虑,我们将切糕视作一个长 P P 、宽、高 R R 的长方体点阵.我们将位于第层中第 x x 行、第列上 (1≤x≤P,1≤y≤Q,1≤z≤R) ( 1 ≤ x ≤ P , 1 ≤ y ≤ Q , 1 ≤ z ≤ R ) 的点称为 (x,y,z) ( x , y , z ) ,它有一个非负的不和谐值 v(x,y,z) v ( x , y , z ) .一个合法的切面满足以下两个条件:
1.与每个纵轴(一共有 P∗Q P ∗ Q 个纵轴)有且仅有一个交点.即切面是一个函数 f(x,y) f ( x , y ) ,对于所有 1≤x≤P,1≤y≤Q 1 ≤ x ≤ P , 1 ≤ y ≤ Q ,我们需指定一个切割点 f(x,y) f ( x , y ) ,且 1≤f(x,y)≤R 1 ≤ f ( x , y ) ≤ R .
2.切面需要满足一定的光滑性要求,即相邻纵轴上的切割点不能相距太远.对于所有的 1≤x,x′≤P 1 ≤ x , x ′ ≤ P 和 1≤y,y′≤Q 1 ≤ y , y ′ ≤ Q ,若 |x−x′|+|y−y′|=1 | x − x ′ | + | y − y ′ | = 1 ,则 |f(x,y)−f(x′,y′)|≤D | f ( x , y ) − f ( x ′ , y ′ ) | ≤ D ,其中 D D 是给定的一个非负整数.
可能有许多切面满足上面的条件,小A希望找出总的切割点上的不和谐值最小的那个,即 ∑v(x,y,f(x,y)) ∑ v ( x , y , f ( x , y ) ) 最小
输入描述 Input Description
第一行是三个正整数
P,Q,
P
,
Q
,
R,表示切糕的长
P
P
、 宽、高
R
R
.第二行有一个非负整数,表示光滑性要求.接下来是R个P行Q列的矩阵,第
z
z
个 矩阵的第行第
y
y
列是.
输出描述 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;
}