HDU 1875 畅通工程再续

畅通工程再续

Time Limit: 1000ms

Memory Limit: 32768KB

This problem will be judged on HDU. Original ID: 1875
64-bit integer IO format: %I64d      Java class name: Main

 

 

相信大家都听说一个“百岛湖”的地方吧，百岛湖的居民生活在不同的小岛中，当他们想去其他的小岛时都要通过划小船来实现。现在政府决定大力发展百岛湖，发展首先要解决的问题当然是交通问题，政府决定实现百岛湖的全畅通！经过考察小组RPRush对百岛湖的情况充分了解后，决定在符合条件的小岛间建上桥，所谓符合条件，就是2个小岛之间的距离不能小于10米，也不能大于1000米。当然，为了节省资金，只要求实现任意2个小岛之间有路通即可。其中桥的价格为 100元/米。

 

Input

输入包括多组数据。输入首先包括一个整数T(T <= 200)，代表有T组数据。
每组数据首先是一个整数C(C <= 100),代表小岛的个数，接下来是C组坐标，代表每个小岛的坐标，这些坐标都是 0 <= x, y <= 1000的整数。

 

Output

每组输入数据输出一行，代表建桥的最小花费，结果保留一位小数。如果无法实现工程以达到全部畅通，输出”oh!”.

 

Sample Input

2
2
10 10
20 20
3
1 1
2 2
1000 1000

Sample Output

1414.2
oh!

Source

2008浙大研究生复试热身赛（2）——全真模拟

 

解题：最小生成树。这题比较适合用prim算法做。我处心积虑的用优先队列优化，终于变得更慢了。

 

这题比较好的体现了prim算法的特点。以前一直不懂prim，后来学了Dijkstra算法后，有点懂prim了，今天这道题目，让我更懂Prim了。

 

prim求这题的优点是它生成树的过程一直是连通的。一旦发现不合适的边，便oh!

 

还有预处理，我们是假设加入了1号顶点。故要done[1] = true;

 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <climits>
#include <vector>
#include <queue>
#include <cstdlib>
#include <string>
#include <set>
#include <stack>
#define LL long long
#define pii pair<int,int>
#define pdi pair<double,int>
#define INF 999999
using namespace std;
const int maxn = 110;
int x[maxn],y[maxn],n;
double d[maxn],mp[maxn][maxn];
bool done[maxn];
priority_queue< pdi,vector< pdi >, greater< pdi > >q;
double calc(int i,int j){
    double tmp = (x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);
    return sqrt(tmp);
}
void prim(){
    while(!q.empty()) q.pop();
    for(int i = 1; i <= n; i++){
        d[i] = mp[1][i];
        done[i] = false;
        q.push(make_pair(d[i],i));
    }
    done[1] = true;//注意此处
    double sum = 0;
    bool flag = true;
    while(!q.empty()){
        int u = q.top().second;
        double w = q.top().first;
        q.pop();
        if(done[u]) continue;
        done[u] = true;
        if(w >= INF){flag = false;break;}
        sum += w;
        for(int i = 1; i <= n; i++){
            if(!done[i] && d[i] > mp[u][i]){
                d[i] = mp[u][i];
                q.push(make_pair(d[i],i));
            }
        }
    }
    if(flag) printf("%.1f\n",sum*100);
    else puts("oh!");
}
int main() {
    int t,i,j;
    double tmp;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        for(i = 1; i <= n; i++)
            scanf("%d %d",x+i,y+i);
        for(i = 1; i <= n; i++){
            for(j = 1; j <= n; j++){
                tmp = calc(i,j);
                if(tmp < 10 || tmp > 1000) mp[i][j] = mp[j][i] = INF;
                else mp[i][j] = mp[j][i] = tmp;
            }
        }
        prim();
    }
    return 0;
}