P1491 集合位置 次短路

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

题干:

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

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

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

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

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

输出格式:

只有一行包含一个数,为第二最短路线的距离(保留两位小数),如果存在多条第一短路径,则答案就是第一最短路径的长度;如果不存在第二最短路径,输出-1

代码:

#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 = 1 << 30;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
    char c;
    bool op = 0;
    while(c = getchar(), c < '0' || c > '9')
        if(c == '-') op = 1;
    x = c - '0';
    while(c = getchar(), c >= '0' && c <= '9')
        x = x * 10 + c - '0';
    if(op) x = -x;
}
template <class T>
void write(T x)
{
    if(x < 0) putchar('-'), x = -x;
    if(x >= 10) write(x / 10);
    putchar('0' + x % 10);
}
const int N = 1e5 + 5;
struct node
{
    int l,r,nxt;
    db w;
}a[10 * N];
int lst[N],len = 0;
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,1,n)
    {
        dis[i] = 2147483647;
    }
    dis[1] = 0;
    qu.push(mp(dis[1],1));
    while(!qu.empty())
    {
        pr u = qu.top();
        qu.pop();
        int x = u.second;
        if(vis[x])
        continue;
        vis[x] = 1;
//        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,1,n)
    {
        dis[i] = 2147483647;
    }
    dis[1] = 0;
    qu.push(mp(dis[1],1));
    while(!qu.empty())
    {
        pr u = qu.top();
        qu.pop();
        int x = u.second;
        if(vis[x])
        continue;
        vis[x] = 1;
        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,1,n)
    read(px[i]),read(py[i]);
    duke(i,1,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 = 2147483647;
    int l,r = n;
    while(r != 1)
    {
        l = r;
        r = pre[l];
//        cout<<l<<" "<<r<<endl;
        minn = min(DIJ(l,r),minn);
    }
    if(minn == 2147483647)
    printf("-1
");
    else
    printf("%.2lf
",minn);
    return 0;
}
原文地址:https://www.cnblogs.com/DukeLv/p/10401172.html