7-50 畅通工程之局部最小花费问题 (35 分)

7-50 畅通工程之局部最小花费问题 (35 分)

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

输入格式:

输入的第一行给出村庄数目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

该题其实为求解最小生成树的题，利用Prim算法，在构造邻接矩阵的时候需要把已经建好的路的权值赋为0即可。 

#include<iostream>
using namespace std;
#define MAXN 110
#define INF 0xffffff

int AdjMatrix[MAXN][MAXN];
int lowcost[MAXN];          //lowcost[i] = -1 表示点i在v集中
int closest[MAXN];
int sumCost = 0;

void prim(int v,int n)
{
    int mincost;
    int k;
    for(int i = 1;i<=n;i++)
    {
        lowcost[i] = AdjMatrix[v][i];
        closest[i] = v;
    }
    lowcost[v] = -1;    //先将起始点加入v集
    for(int i = 0;i<n-1;i++)
    {
        mincost = INF;
        for(int j = 1;j<=n;j++)    //寻找最小权边
        {
            if(lowcost[j]!=-1&&lowcost[j]<mincost)
            {
                mincost = lowcost[j];
                k = j;
            }
        }
        sumCost+=mincost;
        lowcost[k] = -1;         //将k新加入到v集
        for(int j = 1;j<=n;j++)     //更新最短边集
        {
            if(AdjMatrix[k][j]<lowcost[j]&&j!=k)
            {
                lowcost[j] = AdjMatrix[k][j];
                closest[j] = k;
            }
        }
    }
}

int main()
{
    int n;
    cin>>n;
    for(int i = 0;i<MAXN;i++)
        for(int j = 0;j<MAXN;j++)
            AdjMatrix[i][j] = INF;       //邻接矩阵初始化
    for(int i = 0;i<n*(n-1)/2;i++)
    {
        int x,y,cost,flag;
        cin>>x>>y>>cost>>flag;
        if(flag)
            AdjMatrix[x][y] = AdjMatrix[y][x] = 0;
        else
            AdjMatrix[x][y] = AdjMatrix[y][x] = cost;
    }
    prim(1,n);
    cout<<sumCost<<endl;
    return 0;
}