一、题目
二、解法
其实题目就提示你这是一道最小割的题目,首先可以把没有限制的图建出来。
- ( x , y ) (x,y) (x,y)上的点看成 [ 0 , r ] [0,r] [0,r]这 r + 1 r+1 r+1个点,每两个相邻点连边,容量为给定的值。
- ( x , y , 0 ) (x,y,0) (x,y,0)连源点, ( x , y , r ) (x,y,r) (x,y,r)连汇点,容量为 i n f inf inf,表示这条边砍不掉。
发现这样建图答案就是
(
x
,
y
)
(x,y)
(x,y)中最小值的求和,我们考虑把相邻四个点距离差至多为
d
d
d的限制在图上体现出来,大概是这个样子的:
我们把
(
x
,
y
,
k
)
(x,y,k)
(x,y,k)每一个先周围的
(
t
x
,
t
y
,
k
−
d
)
(tx,ty,k-d)
(tx,ty,k−d)的连一条容量为
i
n
f
inf
inf的边,这样操作我们就可以把
k
k
k的选择范围固定在
[
k
−
d
,
r
]
[k-d,r]
[k−d,r]之中,然后对面的还能反向选择,所以就有了一个长度为
2
d
2d
2d的选择区间,在建好的图上跑最大流,就是答案。
#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
const int MAXN = 100005;
using namespace std;
int read()
{
int num=0,flag=1;char c;
while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
return num*flag;
}
int n,m,r,d,S,T,tot=1,ans,f[MAXN],cur[MAXN],dis[MAXN];
int dx[4]={1,-1},dy[4]={0,0,1,-1},a[41][41][41];
queue<int> q;
struct edge
{
int v,c,next;
}e[MAXN*10];
int cal(int x,int y,int z)
{
return z*n*m+(x-1)*m+y;
}
void add_edge(int u,int v,int c)
{
e[++tot]=edge{v,c,f[u]},f[u]=tot;
e[++tot]=edge{u,0,f[v]},f[v]=tot;
}
int bfs()
{
memset(dis,0,sizeof dis);
dis[S]=1;
q.push(S);
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=f[u];i;i=e[i].next)
{
int v=e[i].v;
if(e[i].c>0 && !dis[v])
{
dis[v]=dis[u]+1;
q.push(v);
}
}
}
if(!dis[T]) return 0;
return 1;
}
int dfs(int u,int ept)
{
if(u==T) return ept;
int flow=0,tmp=0;
for(int &i=cur[u];i;i=e[i].next)
{
int v=e[i].v;
if(dis[u]+1==dis[v] && e[i].c>0)
{
tmp=dfs(v,min(e[i].c,ept));
if(!tmp) continue;
ept-=tmp;
e[i].c-=tmp;
e[i^1].c+=tmp;
flow+=tmp;
if(!ept) break;
}
}
return flow;
}
int main()
{
n=read();m=read();r=read();
d=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
add_edge(1,cal(i,j,0)+1,inf);
add_edge(cal(i,j,r)+1,cal(n,m,r)+2,inf);
}
for(int k=1;k<=r;k++)
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
add_edge(cal(i,j,k-1)+1,cal(i,j,k)+1,read());
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
for(int k=d;k<=r;k++)
{
for(int l=0;l<4;l++)
{
int tx=i+dx[l],ty=j+dy[l];
if(tx>=1 && tx<=n && ty>=1 && ty<=m)
add_edge(cal(i,j,k)+1,cal(tx,ty,k-d)+1,inf);
}
}
S=1;T=cal(n,m,r)+2;
while(bfs())
{
for(int i=1;i<=cal(n,m,r)+2;i++)
cur[i]=f[i];
ans+=dfs(S,inf);
}
printf("%d\n",ans);
}