2132: 圈地计划 最小割

做完了happiness这道题就很简单啦,只是把不同换成了相同,那么只要把相邻的点的A,B交换一下就好啦~

#include<iostream>
#include<cstdio>
#include<cstring>
#define inf 1000000007
using namespace std;
int n,m,cnt=1,sum,S,T;
int head[10005],dis[10005],cur[10005],q[10005];
int next[120005],list[120005],key[120005];
int a[105][105],b[105][105],cc[105][105];
int dx[4]={-1,0,1,0};
int dy[4]={0,-1,0,1};
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void insert(int x,int y,int z)
{
    next[++cnt]=head[x];
    head[x]=cnt;
    list[cnt]=y;
    key[cnt]=z;
}
inline int c(int i,int j)
{
    return (i-1)*m+j;
}
inline bool BFS()
{
    memset(dis,-1,sizeof(dis));
    dis[S]=1; q[1]=S; 
    int t=0,w=1,x;
    while (t<w)
    {
        x=q[++t];
        for (int i=head[x];i;i=next[i])
            if (key[i]&&dis[list[i]]==-1)
                dis[list[i]]=dis[x]+1,q[++w]=list[i];
    }
    return dis[T]!=-1;
}
int find(int x,int flow)
{
    if (x==T) return flow;
    int w,used=0;
    for (int i=cur[x];i;i=next[i])
        if (key[i]&&dis[list[i]]==dis[x]+1)
        {
            w=find(list[i],min(flow-used,key[i]));
            key[i]-=w; key[i^1]+=w; used+=w;
            if (key[i]) cur[x]=i;
            if (used==flow) return flow;
        }
    if (!used) dis[x]=-1;
    return used;
}
inline int dinic()
{
    int tmp=0;
    while (BFS()) 
    {
        for (int i=S;i<=T;i++) cur[i]=head[i];
        tmp+=find(S,inf);
    }
    return tmp;
}
int main()
{
    n=read(); m=read(); S=0; T=n*m+1;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            a[i][j]=read(),sum+=a[i][j];
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            b[i][j]=read(),sum+=b[i][j];
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            cc[i][j]=read();
    for (int i=1;i<=n;i++)
    {
        for (int j=1+(i%2==0);j<=m;j+=2)
        {
            swap(a[i][j],b[i][j]);
            insert(S,c(i,j),a[i][j]);
            insert(c(i,j),S,0);
            insert(c(i,j),T,b[i][j]);
            insert(T,c(i,j),0);
            for (int k=0;k<4;k++)
            {
                int x=i+dx[k],y=j+dy[k];
                if (x<1||x>n||y<1||y>m) continue;
                insert(c(i,j),c(x,y),cc[i][j]+cc[x][y]);
                insert(c(x,y),c(i,j),cc[i][j]+cc[x][y]);
                sum+=cc[i][j]+cc[x][y];
            }
        }
        for (int j=1+(i%2==1);j<=m;j+=2)
        {
            insert(S,c(i,j),a[i][j]);
            insert(c(i,j),S,0);
            insert(c(i,j),T,b[i][j]);
            insert(T,c(i,j),0);
        }
    }
    cout << sum-dinic();
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值