【BZOJ1001】狼抓兔子，平面图转对偶图（从最小割到最短路）

Time:2016.08.18
Author:xiaoyimi
转载注明出处谢谢

---

传送门
思路：
显然这个题是求这个图的最小割

在这个n*m个点,2*(n*(m-1)+(n-1)m+(n-1)(m-1))条边的图中直接跑最大流的结果是

那怎么办呢

关于平面图与对偶图，有一个比较牛逼的性质

如果让对偶图的每条边的权值对应平面图中的每条权值，那么对偶图中的最短路就是平面图中的最小割

啥，你不知道什么是平面图，什么是对偶图……

没事，之前的时候我也不知道……

学习一个

总之这个性质非常牛逼就对了

你跑一遍SPFA最短路的时间总比你跑dinic/isap的时间要少的多了吧

注意：
建图略显蛋疼……
discuss里提示还要判断n=1/m=1的情况
代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#define M 1000002
using namespace std;
int n,m,s,t,cnt,tot;
int dis[M<<1],first[M<<1];
int a[1001][1001],b[1001][1001],c[1001][1001],id[1001][1001][2];
bool vis[M<<1];
queue<int> q;
struct edge{
    int v,w,next;
}e[M*6];
int in()
{
    int t=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') t=(t<<1)+(t<<3)+ch-48,ch=getchar();
    return t;
}
void add(int x,int y,int z)
{
    e[++tot]=(edge){y,z,first[x]};first[x]=tot;
    e[++tot]=(edge){x,z,first[y]};first[y]=tot;
}
main()
{
    n=in();m=in();
    t=n*m+1;
    for (int i=1;i<=n;++i)
        for (int j=1;j<m;++j)
            a[i][j]=in();
    for (int i=1;i<n;++i)
        for (int j=1;j<=m;++j)
            b[i][j]=in();
    if (n==1)
    {
        a[0][0]=1<<30;
        for (int i=1;i<m;i++) a[0][0]=min(a[0][0],a[1][i]);
        printf("%d\n",a[0][0]);
        return 0;
    }
    if (m==1)
    {
        b[0][0]=1<<30;
        for (int i=1;i<n;i++) b[0][0]=min(b[0][0],b[i][1]);
        printf("%d\n",b[0][0]);
        return 0;
    }
    for (int i=1;i<n;++i)
        for (int j=1;j<m;++j)
            c[i][j]=in(),
            id[i][j][0]=++cnt,
            id[i][j][1]=++cnt,
            add(id[i][j][0],id[i][j][1],c[i][j]);   
    for (int i=1;i<=n-1;++i)
        add(s,id[i][1][0],b[i][1]),
        add(id[i][m-1][1],t,b[i][m]);
    for (int i=1;i<n;++i)
        for (int j=1;j<m-1;++j)
            add(id[i][j][1],id[i][j+1][0],b[i][j+1]);
    for (int i=1;i<=m-1;++i)
        add(s,id[n-1][i][0],a[n][i]),
        add(id[1][i][1],t,a[1][i]);
    for (int i=1;i<n-1;++i)
        for (int j=1;j<m;++j)
            add(id[i][j][0],id[i+1][j][1],a[i+1][j]);
    q.push(s);
    memset(dis,63,sizeof(dis));
    dis[s]=0;
    for (;!q.empty();q.pop())
    {
        int x=q.front();
        vis[x]=1;
        for (int i=first[x];i;i=e[i].next)
            if (e[i].w+dis[x]<dis[e[i].v])
            {
                dis[e[i].v]=dis[x]+e[i].w;
                if (!vis[e[i].v])
                    vis[e[i].v]=1,
                    q.push(e[i].v);
            }
        vis[x]=0;
    }
    printf("%d\n",dis[t]);
}