【解题报告】 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 data492 655 65 40 551 65100 200103 353516 198400 300200 250=91707.310 0 =0.0210 1020 20=1414.231 12 21000 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 <= 200while(NB--){int Ncity;scanf("%d",&Ncity); // Ncity <= 100for (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;}