【UVA 10816】 Travel in Desert （最小瓶颈树+最短路）

【题意】

　　有n个绿洲， m条道路，每条路上有一个温度，和一个路程长度，从绿洲s到绿洲t，求一条道路的最高温度尽量小， 如果有多条， 选一条总路程最短的。

Input
Input consists of several test cases. Your program must process all of them.
The first line contains two integers N and E (1 ≤ N ≤ 100; 1 ≤ E ≤ 10000) where N represents the
number of oasis and E represents the number of paths between them. Next line contains two distinct
integers S and T (1 ≤ S, T ≤ N) representing the starting point and the destination respectively. The
following E lines are the information the group gathered. Each line contains 2 integers X, Y and 2 real
numbers R and D (1 ≤ X, Y ≤ N; 20 ≤ R ≤ 50; 0 < D ≤ 40). It means there is a path between X and
Y , with length D km and highest temperature RoC. Each real number has exactly one digit after the
decimal point. There might be more than one path between a pair of oases.
Output
Print two lines for each test case. The first line should give the route you find, and the second should
contain its length and maximum temperature.
Sample Input
6 9
1 6
1 2 37.1 10.2
2 3 40.5 20.7
3 4 42.8 19.0
3 1 38.3 15.8
4 5 39.7 11.1
6 3 36.0 22.5
5 6 43.9 10.2
2 6 44.2 15.2
4 6 34.2 17.4
Sample Output
1 3 6
38.3 38.3

【分析】

　　我们可以先求出最大温度的最小值，然后把小于等于这个温度的边加进图中跑最短路。

　　最短路就不说了，现在就是要求最小瓶颈路。

　　最小瓶颈路有两个方法，

　　1、二分+BFS

　　　　二分之后沿着小于等于这个温度的边走，只需判断能否走到终点，所以是mlogn的。

　　2、

　　　　但其实可以nlogn把图上所有两点的最小瓶颈路求出来，就是求出最小瓶颈树，那么两点之间的唯一路径就是他们的最小瓶颈路。

　　　　而最小生成树就是一个最小瓶颈树。

　　[其实这个，我也不是很会证明的说- -谁能告诉我- -]

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<queue>
 #include<cmath>
 using namespace std;
 #define Maxn 110
 #define Maxm 10010
 #define INF 0xfffffff

 int n,m,st,ed;

 struct node
 {
     int x,y,c,d;
     int next;
 }tt[Maxm],t[Maxm*];
 int len,first[Maxn];

 bool cmp(node x,node y) {return x.c<y.c;}
 // double mymax(double x,double y) {return x>y?x:y;}
 int mymax(int x,int y) {return x>y?x:y;}

 int fa[Maxn];
 int ffa(int x)
 {
     if(fa[x]!=x) fa[x]=ffa(fa[x]);
     return fa[x];
 }

 void ins(int x,int y,int c)
 {
     t[++len].x=x;t[len].y=y;t[len].c=c;
     t[len].next=first[x];first[x]=len;
 }

 int mx[Maxn];
 void dfs(int x,int f)
 {
     for(int i=first[x];i;i=t[i].next) if(t[i].y!=f)
     {
         int y=t[i].y;
         mx[y]=mymax(mx[x],t[i].c);
         dfs(y,x);
     }
 }

 queue<int > q;
 int pre[Maxn],dis[Maxn];
 bool inq[Maxn];
 void spfa()
 {
     while(!q.empty()) q.pop();
     // for(int i=1;i<=n;i++) dis[i]=INF;
     memset(dis,,sizeof(dis));
     memset(inq,,sizeof(inq));
     q.push(ed);inq[ed]=;dis[ed]=;
     while(!q.empty())
     {
         int x=q.front();
         for(int i=first[x];i;i=t[i].next)
         {
             int y=t[i].y;
             if(dis[y]>dis[x]+t[i].c)
             {
                 dis[y]=dis[x]+t[i].c;
                 pre[y]=x;
                 if(!inq[y])
                 {
                     q.push(y);
                     inq[y]=;
                 }
             }
         }
         inq[x]=;q.pop();
     }
     if(dis[st]>=INF-) return;
     int now=st;
     while(now!=ed)
     {
         printf("%d ",now);
         now=pre[now];
     }
     printf("%d\n",ed);
     printf("%.1lf %.1lf\n",dis[st]*1.0/,mx[st]*1.0/);
 }

 int main()
 {
     while(scanf("%d%d",&n,&m)!=EOF)
     {
         scanf("%d%d",&st,&ed);
         for(int i=;i<=m;i++)
         {
             double c,d;
             scanf("%d%d%lf%lf",&tt[i].x,&tt[i].y,&c,&d);
             tt[i].c=(int)round(c*);tt[i].d=(int)round(d*);
         }
         sort(tt+,tt++m,cmp);
         for(int i=;i<=n;i++) fa[i]=i;
         int cnt=;
         memset(first,,sizeof(first));
         len=;
         for(int i=;i<=m;i++)
         {
             if(ffa(tt[i].x)!=ffa(tt[i].y))
             {
                 fa[ffa(tt[i].x)]=ffa(tt[i].y);
                 cnt++;
                 ins(tt[i].x,tt[i].y,tt[i].c);
                 ins(tt[i].y,tt[i].x,tt[i].c);
             }
             if(cnt==n-) break;
         }
         mx[ed]=;
         dfs(ed,);
         len=;
         memset(first,,sizeof(first));
         for(int i=;i<=m;i++) if(tt[i].c<=mx[st])
         {
             ins(tt[i].x,tt[i].y,tt[i].d);
             ins(tt[i].y,tt[i].x,tt[i].d);
         }
         spfa();
     }
     return ;
 }

2016-11-01 15:57:34