继续畅通工程
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 20428 Accepted Submission(s): 8801
Problem Description
省政府“畅通工程”的目标是使全省任何两个村庄间都可以实现公路交通(但不一定有直接的公路相连,只要能间接通过公路可达即可)。现得到城镇道路统计表,表中列出了任意两城镇间修建道路的费用,以及该道路是否已经修通的状态。现请你编写程序,计算出全省畅通需要的最低成本。
Input
测试输入包含若干测试用例。每个测试用例的第1行给出村庄数目N ( 1< N < 100 );随后的 N(N-1)/2 行对应村庄间道路的成本及修建状态,每行给4个正整数,分别是两个村庄的编号(从1编号到N),此两村庄间道路的成本,以及修建状态:1表示已建,0表示未建。
当N为0时输入结束。
Output
每个测试用例的输出占一行,输出全省畅通需要的最低成本。
Sample Input
3
1 2 1 0
1 3 2 0
2 3 4 0
3
1 2 1 0
1 3 2 0
2 3 4 1
3
1 2 1 0
1 3 2 1
2 3 4 1
0
Sample Output
3
1
0
Author
ZJU
Source
浙大计算机研究生复试上机考试-2008年
其实就是找最小生成树咯。但是有一些变化,因为有些路已经修好,另一些没修好,那么修好的路不管三七二十一直接加入到我们要的树当中,因为他不要钱肯定是更优的,在这基础上我们再去跑kruskal把剩下的边试一下就阔以了。
#include "iostream"
#include "cstring"
#include "algorithm"
#include "cstdio"
#include "string.h"
using namespace std;
int pre[10005];
typedef struct node
{
int from;
int to;
int cost;
int status;
}node;
node mmap[10005];
bool cmp(node a,node b)
{
return a.cost<b.cost;
}
void initial()
{
for(int i=1;i<=10000;i++)
{
pre[i]=i;
}
}
int find(int now)
{
if(pre[now]!=now)
pre[now]=find(pre[now]);
return pre[now];
}
void join(int a,int b)
{
int p=find(a);
int q=find(b);
if(p!=q)
{
pre[p]=q;
}
}
int main()
{
int n;
while(scanf("%d",&n)&&n!=0)
{
initial();
for(int i=1;i<=(n*(n-1))/2;i++)
{
scanf("%d%d%d%d",&mmap[i].from,&mmap[i].to,&mmap[i].cost,&mmap[i].status);
if(mmap[i].status==1)
{
join(mmap[i].from,mmap[i].to);
//pre[mmap[i].from]=mmap[i].to;
}
}
//kruskal
sort(mmap+1,mmap+1+(n*(n-1))/2,cmp);
int ans=0;
for(int i=1;i<=(n*(n-1))/2;i++)
{
if(find(mmap[i].from)!=find(mmap[i].to))
{
join(mmap[i].from,mmap[i].to);
ans+=mmap[i].cost;
}
}
printf("%d\n",ans);
}
}
然后又用prim写了一波。prim写的话就把已经修好的路的权当作0这样,再跑prim就阔以了
#include "iostream"
#include "cstring"
#include "string.h"
#include "cstdio"
#define INF 9999999
using namespace std;
int map[105][105];
bool vis[105];
double N;
int low[105];
int prim()
{
int res=0,pos,min;
for(int i=1;i<=N;i++)
low[i]=map[1][i];
vis[1]=true;
for(int i=1;i<=N-1;i++)
{
min=INF;
pos=-1;
for(int j=1;j<=N;j++)
{
if(!vis[j]&&min>low[j])
{
min=low[j];
pos=j;
}
}
if(pos==-1)
return -1;
vis[pos]=true;
res=res+min;
for(int j=1;j<=N;j++)
{
if(!vis[j]&&low[j]>map[pos][j])
low[j]=map[pos][j];
}
}
return res;
}
int main()
{
while(scanf("%lf",&N)!=EOF)
{
if(N==0)
return 0;
for(int m=1;m<=104;m++)
for(int n=1;n<=104;n++)
map[m][n]=INF;
for(int m=1;m<=104;m++)
{
low[m]=INF;
map[m][m]=0;
}
memset(vis,false,sizeof(vis));
for(int i=1;i<=(N*(N-1))/2;i++)
{
int u,v,w,s;
scanf("%d%d%d%d",&u,&v,&w,&s);
if(s==1)
map[u][v]=map[v][u]=0;
else
map[u][v]=map[v][u]=w;
}
cout<<prim()<<endl;
}
}