BZOJ 1001 狼抓兔子 st平面图最大流 转化为 对偶图求最短路

/*
2018/10/16
BZOJ 1001
st平面图最大流，转化为对偶图跑最短路
顾业鸣
注意n==1||m==1的特判
*/
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<set>
using namespace std;
typedef long long ll;
struct Edge
{
    int u,v,w,nxt;
}edge[7000000];
int cnt,head[2000000],vis[2000000];
ll dis[2000000];
void addedge(int u,int v,int w)
{
    edge[cnt].u=u;
    edge[cnt].v=v;
    edge[cnt].w=w;
    edge[cnt].nxt=head[u];
    head[u]=cnt++;
}
ll dijkstra(int s,int t)
{
    dis[s]=0;
    set<pair<ll,int> > st;
    st.insert(make_pair(0ll,s));
    while(!st.empty()){
        int x=st.begin()->second;
        ll w=st.begin()->first;
        st.erase(st.begin());
        if(vis[x]==1) continue;
        //printf("%d %lld\n",x,w);
        vis[x]=1;
        for(int i=head[x];i!=-1;i=edge[i].nxt){
            if(vis[edge[i].v]==0&&dis[x]+edge[i].w<dis[edge[i].v]){
                dis[edge[i].v]=dis[x]+edge[i].w;
                st.insert(make_pair(dis[edge[i].v],edge[i].v));
            }
        }
    }
    return dis[t];
}
int main(void)
{
    int n,m,u,v,w;
    for(int i=0;i<2000000;i++){
        head[i]=-1;
        dis[i]=(1ll<<63)-1;
    }
    scanf("%d%d",&n,&m);
    if(n==1||m==1){
        int ans=100000;
        int num=max(n,m);
        for(int i=1;i<num;i++){
            scanf("%d",&w);
            ans=min(ans,w);
        }
        return 0*printf("%d\n",ans);
    }
    for(int i=0;i<n;i++){
        for(int j=1;j<m;j++){
            scanf("%d",&w);
            if(i==0){
                addedge(0,j,w);
                addedge(j,0,w);
                //printf("##1 %d %d\n",0,j);
            }
            else if(i==n-1){
                addedge(2*(i-1)*(m-1)+m-1+j,(n-1)*(m-1)*2+1,w);
                addedge((n-1)*(m-1)*2+1,2*(i-1)*(m-1)+m-1+j,w);
                //printf("##2 %d %d\n",2*(i-1)*(m-1)+m-1+j,(n-1)*(m-1)*2+1);
            }
            else{
                addedge(2*(i-1)*(m-1)+m-1+j,i*2*(m-1)+j,w);
                addedge(i*2*(m-1)+j,2*(i-1)*(m-1)+m-1+j,w);
                //printf("##3 %d %d\n",2*(i-1)*(m-1)+m-1+j,i*2*(m-1)+j);
            }
        }
    }
    for(int i=0;i<n-1;i++){
        for(int j=1;j<=m;j++){
            scanf("%d",&w);
            if(j==1){
                addedge(i*2*(m-1)+m-1+j,(n-1)*(m-1)*2+1,w);
                addedge((n-1)*(m-1)*2+1,i*2*(m-1)+m-1+j,w);
                //printf("$$1 %d %d\n",i*2*(m-1)+m-1+j,(n-1)*(m-1)*2+1);
            }
            else if(j==m){
                addedge(i*2*(m-1)+j-1,0,w);
                addedge(0,i*2*(m-1)+j-1,w);
                //printf("$$2 %d %d\n",i*2*(m-1)+j-1,0);
            }
            else{
                addedge(i*2*(m-1)+m-1+j,i*2*(m-1)+j-1,w);
                addedge(i*2*(m-1)+j-1,i*2*(m-1)+m-1+j,w);
                //printf("$$3 %d %d\n",i*2*(m-1)+m-1+j,i*2*(m-1)+j-1);
            }
        }
    }
    for(int i=0;i<n-1;i++){
        for(int j=1;j<m;j++){
            scanf("%d",&w);
            addedge(i*2*(m-1)+j,i*2*(m-1)+m-1+j,w);
            addedge(i*2*(m-1)+m-1+j,i*2*(m-1)+j,w);
            //printf("!!  %d %d\n",i*2*(m-1)+j,i*2*(m-1)+m-1+j);
        }
    }
    printf("%lld\n",dijkstra(0,(n-1)*(m-1)*2+1));
    return 0;
}