【解题报告】 HDU 1875 畅通工程再续 Kruskal最小生成树 一点关于浮点型在计算机中储存的分析

哎原以为一次ac，结果再在double比较大小的时候砸了，原来比较两个double的大小（不算等于的情况）则可以直接使用<=, <, >, >=比较，但是要比较相等的情况则必须两个相减 再与一个“正0”（0.00000001）或“负0”（-0.00000001）比较即可。
原因就是 : double型的有效数字是15～16位，数值范围-1.7×10^-308~1.7×10^308。
a是一个纯小数，它的位数就是有效位数,而k表示这个浮点数的放大因子，k可以范围很大，使得这个数变动范围很大. 具体原因参加这里 Double附上一个来自维基百科的说明图片

// HDU 1875 畅通工程再续
// 并查集 -> Kruskal最小生成树
// SHIT!在判断两个浮点型的时候又脑残了，想多了。直接比较就好了，要判相等再相减小于某个精度即可
// N 个城市如果连通
// 最少拥有 N-1 条道路, 最多拥有 N(N-1)/2条道路
// 前提是任何两个城市之间最多一条直达通道
//
/* test data
4
9
2 6
55 65 
40 5
51 65
100 200
103 353
516 198
400 300
200 250
			=91707.3

1
0 0 
			=0.0
2
10 10
20 20
			=1414.2
3
1 1
2 2
1000 1000
			= oh!
*/
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;

const double MIND = 10.000;
const double MAXD = 1000.000;
const int MAXV = 101;
int x[MAXV],y[MAXV];
struct Via{
	int a,b;double len;
}via[MAXV*MAXV]; // 按路径长度排序,距离为 via[i].len 的两个小岛为 小岛via[i].a 和 小岛via[i].b 

int father[MAXV];
int num[MAXV];

bool com(Via a, Via b){
	return a.len < b.len ;
}

double Plinelength(int x1,int y1,int x2,int y2){ //计算两点距离
	return sqrt((x1-x2) * (x1-x2) * 1.0 + (y1-y2) * (y1-y2));
}

int find(int x){
	if (father[x] == x) return x;
	return father[x] = find(father[x]);
}

void combination(int a, int b){
	father[b] = a;
	num[a] += num[b];
}

void MST(int n_via,int Ncity){  // 并查集套路
	double res=0;
	for (int i = 0; i < n_via; i++){
		int ra = find(via[i].a);
		int rb = find(via[i].b);
		if (ra != rb){
			combination(ra,rb);
			res += via[i].len;
			if (num[ra] >= Ncity){printf("%.1lf\n",res*100);return ;}
		}
	}
	cout << "oh!" << endl;
}

int buildt(int Ncity){ // 建图
	int k = 0;
	for (int i = 0; i < Ncity - 1; i++)
		for (int j = i + 1; j < Ncity; j++){
			double len = Plinelength(x[i],y[i],x[j],y[j]);
			if (len >= MIND && len <= MAXD){ // 满足要求则算可以连通
				via[k].a = i;	via[k].b = j;
				via[k++].len = len;
			}
		}
	return k; //返回 via 条数
}

int main(){
	freopen("in.txt","r",stdin);
	int NB;cin>>NB; // NB <= 200
	while(NB--){
		int Ncity;
		scanf("%d",&Ncity); // Ncity <= 100
		for (int i = 0 ; i < Ncity; i++){
			scanf("%d%d",&x[i],&y[i]);
			father[i] = i;num[i] = 1;
		}
		if (Ncity==1){cout <<"0.0"<<endl;continue;}
		int k_via = buildt(Ncity);
		sort(&via[0],&via[k_via],com); // 按桥长度排序
		MST(k_via,Ncity); //最小生成树 基于贪心思想
	}
	return 0;
}