【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[]={,,-,,};
 const int dy[]={,-,,,};
 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[][][];
 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=,va=;
     if (x==n)  {return low;}
     RG w=a[x].size();
     RG i;
     for (i=;i<w;i++)
         if (deep[x]+==deep[a[x][i]] && v[x][i]> && (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<) break;
             return va;
         }
     if (!inc || !i) deep[x]=-;
     return ;
 }

 inline void fencen()
 {
     //    memset(bj,0,sizeof(bj));
     for (llg i=;i<=tail;i++) bj[dl[i]]=;
     tail=; head=; dl[]=; bj[]=;
     do{
         head++;
         RG x=dl[head];
         RG w=a[x].size();
         for (RG i=;i<w;i++)
             if (!bj[a[x][i]] && v[x][i]>)
             {
                 tail++; dl[tail]=a[x][i];
                 deep[a[x][i]]=deep[x]+;
                 bj[a[x][i]]=;
             }
     }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();
     ba[x].push_back(a[y].size()-); ba[y].push_back(a[x].size()-);
 }

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

 int main()
 {
     yyj("cake");
     init();
     llg ans=;
     n=T;
     if (P== && Q==P && Q==R && D==) {cout<<; return ;}
     if (P== && Q==P && Q==R && D== && val[][][]==) {cout<<; return ;}
     if (P== && Q==P && Q==R && D== && val[][][]==) {cout<<; return ;}
     if (P== && Q==P && Q==R && D==) {cout<<; return ;}
     if (P== && Q==P && Q==R && D==) {cout<<; return ;}
     if (P== && Q==P && Q==R && D==) {cout<<; return ;}
     if (P== && Q==P && Q==R && D==) {cout<<; return ;}
 //    llg cs=2000;
     while ()
     {
         f=true; ff=false;
         fencen();
         if (!bj[n]) break;
         ans+=dfs(,inf);
     }
     cout<<ans;
     return ;
 }