uva 10816（二分，最短路）

题意：

在沙漠中找一条道路，满足这样的条件： 使最高温度尽量低，如果有多条输出路程最短的一条。

分析： 只要二分最高温度，只留下满足温度限制的道路， 再求最短路即可。在该题的数据规模下，这样做的复杂度完全可以满足。

代码：

#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;
}