Luogu P1491 集合位置 (spfa--次短路)

题目描述

每次有大的活动，大家都要在一起“聚一聚”，不管是去好乐迪，还是避风塘，或者汤姆熊，大家都要玩的痛快。还记得心语和花儿在跳舞机上的激情与释放，还记得草草的投篮技艺是如此的高超，还记得狗狗的枪法永远是'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。

输入输出样例

输入样例#1： 复制

3 3
0 0
1 1
0 2
1 2
1 3
2 3

输出样例#1： 复制

2.83

说明

各个测试点1s

---

大概就是求次短路
好像可以用A*来求k短路，不过我不会。。。
那就换一种思路
既然不是最短路，那么它的路径就一定不能和最短路完全重合
也就是说，至少会有一条边与最短路不同
但是我们不知道哪条/哪些边会与最短路不同
那么我们在第一次spfa时记录下来路径（前驱结点）
然后我们就枚举边，将最短路里的每条边都删去一次之后再跑一边spfa求最短路
这个题显然两个点之间只会有一条边相连，所以这个题中枚举边等价于枚举最短路上的相邻两点
这样求出次短路
而且最多也就只有200个点，所以这种方法是行得通的！
代码(跑的比较慢，空间占的也比较多，有空可能会优化一下。。。)：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <queue>
#include <vector>
#define For(i,l,r) for(int i=l;i<=r;++i)
using namespace std;
int read()
{
    char c;
    bool t=0;
    int a=0;
    while((c=getchar())==' '||c=='\n'||c=='\r');
    if(c=='-')
    {
        t=1;
        c=getchar();
    }
    while(isdigit(c))
    {
        a*=10;
        a+=(c-'0');
        c=getchar();
    }
    return a*(t?-1:1);
}
struct line{
    int to,next;
    double v;
}edge[80010];
struct poi{
    int x,y;
}node[201];
int last[201],pre[201],n,m;
double dis[201];
bool vis[201];
double dist(int x,int y)
{
    return sqrt((node[x].x-node[y].x)*(node[x].x-node[y].x)+(node[x].y-node[y].y)*(node[x].y-node[y].y));
}
void add(int from,int to,int i)
{
    edge[i].next=last[from];
    last[from]=i;
    edge[i].to=to;
    edge[i].v=dist(from,to);
}
void spfa()
{
    int tx;
    memset(dis,127,sizeof dis);
    queue<int> q;
    q.push(1);dis[1]=0;vis[1]=1;
    while(!q.empty())
    {
        tx=last[q.front()];
        while(tx)
        {
            if(dis[q.front()]+edge[tx].v<dis[edge[tx].to])
            {
                dis[edge[tx].to]=dis[q.front()]+edge[tx].v;
                pre[edge[tx].to]=q.front();
                if(!vis[edge[tx].to])
                {
                    q.push(edge[tx].to);
                    vis[edge[tx].to]=1;
                }
            }
            tx=edge[tx].next;
        }
        vis[q.front()]=0;
        q.pop();
    }
    return;
}
void spfa2(int x,int y)
{
    int tx;
    memset(dis,127,sizeof dis);
    queue<int> q;
    q.push(1);dis[1]=0;vis[1]=1;
    while(!q.empty())
    {
        tx=last[q.front()];
        while(tx)
        {
            if((q.front()==x&&edge[tx].to==y)||(q.front()==y&&edge[tx].to==x))
            {
                tx=edge[tx].next;
                continue;
            }
            if(dis[q.front()]+edge[tx].v<dis[edge[tx].to])
            {
                dis[edge[tx].to]=dis[q.front()]+edge[tx].v;
                if(!vis[edge[tx].to])
                {
                    q.push(edge[tx].to);
                    vis[edge[tx].to]=1;
                }
            }
            tx=edge[tx].next;
        }
        vis[q.front()]=0;
        q.pop();
    }
    return;
}
int main()
{
    int tx,ty,now;
    double ans=0,tmp;
    n=read();m=read();now=n;
    m+=m;
    For(i,1,n)
    {
        node[i].x=read();
        node[i].y=read();
    }
    For(i,1,m)
    {
        tx=read();
        ty=read();
        add(tx,ty,i++);
        add(ty,tx,i);
    }
    spfa();
    while(now!=1)
    {
        spfa2(now,pre[now]);
        if(dis[n]<ans||!ans)
         ans=dis[n];
        now=pre[now];
    }
    printf("%.2f",ans);
    return 0;
}