BZOJ 3144: [Hnoi2013]切糕

3144: [Hnoi2013]切糕

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1495  Solved: 819
[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

Sample Output

6

HINT

最佳切面的f为f(1,1)=f(2,1)=2,f(1,2)=f(2,2)=1

Source

分析：

限制距离模型：

可以把题意看成有p*q个网格，每个网格有r个点，要在每个格子中选一个点，并且相邻的格子选择的点距离不超过d，求最小代价...

考虑如果没有距离限制怎么建图...把每个格子拆成r个点，串成一条链和ST相连，求最小割就是答案...

现在有了距离限制...怎么办？？最常用的限制方法就是添加容量为+∞的边...

我们把(i,j,k)向(i',j',k-d)(相邻的格子)连边...这个正确性画一下图YY一下就好...

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
//by NeighThorn
#define inf 0x3f3f3f3f
using namespace std;

const int maxn=100000+5,maxm=1000000+5;

int n,m,h,d,S,T,cnt,hd[maxn],fl[maxm],to[maxm],nxt[maxm],pos[maxn];

int mv[4][2]={0,1,0,-1,1,0,-1,0};

inline bool bfs(void){
    memset(pos,-1,sizeof(pos));
    int head=0,tail=0,q[maxn];
    q[0]=S,pos[S]=0;
    while(head<=tail){
        int top=q[head++];
        for(int i=hd[top];i!=-1;i=nxt[i])
            if(pos[to[i]]==-1&&fl[i])
                pos[to[i]]=pos[top]+1,q[++tail]=to[i];
    }
    return pos[T]!=-1;
}

inline int find(int v,int f){
    if(v==T)
        return f;
    int res=0,t;
    for(int i=hd[v];i!=-1&&f>res;i=nxt[i])
        if(pos[to[i]]==pos[v]+1&&fl[i])
            t=find(to[i],min(fl[i],f-res)),res+=t,fl[i]-=t,fl[i^1]+=t;
    if(!res)
        pos[v]=-1;
    return res;
}

inline int dinic(void){
    int res=0,t;
    while(bfs())
        while(t=find(S,inf))
            res+=t;
    return res;
}

inline void add(int s,int x,int y){
    fl[cnt]=s;to[cnt]=y;nxt[cnt]=hd[x];hd[x]=cnt++;
    fl[cnt]=0;to[cnt]=x;nxt[cnt]=hd[y];hd[y]=cnt++;
}

signed main(void){
    // freopen("in.txt","r",stdin);
    memset(hd,-1,sizeof(hd));
    scanf("%d%d%d%d",&n,&m,&h,&d);
    S=0,T=n*m*h+1;
    for(int k=1,x;k<=h;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++){
                scanf("%d",&x);
                if(k==1)
                    add(x,S,((i-1)*m+j-1)*h+k);
                else
                    add(x,((i-1)*m+j-1)*h+k-1,((i-1)*m+j-1)*h+k);
                if(k==h)
                    add(inf,((i-1)*m+j-1)*h+k,T);
            }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            for(int k=d+1;k<=h;k++)
                for(int t=0;t<4;t++){
                    int x=i+mv[t][0],y=j+mv[t][1];
                    if(x>=1&&x<=n&&y>=1&&y<=m)
                        add(inf,((i-1)*m+j-1)*h+k,((x-1)*m+y-1)*h+k-d);
                }
    printf("%d
",dinic());
    return 0;   
}//Cap ou pas cap. Cap.

---

By NeighThorn