畅通工程 HDU - 1863（图论，并查集，Kruskal算法）

Problem

省政府“畅通工程”的目标是使全省任何两个村庄间都可以实现公路交通（但不一定有直接的公路相连，只要能间接通过公路可达即可）。经过调查评估，得到的统计表中列出了有可能建设公路的若干条道路的成本。现请你编写程序，计算出全省畅通需要的最低成本。   

Input 

测试输入包含若干测试用例。每个测试用例的第1行给出评估的道路条数 N、村庄数目M ( < 100 )；随后的 N 
行对应村庄间道路的成本，每行给出一对正整数，分别是两个村庄的编号，以及此两村庄间道路的成本（也是正整数）。为简单起见，村庄从1到M编号。当N为0时，全部输入结束，相应的结果不要输出。

Output

对每个测试用例，在1行里输出全省畅通需要的最低成本。若统计数据不足以保证畅通，则输出“?”。 

Sample Input

3 3
1 2 1
1 3 2
2 3 4
1 3
2 3 2
0 100

Sample Output

3
?

Kruskal思想（求最小生成树）：

有两个集合，一个U（初始状态为空集），一个V（初始状态为边的全集），每次从集合V中选出最短的且不与集合U中现有的边形成回路的边，加入到集合U中。即可得到最小生成树。

如何判断会不会形成回路？

对一条边的两个端点执行unite(x,y)，如果x，y的根源父节点是同一个，说明会形成回路（因为此时x与y也连成边）。

如果x，y的根源父节点不是同一个，则加入集合U中不会形成回路，所以这条边可以被加到集合U中。

Kruskal算法中，每次加入一条不会与现有边形成回路的边，当边的数目=点的数目-1，则形成一个最小生成树，说明存在一条任意两点之间可达的路径。

#include <iostream>

using namespace std;
struct edge
{
    int a;
    int b;
    int value;
};
edge road[5000];
int village[105];
int flag;
void init(int n)
{
    for(int i=1;i<=n;i++)
    {
        village[i]=i;
    }
}
int cmp(const void*x,const void*y)//qsort函数中cmp方法的重载（结构体）
{
    return (*(struct edge*)x).value-(*(struct edge*)y).value;
}
int Find(int x)
{
    if(village[x]==x)
        return x;
    else return village[x]=Find(village[x]);
}
void unite(int x,int y)
{
    x=Find(x);
    y=Find(y);
    if(x!=y)
    {
        village[y]=x;
        flag=1;    //flag=1说明x，y不会形成回路。
    }

}
int main()
{
    int n,m,cnt,sum;
    while((cin>>n>>m)&&n)
    {
        init(n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d%d",&road[i].a,&road[i].b,&road[i].value);
        }
        qsort(road,n,sizeof(road[0]),cmp);//按边的权值从小到大进行排序
        sum=0,cnt=0;
        for(int i=1;i<=n;i++)
        {
            flag=0;
            unite(road[i].a,road[i].b);
            if(flag)
            {
                sum+=road[i].value;    //能被加到集合U中，更新总长度
                cnt++;    //更新边的数目
            }
        }
        if(cnt==m-1)    //若边的数目刚好等于点的数目-1，说明最小生成树能构成
            cout<<sum<<endl;
        else
            cout<<"?"<<endl;    //否则，无法生成最小生成树
    }
    return 0;
}

变形

畅通工程再续 HDU - 1875

相信大家都听说一个“百岛湖”的地方吧，百岛湖的居民生活在不同的小岛中，当他们想去其他的小岛时都要通过划小船来实现。现在政府决定大力发展百岛湖，发展首先要解决的问题当然是交通问题，政府决定实现百岛湖的全畅通！经过考察小组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!

这道题目和前面那题基本一样，但是由于题目没有给出每两个岛屿之间的距离，所以我们需要自己根据坐标算出距离。

#include <iostream>
#include <math.h>
#include<algorithm>
using namespace std;
struct edge
{
    int a;//第a个小岛
    int b;//第b个小岛
    double dist;//之间的距离
};
struct node
{
    int x;//小岛横坐标
    int y;//小岛纵坐标
    int par;//小岛父节点
};
edge road[5000];
node island[105];
int flag;
void init(int n)
{
    for(int i=0;i<n;i++)
    {
        island[i].par=i;//父节点假设是自己
    }
}
int Find(int x)
{
    if(island[x].par==x)
    return x;
    else return island[x].par=Find(island[x].par);

}
void unite(int x,int y)
{
    x=Find(x);
    y=Find(y);
    if(x!=y)
    {
        island[y].par=x;
        flag=1;
    }
}
double solvedist(int x1,int x2,int y1,int y2)
{
    return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
}
/*int cmp(const void *a,const void*b)
{
    return (*(struct edge*)a).dist-(*(struct edge*)b).dist;
}*/
int cmp(edge a,edge b)
{
    return a.dist<b.dist;
}
int main()
{
    int t,c,cnt,k;
    double dist,sum;
    cin>>t;
    while(t--)
    {
        cin>>c;
        for(int i=0;i<c;i++)
        {
            cin>>island[i].x>>island[i].y;
        }
        init(c);
        cnt=0;
        for(int i=0;i<c-1;i++)    //每两个岛屿两两求出它们之间的距离，若超出距离范围则直接略过这条路线，符合条件的加入road
        {
            for(int j=i+1;j<c;j++)
            {
                dist=solvedist(island[i].x,island[j].x,island[i].y,island[j].y);
                if(dist>=10&&dist<=1000)
                {
                    road[cnt].a=i;
                    road[cnt].b=j;
                    road[cnt].dist=dist;
                    cnt++;
                }
            }
        }
    //    qsort(road,cnt,sizeof(road[0]),cmp);
        sort(road,road+cnt,cmp);
        sum=0;k=0;
        for(int i=0;i<cnt;i++)
        {
            flag=0;
            unite(road[i].a,road[i].b);
            if(flag)
            {
                sum+=road[i].dist;
                k++;
            }
        }
        if(k==c-1)
        {
            printf("%.1lf\n",100*sum);
        }
        else printf("oh!\n");
    }
    return 0;
}