【洛谷】P1491 集合位置【次短路】两种解法

名欺！

已于 2023-04-11 22:54:25 修改

阅读量199

点赞数 2

文章标签：算法 c++ 图论

于 2023-04-11 22:52:53 首次发布

本文链接：https://blog.csdn.net/qq_52156689/article/details/130095569

版权

集合位置

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

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

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

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

输入格式

第一行是两个整数 $\le n \le 200)$ 和 $m$ ，表示一共有 $n$ 个点和 $m$ 条路，以下 $n$ 行每行两个数 $x_i$ ， $y_i$ ， $\le x_i,y_i \le 500),$ 代表第 $i$ 个点的坐标，再往下的 $m$ 行每行两个整数 $p_j$ ， $q_j,(1 \le p_j,q_j \le n)$ ，表示两个点相通。

输出格式

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

样例 1

样例输入 1

样例输出 1

2.83

不标准板子（枚举边的端点更新答案）ans = min(ans , d1[u]+dis(u->v1(v1!=v)) + d2[v])

//这题不可以重复跑，整体思路，不删边（删边是次短路的板子）。我们先1->n跑一遍最短路d1[N]存正着跑最短路，
//再n->1 跑一遍最短路,d2[N]存正着跑最短路
//再存一下最短路的路径 ，遍历路径上的每条边（u->v）答案ans = min(ans , d1[u]+dis(u->v1(v1!=v)) + d2[v]),当然在
//d2[v]的路径中可能会出现d1[u]经过的点(这题不能有重复的点)。所以在每次比较ans时，我们去暴力判断一下 
//这两条路径中是否有重复的点，所以在跑最短路是随便记录两个方向的路径fr1[N] , fr2[N],如果fr1[u]->1路径中
//有n->fr2[v]经过的点,此时d1[u]+d1[u]+dis( u->v1(v1!=v) )  + d2[v]不能去更新ans。 

#include<bits/stdc++.h>
#define int long long
using namespace std; 
const int N = 6e5+10;
const double INF = 1e15;
double x[N] ,y[N];
int n , m;
int h[N] , ne[N] , idx , e[N];
double val[N];
double clac(int p  ,int q){
	return sqrt((x[p]-x[q])*(x[p]-x[q])+(y[p]-y[q])*(y[p]-y[q]));
}
struct node{
	int u ;
	double dis ;
	bool operator <(const node &A)const{
	return dis > A.dis;
	}
};
int st[N];
double d[N] ,d1[N] , d2[N];
int fr1[N] , fr2[N];
int findfr(int u  , int v){ //判断两条路径中有没有相交 
   unordered_map <int , int> mp;
		for(int j = u; fr1[j]!=0 ; j=fr1[j]){
			mp[j]++;
		}
		for(int j = v; fr2[j]!=0 ; j=fr2[j]){
			if(mp[j]) return 1;
			mp[j]++;
		}
	return 0;
}
void dijstra(int s){
	for(int i = 0; i <= n ; i++ ) d[i] = INF , st[i]= 0;
	
	priority_queue<node> q;
	d[s] = 0;
	q.push((node){s ,d[s]});
	while(!q.empty()){
		node top = q.top();
		q.pop();
		int u = top.u ;
		double dis =top.dis;
		if(st[u])continue;
		st[u]=1;
		for(int i = h[u] ; i!=-1 ; i=ne[i]){
			int j = e[i];
			
			if(d[j] > val[i]+ dis){
				d[j] = dis+val[i];
				q.push((node){j , d[j]});
				if(s == 1) fr1[j] = u;///记录路径 
				else fr2[j] = u;
			}
		}
	}
	
}
void add(int a ,int b , double c){
	e[idx] = b;
	val[idx] = c;
	ne[idx] = h[a];
	h[a] = idx++;
}
signed main(){
	std::ios::sync_with_stdio(false);
	std::cin.tie(0);
	cin >> n >> m;
	for(int i = 1 ; i <= n ;i++){
		cin >>x[i] >> y[i];
	}
	memset(h , -1 ,sizeof h);
	for(int i = 1; i <= m ;i++){
		int p , q;
		cin >> p >>q;
	    add(p , q , clac(p ,q));
		add(q,p , clac(p,q)); 
	}
	dijstra(1);
	for(int i =1 ;i <= n ;i++) d1[i] = d[i];
	dijstra(n);
	for(int i = 1; i <= n ; i++) d2[i] = d[i];
	//存图，跑图 
	
	double ans = INF;
	for(int j = n; fr1[j]!=0 ; j=fr1[j]){ // 遍历每条最短路径上的边 
	    int u = fr1[j];
	    for(int i = h[u] ; i!=-1 ; i= ne[i]){
	    	int v = e[i];
	    	if(v != j && !findfr(u , v)){  
	    		ans = min(ans , val[i]+d1[u] + d2[v]);
	    	}
	    }
	}
	if(ans!=INF)
	cout <<setiosflags(ios::fixed) << setprecision(2) <<ans <<endl;
	else cout << -1 <<endl; 
	return 0;
}

