BZOJ 1150 数据备份

1150: [CTSC2007]数据备份Backup
Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 1480 Solved: 592
[Submit][Status][Discuss]
Description

　　你在一家 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 的网络电缆，满足距离之和最小的要求。
Input

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

　　输出应由一个正整数组成，给出将2K个相异的办公楼连成k对所需的网络电缆的最小总长度。
Sample Input
5 2
1
3
4
6
12
Sample Output
4

---

【分析】
这道题比较6
翻译一下题目：数轴上有n个点，选出其中的k对点，每个点只能选一次，使得每对点的距离总和最小。
首先对数据差分。求出线段的属性:左边线段下标，右边线段下标，和该线段长度。并将它们push进一个小根堆！左线段与右线段下标用两个数组pre和next来储存，线段原下标和长度一起push进堆里。

下面就非常精妙了
加入一个线段后，它两边的线段就不能再加进去了。但是这样裸着贪心答案是错的。举个反例： 100 2 1 2 100 选出两对点，如果直接贪心选的是1和100，但显然2和2更优。那么我们将1push进去之后，将2 1 2合并（改动该线段的左右线段下标等等，具体见代码）,ans+=1。然后再push一个长度：2+2-1。
如果在之后的贪心中，2+2-1更优，那么我们就选出2+2-1，并将ans+=2+2-1
这样的话就相当于选了左右两个线段，而不选两个线段中间的线段。
操作：合并之后把线段x的左线段l,右线段r的长度修改为inf。在大循环中如果发现该线段当前长度与初长度不符，则直接弹出。
具体会用到make_pair等乱七八糟的东西，用于把两个属性同步放入堆中。

（感谢学姐的6到爆炸的博客）

---

【代码】

//BZOJ 1150 数据备份 
#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<map>
#define inf 1000000001
#define mp make_pair
#define pa pair<int,int>
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
priority_queue <pa,vector<pa>,greater<pa> > q;
const int T=100002;
int len[T],pre[T],next[T];
int b[T],a[T],n,m,ans;
int main()
{
    int i,j,k;
    scanf("%d%d",&n,&m);
    scanf("%d",&b[1]);
    fo(i,2,n)
    {
        scanf("%d",&b[i]);
        len[i]=b[i]-b[i-1];
        pre[i]=i-1;
        next[i]=i+1;
        q.push(mp(len[i],i));
    }
    pre[2]=0,next[n]=0;   //处理边界
    fo(i,1,m)
    {
        while(q.top().first!=len[q.top().second])
          q.pop();  //如果该点已经被删除，弹出 
        int x=q.top().second,l=pre[x],r=next[x];
        ans+=len[x];
        q.pop();
        next[x]=next[r];
        pre[next[r]]=x;
        pre[x]=pre[l];
        next[pre[l]]=x;   //合并三条线段
        if(l && r) len[x]=min(inf,len[l]+len[r]-len[x]);
        else len[x]=inf;
        len[l]=len[r]=inf;  //删除
        q.push(mp(len[x],x));
    }
    printf("%d\n",ans);
    return 0;
}