[BZOJ3144][HNOI2013]切糕（最小割）

题目：

我是超链接

题解：

这么经典的题目我现在才做？
首先和s，t相连的就是每个格子都要有40个点，起码要在40个点里选择一个（最小割）最小代价
那么这个限制d怎么办呢？相邻的格子如果某一个选择了h，那么另一个选择的范围只能是[h-d,h+d]，那么相连上相邻格子的h-d，h+d点，那么必须要选择下一个格子的[h-d,h+d]，不然S还是能到T。

因为x和y相连，y也会和x相连，那么只需要连接相邻部分的z，z-d就行了，另一个部分会在y那里连
所以这道题我应该是做过？

代码：

#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#define INF 1e9
using namespace std;
const int N=64005;
const int M=N*20;
int tot=-1,point[N],nxt[M],v[M],remind[M],a[45][45][45],f[45][45][45],dis[N],cur[N];
int C[4][2]={{1,0},{-1,0},{0,1},{0,-1}};
void addline(int x,int y,int l)
{
    ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; remind[tot]=l;
    ++tot; nxt[tot]=point[y]; point[y]=tot; v[tot]=x; remind[tot]=0;
}
bool bfs(int s,int t)
{
    for (int i=s;i<=t;i++) cur[i]=point[i];
    queue<int> q;
    q.push(s); memset(dis,0x7f,sizeof(dis));
    dis[s]=0;
    while (!q.empty())
    {
        int now=q.front(); q.pop();
        for (int i=point[now];i!=-1;i=nxt[i])
            if (remind[i]&&dis[v[i]]>dis[now]+1)
            {
                dis[v[i]]=dis[now]+1;
                q.push(v[i]);
            }
    }
    return dis[t]<INF;
}
int dfs(int now,int t,int limit)
{
    if (now==t || !limit) return limit;
    int flow=0,f;
    for (int i=cur[now];i!=-1;i=nxt[i])
    {
        cur[now]=i;
        if (dis[v[i]]==dis[now]+1 && remind[i])
        {
            f=dfs(v[i],t,min(limit,remind[i]));
            remind[i]-=f; remind[i^1]+=f;
            flow+=f; limit-=f;
            if (!limit) return flow;
        }
    }
    return flow;
}
int dinic(int s,int t)
{
    int ans=0;
    while (bfs(s,t)) ans+=dfs(s,t,INF);
    return ans;
}
int main()
{
    memset(point,-1,sizeof(point));
    int n,m,r,D;scanf("%d%d%d%d",&n,&m,&r,&D);
    int num=0;
    for (int i=1;i<=r;i++)
       for (int j=1;j<=n;j++)
            for (int k=1;k<=m;k++) scanf("%d",&f[j][k][i]),a[j][k][i]=++num;
    int s=0,t=n*m*r+1;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            addline(s,a[i][j][1],f[i][j][1]);
            for (int k=2;k<=r;k++)
                addline(a[i][j][k-1],a[i][j][k],f[i][j][k]);
            addline(a[i][j][r],t,INF);
        }
    for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) for (int k=D;k<=r;k++)
        for (int t=0;t<4;t++)
        {
            int xx=i+C[t][0],yy=j+C[t][1];
            if (xx<=0 || yy<=0 || xx>n || yy>m) continue;
            addline(a[i][j][k],a[xx][yy][k-D],INF);
        }
    printf("%d",dinic(s,t));
}