BZOJ 1412 狼和羊的故事

首先，题目目的就是为了分割狼群和羊群，即建立超级源和超级汇求最小割从而转化成用网络流来处理。

如果没有空地，那么就是简单的二分图最大匹配，但是题中有空地的出现，所以需要在点与点之间建立双向边（不算后向弧），这样才能满足题意（我一开始挂到了这里）

理解透了还是很简单的

代码付上

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <queue>
#include <iostream>
using namespace std;
#define N 105
#define T 10005
#define S 0
int head[N*N],dep[N*N],cnt,n,m;
struct node
{
    int to,next,val;
}e[N*N*N];
inline void add(int x,int y,int z)
{
    e[cnt].to=y;
    e[cnt].val=z;
    e[cnt].next=head[x];
    head[x]=cnt++;
    return ;
}
inline void insert(int x,int y,int z)
{
    add(x,y,z);
    add(y,x,0);
    return ;
}
int bfs()
{
    memset(dep,-1,sizeof(dep));
    queue <int >q;
    q.push(S);
    dep[S]=1;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=head[x];i!=-1;i=e[i].next)
        {
            int to1=e[i].to;
            if(dep[to1]==-1&&e[i].val)
            {
                q.push(to1);
                dep[to1]=dep[x]+1;
            }
        }
    }
    return dep[T]==-1?0:1;
}
int dfs(int x,int maxf)
{
    if(!maxf)return 0;
    if(x==T)return maxf;
    int tflow=maxf,nowf;
    for(int i=head[x];i!=-1;i=e[i].next)
    {
        int to1=e[i].to;
        if(dep[to1]==dep[x]+1&&e[i].val&&dep[x]!=-1)
        {
            nowf=dfs(to1,min(e[i].val,tflow));
            if(!nowf)
            {
                dep[to1]=-1;
                continue;
            }
            tflow-=nowf;
            e[i].val-=nowf,e[i^1].val+=nowf;
            if(!tflow)break;
        }
    }
    dep[x]=-1;
    return maxf-tflow;
}
int map[N][N];
int dx[4]={0,1,0,-1};
int dy[4]={1,0,-1,0};
int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            int x;
            scanf("%d",&x);
            map[i][j]=x;
            if(x==2)insert((i-1)*m+j,T,1<<30);
            if(x==1)insert(S,(i-1)*m+j,1<<30);
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            for(int k=0;k<4;k++)
            {
                int ty=dx[k]+i,tx=dy[k]+j;
                if(ty==0||ty==n+1)continue;
                if(tx==0||tx==m+1)continue;
                int f=(i-1)*m+j,t=(ty-1)*m+tx;
                insert(f,t,1);
            }
        }
    }
    int ans=0;
    while(bfs())
    {
        ans+=dfs(S,1<<30);
    }
    printf("%d
",ans);
    return 0;
}