题面链接:https://codeforces.com/problemset/problem/1245/D

题意大概是给你一些城市的坐标,可以在城市中建立发电站,也可以让某个城市和已经建好发电站的城市连接,保证在这两种操作下使得所有的城市供电,在城市建发电站需要花费Ci,城市a和城市b连接需要花费(|Xa-Xb| + |Ya-Yb| )*(Ka+Kb),求最小的花费让所有城市通电。

思路:首先建立一个源点,连接各个城市,边权就是建发电站费用Ci,若选择此边则表示是在该城市建了发电站。

剩下的城市互相按照题意要求建边,建完全图。

最终模型就抽象成求最小生成树,用克鲁斯卡尔算法跑一下即可

AC代码:

#include<iostream>
#include<vector>
#include<set>
#include<cstdio>
#include<algorithm>
#define maxn 2005
using namespace std;
typedef long long ll;
struct node{
int from,to;
ll cost;
}e[maxn*maxn];//存储边的信息
int father[maxn];
int cnt,N,city;
vector<int> sta;//存建发电站的城市
vector<int> citya,cityb;//存储a b城市表示ab之间有电线
bool cmp(node a,node b){//比较函数
return a.cost < b.cost ;
}
void init(){//初始化
cnt = 0;
for(int i = 0;i<=N;i++){
father[i] = i;
}
}
int find(int x){//寻找根
if(x == father[x]) return x;
return father[x] = find(father[x]);
}
bool same(int x,int y){//判断是否在一个集合
return find(x)==find(y);
}
void unionSet(int x,int y){ //并查集合并
int u = find(x),v = find(y);
if(u==v) return ;
father[u] = v;
}
long long kruskal(){ //求最小生成树
ll res = 0;
std::sort(e,e+cnt,cmp);
for(int i = 0;i<cnt;i++){
if(same(e[i].from,e[i].to )) continue;
unionSet(e[i].from ,e[i].to );
if(e[i].from == 0 || e[i].to == 0){//判一下是否在城市建了发电站(利用源点判断)
if(e[i].from == 0) sta.push_back(e[i].to );
else sta.push_back(e[i].from);
}
else{//如果不是,存一些a,b两座连接的城市
citya.push_back(e[i].from ),cityb.push_back(e[i].to);
}
res+=e[i].cost ;
}
return res;
}
void addedge(int u,int v,ll C){ //建边操作
e[cnt].from = u,e[cnt].to = v,e[cnt].cost = C;
cnt++;
}
int main(){
cin>>N;
init();
vector<long long> vx,vy;
vector<long long> c,k;
for(int i = 0;i<N;i++){
int x,y;cin>>x>>y;
vx.push_back(x),vy.push_back(y);
}
for(int i = 0;i<N;i++){
ll tc;cin>>tc;
c.push_back(tc);
}
for(int i = 0;i<N;i++){
ll tk;cin>>tk;
k.push_back(tk);
}
for(int i = 0;i<N;i++){//建立源点,和各个城市连接
addedge(0,i+1,c[i]);
}
for(int i = 0;i<N;i++){//城市和城市之间建完全图
for(int j = i+1;j<N;j++){
ll c = (abs(vx[i]-vx[j])+abs(vy[i]-vy[j]))*(k[i]+k[j]);
addedge(i+1,j+1,c);
}
}
ll ans = kruskal();
cout<<ans<<endl;
cout<<sta.size()<<endl;
for(int i = 0;i<sta.size();i++){
if(i==0) cout<<sta[i];
else cout<<" "<<sta[i];
}
cout<<endl;
cout<<citya.size()<<endl;
for(int i = 0;i<citya.size() ;i++){
cout<<citya[i]<<" "<<cityb[i]<<endl;
}
return 0;
}

最新文章

  1. Alpha阶段总结
  2. 发邮件 和 excel导出中文文件名
  3. 混合模式程序集是针对“v2.0.50727”版的运行时生成的,在没有配置其他信息的情况下,无法在 4.0 运行时中加载该程序集。
  4. Delphi又要换东家了
  5. C++常量(C++数值常量、字符串常量、符号常量)
  6. jsp片段
  7. IIS 发布网站到外网
  8. java.lang.NoSuchMethodError: org.springframework.beans.factory.annotation.InjectionMetadata.&lt;init&gt;(L
  9. Android4.2.2由于越来越多的物理按键(frameworks)
  10. 关于PHPAPI ZEND_API TSRM_API宏的定义
  11. Python之路:堡垒机实例以及数据库操作
  12. 支付宝WAP支付接口开发
  13. 一个基于php+mysql的外卖订餐网站(带源码)
  14. Hibernate初体验及简单错误排除
  15. 2019-04-15 Python之利用matplotlib和numpy的简单绘图
  16. 下拉框click事件与搜索框blur事件的爱恨纠葛
  17. Java 中变量初始化、子类和父类构造器调用的顺序
  18. python numpy 三行代码打乱训练数据
  19. Windows7下安装redmine-3.4.6
  20. 学生&amp;部门智能匹配程序

热门文章

  1. Android 基础知识 -- Linux环境搭建
  2. C#常规TcpListener
  3. cssflex兼容(横向顶边对齐)
  4. windows上快速新建1g的文件
  5. 解决pjax重复绑定
  6. 01:认识QT
  7. Python3标准库:difflib差异计算工具
  8. Mysql部分存储引擎介绍
  9. 使用 C++11 编写可复用多线程任务池
  10. 清理 /dev/vda1 系统磁盘