[Codeforces 1245D] Shichikuji and Power Grid (最小生成树)

[Codeforces 1245D] Shichikuji and Power Grid (最小生成树)

题面

有n个城市,坐标为((x_i,y_i)),还有两个系数(c_i,k_i).在每个城市建立发电站需要费用(c_i).如果不建立发电站，要让城市通电，就需要与有发电站的城市连通。i与j之间连一条无向的边的费用是((k_i+k_j))*两个城市之间的曼哈顿距离.求让每个城市都通电的最小费用,并输出任意一个方案。

分析

把选每个点的代价转成虚拟原点到这个点的边,这个套路很常见，但在最小生成树题里还是第一次见到。

城市之间两两连边，边权按题目里提到的计算。然后建立一个虚拟源点，向每个城市(i)连边，边权为(c_i).对整个图跑一个最小生成树即可.

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define maxn 2000
using namespace std;
typedef long long ll;
int n; 
struct node{
	int x;
	int y;
	int c;
	int k;
}a[maxn+5];
inline ll dist(ll x1,ll y1,ll x2,ll y2){
	return abs(x1-x2)+abs(y1-y2);
} 

struct edge{
	int from;
	int to;
	ll len;
	int next;
}E[maxn*maxn+5];
int head[maxn+5];
int sz=0;
void add_edge(int u,int v,ll w){
	sz++;
	E[sz].from=u;
	E[sz].to=v;
	E[sz].len=w;
	E[sz].next=head[u];
	head[u]=sz;
} 
bool cmp(edge p,edge q){
	return p.len<q.len;
}

int fa[maxn+5];
int find(int x){
	if(fa[x]==x) return x;
	else return fa[x]=find(fa[x]);
}

vector< pair<int,int> >res;
ll kruskal(){
	ll ans=0;
	for(int i=1;i<=n+1;i++) fa[i]=i;
	sort(E+1,E+1+sz,cmp);
	for(int i=1;i<=sz;i++){
		int u=E[i].from,v=E[i].to;
		int fu=find(u),fv=find(v);
		if(fu!=fv){
			ans+=E[i].len;
			fa[fu]=fv;
			res.push_back(make_pair(E[i].from,E[i].to));
		} 
	}
	return ans;
}

vector<int>poi;//poi~ 
vector< pair<int,int> >edg;
int main(){
	int st;
	scanf("%d",&n);
	st=n+1;
	for(int i=1;i<=n;i++) scanf("%d %d",&a[i].x,&a[i].y);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i].c);
		add_edge(st,i,a[i].c);
	}
	for(int i=1;i<=n;i++) scanf("%d",&a[i].k);
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			add_edge(i,j,(a[i].k+a[j].k)*dist(a[i].x,a[i].y,a[j].x,a[j].y));
		}
	}
	printf("%I64d
",kruskal());
	for(auto ed : res){
		if(ed.first==st) poi.push_back(ed.second);
		else edg.push_back(ed);
	}
	printf("%d
",poi.size());
	for(int x : poi) printf("%d ",x);
	printf("
");
	printf("%d
",edg.size());
	for(auto x : edg) printf("%d %d
",x.first,x.second); 
}