狼抓兔子

不考虑时间,裸的dinic。注意无向图反向弧直接w。

我们发现图是平面图,平面图最小割——对偶图最短路。
对偶图的构造——
对偶图的性质:
1. G*的点数等于G的面数,G*的边数等于G的边数
2. G*中的环对应G中的割
∴可应用此性质求最小割!!!
应用对偶图求最小割
步骤:
1. 在平面图中的起点和终点之间加一条边
2. 将平面图(即原来的图)转化为对偶图,记起点和终点构成的边组成的面为源点,无界面为汇点,每条边的权等于原图中对应边的权
3. 删去源点和汇点之间的边
4. 在对偶图上求最短路即为原图的最大流(或最小割)

#include<bits/stdc++.h>
using namespace std;

const int MAXN=1e7+5;
const int INF=1e9+7;

struct edge{
    int to,next,w;
}e[MAXN<<1];

int n,m,s,t;

int head[MAXN],cur[MAXN],cnt=1;

inline void add(int u,int v,int w){
    e[++cnt]=(edge){v,head[u],w},head[u]=cnt;
    e[++cnt]=(edge){u,head[v],w},head[v]=cnt;
}

int dep[MAXN];
queue<int>q;
bool bfs(int x){
    memset(dep,0,sizeof(dep));
    q.push(x);dep[x]=1;
    while(q.size()){
        int u=q.front();q.pop();
        for(int i=head[u];i;i=e[i].next){
            int v=e[i].to,w=e[i].w;
            if(e[i].w&&!dep[v]){
                dep[v]=dep[u]+1;
                q.push(v);
            }
        }   
    }
    if(!dep[t])return 0;
    return 1;
}

int dfs(int u,int flow){
    if(flow==0||u==t)return flow;
    for(int &i=cur[u];i;i=e[i].next){
        int v=e[i].to,w=e[i].w;
        if(dep[v]==dep[u]+1&&w){
            int tem=dfs(v,min(flow,w));
            if(tem){
                e[i].w-=tem;
                e[i^1].w+=tem;
                return tem;
            }
        }
    }
    return 0;
} 

int dinic(){
    int ans=0;
    while(bfs(s)){
        for(int i=s;i<=t;i++)cur[i]=head[i];
        while(int d=dfs(s,INF))ans+=d;
    }
    return ans;
} 

int main(){
    scanf("%d%d",&n,&m);
    s=0,t=n*m+1;
    for(int i=1;i<=n;i++){
        int tem;
        for(int j=1;j<m;j++){
            scanf("%d",&tem);
            int dian=(i-1)*m+j;
            add(dian,dian+1,tem);
        }
    }
    for(int i=1;i<n;i++){
        int tem;
        for(int j=1;j<=m;j++){  
            scanf("%d",&tem);
            int dian=(i-1)*m+j;
            add(dian,dian+m,tem);
        }
    }
    for(int i=1;i<n;i++){
        int tem;
        for(int j=1;j<m;j++){
            scanf("%d",&tem);
            int dian=(i-1)*m+j;
            add(dian,dian+m+1,tem);
        }
    }
    add(s,1,INF);
    add(n*m,t,INF);
    printf("%d\n",dinic());
    return 0;
} 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值