分层图最短路【bzoj2834】: 回家的路

分层图最短路【bzoj2834】: 回家的路

题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=2834

这道题难在建边。

自己写的时候想到了两边排序这些细节,但是思路出错在跑最短路的时候。

前几天做的水壶也是难在建边。一想那个题我就开始敲广搜,可是要比自己想的难。之后冷静了一下想到了排序建边。建完边之后就可以开始跑最短路了。

这道题确实加深了我对分层图最短路的理解。

因为当时自己敲的时候还是设了两维的状态,第二维我设的是方向,因为只会横着走或者竖着走,只要在转弯的时候枚举两种情况就可以了。并没有证出是错的,但是太难写了。

瞟了眼题解发现原来我们可以拆点建边,只要在自己的等价点位置连一条边权为1的边表示转弯就可以了。

妙啊。

原来自己对分层图这东西还是一知半解,继续刷吧。

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
const int wx=2000017;
inline int read(){
    int sum=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){sum=(sum<<1)+(sum<<3)+ch-'0';ch=getchar();}
    return sum*f;
}
int head[wx],dis[wx],vis[wx];
int n,m,sx,sy,tx,ty,num,S,T; 
struct e{
    int nxt,to,dis;
}edge[wx*2];
void add(int from,int to,int dis){
    edge[++num].nxt=head[from];
    edge[num].to=to;
    edge[num].dis=dis;
    head[from]=num;
}
int dx[]={0,1,0,-1,0};
int dy[]={0,0,1,0,-1};
int Abs(int a){
    if(a<0)return -a;
    return a;
}
struct node{
    int x,y,id;
}t[wx];
queue<int > q;
bool cmp1(node a,node b){
    if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}
bool cmp2(node a,node b){
    if(a.y==b.y)return a.x<b.x;
    return a.y<b.y;
}
void SPFA(int s){
    for(int i=1;i<=m*2+4;i++)dis[i]=0x3f3f3f3f;
    dis[s]=0;vis[s]=1;q.push(s);
    while(q.size()){
        int u=q.front();q.pop();vis[u]=0;
        for(int i=head[u];i;i=edge[i].nxt){
            int v=edge[i].to;
            if(dis[v]>dis[u]+edge[i].dis){
                dis[v]=dis[u]+edge[i].dis;
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
}
int main(){
    n=read();m=read();
    for(int i=1;i<=m+2;i++){
        t[i].x=read();t[i].y=read();t[i].id=i;
    }
    sort(t+1,t+3+m,cmp1);
    for(int i=1;i<m+2;i++){
        if(t[i+1].x==t[i].x)add(t[i].id,t[i+1].id,2*(t[i+1].y-t[i].y)),add(t[i+1].id,t[i].id,2*(t[i+1].y-t[i].y));
    }
    sort(t+1,t+3+m,cmp2);
    for(int i=1;i<m+2;i++){
        if(t[i+1].y==t[i].y)add(t[i].id+m+2,t[i+1].id+m+2,2*(t[i+1].x-t[i].x)),add(t[i+1].id+m+2,t[i].id+m+2,2*(t[i+1].x-t[i].x)); 
    }
    for(int i=1;i<=m;i++)add(i,i+m+2,1),add(i+2+m,i,1);
    add(m+1,m+2+m+1,0);add(m+2+m+1,m+1,0);add(m+2,m+m+4,0);add(m+m+4,m+2,0);
    SPFA(m+1);
    if(dis[m+2]==0x3f3f3f3f)dis[m+2]=-1;
    printf("%d\n",dis[m+2]);
    return 0;
}

猜你喜欢

转载自www.cnblogs.com/wangxiaodai/p/9760580.html