题意:
在沙漠中找一条道路,满足这样的条件: 使最高温度尽量低,如果有多条输出路程最短的一条。
分析: 只要二分最高温度,只留下满足温度限制的道路, 再求最短路即可。在该题的数据规模下,这样做的复杂度完全可以满足。
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
const int maxn= 110 ;
const int INF = 1e9 ;
const double eps = 1e-6 ;
int n , m;
struct HeapNode{
int u ;
double d;
bool operator < (const HeapNode & b)const { return d > b.d ;}
};
struct Edge{
int from , to ;
double dis , deg ;
};
vector<int>G[maxn] ;
Edge E[maxn*maxn*2] ;
int edge_cnt ;
void init()
{
edge_cnt = 0;
for(int i=0 ;i<=n ; i++) G[i].clear();
}
void add_edge(int u ,int v ,double dis , double deg)
{
E[++edge_cnt] = (Edge) { u ,v ,dis , deg} ;
G[u].push_back(edge_cnt);
}
int pre[maxn] ;
bool done[maxn] ;
double d[maxn] ;
void dijkstra(int S, int T , double limit)
{
memset(done , 0 ,sizeof(done));
memset(pre , 0 ,sizeof(pre)) ;
for(int i=0 ;i<=n ; i++) d[i] = INF ;
priority_queue<HeapNode>Q ;
d[S] = 0 ;
Q.push((HeapNode){S,0}) ;
while(!Q.empty())
{
int u = Q.top().u ; Q.pop() ;
if(done[u]) continue ;
done[u] = true ;
if(u == T) return ;
int sz = G[u].size() ;
for(int i=0 ;i<sz ;i++)
{
Edge & e = E[G[u][i]] ;
if(e.deg > limit) continue ;
if(d[e.to] > d[u] + e.dis)
{
d[e.to] = d[u] + e.dis ;
pre[e.to] = u ;
Q.push((HeapNode){e.to , d[e.to]}) ;
}
}
}
}
bool test(int S , int T , double limit)
{
dijkstra(S , T ,limit) ;
if(d[T] - INF < -eps ) return true ;
return false ;
}
int main()
{
int S , T ;
while(scanf("%d%d" ,&n ,&m) == 2)
{
scanf("%d%d" ,&S ,&T) ;
init() ;
int u ,v ;
double dist , degree ;
for(int i=0 ;i<m ;i++)
{
scanf("%d%d%lf%lf" ,&u ,&v ,°ree ,&dist) ;
add_edge(u ,v , dist , degree) ;
add_edge(v ,u , dist , degree) ;
}
double L = 19.9 , H = 50.1 ;
while(H-L > 0.01)
{
double mid = (H+L) / 2.0 ;
if(test(S , T , mid)) H = mid ;
else L = mid ;
}
test(S ,T , H) ;
vector<int>ans ;
for(int p=T;;)
{
ans.push_back(p) ;
if(p == S) break ;
p = pre[p] ;
}
printf("%d" , ans[ans.size()-1]) ;
for(int i=ans.size()-2 ;i>=0 ;i--)
{
printf(" %d" , ans[i]);
}
printf("\n") ;
printf("%.1f %.1f\n" , d[T] , (H+L)/2.0) ;
}
return 0;
}