bzoj3144 [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

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

分析：

tip

stl的queu慢出天际，

以后一定要手写队列

空间不要开太大

一开始的内存4Wkb，慢出天际，
学长的只有3W-，改了一下内存，立马时间砍一半

这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int INF=0x33333333;
const int N=1000010;
const int M=60010;
struct node{
    int x,y,v,nxt;
};
node way[N<<1];
int st[M],tot=-1,cur[M],deep[M],P,Q,R,z[41][41][41],s,t,D;
int q[M],tou,wei;

int get(int x,int y,int z){return (z-1)*P*Q+(x-1)*Q+y;}

void add(int u,int w,int z)
{
    tot++;
    way[tot].x=u;way[tot].y=w;way[tot].v=z;way[tot].nxt=st[u];st[u]=tot;
    tot++;
    way[tot].x=w;way[tot].y=u;way[tot].v=0;way[tot].nxt=st[w];st[w]=tot;
}

int bfs(int s,int t)
{
    memset(deep,-1,sizeof(deep));
    tou=wei=0;
    q[++wei]=s;
    deep[s]=1;
    for (int i=s;i<=t;i++) cur[i]=st[i];
    do
    {
        int r=q[++tou];
        for (int i=st[r];i!=-1;i=way[i].nxt)
            if (way[i].v&&deep[way[i].y]==-1)
            {
                deep[way[i].y]=deep[r]+1;
                q[++wei]=way[i].y;
            }
    }while (tou<wei);
    return deep[t]!=-1;
}

int dfs(int now,int t,int limit)
{
    if (!limit||now==t) return limit;
    int i,f,flow=0;
    for (i=cur[now];i!=-1;i=way[i].nxt)
    {
        cur[now]=i;
        if (deep[way[i].y]==deep[now]+1&&way[i].v&&(f=dfs(way[i].y,t,min(limit,way[i].v))))
        {
            flow+=f;
            limit-=f;
            way[i].v-=f;
            way[i^1].v+=f;
            if (!limit) break;
        }
    }
    return flow;
}

int doit()
{
    int ans=0;
    while (bfs(s,t))
        ans+=dfs(s,t,INF);
    return ans;
}

int main()
{
    memset(st,-1,sizeof(st));
    scanf("%d%d%d",&P,&Q,&R);
    scanf("%d",&D);
    s=0;t=P*Q*(R+1)+1;
    for (int k=1;k<=R;k++)  //
        for (int i=1;i<=P;i++)
            for (int j=1;j<=Q;j++)
            {
                scanf("%d",&z[i][j][k]);
            }  
    for (int i=1;i<=P;i++)
        for (int j=1;j<=Q;j++)
        {
            add(s,get(i,j,1),INF);
            add(get(i,j,R+1),t,INF);
        }     
    for (int k=1;k<=R+1;k++)
        for (int i=1;i<=P;i++)
            for (int j=1;j<=Q;j++)
            {
                if (k!=R+1) add(get(i,j,k),get(i,j,k+1),z[i][j][k]);
                if (k>D)
                {
                    int h=k-D;
                    if (i>1) add(get(i,j,k),get(i-1,j,h),INF);
                    if (j>1) add(get(i,j,k),get(i,j-1,h),INF);
                    if (i<P) add(get(i,j,k),get(i+1,j,h),INF);
                    if (j<Q) add(get(i,j,k),get(i,j+1,h),INF);
                }
            }
    printf("%d",doit());
    return 0;
}