标准板子（删边）

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
#include<cmath>
using namespace std;
const int  N = 210;
const int M = 5e4+10;
const int INF = 1e9+7;
int n , m ;
double d[N];
double x[M] , y[M];
int fre[N];
bool flag[N];
int e[2*M] , idx , ne[M*2] ,h[2*M] ;
double val[2*M];
void add(int a, int b , double c ){
    e[idx] = b;
    val[idx] = c;
    ne[idx] = h[a];
    h[a] = idx++;
}

typedef struct Node{
    int u ;
    double di;
    bool operator < (const Node &A )const{
        return di > A.di;
    }
}node;

void add(int a ,int b ,int c ){
    e[idx] = b;
    ne[idx] = h[a];
    val[idx] = c;
    h[a] = idx++;
}

double calc(double x1 ,double y1 , double x2 ,double y2){
  
    return  sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}

void dijstra(int s, int x , int y){
    priority_queue<Node> q;
    
    for(int i =1 ; i <= n ; i++) {d[i] = INF;flag[i] = false;}
    d[s] = 0;
    node top={1,0};
    
    q.push(top);
    while(!q.empty()){
        node top = q.top();
        int u = top.u ;
        double dis = top.di;q.pop();
       // cout << "u: " << u <<endl;
        if(flag[u]) continue;
        flag[u] =true;
    
        
        for(int i = h[u] ; i !=-1 ; i = ne[i]){
            int j = e[i];
            //cout << "j: " << j <<endl;
            if((j == y&& u == x)||(j==x&&u==y)) continue;
            if(d[j] >dis + val[i]){
                d[j] = dis + val[i];
                if(x==-1&&y==-1) 
                fre[j] = u;
              
                  node root ={j,d[j]};
                  q.push(root);  
                
                
            }
        }
    }
    
}

int main(){
    scanf("%d%d", &n, &m);
    for(int i = 1 ; i  <= n ; i++){
       scanf("%lf%lf",&x[i],&y[i]);
    }
    memset(h, -1 ,sizeof h);
    for(int i =1 ; i <= m ; i++){
        int a , b;
        scanf("%d%d", &a, &b);
        double c = calc(x[a],y[a],x[b],y[b]);
        add(a,b,c);
        add(b,a,c);
    }
    dijstra(1,-1,-1);//不删边跑最短路保存路径。
    double ans = INF;
    for(int i = n ; i != 1 ; i = fre[i]){
        dijstra(1,fre[i],i);//删边跑.
        ans = min(ans,d[n]);
    }
    if(ans==INF)printf("-1\n");
    else
    printf("%.2lf\n",ans);
    return 0;
}

名欺！

关注

2
点赞
踩
0

收藏

觉得还不错? 一键收藏
1
评论
【洛谷】P1491 集合位置【次短路】两种解法

每次有大的活动，大家都要在一起“聚一聚”，不管是去好乐迪，还是避风塘，或者汤姆熊，大家都要玩的痛快。还记得心语和花儿在跳舞机上的激情与释放，还记得草草的投篮技艺是如此的高超，还记得狗狗的枪法永远是’S’……个点的坐标，其中第一个为野猫的出发位置，最后一个为大家的集合位置，并给出哪些位置点是相连的。野猫从出发点到达集合点，总会挑一条最近的路走，如果野猫没找到最近的路，他就会走第二近的路。野猫是公认的“路盲”，野猫自己心里也很清楚，每次都提前出门，但还是经常迟到，这点让大家很是无奈。个点的坐标，再往下的。
复制链接

扫一扫