[BZOJ 1412][ZJOI2009]狼和羊的故事（最小割）

Description

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

Solution

T了两个点，好气哦…先存代码

立一个要学ISAP的flag，题解回来再补

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
int s,t,n,m,a[100][100];
int head[10005],level[10005],cnt=0;
int dx[4]={1,-1,0,0},dy[4]={0,0,1,-1};
int read()
{
    int x=0,f=1;char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-')f=-1;c=getchar();
    }
    while(c>='0'&&c<='9'){
        x=x*10+c-'0';c=getchar();
    }
    return x*f;
}
struct Node
{
    int next,to,cap;
}Edges[500005];
void addedge(int u,int v,int c)
{
    Edges[cnt].next=head[u];
    head[u]=cnt;
    Edges[cnt].to=v;
    Edges[cnt++].cap=c;
}
void insert(int u,int v,int c)
{
    addedge(u,v,c);
    addedge(v,u,0);
}
int pos(int x,int y){return (x-1)*m+y;}
queue<int>q;
bool bfs()
{
    memset(level,-1,sizeof(level));
    q.push(s);level[s]=0;
    while(!q.empty())
    {
        int u=q.front();q.pop();
        for(int i=head[u];~i;i=Edges[i].next)
        {
            int v=Edges[i].to;
            if(level[v]==-1&&Edges[i].cap)
            level[v]=level[u]+1,q.push(v);
        }
    }
    if(level[t]==-1)return false;
    return true;
} 
int dfs(int u,int f)
{
    if(u==t)return f;
    int flow=0,d;
    for(int i=head[u];~i&&flow<f;i=Edges[i].next)
    {
        int v=Edges[i].to;
        if(level[v]==level[u]+1&&Edges[i].cap)
        {
            d=dfs(v,min(f-flow,Edges[i].cap));
            flow+=d;
            Edges[i].cap-=d;
            Edges[i^1].cap+=d;
        }
    }
    if(!flow)level[u]=-1;
    return flow;
}
int dinic()
{
    int res=0,d;
    while(bfs())
    {
        while(d=dfs(s,INF))
        res+=d;
    }
    return res;
}
int main()
{
    memset(head,-1,sizeof(head));
    n=read(),m=read();
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    a[i][j]=read();
    s=0,t=n*m+1;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    {
        if(a[i][j]!=2)
        for(int k=0;k<4;k++)
        {
            int x=i+dx[k],y=j+dy[k];
            if(x>0&&x<=n&&y>0&&y<=m&&(a[i][j]||a[x][y]!=1))
            insert(pos(i,j),pos(x,y),1);
        }
        if(a[i][j]==1)
        insert(s,pos(i,j),INF);
        else if(a[i][j]==2)
        insert(pos(i,j),t,INF);
    }
    printf("%d
",dinic());
    return 0;
}