【bzoj1001】[BeiJing2006]狼抓兔子

最大流

数据有些鬼。。开小容易RE，开大就MLE。。

存双向边

然后直接跑dinic

注意读入有些复杂= =

速度还是跟得上的

Problem: 1001

User: YJY

Language: C++

Result: Accepted

Time:2004 ms

Memory:106744 kb

#include<algorithm>

#include<iostream>

#include<cstdlib>

#include<cstring>

#include<cstdio>

#include<vector>

#include<cmath>

#include<stack>

#include<queue>

using namespace std;

 

#define INF 0x7fffffff

#define MAXN 3000010

 

int n,m;

int ans,k;

 

struct Node

{

    int to,next,v;

}e[MAXN<<1];

 

int head[MAXN];

int dis[MAXN],q[MAXN];

 

void link(int u,int v,int w)

{

    k++;

    e[k].to=v;

    e[k].v=w;

    e[k].next=head[u];

    head[u]=k;

}

 

bool bfs()

{

    queue<int>Q;

    int now;

    memset(dis,-1,sizeof(dis));

    Q.push(1);

    dis[1]=0;

    while (!Q.empty())

    {  

        now=Q.front();

        Q.pop();

        for (int i=head[now];i;i=e[i].next)       

        {

            if (e[i].v && dis[e[i].to]<0)

            {

                Q.push(e[i].to);

                dis[e[i].to]=dis[now]+1;                

            }

        }

    }

    if (dis[n*m]==-1)

        return 0;

    return 1;

}

 

int dfs(int x,int f)

{

    int w,used=0;

    if (x==n*m)

        return f;

    for (int i=head[x];i;i=e[i].next)

    {

        if (e[i].v && dis[e[i].to]==dis[x]+1)

        {

            w=f-used;

            w=dfs(e[i].to,min(w,e[i].v));

            e[i].v-=w;

            e[i+1].v+=w;

            used+=w;

            if (used==f)

                return f;

        }

    }

    if (!used)

        dis[x]=-1;

    return used;

}

 

void dinic()

{

    while (bfs())

        ans+=dfs(1,INF);

}

 

int main()

{

    scanf("%d%d",&n,&m);

    int x;

    for(int i=1;i<=n;i++)

        for(int j=1;j<m;j++)

        {

            scanf("%d",&x);

            link(m*(i-1)+j,m*(i-1)+j+1,x);

            link(m*(i-1)+j+1,m*(i-1)+j,x);

        }

    for(int i=1;i<n;i++)

        for(int j=1;j<=m;j++)

        {

            scanf("%d",&x);

            link(m*(i-1)+j,m*(i)+j,x);

            link(m*(i)+j,m*(i-1)+j,x);

        }

    for(int i=1;i<n;i++)

        for(int j=1;j<m;j++)

        {

            scanf("%d",&x);

            link(m*(i-1)+j,m*(i)+j+1,x);

            link(m*(i)+j+1,m*(i-1)+j,x);

        }

    dinic();

    printf("%d",ans);

    return 0;

}