BZOJ 4152 [AMPPZ2014]The Captain 最短路题解

Time Limit: 20 Sec Memory Limit: 256 MB
Submit: 1108 Solved: 433

---

Description

给定平面上的n个点，定义(x1,y1)到(x2,y2)的费用为min(|x1-x2|,|y1-y2|)，求从1号点走到n号点的最小费用。

---

Input

第一行包含一个正整数n(2<=n<=200000)，表示点数。

接下来n行，每行包含两个整数x[i],yi，依次表示每个点的坐标。

---

Output

一个整数，即最小费用。

---

Sample Input

5

2 2

1 1

4 5

7 1

6 7

---

Sample Output

2

---

题解：

显然这道题就是一个最短路的求解，我们考虑这样一张图的几个点

我们会发现如果我们对A，B，C点都按x排好序了，那么我们从A到C的方案中，直接从A到C不会比从A到再到C的优，因为无论如何，我已经排好了序，那么最差的情况我也是每次都选择走x，这样就和直接过去走x的代价是等价的，如果我每次都走y，就和直接走y过去的代价是等价的，所以我既可以选择走x，也可以选择走y，所以对于x排序，然后相邻的建边（之所以要相邻，从a建到b，再从b建到c的原因已经说了，因为这样不会比从A直接到C差），对于y也排序，按照y的顺序建边即可

提示：①这道题卡SPFA，所以我用的Dijkstra的堆优化过的②这道题爆int，要用long long

---

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<queue>
#define ppap pair<long long,int>
const int MAXN=200000*2;
using namespace std;
int readin(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int head[MAXN],tail,n;
long long dis[MAXN];
bool vis[MAXN];
struct Point{int x,y,id;}p[MAXN];
struct Line{int from,to,num,nxt;}line[MAXN*2+10];
bool cmp(const Point &A,const Point &B){return A.x<B.x;}
bool cmpp(const Point &A,const Point &B){return A.y<B.y;}
void add_line(int from,int to,int num){tail++;line[tail].from=from;line[tail].to=to;line[tail].nxt=head[from];line[tail].num=num;head[from]=tail;}
void add(int a,int b,int num){add_line(a,b,num);add_line(b,a,num);}
void build(){
    sort(p+1,p+n+1,cmp);
    for(register int i=1;i<=n-1;i++) 
    add(p[i].id,p[i+1].id,p[i+1].x-p[i].x);
    sort(p+1,p+n+1,cmpp);
    for(register int i=1;i<=n-1;i++) 
    add(p[i].id,p[i+1].id,p[i+1].y-p[i].y);
}
void SPFA(){
    priority_queue<ppap,vector<ppap>,greater<ppap> >q;
    for(register int i=1;i<=n;i++) dis[i]=0x7fffffff;
    dis[1]=0;
    q.push(make_pair(0,1));
    while(!q.empty()){
        int u=q.top().second;q.pop();
        if(vis[u])continue; vis[u]=true;
        for(register int i=head[u];i;i=line[i].nxt){
            int v=line[i].to;
            if(dis[v]>dis[u]+line[i].num){
                dis[v]=dis[u]+line[i].num;
                q.push(make_pair(dis[v],v));
            }
        }
    }
    printf("%lld\n",dis[n]);
}
int main(){
    //freopen(".txt","r",stdin);
    //freopen(".out","w",stdout);
    n=readin();
    for(register int i=1;i<=n;i++){p[i].x=readin();p[i].y=readin();p[i].id=i;}
    build();
    SPFA();
    return 0;
}