1835: [ZJOI2010]base 基站选址

题目描述

有N个村庄坐落在一条直线上,第i(i>1)个村庄距离第1个村庄的距离为Di。需要在这些村庄中建立不超过K个通讯基站,在第i个村庄建立基站的费用为Ci。如果在距离第i个村庄不超过Si的范围内建立了一个通讯基站,那么就成它被覆盖了。如果第i个村庄没有被覆盖,则需要向他们补偿,费用为Wi。现在的问题是,选择基站的位置,使得总费用最小。 输入数据 (base.in) 输入文件的第一行包含两个整数N,K,含义如上所述。 第二行包含N-1个整数,分别表示D2,D3,…,DN ,这N-1个数是递增的。 第三行包含N个整数,表示C1,C2,…CN。 第四行包含N个整数,表示S1,S2,…,SN。 第五行包含N个整数,表示W1,W2,…,WN。

输入输出格式

输入格式:

输入文件的第一行包含两个整数N,K,含义如上所述。

第二行包含N-1个整数,分别表示D2,D3,…,DN ,这N-1个数是递增的。

第三行包含N个整数,表示C1,C2,…CN。

第四行包含N个整数,表示S1,S2,…,SN。

第五行包含N个整数,表示W1,W2,…,WN。

输出格式:

输出文件中仅包含一个整数,表示最小的总费用。

输入输出样例

输入样例#1:

3 2
1 2
2 3 2
1 1 0
10 20 30
输出样例#1:

4

说明

40%的数据中,N<=500;

100%的数据中,K<=N,K<=100,N<=20,000,Di<=1000000000,Ci<=10000,Si<=1000000000,Wi<=10000。


到洛谷偷了题面

理解了题解之后还是比较好写的

显然f[i][j]表示前i个村子建了j个基站且第i个村子建了基站

f[i][j]=c[i]+min{f[k][j-1]+cost(k,j)}

cost(k,j)表示k有一个基站,j有一个基站,k..j的补偿代价

关键就是快速计算这个东西了

线段树优化,就是用线段树区间min来logn获得转移来的状态中最小值吧

j这一维显然可以滚掉

想办法让线段树每个点表示了选这个点作为转移点时的代价

先把f[][j-1]建树,然后处理cost的问题

当i-->i+1时,发现左端点不变,右段点右移了,那么哪些刚好最远i位置可以覆盖到的点就可能要补偿了

所以对于点x,通过二分计算st[x]和ed[x]为x最左和最右到哪,然后用链表记录ed[x]为某个值的点有哪些,

对于ed[x]=i的点线段树[1,st[x]-1]区间加w[x],因为这些点右面不能被覆盖,左面再不能的话就要补偿了

复杂度 k*n*logn,区间加n次,区间min也有n次

注意:

1.j==1的时候O(n)特判就行了

2.n++ k++后 d[n]=w[n]=INF c[n]=0,f[n]就是最优解了

//
//  main.cpp
//  bzoj1835
//
//  Created by Candy on 2017/1/8.
//  Copyright © 2017年 Candy. All rights reserved.
//

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define lc o<<1
#define rc o<<1|1
#define m ((l+r)>>1)
#define lson o<<1,l,m
#define rson o<<1|1,m+1,r
,INF=1e9+;
inline int read(){
    ,f=;
    ;c=getchar();}
    +c-';c=getchar();}
    return x*f;
}
int n,k,d[N],c[N],s[N],w[N];
int st[N],ed[N],f[N];

struct edge{
    int v,ne;
}e[N];
int h[N],cnt;
inline void ins(int u,int v){
    cnt++;
    e[cnt].v=v;e[cnt].ne=h[u];h[u]=cnt;
}

struct node{
    int mn,tag;
}t[N<<];
inline void merge(int o){t[o].mn=min(t[lc].mn,t[rc].mn);}
inline void paint(int o,int d){
    t[o].tag+=d;
    t[o].mn+=d;
}
inline void pushDown(int o){
    if(t[o].tag){
        paint(lc,t[o].tag);
        paint(rc,t[o].tag);
        t[o].tag=;
    }
}
void build(int o,int l,int r){
    t[o].tag=;
    if(l==r) t[o].mn=f[l];
    else{
        build(lson);
        build(rson);
        merge(o);
    }
}
void segAdd(int o,int l,int r,int ql,int qr,int d){
    if(ql>qr) return;
    if(ql<=l&&r<=qr) paint(o,d);
    else{
        pushDown(o);
        if(ql<=m) segAdd(lson,ql,qr,d);
        if(m<qr) segAdd(rson,ql,qr,d);
        merge(o);
    }
}
int segQue(int o,int l,int r,int ql,int qr){
    ;
    if(ql<=l&&r<=qr) return t[o].mn;
    else{
        pushDown(o);
        int mn=INF;
        if(ql<=m) mn=min(mn,segQue(lson,ql,qr));
        if(m<qr) mn=min(mn,segQue(rson,ql,qr));
        return mn;
    }
}

