天天看點

BZOJ1001: [BeiJing2006]狼抓兔子(洛谷P4001)平面圖最小割->對偶圖最短路

平面圖最小割->對偶圖最短路

BZOJ題目傳送門

洛谷題目傳送門

比較經典的一道平面圖最小割轉對偶圖最短路的題。

推薦這個

我的對偶圖示号方法:

BZOJ1001: [BeiJing2006]狼抓兔子(洛谷P4001)平面圖最小割->對偶圖最短路

注意特判n==1||m==1的情況

代碼:

#include<queue>
#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1005
#define K N*N<<1
#define M N*N*3<<1
using namespace std;
struct edge{ int next,to,dis; }ed[M];
int n,m,k,h[K],s,t,dis[K];
bool f[K];
queue <int> que;//卡了十幾發空間都沒過,就直接搞STL了
//讀優不知道哪裡出了問題,就沒打了
inline void addedge(int x,int y,int z){
    ed[++k].next=h[x],ed[k].to=y,ed[k].dis=z,h[x]=k;
}
inline int spf(){//spfa
    memset(dis,,sizeof(dis));
    memset(f,false,sizeof(f));
    while (!que.empty()) que.pop();
    que.push(s),dis[s]=;
    while (!que.empty()){
        int x=que.front(); que.pop(); f[x]=false;
        for (int i=h[x],v;i;i=ed[i].next)
            if (dis[x]+ed[i].dis<dis[v=ed[i].to]){
                dis[v]=dis[x]+ed[i].dis;
                if (!f[v]) que.push(v),f[v]=true;
            }
    }
    return dis[t];
}
inline int id(int x,int y){//對偶圖的标号
    return (x-)*(m-)+y<<;
}
inline void calc(){
    int x,ans=;
    while (~scanf("%d",&x)) ans=min(ans,x);
    printf("%d\n",ans);
}
int main(){
    scanf("%d%d",&n,&m),s=,t=((n-)*(m-)<<)+;
    if (n==||m==) return calc(),;//一定要特判
    for (int i=;i<=n;i++)
        for (int j=;j<m;j++){
            int x,u=id(i,j),v=id(i-,j)-; scanf("%d",&x);
            u=min(u,t),v=max(v,s),addedge(u,v,x),addedge(v,u,x);
        }
    for (int i=;i<n;i++)
        for (int j=;j<=m;j++){
            int x,u=id(i,j)-,v=u-; scanf("%d",&x);
            if (j==) v=t; if (j==m) u=s;
            addedge(u,v,x),addedge(v,u,x); 
        }
    for (int i=;i<n;i++)
        for (int j=;j<m;j++){
            int x,u=id(i,j),v=u-; scanf("%d",&x);
            addedge(u,v,x),addedge(v,u,x);
        }
    return printf("%d\n",spf()),;
}