BZOJ 1001 狼抓兔子

第一眼感觉终于看到了一道自己会做的题,不就是最小割裸题嘛简单。然后,发现正解是传说中的最小割转最短路。。

比较玄学。

//Twenty
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
using namespace std;
const int maxn=1000+299;
const int N=2*1000*1000+29;
int s,t,n,m,hl[maxn][maxn],sl[maxn][maxn],xl[maxn][maxn],tot;
int fir[N],nxt[N*4],to[N*4],val[N*4],ecnt,vis[N],dis[N],ans,mi1=1e9,mi2=1e9;
void add(int u,int v,int w) {
    nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v; val[ecnt]=w;
    nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u; val[ecnt]=w;
}
queue<int>que;
void spfa() {
    while(!que.empty()) que.pop();
    memset(vis,0,sizeof(vis));
    memset(dis,127,sizeof(dis));
    dis[s]=0; vis[s]=1;
    que.push(s);
    while(!que.empty()) {
        int now=que.front();
        que.pop(); vis[now]=0;
        for(int i=fir[now];i;i=nxt[i]) {
            if(dis[to[i]]>dis[now]+val[i]) {
                dis[to[i]]=dis[now]+val[i];
                if(!vis[to[i]]) {
                    vis[to[i]]=1;
                    que.push(to[i]);
                } 
             }
           }
    }
}
int main() {
    //freopen(".in","r",stdin);
    //freopen(".out","w",stdout);
    while(scanf("%d%d",&n,&m)==2){
    ecnt=0; memset(fir,0,sizeof(fir));
    for(int i=1;i<=n;i++) for(int j=1;j<m;j++) {
    scanf("%d",&hl[i][j]); 
    mi1=min(mi1,hl[i][j]); 
    }
    for(int i=1;i<n;i++) for(int j=1;j<=m;j++) {
    scanf("%d",&sl[i][j]); 
    mi2=min(mi2,sl[i][j]); 
    }
    for(int i=1;i<n;i++) for(int j=1;j<m;j++) scanf("%d",&xl[i][j]);
    s=0; t=(n-1)*(m-1)*2+1;
    for(int i=1;i<n;i++) {
        for(int j=1;j<m;j++) {
            tot++;
            if(j==1) 
                 add(tot,s,sl[i][j]);
            if(i==n-1) 
                add(tot,s,hl[i+1][j]);
            add(tot,tot+1,xl[i][j]);
            tot++;
            if(i==1) 
                add(tot,t,hl[i][j]);
            if(j==m-1) 
                add(tot,t,sl[i][j+1]);
            if(i!=1) 
                add(tot,(i-2)*(m-1)*2+(j-1)*2+1,hl[i][j]);            
            if(j!=m-1) 
                add(tot,tot+1,sl[i][j+1]);
        }
    }
    spfa();
    if(n==1||m==1)  ans=min(mi1,mi2);
    else ans=dis[t];
    printf("%d
",ans);
    }
    return 0;
}
View Code
原文地址:https://www.cnblogs.com/Achenchen/p/7516738.html