博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ1150[CTSC2007]数据备份Backup——模拟费用流+堆+链表
阅读量:6976 次
发布时间:2019-06-27

本文共 2540 字,大约阅读时间需要 8 分钟。

题目描述

 

  你在一家 IT 公司为大型写字楼或办公楼(offices)的计算机数据做备份。然而数据备份的工作是枯燥乏味
的,因此你想设计一个系统让不同的办公楼彼此之间互相备份,而你则坐在家中尽享计算机游戏的乐趣。已知办公
楼都位于同一条街上。你决定给这些办公楼配对(两个一组)。每一对办公楼可以通过在这两个建筑物之间铺设网
络电缆使得它们可以互相备份。然而,网络电缆的费用很高。当地电信公司仅能为你提供 K 条网络电缆,这意味
着你仅能为 K 对办公楼(或总计2K个办公楼)安排备份。任一个办公楼都属于唯一的配对组(换句话说,这 2K 
个办公楼一定是相异的)。此外,电信公司需按网络电缆的长度(公里数)收费。因而,你需要选择这 K 对办公
楼使得电缆的总长度尽可能短。换句话说,你需要选择这 K 对办公楼,使得每一对办公楼之间的距离之和(总距
离)尽可能小。下面给出一个示例,假定你有 5 个客户,其办公楼都在一条街上,如下图所示。这 5 个办公楼分
别位于距离大街起点 1km, 3km, 4km, 6km 和 12km 处。电信公司仅为你提供 K=2 条电缆。
  上例中最好的配对方案是将第 1 个和第 2 个办公楼相连,第 3 个和第 4 个办公楼相连。这样可按要求使用
 K=2 条电缆。第 1 条电缆的长度是 3km-1km=2km ,第 2 条电缆的长度是 6km-4km=2km。这种配对方案需要总长
 4km 的网络电缆,满足距离之和最小的要求。

输入

 

第一行包含整数n和k
其中n(2≤n≤100000)表示办公楼的数目,k(1≤k≤n/2)表示可利用的网络电缆的数目。
接下来的n行每行仅包含一个整数(0≤s≤1000000000),表示每个办公楼到大街起点处的距离。
这些整数将按照从小到大的顺序依次出现。

输出

输出应由一个正整数组成,给出将2K个相异的办公楼连成k对所需的网络电缆的最小总长度。

样例输入

5 2
1
3
4
6
12

样例输出

4
 
首先选择的每对点一定是相邻的点,将相邻坐标相减即可得到$n-1$个可选方案,我们需要在这些方案中选择$k$个。因为每个点只能在最多一个方案中,所以相邻的方案不能同时选择。这个问题显然可以用费用流来解决:将$n$个点中每个点分成两个点$a_{i},b_{i}$,源点连向$a_{i}$,$b_{i}$连向汇点,$a_{i}$连向$b_{i-1}$和$b_{i+1}$。这样跑一遍费用流既可以得到答案二倍大小的最小费用(因为在建模中$i$匹配$i+1$和$i+1$匹配$i$被算了两次),我们考虑这个图的特殊性质:每个点只向它两边的两个点连边(即不考虑源点与汇点连的边,剩余边是两条交错的链)。那么当选了一种方案$i$进行增广并建反向边之后,如果在之后增广到$i$相邻的一个方案(假设是$i-1$),那么根据费用流每次选择最短路的原则一定会走$i$方案的反向边然后再走$i+1$方案。这样我们就可以模拟这个过程,用堆来维护当前代价最小的方案,每次取出最优方案后(假设最优方案为$b$,相邻方案为$a,c$),将最优方案与相邻方案取出,并在原位置放入代价为$a+c-b$的新方案,这样如果选了新方案就说明实际选择了$a$和$c$而放弃了$b$,也就实现了费用流的反悔过程。而这样每次取出堆顶时也保证了多选了一个方案(不过这里注意两端的处理)。而相邻关系用双向链表维护即可。
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define ll long long#define pr pair
using namespace std;int n,k;priority_queue< pr,vector
,greater
>q;int cnt;int pre[600010];int suf[600010];int a[600010];int vis[300010];ll s[300010];int t[300010];int v[600010];ll ans;int main(){ scanf("%d%d",&n,&k); for(int i=1;i<=n;i++) { scanf("%d",&a[i]); if(i!=1) { v[i]=a[i]-a[i-1]; } } for(int i=2;i<=n;i++) { cnt++; s[cnt]=1ll*v[i]; t[cnt]=cnt; pre[cnt]=cnt-1; suf[cnt]=cnt+1; q.push(make_pair(s[cnt],t[cnt])); } pre[1]=0; suf[n-1]=0; s[0]=1<<30; while(k) { int now=q.top().second; ll res=q.top().first; q.pop(); if(!vis[now]) { ans+=res; vis[now]=1; vis[pre[now]]=1; vis[suf[now]]=1; cnt++; t[cnt]=cnt; s[cnt]=s[pre[now]]+s[suf[now]]-res; pre[cnt]=pre[pre[now]]; suf[cnt]=suf[suf[now]]; suf[pre[pre[now]]]=cnt; pre[suf[suf[now]]]=cnt; q.push(make_pair(s[cnt],t[cnt])); k--; } } printf("%lld",ans);}

转载于:https://www.cnblogs.com/Khada-Jhin/p/10421588.html

你可能感兴趣的文章
Wall Treatment
查看>>
iOS开发的一些基础知识
查看>>
(转)Maven常用命令
查看>>
iOS中使用UIWebView与JS进行交互
查看>>
ios之sqllite3简单使用
查看>>
IOS实现自动循环滚动广告--ScrollView的优化和封装
查看>>
IOS 数据存储之 Core Data详解
查看>>
阿里云RDS读写分离数据查询延迟解决
查看>>
微信扫码支付asp.net(C#)实现步骤
查看>>
How to create ISO on macOS
查看>>
python中动态导入模块
查看>>
Unity UGUI Button 无法点击问题一例
查看>>
2019.01.19-2018年6月NEYC集训counting
查看>>
codeforces 801B Valued Keys
查看>>
C语言不容易识别的坑
查看>>
Redis集群
查看>>
2016-8-15
查看>>
【网站分享】文档转换网站分享
查看>>
linux配置本地源
查看>>
爬虫_电影天堂 热映电影(xpath)
查看>>