bzoj 3144: [Hnoi2013]切糕(最小割)
来源:互联网 发布:实用数据再分析法 pdf 编辑:程序博客网 时间:2024/06/08 18:08
3144: [Hnoi2013]切糕
Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1393 Solved: 768
[Submit][Status][Discuss]
Description
Input
第一行是三个正整数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)。
100%的数据满足P,Q,R≤40,0≤D≤R,且给出的所有的不和谐值不超过1000。
Output
仅包含一个整数,表示在合法基础上最小的总不和谐值。
Sample Input
2 2 2
1
6 1
6 1
2 6
2 6
1
6 1
6 1
2 6
2 6
Sample Output
6
HINT
最佳切面的f为f(1,1)=f(2,1)=2,f(1,2)=f(2,2)=1
Source
题解:最小割
对于z轴,每一层建p*q个虚点。
s->(i,j,1) inf
(i,j,k)->(i,j,k+1) v[i][j][k]
(i,j,k)->(i-1,j,k-d) inf
(i,j,k)->(i,j+1,k-d) inf
(i,j,k)->(i+1,j,k-d) inf
(i,j,k)->(i,j-1,k-d) inf
(i,j,z)->t v[i][j][z]
#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<queue>#define N 200003#define inf 1000000000using namespace std;int n,m,h,d;int point[N],next[N],v[N],remain[N],tot;int deep[N],num[N],cur[N],last[N];int pos[43][43][43],v1[43][43][43];int xx[10]={0,1,0,-1},yy[10]={1,0,-1,0};void add(int x,int y,int z){//cout<<x<<" "<<y<<" "<<z<<endl;tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z;tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0;}int addflow(int s,int t){int now=t; int ans=inf;while (now!=s){ans=min(ans,remain[last[now]]);now=v[last[now]^1];}now=t;while (now!=s){remain[last[now]]-=ans;remain[last[now]^1]+=ans;now=v[last[now]^1];}return ans;}void bfs(int s,int t){for (int i=s;i<=t;i++) deep[i]=t;queue<int> p; p.push(t);deep[t]=0;while(!p.empty()){int now=p.front(); p.pop();for (int i=point[now];i!=-1;i=next[i]) if (deep[v[i]]==t&&remain[i^1]) deep[v[i]]=deep[now]+1,p.push(v[i]);}}int isap(int s,int t){bfs(s,t);for (int i=s;i<=t;i++) cur[i]=point[i];for (int i=s;i<=t;i++) num[deep[i]]++;int now=s; int ans=0;while (deep[s]<t){if (now==t){ans+=addflow(s,t);now=s;}bool f=false;for (int i=cur[now];i!=-1;i=next[i]) if (deep[now]==deep[v[i]]+1&&remain[i]) { cur[now]=i; last[v[i]]=i; now=v[i]; f=true; break; }if (!f){int minn=t;for (int i=point[now];i!=-1;i=next[i]) if (remain[i]) minn=min(minn,deep[v[i]]);if (!--num[deep[now]]) break;deep[now]=minn+1;num[deep[now]]++;cur[now]=point[now];if (now!=s) now=v[last[now]^1];}}return ans;}int main(){freopen("a.in","r",stdin);tot=-1;memset(point,-1,sizeof(point));memset(next,-1,sizeof(next));scanf("%d%d%d%d",&n,&m,&h,&d);int s=0; int t=n*m*h+1; int cnt=0;for (int k=1;k<=h;k++){for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) { scanf("%d",&v1[i][j][k]); pos[i][j][k]=++cnt; if (k==1) add(s,cnt,inf); else add(pos[i][j][k-1],cnt,v1[i][j][k-1]); if (k==h) add(cnt,t,v1[i][j][k]); if (k-d>0) { for (int t=0;t<4;t++) { int nowx=i+xx[t]; int nowy=j+yy[t]; if (nowx>0&&nowy>0&&nowx<=n&&nowy<=m) add(cnt,pos[nowx][nowy][k-d],inf); } } }}printf("%d\n",isap(s,t));}
0 0
- bzoj 3144: [Hnoi2013]切糕(最小割)
- [BZOJ 3144][Hnoi2013]切糕(最小割)
- bzoj 3144: [Hnoi2013]切糕 最小割
- BZOJ 3144: [Hnoi2013]切糕 最小割
- [BZOJ]3144: [Hnoi2013]切糕 最小割
- [最小割] BZOJ 3144: [Hnoi2013]切糕
- 3144: [Hnoi2013]切糕 思路题 最小割
- [BZOJ3144][Hnoi2013]切糕(最小割)
- [BZOJ3144][HNOI2013]切糕(最小割)
- bzoj3144 [Hnoi2013]切糕(最小割)
- 【bzoj3144】【HNOI2013】【切糕】【最小割】
- 【BZOJ3144】[Hnoi2013]切糕【最小割】
- [BZOJ3144][HNOI2013]切糕 最小割
- BZOJ 3144: [Hnoi2013]切糕
- 【BZOJ 3144】 [Hnoi2013]切糕
- BZOJ 3144 HNOI 2013 切糕 最小割
- bzoj 3144 切糕 最小割 解题报告
- bzoj3144: [Hnoi2013]切糕 DINIC最小割
- 下载过程中修改文件名
- Oracle删除数据如何恢复数据
- 什么是NoSql
- Generating all permutations of a given string(JAVA)
- 无刷新提交表单(php+js+ajax简单实例 )
- bzoj 3144: [Hnoi2013]切糕(最小割)
- crypto-js,require,nodejs,加密使用手册
- C++ char*,const char*,string的相互转换
- Java JDBC事务与JTA分布式事务
- 性能优化(一):迭代DOM真的很慢吗?
- el表达式的function标签
- JAVA的节点流和处理流以及流的关闭顺序
- C/C++求职宝典21个重点笔记(容易混淆的东西)
- js中几种实用的跨域方法原理详解