做这道题的动机就是想练习一下堆的应用,顺便补一下好久没看的图论算法。

Dijkstra算法概述

//从0出发的单源最短路

dis[][] = {INF}
ReadMap(dis);
for i = 0 -> n - 1
    d[i] = dis[0][i]
while u = GetNearest(1 .. n - 1, !been[])
    been[u] = 1
    for_each edge from u
        d[edge.v] = min(d[edge.v], d[u] + dis[u][edge.v])

上述算法遍历所有节点,每次 GetNearest() 循环一次,并遍历了所有边,算法复杂度 O(V2+E) = O(V2)

其中 GetNearest() 总取未去过的点中 d[] 最小的点,可以用小根堆维护 d[] 数组优化。

堆优化算法概述

//从0出发的单源最短路

dis[][] = {INF}
ReadMap(dis);
for i = 0 -> n - 1
    d[i] = dis[0][i]
heap = BuildHeap(d[], n - 1)
while u = heap.pop()
    been[u] = 1
    for_each edge from u
        if IsBetterDist()
            d[edge.v] = d[u] + dis[u][edge.v]
            heap.up(heap.IndexOf(edge.v))

建堆 O(VlogV),n次pop() O(VlogV),遍历边并更新堆 O(ElogV),优化算法复杂度 O(2VlogV+ElogV) = O((V+E)logV)

其他心得

* 几何类最短路一定要判重点

//POJ 2502
//Dijkstra 堆优化试水
//飘忽不定的迷离的数据
//AC 2016-10-17

#define HEAP_OPTIMIZE

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#define MAXN (200 + 2) + 10
#define SUB_SPEED (40000.0 / 60.0)
#define WALK_SPEED (10000.0 / 60.0)
#define INF 0x7f7f7f7f
using namespace std;

inline double sqr(int x){
    return x * x;
}

struct point {
    int x, y;
    point(){}
    point(int X, int Y): x(X), y(Y){}
    inline double norm(){
        return sqrt(sqr(x) + sqr(y));
    }
    friend point operator - (const point &p1, const point &p2){
        return point(p1.x - p2.x, p1.y - p2.y);
    }
    friend bool operator == (const point &p1, const point &p2){
        return (p1.x == p2.x) && (p1.y == p2.y);
    }
}pt[MAXN];

int n = 0;

int IndexOf(const point &p){
    for (int i = 0; i < n; i++){
        if (pt[i] == p)
            return i;
    }
    return n;
}

double dis[MAXN][MAXN], d[MAXN];
bool been[MAXN];

struct BHeap{
    int heap[MAXN], n, index[MAXN];
    BHeap(): n(0){}
    BHeap(int N): n(0){
        for (int i = 1; i <= N; i++)
            push(i);
    }
    void down(int i){
        for (int j = 2 * i; j <= n; j <<= 1){
            j += j < n && d[heap[j]] > d[heap[j + 1]];
            if (d[heap[j]] < d[heap[i]]){
                swap(index[heap[i]], index[heap[j]]);
                swap(heap[i], heap[j]);
                i = j;
            }
            else break;
        }
    }
    void up(int i){
        for (int j = i/2; j > 0; j >>= 1){
            if (d[heap[j]] > d[heap[i]]){
                swap(index[heap[i]], index[heap[j]]);
                swap(heap[i], heap[j]);
                i = j;
            }
            else break;
        }
    }
    void push(int a){
        heap[++n] = a;
        index[a] = n;
        up(n);
    }
    int pop(){
        if (!n) return 0;
        swap(index[heap[1]], index[heap[n]]);
        swap(heap[1], heap[n--]);
        down(1);
        return heap[n + 1];
    }
}heap;

int main(){
    freopen("fin.c", "r", stdin);
    for (int i = 0; i < MAXN; i++)
        for(int j = 0; j < MAXN; j++)
			if (i == j) dis[i][j]=0;
			else dis[i][j] = INF;
    scanf("%d%d%d%d", &pt[0].x, &pt[0].y, &(pt[1].x), &pt[1].y);
    n = 2;
    dis[0][1] = dis[1][0] = (pt[0] - pt[1]).norm() / WALK_SPEED;
    point p;
    while (~scanf("%d%d", &p.x, &p.y)){
        int iter = IndexOf(p), olditer = - 1;
        if (iter == n){
            pt[n++] = p;
            for (int i = 0; i < n; i++){
                dis[i][n - 1] = dis[n - 1][i] = (pt[n - 1] - pt[i]).norm() / WALK_SPEED;
            }
        }
        olditer = iter;
        while (scanf("%d%d", &p.x, &p.y), ~p.x && ~p.y){
            int iter = IndexOf(p);
            if (iter == n){
                pt[n++] = p;
                for (int i = 0; i < n; i++){
                    if (i != olditer)
                        dis[i][n - 1] = dis[n - 1][i] = (pt[n - 1] - pt[i]).norm() / WALK_SPEED;
                    else
                        dis[i][n - 1] = dis[n - 1][i] = (pt[n - 1] - pt[i]).norm() / SUB_SPEED;
                }
            }
            else dis[olditer][iter] = dis[iter][olditer] = (pt[iter] - pt[olditer]).norm() / SUB_SPEED;
            olditer = iter;
        }
    }
    for (int i = 1; i < n; i++)
        d[i] = dis[0][i];
    #ifdef HEAP_OPTIMIZE
        heap = BHeap(n - 1);
    #endif
    while (
    #ifdef HEAP_OPTIMIZE
        int v = heap.pop()
    #else
        1
    #endif
    ){
        #ifndef HEAP_OPTIMIZE
	        double minn = INF;
	        int v = - 1;
	        for (int i = 1; i < n; i ++)
	            if ((!been[i])&&(d[i] < minn)){
	                minn = d[i];
	                v = i;
	            }
	        if (!~v) break;
        #endif
        been[v] = 1;
        for (int i = 1; i < n; i++){
            if (!been[i]){
                if (d[v] + dis[v][i] < d[i]){
                    d[i] = d[v] + dis[v][i];
                    #ifdef HEAP_OPTIMIZE
                        heap.up(heap.index[i]);
                    #endif
                }
            }
        }
    }
    printf("%.f\n", d[1]);
}

