2-2 畅通工程之局部最小花费问题

2-2 畅通工程之局部最小花费问题

分数 30
作者 DS课程组
单位 浙江大学

某地区经过对城镇交通状况的调查，得到现有城镇间快速道路的统计数据，并提出“畅通工程”的目标：使整个地区任何两个城镇间都可以实现快速交通（但不一定有直接的快速道路相连，只要互相间接通过快速路可达即可）。现得到城镇道路统计表，表中列出了任意两城镇间修建快速路的费用，以及该道路是否已经修通的状态。现请你编写程序，计算出全地区畅通需要的最低成本。

输入格式:

输入的第一行给出村庄数目N (1≤N≤100)；随后的N(N−1)/2行对应村庄间道路的成本及修建状态：每行给出4个正整数，分别是两个村庄的编号（从1编号到N），此两村庄间道路的成本，以及修建状态 — 1表示已建，0表示未建。

输出格式:

输出全省畅通需要的最低成本。

输入样例:

4
1 2 1 1
1 3 4 0
1 4 1 1
2 3 3 0
2 4 2 1
3 4 5 0

输出样例:

3

代码长度限制
16 KB
时间限制
400 ms
内存限制
64 MB
C++ (g++)

思路：

由题可知让我们求最小生成树。
计算出全地区畅通需要的最低成本，所以修好的路就不用再花费了，输出的就是没连通路的花费，所以查找没有再那个集合的村子然后找到最小的权重值和，即可(最小生成树）。

AC代码：

法一：Prim算法

#include<bits/stdc++.h>
using namespace std;
#define inf 999999
int n,e[101][101],vis[101],dis[101];
int main()
{
    int a,b,price,statue;//村庄a和b，成本，修没修
    cin>>n;
    for(int i=0;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
            e[i][j]=e[j][i]=inf;//先把值设置为最大
        e[i][i]=0;
    }//初始化
    for(int i=1;i<=n*(n-1)/2;i++)
    {
        cin>>a>>b>>price>>statue;
        if(statue) e[a][b]=e[b][a]=0;//如果已经创建，则赋值为0
        else e[a][b]=e[b][a]=price;//存储
    }
    vis[1]=1;//判断是否访问过
    for(int i=0;i<=n;i++)
        dis[i]=e[1][i];//存值普利姆算法
    int count=1;//判断是否访问完
    int sum=0;//输出最小值
    while(count<n){
        int min=inf,d=0;
        for(int i=1;i<=n;i++)
            if(vis[i]==0&&dis[i]<min)
                min=dis[i],d=i;//更新最小
        vis[d]=1;count++;sum+=min;//更新值
        for(int i=0;i<=n;i++)
            if(vis[i]==0&&dis[i]>e[d][i])
                dis[i]=e[d][i];//更新最小值
    }
    cout<<sum;
}

法二：kruskal算法

#include<bits/stdc++.h>
using namespace std;
int n;
typedef struct{
    int from;//起点
    int next;//终点
    int w;//权值
    int f;//状态
}node;
int father[1000];
node edge[10000];
struct rule{
    bool operator()(node a,node b)
    {
        if(a.f!=b.f) return a.f>b.f;
        else{
            return a.w<b.w;
        }
    }
};
int find(int x)
{
    if(father[x]==x) return x;
    else{
        return father[x]=find(father[x]);
    }
}
int main()
{
    cin>>n;
    for(int i=0;i<n*(n-1)/2;i++)
    {
        cin>>edge[i].from>>edge[i].next>>edge[i].w>>edge[i].f;
    }
    for(int i=0;i<=n;i++) father[i]=i;
    //开始排序
    sort(edge,edge+n,rule());
    int sum=0;//边数
    //将已经修成的先放到一个集合中，类似于贪心算法
    for(int i=0;i<n*(n-1)/2;i++){
        int x=find(edge[i].from);
        int y=find(edge[i].next);
        if(x!=y&&edge[i].f==1){
            father[x]=y;
            sum++;
        }
    }
    int ans=0;
    for(int i=0;i<n*(n-1)/2;i++){
        int x=find(edge[i].from);
        int y=find(edge[i].next);
        if(x!=y&&edge[i].f==0){
            sum++;
            father[x]=y;
            ans+=edge[i].w;
            if(sum==n-1) break;
        }
    }
    cout<<ans;
    return 0;
}