P1491 集合位置 次短路

这个题是一个次短路的裸题，就是把最短路路径求出来之后依次删边，然后跑最短路，在这些情况里取最小值就行了。

题干：

每次有大的活动，大家都要在一起“聚一聚”，不管是去好乐迪，还是避风塘，或者汤姆熊，大家都要玩的痛快。还记得心语和花儿在跳舞机上的激情与释放，还记得草草的投篮技艺是如此的高超，还记得狗狗的枪法永远是'S'……还有不能忘了，胖子的歌声永远是让我们惊叫的！！

今天是野猫的生日，所以想到这些也正常，只是因为是上学日，没法一起去玩了。但回忆一下那时的甜蜜总是一种幸福嘛。。。

但是每次集合的时候都会出现问题！野猫是公认的“路盲”，野猫自己心里也很清楚，每次都提前出门，但还是经常迟到，这点让大家很是无奈。后来，野猫在每次出门前，都会向花儿咨询一下路径，根据已知的路径中，总算能按时到了。

现在提出这样的一个问题：给出n个点的坐标，其中第一个为野猫的出发位置，最后一个为大家的集合位置，并给出哪些位置点是相连的。野猫从出发点到达集合点，总会挑一条最近的路走，如果野猫没找到最近的路，他就会走第二近的路。请帮野猫求一下这条第二最短路径长度。
输入输出格式
输入格式：

第一行是两个整数n(<=n<=)和m，表示一共有n个点和m条路，以下n行每行两个数xi，yi，(-<=xi,yi<=),代表第i个点的坐标，再往下的m行每行两个整数pj，qj,(<=pj,qj<=n)，表示两个点相通。

输出格式：

只有一行包含一个数，为第二最短路线的距离（保留两位小数），如果存在多条第一短路径，则答案就是第一最短路径的长度；如果不存在第二最短路径，输出－。

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(register int i = a;i <= n;++i)
#define lv(i,a,n) for(register int i = a;i >= n;--i)
#define clean(a) memset(a,0,sizeof(a))
#define mp make_pair
#define pr pair<db,int>
const int INF =  << ;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
    char c;
    bool op = ;
    while(c = getchar(), c < '' || c > '')
        if(c == '-') op = ;
    x = c - '';
    while(c = getchar(), c >= '' && c <= '')
        x = x *  + c - '';
    if(op) x = -x;
}
template <class T>
void write(T x)
{
    if(x < ) putchar('-'), x = -x;
    if(x >= ) write(x / );
    putchar('' + x % );
}
const int N = 1e5 + ;
struct node
{
    int l,r,nxt;
    db w;
}a[ * N];
int lst[N],len = ;
int n,m,px[N],py[N];
db dis[N];
int vis[N],pre[N];
void add(int x,int y,db w)
{
    a[++len].l = x;
    a[len].r = y;
    a[len].w = w;
    a[len].nxt = lst[x];
    lst[x] = len;
}
db dist(int x,int y)
{
    return sqrt((db)(px[x] - px[y]) * (px[x] - px[y]) + (db)(py[x] - py[y]) * (py[x] - py[y]));
}
priority_queue <pr,vector<pr>,greater<pr> > qu;
void dij()
{
    clean(vis);
    duke(i,,n)
    {
        dis[i] = ;
    }
    dis[] = ;
    qu.push(mp(dis[],));
    while(!qu.empty())
    {
        pr u = qu.top();
        qu.pop();
        int x = u.second;
        if(vis[x])
        continue;
        vis[x] = ;
//        cout<<x<<endl;
        for(int k = lst[x];k;k = a[k].nxt)
        {
            int y = a[k].r;
//            cout<<x<<" "<<y<<" "<<dis[y]<<" "<<dis[x]<<" "<<a[k].w<<endl;
            if(dis[y] > dis[x] + a[k].w)
            {
                dis[y] = dis[x] + a[k].w;
                pre[y] = x;
                qu.push(mp(dis[y],y));
            }
        }
    }
//    cout<<dis[n]<<endl;
}
db DIJ(int l,int r)
{
    clean(vis);
    duke(i,,n)
    {
        dis[i] = ;
    }
    dis[] = ;
    qu.push(mp(dis[],));
    while(!qu.empty())
    {
        pr u = qu.top();
        qu.pop();
        int x = u.second;
        if(vis[x])
        continue;
        vis[x] = ;
        for(int k = lst[x];k;k = a[k].nxt)
        {
            int y = a[k].r;
            if((x == l && y == r) || (y == l && x == r))
            {
//                cout<<"ok"<<endl;
                continue;
            }
            if(dis[y] > dis[x] + a[k].w)
            {
                dis[y] = dis[x] + a[k].w;
                pre[y] = x;
                qu.push(mp(dis[y],y));
            }
        }
    }
    return dis[n];
}
int main()
{
    read(n);read(m);
    duke(i,,n)
    read(px[i]),read(py[i]);
    duke(i,,m)
    {
        int x,y;
        read(x);read(y);
//        cout<<dist(x,y)<<" "<<px[x]<<" "<<py[x]<<" "<<px[y]<<" "<<py[y]<<endl;
        add(x,y,dist(x,y));
        add(y,x,dist(x,y));
    }
    dij();
    db minn = ;
    int l,r = n;
    while(r != )
    {
        l = r;
        r = pre[l];
//        cout<<l<<" "<<r<<endl;
        minn = min(DIJ(l,r),minn);
    }
    if(minn == )
    printf("-1\n");
    else
    printf("%.2lf\n",minn);
    return ;
}