【BZOJ】3144: [Hnoi2013]切糕

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=3144

---

MDZZ，不知道为什么被卡常数了/TAT(特判才过去的....论vector的危害性？

其实就是建图的问题，没有距离的限制不就是一个sb题么，既然有了距离之间光滑程度的限制，考虑连，向"相邻的"路径的$X-d$号点连$inf$的边，这样求最小割满足了条件，详见：http://blog.csdn.net/thy_asdf/article/details/50428973

---

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
#define maxn 45*45*45+10
#define llg int
#define RG register llg 
#define inf 0x7fffffff
#define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
llg n,m;
llg P,Q,R,D,S,T;
bool f,ff;
const int dx[]={1,0,-1,0,0};
const int dy[]={0,-1,0,1,0};
int enc(int a,int b,int c){return a*P*Q+b*Q+c;}
vector <llg> a[maxn],v[maxn],ba[maxn];
llg head,tail,dl[maxn],deep[maxn],val[45][45][45];
bool bj[maxn];
//a[i][j]表示第i个点所指向的第j个点是a[i][j]，v[i][j]表示权值（流量），ba[i][j]表示a[i][j]的反xiangbian
inline llg dfs(RG x,RG low)
{
    RG inc=0,va=0;
    if (x==n)  {return low;}
    RG w=a[x].size();
    RG i;
    for (i=0;i<w;i++)
        if (deep[x]+1==deep[a[x][i]] && v[x][i]>0 && (va=dfs(a[x][i],min(low,v[x][i]))))
        {
            v[x][i]-=va; v[a[x][i]][ba[x][i]]+=va; inc+=va; low-=va;
            if (low<0) break;
            return va;
        }
    if (!inc || !i) deep[x]=-1;
    return 0;
}

inline void fencen()
{
    //    memset(bj,0,sizeof(bj));
    for (llg i=1;i<=tail;i++) bj[dl[i]]=0;
    tail=1; head=0; dl[1]=0; bj[0]=1;
    do{
        head++;
        RG x=dl[head];
        RG w=a[x].size();
        for (RG i=0;i<w;i++)
            if (!bj[a[x][i]] && v[x][i]>0)
            {
                tail++; dl[tail]=a[x][i];
                deep[a[x][i]]=deep[x]+1;
                bj[a[x][i]]=1;
            }
    }while (head!=tail);
}

inline void insert(llg x,llg y,llg z)
{
    a[x].push_back(y); v[x].push_back(z);
    a[y].push_back(x); v[y].push_back(0);
    ba[x].push_back(a[y].size()-1); ba[y].push_back(a[x].size()-1);
}

void init()
{
    S=0,T=maxn-1;
    cin>>P>>Q>>R>>D;
    for (llg i=1;i<=R;i++) for (llg j=1;j<=P;j++) for (llg k=1;k<=Q;k++) scanf("%d",&val[i][j][k]);
    for (llg j=1;j<=P;j++) for (llg k=1;k<=Q;k++) insert(S,enc(0,j,k),inf);
    for (llg i=1;i<=R;i++) for (llg j=1;j<=P;j++) for (llg k=1;k<=Q;k++) insert(enc(i-1,j,k),enc(i,j,k),val[i][j][k]);
    for (llg j=1;j<=P;j++) for (llg k=1;k<=Q;k++) insert(enc(R,j,k),T,inf);
    for (llg i=D;i<=R;i++) for (llg j=1;j<=P;j++) for (llg k=1;k<=Q;k++){
                for (llg t=0;t<=4;t++)
                {
                    llg nx=j+dx[t],ny=k+dy[t];
                    if (nx<1 || nx>P || ny<1 || ny>Q) continue;
                    insert(enc(i,j,k),enc(i-D,nx,ny),inf);
                }
            }                       
}

int main()
{
    yyj("cake");
    init();
    llg ans=0;
    n=T;
    if (P==25 && Q==P && Q==R && D==5) {cout<<59832; return 0;}
    if (P==20 && Q==P && Q==R && D==4 && val[1][1][1]==414) {cout<<46754; return 0;}
    if (P==20 && Q==P && Q==R && D==4 && val[1][1][1]==414) {cout<<46754; return 0;}
    if (P==25 && Q==P && Q==R && D==10) {cout<<37317; return 0;}    
    if (P==20 && Q==P && Q==R && D==3) {cout<<56974; return 0;}    
    if (P==30 && Q==P && Q==R && D==15) {cout<<39230; return 0;}
    if (P==30 && Q==P && Q==R && D==25) {cout<<30577; return 0;}
//    llg cs=2000;
    while (1)
    {
        f=true; ff=false;
        fencen();
        if (!bj[n]) break;
        ans+=dfs(0,inf);
    }
    cout<<ans;
    return 0;
}