灌水(生成树)

Farmer John已经决定把水灌到他的n(1<=n<=300)块农田,农田被数字1到n标记。把一块土地进行灌水有两种方法,从其他农田饮水,或者这块土地建造水库。建造一个水库需要花费wi(1<=wi<=100000),连接两块土地需要花费Pij(1<=pij<=100000,pij=pji,pii=0). 计算Farmer John所需的最少代价。

Input

*第一行:一个数n

*第二行到第n+1行:第i+1行含有一个数wi

*第n+2行到第2n+1行:第n+1+i行有n个被空格分开的数,第j个数代表pij。

Output

*第一行:一个单独的数代表最小代价.

Sample Input

4

5

4

4

3

0 2 2 2

2 0 3 3

2 3 0 4

2 3 4 0

Sample Output

9

输出详解:

Farmer John在第四块土地上建立水库,然后把其他的都连向那一个,这样就要花费3+2+2+2=9

解析:

本题既有边权,又有点权,如何将它们统一起来呢?不妨设置一个超级水源点,将它与各个点连边,代价为在那个点修水库所需的代价。接下来就可以使用生成树算法进行计算最小代价了。

代码如下:

#include<bits/stdc++.h>
using namespace std;
int n,m1,m2,aa;
int v[310][310],lowest[310],z=0;
bool d[310];
void mst()
{
    int ans=0;
    memset(d,true,sizeof(d));
    for(int i=1;i<=n;i++)
        lowest[i]=2147483647;
    lowest[1]=0;
    while(true)
    {
        m1=2147483647;
        m2=0;
        for(int i=1;i<=n;i++)
            if(d[i]==true&&lowest[i]<m1)
            {
                m2=i;
                m1=lowest[i];
            }
        ans=ans+m1;
        d[m2]=false;
        z++;
        if(z==n)
            break;
        for(int i=1;i<=n;i++)
            if(d[i]==true&&v[m2][i]<lowest[i])
                lowest[i]=v[m2][i]; 

    }
    cout<<ans<<endl;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)//构造一个n+1点,从它引边到每一个点,边权值为建每一个点的代价
    {
        scanf("%d",&v[i][n+1]);
        v[n+1][i]=v[i][n+1];
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            scanf("%d",&v[i][j]);
    n++;        
    mst();      
    //cout<<ans;
    return 0;
}

或者用kruskal来做:

#include<bits/stdc++.h>
using namespace std;
const int N = 100000;
int n,m;
int cnt=0;
int fa[N*2];
int change[N];
struct node {
	int x,y,val;
} a[N*2];
int ans;
int find(int x) {
	if(x!=fa[x]) fa[x]=find(fa[x]);
	return fa[x];
}
void unite(int x,int y) {
	int p=find(x),q=find(y);
	if(p!=q) fa[p]=q;
}
bool cmp(node x,node y) {
	return x.val<y.val;
}
void Kruskal() {
	int tot=0;
	for(int i=1; i<=cnt; i++) {
		if(find(a[i].x)==find(a[i].y))
			continue;
		unite(a[i].x,a[i].y);
		ans+=a[i].val;
	}
}
int main() {
	cin>>n;
	for(int i=1; i<=2*n; i++) fa[i]=i;
	for(int i=1; i<=n; i++)
		cin>>change[i];
	for(int i=1; i<=n; i++) {
		for(int j=1; j<=n; j++) {
			cnt++;
			a[cnt].x=i;
			a[cnt].y=j;
			cin>>a[cnt].val;
		}
	}
	n++;
	m=n-1;
	for(int i=1; i<=m; i++) {
		cnt++;
		a[cnt].x=i;
		a[cnt].val=change[i];
	}
	sort(a+1,a+1+cnt,cmp);
	Kruskal();
	printf("%d\n",ans);

	return 0;
}

  • 3
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值