zoukankan      html  css  js  c++  java
  • BZOJ 1412: [ZJOI2009]狼和羊的故事【网络流】

    Description

    “狼爱上羊啊爱的疯狂,谁让他们真爱了一场;狼爱上羊啊并不荒唐,他们说有爱就有方向......” Orez听到这首歌,心想:狼和羊如此和谐,为什么不尝试羊狼合养呢?说干就干! Orez的羊狼圈可以看作一个n*m个矩阵格子,这个矩阵的边缘已经装上了篱笆。可是Drake很快发现狼再怎么也是狼,它们总是对羊垂涎三尺,那首歌只不过是一个动人的传说而已。所以Orez决定在羊狼圈中再加入一些篱笆,还是要将羊狼分开来养。 通过仔细观察,Orez发现狼和羊都有属于自己领地,若狼和羊们不能呆在自己的领地,那它们就会变得非常暴躁,不利于他们的成长。 Orez想要添加篱笆的尽可能的短。当然这个篱笆首先得保证不能改变狼羊的所属领地,再就是篱笆必须修筑完整,也就是说必须修建在单位格子的边界上并且不能只修建一部分。

    Input

    文件的第一行包含两个整数n和m。接下来n行每行m个整数,1表示该格子属于狼的领地,2表示属于羊的领地,0表示该格子不是任何一只动物的领地。

    Output

    文件中仅包含一个整数ans,代表篱笆的最短长度。

    Sample Input

    2 2
    2 2 
    1 1 

    Sample Output

    2

    数据范围
    10%的数据  n,m≤3
    30%的数据   n,m≤20
    100%的数据  n,m≤100

    思路:一开始没明白题意。。。。以为在1或者2的边上围一圈就好了TUT

    明白过来题意后发现这货就是个最小割嘛,这样就使源点和汇点不连通了

    //bzoj1412

    #include <stdio.h>

    #include <string.h>

    #include <algorithm>

    #include <queue>

    #define maxn 300000

    #define inf 0x3f3f3f3f

    using namespace std;

    const int dx[10]={0,0,0,1,-1};

    const int dy[10]={0,1,-1,0,0};

    int head[maxn],next[maxn],point[maxn],now;

    int flow[maxn],dist[maxn];

    void add(int x,int y,int v)

    {

        next[++now]=head[x];

        head[x]=now;

        flow[now]=v;

        point[now]=y;

        next[++now]=head[y];

        head[y]=now;

        flow[now]=0;

        point[now]=x;

    }

    int bfs(int s,int t)

    {

        queue<int>q;

        q.push(s);

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

        dist[s]=0;

        while(!q.empty())

        {

            int k=q.front();

            q.pop();

            for(int i=head[k];i;i=next[i])if(flow[i]&&dist[point[i]]==-1)

            {

                int u=point[i];

                dist[u]=dist[k]+1;

                q.push(u);

            }

        }

        return dist[t]!=-1;

    }

    int dfs(int s,int d,int t)

    {

        if(s==t)return d;

        int res=0;

        for(int i=head[s];i&&res<d;i=next[i])

        {

            int u=point[i];

            if(flow[i]&& dist[u]==dist[s]+1)

            {

                int dd=dfs(u,min(flow[i],d-res),t);

                if(dd)

                {

                    flow[i]-=dd;

                    flow[((i-1)^1)+1]+=d;

                    res+=dd;

                }

            }

        }

        if(res==0)dist[s]=-1;

        return res;

    }

    int main()

    {

        int n,m,s=maxn-2,t=maxn-3,x,ans=0;

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

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

        {

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

            {

                scanf("%d",&x);

                if(x==1)add(s,i*m+j,inf);

                if(x==2)add(i*m+j,t,inf);

                for(int k=1;k<=4;k++)

                {

                    int x=i+dx[k],y=j+dy[k];

                    if(1<=x&&x<=n&&1<=y&&y<=m)

                    {

                        add(i*m+j,x*m+y,1);

                    }

                }

            }

        }

        while(bfs(s,t))ans+=dfs(s,inf,t);

        printf("%d ",ans);

        return 0;

    }

  • 相关阅读:
    2015.2.27 UltraEdit中显示XML结构
    2015.1.31 DataGridView自动滚动到某行
    2015.1.15 利用函数实现将一行记录拆分成多行记录 (多年想要的效果)
    2015.1.15 利用Oracle函数返回表结果 重大技术进步!
    2015.1.15 利用Oracle函数插入表结构 Bulk collect into 不用循环,简洁高效
    2015.1.8 Left join 左连接
    2015.1.10 解决DataGridView SelectionChanged事件自动触发问题
    delphi 遍历窗口
    delphi 访问 protected 属性 哈哈
    clientdataset 读取excel 如果excel 文件不存在的时候 相应的gird 会不显示数据, 鼠标掠过 gird 格子 才会显示数据。 这是一个bug 哈哈
  • 原文地址:https://www.cnblogs.com/philippica/p/4094333.html
Copyright © 2011-2022 走看看