poj 2502 Subway

题目链接：

　　http://poj.org/problem?id=2502

题目大意：

　　一个学生去上学，从家到学校，可以有若干个地铁路线，每个地铁路线有若干站，给出步行和地铁的速度，问：最短用多长时间从家到达学校？

解题思路：

　　有地铁的两点建立地铁路线，没有的步行，求最短路，但是有一点坑的是地铁线有可能是曲线，也就是说从a站到b站再到c站的距离和大于a站直接到c站的距离

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <cmath>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;

#define maxn 310
#define INF 0x3f3f3f3f
#define walk (10000.0/60)
#define sub (40000.0/60)
double map[maxn][maxn], dist[maxn];
int k = 2;

void init ();
void dijkstra ();

int main ()
{
    int i, j, s;
    double x[maxn], y[maxn];

    init ();
    while (scanf ("%lf %lf %lf %lf", &x[0], &y[0], &x[1], &y[1]) != EOF)
    {
        while (scanf ("%lf %lf", &x[k], &y[k]) != EOF)
        {
            s = k;
            while (x[k]!=-1 || y[k] != -1)
            {
                k ++;
                scanf ("%lf %lf", &x[k], &y[k]);
            }

            for (i=s+1; i<k; i++)
                    map[i-1][i] = map[i][i-1] = sqrt((x[i] - x[i-1])*(x[i] - x[i-1]) + (y[i] - y[i-1])*(y[i] - y[i-1])) / sub;

        }
        for (i=0; i<k; i++)
            for (j=i+1; j<k; j++)
                map[i][j] = map[j][i] = min (map[j][i], sqrt((x[i] - x[j])*(x[i] - x[j]) + (y[i] - y[j])*(y[i] - y[j])) / walk);
        dijkstra ();
        printf ("%.0f
", dist[1]);
    }
    return 0;
}

void init ()
{
    for (int i=0; i<maxn; i++)
        for (int j=0; j<maxn; j++)
            if (i == j)
                map[i][j] = 0;
            else
                map[i][j] = INF;
}
void dijkstra ()
{
    bool vis[maxn];
    int temp, index, i, j;
    memset (vis, false, sizeof(vis));
    vis[0] = 1;
    for (i=0; i<k; i++)
        dist[i] = map[0][i];
    for (i=1; i<k; i++)
    {
        temp = INF;
        for (j=0; j<k; j++)
            if (!vis[j] && temp > dist[j])
            {
                temp = dist[j];
                index = j;
            }

        vis[index] = true;

        for (j=0; j<k; j++)
            if (!vis[j])
                dist[j] = min(dist[j], dist[index]+map[index][j]);
    }
}