bzoj1001 [BeiJing2006]狼抓兔子 对偶图

Description

---

给定一个网格图和边权求(1,1)到(n,m)的最小割

Solution

---

显然直接最大流是会T的，这里需要模型转换一下
详见这篇论文
转成对偶图跑spfa即可
注意只有一行和只有一列的情况要特判一下

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

const int INF=0x3f3f3f3f;
const int N=4000005;
const int E=8000005;

struct edge {int y,w,next;} e[E];

std:: queue <int> que;

int dis[N];
int ls[N],edCnt;

bool vis[N];

void add_edge(int x,int y,int w) {
    e[++edCnt]=(edge) {y,w,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge) {x,w,ls[y]}; ls[y]=edCnt;
}

void spfa(int st) {
    fill(dis,63); dis[st]=0;
    que.push(st);
    while (!que.empty()) {
        int now=que.front(); que.pop();
        for (int i=ls[now];i;i=e[i].next) {
            if (dis[now]+e[i].w>=dis[e[i].y]) continue;
            dis[e[i].y]=dis[now]+e[i].w;
            if (vis[e[i].y]) continue;
            vis[e[i].y]=1;
            que.push(e[i].y);
        }
        vis[now]=0;
    }
}

void build(int n,int m) {
    int s=0,t=(n-1)*(m-1)*2+1;
    rep(j,1,m-1) {
        int w; scanf("%d",&w);
        add_edge(j*2,t,w);
    }
    rep(i,2,n-1) rep(j,1,m-1) {
        int w; scanf("%d",&w);
        add_edge(2*((i-1)*(m-1)+j),2*((i-2)*(m-1)+j)-1,w);
    }
    rep(j,1,m-1) {
        int w; scanf("%d",&w);
        add_edge(s,2*((n-2)*(m-1)+j)-1,w);
    }
    rep(i,1,n-1) {
        int w; scanf("%d",&w);
        add_edge(s,2*((i-1)*(m-1))+1,w);
        rep(j,2,m-1) {
            scanf("%d",&w);
            add_edge(2*((i-1)*(m-1)+j)-1,2*((i-1)*(m-1)+j-1),w);
        }
        scanf("%d",&w);
        add_edge(2*((i-1)*(m-1)+m-1),t,w);
    }
    rep(i,1,n-1) rep(j,1,m-1) {
        int w; scanf("%d",&w);
        add_edge(2*((i-1)*(m-1)+j)-1,2*((i-1)*(m-1)+j),w);
    }
}

int main(void) {
    int n,m; scanf("%d%d",&n,&m);
    if (n==1&&m==1) {
        puts("0");
        return 0;
    }
    int ans=INF;
    if (n==1) {
        rep(i,1,m-1) {
            int x; scanf("%d",&x);
            ans=std:: min(ans,x);
        }
        printf("%d\n", ans);
        return 0;
    }
    if (m==1) {
        rep(i,1,n-1) {
            int x; scanf("%d",&x);
            ans=std:: min(ans,x);
        }
        printf("%d\n", ans);
        return 0;
    }
    build(n,m);
    spfa(0);
    printf("%d\n", dis[(n-1)*(m-1)*2+1]);
    return 0;
}