hdu 3046 Pleasant sheep and big big wolf (最小割 dinic )

题目链接:哆啦A梦传送门

题意:给出n*m方格,1表示羊格,2表示狼格,0表示空格。现在我们要将羊群与狼群分隔开,我们可以在相邻格的边界上建上栅栏,问:最少需要多少栅栏,就能把羊群与狼群分隔开?

 

题解:最小割+dinic。

我们将源点到狼点建边,权值为INF,羊点到汇点建边,权值为INF。若两格有公共边就连一条容量为1的无向边。

因为题目不能让羊群和狼群相通,故要我们求最小割,使得拿掉这些边(此题就是建上栅栏),就能让羊群与狼群不能相通。

这里与 hdu1565有所不同,这里是源点到狼点(羊点到汇点)权值为INF,而hdu1565是源点到奇数点权值为所在点的值。为什么呢?因为我们要把狼群与羊群给分割开,也就是要能割掉它们的边,故我们设源点到狼点的权值为INF,中间点之间的权值不为INF。

#include<bits/stdc++.h>

using namespace std;

const int N=400000;
const int INF=0x3f3f3f3f;

struct edge{
    int from,to,w,next;
}e[N];
int cnt;

int a[300][300],dep[N],head[N];
queue<int> que;

int dir[4][2]={{0,1},{0,-1},{-1,0},{1,0}};

void init(){
    memset(head,-1,sizeof(head));
    cnt=0;
}
void add(int u,int v,int w)
{
    e[cnt].from=u;
    e[cnt].to=v;
    e[cnt].w=w;
    e[cnt].next=head[u];
    head[u]=cnt++;
}

int bfs(int s,int t)
{
    while(!que.empty()) que.pop();

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

    dep[s]=0;
    que.push(s);

    while(!que.empty()){
        int u=que.front();
        que.pop();
        if(u==t) return 1;
            for(int i=head[u];i!=-1;i=e[i].next){
                int w=e[i].w;
                int v=e[i].to;
            if(w&&dep[v]==-1){
                dep[v]=dep[u]+1;
                que.push(v);
            }
        }

    }
    return 0;
}

int dfs(int u,int mi,int t){
    if(u==t) return mi;
    int tmp;
    int ret=0;
//    for(int v=1;v<=t;v++){
        for(int i=head[u];i!=-1;i=e[i].next){
            int w=e[i].w;
            int v=e[i].to;
        if(w&&dep[v]==dep[u]+1){
           tmp=dfs(v,min(mi-ret,w),t) ;
           e[i].w-=tmp;
            e[i^1].w+=tmp;
            ret+=tmp;
            if(ret==mi) return ret;
        }
    }
    return ret;
}

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

    }
    return ans;
}
int main()
{

    int n,m;
    int T=0;
    while(~scanf("%d%d",&n,&m))
    {

        init();
        int sum=0;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
            scanf("%d",&a[i][j]),sum+=a[i][j];

        int s=0,t=n*m+1;
        int cur;
        for(int i=1;i<=n;i++) ///建图
        {
            for(int j=1;j<=m;j++){
                    cur=(i-1)*m+j;
                if(a[i][j]==1){
                    add(cur,t,INF);///给羊格到汇点建边
                    add(t,cur,0); ///反向边
                }
                if(a[i][j]==2){ ///给源点到狼格建边
                    add(s,cur,INF);
                    add(cur,s,0);///反向边
                }

                for(int k=0;k<4;k++){ ///给相邻格建边
                    int xx=i+dir[k][0];
                    int yy=j+dir[k][1];
                    if(xx>0&&xx<=n&&yy>0&&yy<=m){
                        add(cur,(xx-1)*m+yy,1);
                        add((xx-1)*m+yy,cur,0);///反向边
                    }
                }
            }
        }
        printf("Case %d:\n",++T);
        printf("%d\n",dinic(s,t));
    }
    return 0;
}


wAAACH5BAEKAAAALAAAAAABAAEAAAICRAEAOw==

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值