洛谷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。

输入输出样例

输入样例#1： 

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

输出样例#1： 

2.83

说明

各个测试点1s

注：此题有毒。。。

不能用 k短路。。。

即：不能走过一条边后又走回来。。。

坑。。。

所以，spfa，再每次删去最短路上一条边，跑spfa，求最小值。

附代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<cmath>
#include<cstring>
#define MAXN 210
#define MAXM 40010
#define MAX 1e9
using namespace std;
int n,m,s,t,k,c=2;
int head[MAXN],fa[MAXN];
double ans,path[MAXN];
bool flag=true,vis[MAXN];
struct node1{
    int next,to;
    double w;
}a[MAXM<<1];
struct node2{
    double x,y;
}b[MAXN];
inline int read(){
	int date=0,w=1;char c=0;
	while(c<'0'||c>'9'){if(c=='-')w=-1;c=getchar();}
	while(c>='0'&&c<='9'){date=date*10+c-'0';c=getchar();}
	return date*w;
}
inline double dis(int i,int j){
    return sqrt((b[i].x-b[j].x)*(b[i].x-b[j].x)+(b[i].y-b[j].y)*(b[i].y-b[j].y));
}
inline int relax(int u,int v,double w){
    if(path[v]>path[u]+w){
        path[v]=path[u]+w;
        if(flag)fa[v]=u;
        return 1;
    }
    return 0;
}
inline void add(int u,int v,double w){
    a[c].to=v;a[c].w=w;
    a[c].next=head[u];
    head[u]=c++;
    a[c].to=u;a[c].w=w;
    a[c].next=head[v];
    head[v]=c++;
}
void spfa(){
	int u,v;
	queue<int> q;
	for(int i=1;i<=n;i++){path[i]=MAX;vis[i]=false;if(flag)fa[i]=0;}
	path[s]=0;
	vis[s]=true;
	q.push(s);
	while(!q.empty()){
		u=q.front();
		q.pop();
		vis[u]=false;
		for(int i=head[u];i;i=a[i].next){
			v=a[i].to;
			if(relax(u,v,a[i].w)&&!vis[v]){
				vis[v]=true;
				q.push(v);
			}
		}
	}
}
int main(){
    int u,v;
    double w;
    n=read();m=read();
    k=2;s=1;t=n;
    for(int i=1;i<=n;i++){b[i].x=read();b[i].y=read();}
    for(int i=1;i<=m;i++){
        u=read();v=read();w=dis(u,v);
        add(u,v,w);
    }
    spfa();
    flag=false;ans=MAX;
    if(path[t]==MAX){
        printf("-1\n");
        return 0;
    }
    for(int i=t;i!=s;i=fa[i])
    for(int j=head[i];j;j=a[j].next){
        int v=a[j].to;
        if(v==fa[i]){
            double w=a[j].w;
            a[j].w=a[j^1].w=MAX;
            spfa();
            a[j].w=a[j^1].w=w;
            ans=min(ans,path[t]);
            break;
        }
    }
    if(ans!=MAX)printf("%.2lf\n",ans);
    else printf("-1\n");
	return 0;
}