小白算法练习 安慰奶牛 最小生成树 lanqiao prim

问题描述

Farmer John变得非常懒，他不想再继续维护供奶牛之间供通行的道路。道路被用来连接N个牧场，牧场被连续地编号为1到N。每一个牧场都是一个奶牛的家。FJ计划除去P条道路中尽可能多的道路，但是还要保持牧场之间 的连通性。你首先要决定那些道路是需要保留的N-1条道路。第j条双向道路连接了牧场Sj和Ej(1 <= Sj <= N; 1 <= Ej <= N; Sj != Ej)，而且走完它需要Lj的时间。没有两个牧场是被一条以上的道路所连接。奶牛们非常伤心，因为她们的交通系统被削减了。你需要到每一个奶牛的住处去安慰她们。每次你到达第i个牧场的时候(即使你已经到过)，你必须花去Ci的时间和奶牛交谈。你每个晚上都会在同一个牧场(这是供你选择的)过夜，直到奶牛们都从悲伤中缓过神来。在早上 起来和晚上回去睡觉的时候，你都需要和在你睡觉的牧场的奶牛交谈一次。这样你才能完成你的 交谈任务。假设Farmer John采纳了你的建议，请计算出使所有奶牛都被安慰的最少时间。

输入格式

第1行包含两个整数N和P。

接下来N行，每行包含一个整数Ci。

接下来P行，每行包含三个整数Sj, Ej和Lj。

输出格式

输出一个整数, 所需要的总时间(包含和在你所在的牧场的奶牛的两次谈话时间)。

样例输入

5 7
10
10
20
6
30
1 2 5
2 3 5
2 4 12
3 4 17
2 5 15
3 5 6

样例输出

176

数据规模与约定

5 <= N <= 10000，N-1 <= P <= 100000，0 <= Lj <= 1000，1 <= Ci <= 1,000。

Code

#include<iostream>
#include<vector>
#include<map>
using namespace std;
const int Max=10008;
const int INF=0x3f3f3f3f;
struct Node{
	int to;
	int weight;
};
vector<Node>vec[Max];
map<int,int>ma;
int vis[Max]={0};
int dist[Max]={0};
int pre[Max]={0};
void init(int f,int e,int w){
	Node node;
	node.to=e;
	node.weight=w;
	vec[f].push_back(node);
	Node node1;
	node1.to=f;
	node1.weight=w;
	vec[e].push_back(node1);
}
void prim(int d,int b,int u){
	vis[u]=1;
	pre[u]=-1;
	for(int i=1;i<=d-1;i++)
	{
		int Min=INF;
		int k=0;
		for(int j=1;j<=d;j++)
		{
			if(vis[j]==0 && dist[j]!=-1 && dist[j]<Min)
			{
				Min=dist[j];
				k=j;
			}
		}
		vis[k]=1;
		//松弛
		for(int j=0;j<vec[k].size();j++)
		{
			if(vis[vec[k][j].to]==0 && dist[vec[k][j].to]>vec[k][j].weight)
			{
				dist[vec[k][j].to]=vec[k][j].weight;
				pre[vec[k][j].to]=k;
			}
		} 
	}
	return;
}
int main(){
	
//	freopen("E://1001.in","r",stdin);
	//input & init
	int d,b;//点，边
	cin>>d>>b;
	int tree[Max];
	for(int i=1;i<=d;i++)
	{
		cin>>tree[i];
		ma.insert(pair<int,int>(tree[i],i));
		dist[i]=INF;
	}
	int Min=ma.begin()->second;	
	for(int i=1;i<=b;i++)
	{
		int f,e,w;
		cin>>f>>e>>w;
		init(f,e,2*w+tree[f]+tree[e]);
	}
	dist[Min]=0;
	for(int i=0;i<vec[Min].size();i++)
	{
		dist[vec[Min][i].to]=vec[Min][i].weight;
	}
	for(int i=1;i<=d;i++)
	{
		pre[i]=Min;
	}
	//prim
	prim(d,b,Min);
	//output
	int sum=0;	
	for(int i=1;i<=d;i++)
	{
		sum+=dist[i];
	}
	sum+=tree[Min];
	cout<<sum<<endl;
} 
const int Max=10008;
const int INF=0x3f3f3f3f;
struct Node{
	int to;
	int weight;
};
vector<Node>vec[Max];
map<int,int>ma;
int vis[Max]={0};
int dist[Max]={0};
int pre[Max]={0};
void init(int f,int e,int w){
	Node node;
	node.to=e;
	node.weight=w;
	vec[f].push_back(node);
	Node node1;
	node1.to=f;
	node1.weight=w;
	vec[e].push_back(node1);
}
void prim(int d,int b,int u){
	vis[u]=1;
	pre[u]=-1;
	for(int i=1;i<=d-1;i++)
	{
		int Min=INF;
		int k=0;
		for(int j=1;j<=d;j++)
		{
			if(vis[j]==0 && dist[j]!=-1 && dist[j]<Min)
			{
				Min=dist[j];
				k=j;
			}
		}
		vis[k]=1;
		//松弛
		for(int j=0;j<vec[k].size();j++)
		{
			if(vis[vec[k][j].to]==0 && dist[vec[k][j].to]>vec[k][j].weight)
			{
				dist[vec[k][j].to]=vec[k][j].weight;
				pre[vec[k][j].to]=k;
			}
		} 
	}
	return;
}
int main(){
	
//	freopen("E://1001.in","r",stdin);
	//input & init
	int d,b;//点，边
	cin>>d>>b;
	int tree[Max];
	for(int i=1;i<=d;i++)
	{
		cin>>tree[i];
		ma.insert(pair<int,int>(tree[i],i));
		dist[i]=INF;
	}
	int Min=ma.begin()->second;	
	for(int i=1;i<=b;i++)
	{
		int f,e,w;
		cin>>f>>e>>w;
		init(f,e,2*w+tree[f]+tree[e]);
	}
	dist[Min]=0;
	for(int i=0;i<vec[Min].size();i++)
	{
		dist[vec[Min][i].to]=vec[Min][i].weight;
	}
	for(int i=1;i<=d;i++)
	{
		pre[i]=Min;
	}
	//prim
	prim(d,b,Min);
	//output
	int sum=0;	
	for(int i=1;i<=d;i++)
	{
		sum+=dist[i];
	}
	sum+=tree[Min];
	cout<<sum<<endl;
} 

我一开是想写最小生成树的，结果写了dijkstra最短路径，话说这两个代码简直一毛一样。。后来才发现的，xjb乱写写了好几个小时。