void dp(){
    ;
    ;i<=n;i++){
        f[i]=_+c[i];
        for(int k=h[i];k;k=e[k].ne)
            _+=w[e[k].v];
        //printf("f j1 %d\n",f[i]);
    }

    ;j<=k;j++){
        build(,,n);
        ;i<=n;i++){
            f[i]=segQue(,,n,,i-)+c[i];
            for(int k=h[i];k;k=e[k].ne){
                int v=e[k].v;
                segAdd(,,n,,st[v]-,w[v]);
            }
        }
        ans=min(ans,f[n]);
    }
    printf("%d",ans);
}
int main(int argc, const char * argv[]) {
    n=read();k=read();
    ;i<=n;i++) d[i]=read();
    ;i<=n;i++) c[i]=read();
    ;i<=n;i++) s[i]=read();
    ;i<=n;i++) w[i]=read();
    n++;k++;
    d[n]=INF;w[n]=INF;
    ;i<=n;i++){
        st[i]=lower_bound(d+,d++n,d[i]-s[i])-d;
        ed[i]=lower_bound(d+,d++n,d[i]+s[i])-d;
        if(d[ed[i]]-d[i]>s[i]) ed[i]--;
        ins(ed[i],i);
        //printf("sted %d %d %d\n",i,st[i],ed[i]);
    }
    dp();
    ;
}

随机推荐

  1. Intellij IDEA 快捷键整理

    CSDN 2016博客之星评选结果公布      [系列直播]算法与游戏实战技术      "我的2016"主题征文活动 Intellij IDEA 快捷键整理(TonyCody) ...

  2. Mac平台与Windows平台下AndroidStudio增量升级

    Android Studio增量升级什么情况下使用最合适呢? 比如现在的as版本是2.2版本,而你的as版本2.0版本,这个时候点Check For Updates就没有反应了,因为你已经2个有版本没 ...

  3. R----dplyr包介绍学习

    dplyr包:plyr包的替代者,专门面对数据框,将ddplyr转变为更易用的接口 %>%来自dplyr包的管道函数,其作用是将前一步的结果直接传参给下一步的函数,从而省略了中间的赋值步骤,可以 ...

  4. Storm Bolt接口

          Bolt是Topology中数据处理的基本单元,也是Storm针对处理过程的编程单元.Topology中所有的处理都是在这些bolt中完成的. Bolt可以将数据项发送至多个数据流(Str ...

  5. 【poj1740】 A New Stone Game

    http://poj.org/problem?id=1740 (题目链接) 男人八题之一 题意 对于n堆石子,每堆若干个,两人轮流操作,每次操作分两步,第一步从某堆中去掉至少一个,第二步(可省略)把该 ...

  6. JS中setAttribute的兼容性问题(摘自leejersey)

    class和className兼容方法: object.setAttribute("class","content") 在IE8.Chrome.火狐.Opera ...

  7. HBase Coprocessor 剖析与编程实践(转载http://www.cnblogs.com/ventlam/archive/2012/10/30/2747024.html)

    HBase Coprocessor 剖析与编程实践 1.起因(Why HBase  Coprocessor) HBase作为列族数据库最经常被人诟病的特性包括:无法轻易建立“二级索引”,难以执行求和. ...

  8. jst通用删除数组中重复的值和删除字符串中重复的字符

    以下内容属于个人原创,转载请注明出处,非常感谢! 删除数组中重复的值或者删除字符串重复的字符,是我们前端开发人员碰到很多这样的场景.还有求职者在被面试时也会碰到这样的问题!比如:问删除字符串重复的字符 ...

  9. 屏幕旋转时调用PopupWindow update方法更新位置失效的问题及解决方案

       接到一个博友的反馈,在屏幕旋转时调用PopupWindow的update方法失效.使用场景如下:在一个Activity中监听屏幕旋转事件,在Activity主布局文件中有个按钮点击弹出一个Pop ...

  10. Java多线程volatile和synchronized总结

    volatile是轻量级的synchronized,在多处理器(多线程)开发中保证了共享变量的"可见性".可见性表示当一个线程修改了一个共享变量时,另外一个线程能读到这个修改的值. ...