【题解】[BeiJing2008]狼抓兔子—BZOJ1001。

(胡扯时间)今天炒鸡无聊就打算BZOJ开始从第一道题开始顺着打,这样未来一段时间内也就有事干了。结果发现A+B切掉后就遭遇了一个"小小"的瓶颈(真不友好。

好了说题说题。看题第一眼我下意识想打个最短路直接跑过,结果发现没有给兔子总共有多少(那就肯定最短不了了。

然后内心***之后发现还不是平面图,那我赶紧先认真起来好好打(恶)个(补)对偶图吧。

此时我又想到了qy巨佬教育我们的"脑补一下这个题,是不是很简单。"

脑补越强,对偶图越简单。

正解time:

显而易见这是一道赤果果的网络流,那就转化为对偶图然后直接跑个堆优化Dij吧。

在我MLE两遍T了一次以后发现还是得用spfa(qwq。

个人感觉最难的其实还是在平面图->对偶图的时候需要脑补每个面的编号进行加边。

我当时一直不太理解两个面连边是穿过的曾经的那个边就为现在这个边的权值应该如何实现。

后来发现完全可以直接建对偶图而不是平面图->对偶图。(orz。

我的加边代码可能和你们的不太一样(编号问题,所以推荐自己脑补一下自己的代码。

死于快读。!!!!!!!

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cctype>
#include<queue>
#define rg register
using namespace std;
//inline int read(){
//  rg int s=0,f=0;
//  rg char ch=getchar();
//  while(!isdigit(ch)) f|=(ch=='-'),ch=getchar();  
//  while(isdigit(ch)) s=(s<<1)+(s<<3)+(ch^48),ch=getchar();
//  return f?-s:s;
//}
int n,m,s,t,cnt;
const int MAX=6000015;
struct edge{
    int nxt;
    int to;
    int val;
}e[MAX];
int dis[MAX],head[MAX];
bool vis[MAX];
queue<int>q;
inline void add(int u,int v,int w){
    e[++cnt].nxt=head[u];
    e[cnt].to=v;
    e[cnt].val=w;
    head[u]=cnt;
}
int spfa(){
    memset(dis,0x3f,sizeof(dis));
    dis[s]=0;
    q.push(s); 
    while(!q.empty()){
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(rg int i=head[u];i;i=e[i].nxt){
            if(dis[e[i].to]>dis[u]+e[i].val){
                dis[e[i].to]=dis[u]+e[i].val;
                if(!vis[e[i].to]){
                    vis[e[i].to]=1;
                    q.push(e[i].to);
                }
            }
        }
    }
    return dis[t];
}
int main(){
//  n=read(),m=read();
    scanf("%d %d",&n,&m);
    s=0,t=(n-1)*(m-1)*2+1;
    for(rg int i=1;i<=m-1;++i){
//      int w=read();
        int w;
        scanf("%d",&w);
        add(i*2,t,w);
        add(t,i*2,w);
    }
    for(rg int i=2;i<=n-1;++i){
        for(rg int j=1;j<=m-1;++j){
//          int w=read();
            int w;
            scanf("%d",&w);
            add((i-1)*(m-1)*2+j*2,(i-1)*(m-1)*2+j*2-m*2+1,w);
            add((i-1)*(m-1)*2+j*2-m*2+1,(i-1)*(m-1)*2+j*2,w);
        }
    }
    for(rg int i=1;i<=m-1;++i){
//      int w=read();
        int w;
        scanf("%d",&w);
        add((n-2)*2*(m-1)+i*2-1,s,w);
        add(s,(n-2)*2*(m-1)+i*2-1,w);
    }
    for(rg int i=1;i<=n-1;++i){
//      int w=read();
        int w;
        scanf("%d",&w);
        add((i-1)*(m-1)*2+1,s,w);
        add(s,(i-1)*(m-1)*2+1,w);
        for(rg int j=1;j<=m-2;++j){
//          w=read();
            scanf("%d",&w);
            add((i-1)*(m-1)*2+j*2+1,(i-1)*(m-1)*2+j*2,w);
            add((i-1)*(m-1)*2+j*2,(i-1)*(m-1)*2+j*2+1,w);
        }
//      w=read();
        scanf("%d",&w);
        add((i-1)*(m-1)*2+(m-2)*2+2,t,w);
        add(t,(i-1)*(m-1)*2+(m-2)*2+2,w);
    }
    for(rg int i=1;i<=n-1;++i){
        for(rg int j=1;j<=m-1;++j){
//          int w=read();
            int w;
            scanf("%d",&w);
            add((i-1)*(m-1)*2+j*2-1,(i-1)*(m-1)*2+j*2,w);
            add((i-1)*(m-1)*2+j*2,(i-1)*(m-1)*2+j*2-1,w);
        }
    }
    int ans=spfa();
    printf("%d",ans);
    return 0;
}

转载于:https://www.cnblogs.com/Sinyess/p/10603869.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值