BZOJ 3144 网络流最小割 解题报告

3144: [Hnoi2013]切糕

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

[解题报告]
比较详细的题解
http://blog.csdn.net/zarxdy34/article/details/45272055

代码如下:

#include<cstdio>
#include<cstring> 
#include<string>
#include<algorithm>
#include<map>
#include<queue>  
#include<vector>
using namespace std;  
#define inf 0x3f3f3f3f  
#define maxv 52000
#define maxe 520000

int nume=0,head[maxv],e[maxe][3];

void inline adde(int i,int j,int c)  
{  
    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;  
    e[nume++][2]=c; 
    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;  
    e[nume++][2]=0;  
}  

int ss,tt,n,m;  
int vis[maxv],lev[maxv]; 

bool bfs()  
{  
    for(int i=0;i<maxv;i++)  
    vis[i]=lev[i]=0;  
    queue<int>q;  
    q.push(ss);  
    vis[ss]=1;  
    while(!q.empty())  
    {  
        int cur=q.front();  
        q.pop();  
        for(int i=head[cur];i!=-1;i=e[i][1])  
        {  
            int v=e[i][0];  
            if(!vis[v]&&e[i][2]>0)  
            {  
                lev[v]=lev[cur]+1;  
                vis[v]=1;  
                q.push(v);  
            }  
        }  
    }  
    return vis[tt];  
}  
int dfs(int u,int minf)  
{  
    if(u==tt||minf==0) return minf;  
    int sumf=0,f;  
    for(int i=head[u];i!=-1&&minf;i=e[i][1])  
    {  
        int v=e[i][0];  
        if(lev[v]==lev[u]+1&&e[i][2]>0)  
        {  
            f=dfs(v,minf<e[i][2]?minf:e[i][2]);  
            e[i][2]-=f;e[i^1][2]+=f;  
            sumf+=f;minf-=f;  
        }  
    }  
    if(!sumf) lev[u]=-1;  
    return sumf;  
}  
int Dinic()  
{  
    int sum=0;  
    while(bfs()) sum+=dfs(ss,inf);  
    return sum;  
} 

int p,q,r,d;
int mp[50][50][50];
int mov[4][2]={{0,1},{1,0},{0,-1},{-1,0}};

int main()
{
    nume=0;
    memset(head,-1,sizeof(head));
    scanf("%d%d%d%d",&p,&q,&r,&d);
    int cnt=0;ss=++cnt;
    for(int i=1;i<=r+1;++i)
    for(int j=1;j<=p;++j)
    for(int k=1;k<=q;++k) mp[i][j][k]=++cnt;
    tt=++cnt;
    for(int i=1;i<=r;++i)
    for(int j=1;j<=p;++j)
    for(int k=1;k<=q;++k) 
    {
        int x;scanf("%d",&x);
        adde(mp[i][j][k],mp[i+1][j][k],x);
    }
    for(int i=1;i<=p;++i)
    for(int j=1;j<=q;++j)
    {
        adde(ss,mp[1][i][j],inf);
        adde(mp[r+1][i][j],tt,inf);
    }
    for (int i=d+1;i<=r+1;i++)  
    for (int j=1;j<=p;j++)  
    for (int k=1;k<=q;k++)  
    for (int mv=0;mv<4;mv++)  
    if (mp[i][j+mov[mv][0]][k+mov[mv][1]])  
    adde(mp[i][j][k],mp[i-d][j+mov[mv][0]][k+mov[mv][1]],inf);  
    printf("%d\n",Dinic());
    return 0;
}