POJ 2502 - Subway Dijkstra堆优化试水的更多相关文章

  1. POJ 2502 Subway (Dijkstra 最短+建设规划)

    Subway Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 6689   Accepted: 2176 Descriptio ...

  2. Bzoj 2834: 回家的路 dijkstra,堆优化,分层图,最短路

    2834: 回家的路 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 62  Solved: 38[Submit][Status][Discuss] D ...

  3. POJ 1511 - Invitation Cards 邻接表 Dijkstra堆优化

    昨天的题太水了,堆优化跑的不爽,今天换了一个题,1000000个点,1000000条边= = 试一试邻接表 写的过程中遇到了一些问题,由于习惯于把数据结构封装在 struct 里,结果 int [10 ...

  4. Dijkstra堆优化学习

    最短路径例题 今天特地学习了Dijkstra的堆优化(主要是慕名已久). 我们需要一个堆来记录[编号,到编号这个点的最短路径值(当然只是当前的)] 与原来的Dijkstra操作基本一致,主要有以下几点 ...

  5. 【Dijkstra堆优化】洛谷P2243电路维修

    题目背景 Elf 是来自Gliese 星球的少女,由于偶然的原因漂流到了地球上.在她无依无靠的时候,善良的运输队员Mark 和James 收留了她.Elf 很感谢Mark和James,可是一直也没能给 ...

  6. (简单) POJ 2502 Subway,Dijkstra。

    Description You have just moved from a quiet Waterloo neighbourhood to a big, noisy city. Instead of ...

  7. Dijkstra+计算几何 POJ 2502 Subway

    题目传送门 题意:列车上行驶40, 其余走路速度10.问从家到学校的最短时间 分析:关键是建图:相邻站点的速度是40,否则都可以走路10的速度.读入数据也很变态. #include <cstdi ...

  8. POJ 2502 Subway(迪杰斯特拉)

    Subway Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 6692   Accepted: 2177 Descriptio ...

  9. POJ 2502 Subway

    Subway Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 4928   Accepted: 1602 Descriptio ...

随机推荐

  1. Entity Framework 6 Recipes 2nd Edition(11-12)译 -&gt; 定义内置函数

    11-12. 定义内置函数 问题 想要定义一个在eSQL 和LINQ 查询里使用的内置函数. 解决方案 我们要在数据库中使用IsNull 函数,但是EF没有为eSQL 或LINQ发布这个函数. 假设我 ...

  2. C# List 排序各种用法与比较

    下面介绍各种List的sort的用法与比较 首先,我们建一个People的实体,有name.age.sex的属性,我们要排序的字段是年龄age 新建一个实体类 public class People ...

  3. HDU 5398 (动态树)

    Problem GCD Tree 题目大意 n个点的无向完全图,标号1~n,每条边u-->v 的权值为gcd(u,v),求其最大生成树,输出最大边权和. n<=10^5,有多个询问. 解题 ...

  4. Java多线程系列--“JUC锁”06之 Condition条件

    概要 前面对JUC包中的锁的原理进行了介绍,本章会JUC中对与锁经常配合使用的Condition进行介绍,内容包括:Condition介绍Condition函数列表Condition示例转载请注明出处 ...

  5. Android使用SQLite数据库(4)

    读取SQLite数据库中的字符串字段,使用Cursor的getString方法(其他类型的字段也有相应的读取方法): public abstract String getString (int col ...

  6. SQL复制一个表的数据到另一个表

    最近做一个项目,由于客户数据量大,为了不将数据彻底删除,于是将数据移动到历史表,原始表的数据删除.由于技术有限,想不到好的方法,于是写个存储过程 执行,为了防止执行过程中出现异常,执行不完整.用到hI ...

  7. NoSQL精粹(NoSQL Distilled)——序言

    之前说到博客长草的问题,想了想除了很忙特别忙非常忙各种瞎忙忙你妹啊外,主要还是不知道写什么好--到这家公司的两年中从JS到领域驱动到缓存服务器从前端到后端各种折腾,有些东西虽然有所心得,不过既然前人已 ...

  8. oracle数据泵备份(Expdp命令)[转]

      Oracle备份方式主要分为数据泵导出备份.热备份与冷备份三种,今天首先来实践一下数据泵备份与还原.数据泵导出/导入属于逻辑备份,热备份与冷备份都属于物理备份.oracle10g开始推出了数据泵( ...

  9. oracle sqlplus及常用sql语句

    常用sql语句 有需求才有动力 http://blog.csdn.net/yitian20000/article/details/6256716 常用sql语句 创建表空间:create tables ...

  10. Javascript学习笔记1 数论

    1.Javascript不用担心内存的回收与对象的销毁! 2.Javascript有:±infinity.NaN全局变量表示 被0整除的±无穷 和 非数字.undefined和null表示 未定义 和 ...