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

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

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

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

思路
关键词：最小生成树（对于已经存在道路联通的两点，可设置其连接代价为0，不连通的连接代价“无限大”）

#include<stdlib.h>
#include<stdio.h>

/*
 评测结果 时间    结果  得分  题目  编译器     用时（ms）  内存（MB）  用户
2016-04-18 16:57    答案正确    30  08-图7   gcc     41  5   569985011
测试点结果 测试点   结果  得分/满分   用时（ms）  内存（MB）  要点提示
测试点1    答案正确    15/15   1   1   sample换数字，各种回路判断
测试点2    答案正确    2/2     2   1   M<N-1，不可能有生成树
测试点3    答案正确    2/2     1   1   M达到N-1，但是图不连通
测试点4    答案正确    5/5     41  5   最大N和M，连通
测试点5    答案正确    6/6     14  5   最大N和M，不连通
*/
struct data {
    int VerNode,Coast;
};
typedef struct LNode* Vertex;
struct LNode {
    struct data a[1001];
    int Left,Right;
};
void InsertV(Vertex,int ,int,int);
int Prim(Vertex ,int );
Vertex CreatV(int);
int Comp(const void*a,const void*b) {
    struct data x=*(struct data *)a;
    struct data y=*(struct data *)b;
    return x.Coast - y.Coast ;
}

int main() {
    int N;
    scanf("%d",&N);
    Vertex V=CreatV(N+1);

    int from,to,lenth,flag;
    for(int i=0; i<N*(N-1)/2; i++) {
        scanf("%d%d%d%d",&from,&to,&lenth,&flag);
        if(flag==0)
        InsertV(V,from,to,lenth);
        else InsertV(V,from,to,0);
    }
    for(int i=1; i<=N; i++) {
        qsort(V[i].a,V[i].Right-V[i].Left,sizeof(struct data),Comp);
//      printf("%d",i);
//      for(int j=V[i].Left;j<V[i].Right ;j++){
//          printf("{-%d=%d}",V[i].a[j].VerNode,V[i].a[j].Coast  );
//      }printf("\n");
    }

    printf("%d",Prim(V,N));
    return 0;
}

int Prim(Vertex V,int N) {
    int right=0,left=0,Node[1001],Coast=0;
    for(int i=1; i<=N; i++) {
        if(V[i].Right >V[i].Left) {
            if(right) {
                if(V[i].a[0].Coast <V[Node[0]].a[0].Coast) {
                    Node[0]=i;
                }
            } else {
                Node[right++]=i;
            }

        }
    }
    Coast+=V[Node[0]].a[0].Coast ;
    Node[right++]=V[Node[0]].a[0].VerNode;
    V[Node[0]].Left ++;
//  printf("-%d-",right);
    while(right<N) {
        int flag=-1;
        for(int i=0; i<right; i++) {
            if(V[Node[i]].Left<V[Node[i]].Right) {
                if(flag!=-1) {
                    if(V[Node[i]].a[V[Node[i]].Left].Coast <V[Node[flag]].a[V[Node[flag]].Left ].Coast)
                        flag=i;
                } else {
                    flag=i;
                }
            }
        }
        if(-1==flag)break;
        int i;
        for(i=0;i<right;i++){
            if(V[Node[flag]].a[V[Node[flag]].Left].VerNode==Node[i])break;
        }
        if(i==right){
            Node[right++]=V[Node[flag]].a[V[Node[flag]].Left].VerNode;
            Coast+=V[Node[flag]].a[V[Node[flag]].Left].Coast ;
        }
        V[Node[flag]].Left++;
    }
    for(int i=0;i<right;i++){
//      printf("[%d]",Node[i]);
    }
if(right>=N)return Coast;
    else return -1;
}

Vertex CreatV(int num) {
    Vertex V=(Vertex)malloc(sizeof(struct LNode)*num);
    for(int i=0; i<num; i++) {
        V[i].Left =0;
        V[i].Right =0;
    }

    return V;
}
void InsertV(Vertex V,int from,int to,int lenth) {
    V[from].a[V[from].Right].Coast =lenth;
    V[from].a[V[from].Right++].VerNode =to;
    V[to].a[V[to].Right].Coast =lenth;
    V[to].a[V[to].Right++].VerNode =from;
}