BZOJ1835: [ZJOI2010]base 基站选址【线段树优化DP】
Description
有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。
Input
输出文件中仅包含一个整数,表示最小的总费用。
Output
3 2 1 2 2 3 2 1 1 0 10 20 30
Sample Input
4
Sample Output
40%的数据中,N<=500;
100%的数据中,K<=N,K<=100,N<=20,000,Di<=1000000000,Ci<=10000,Si<=1000000000,Wi<=10000
思路
首先一个思路是前i个节点放j个基站
为了方便转移,强制第i个节点上一定会放上一个基站
然后怎么办?
有一个转移是\(dp_{i,k}=dp_{j,k}+calc(j+1,i-1)+c_i\)
然后考虑怎么快速计算对于每个i的\(dp_{j,k}+calc(j+1,i-1)\)
我们令每一个节点可以覆盖到它的基站位置是\([pre_i,nxt_i]\)
那么\(calc(j+1,i-1)=\sum_{p=j+1}^{i-1}[pre_i>j\&\&nxt_i<i]w_i\)
如果把i向右平移,发现受影响的新增节点只有\(nxt_t=i\)的所有节点
那么就按照\(nxt_t\)排序,每次因为我们维护所有的\(dp_{j,k}+calc(j+1,i-1)\)
所以对于所有\([1,pre_t-1]\)的dp值全部要加上\(w_t\)
然后就用线段树维护一下就可以了
注意当k是0的时候初值的处理
然后应该没什么大问题了
是我傻逼写错了线段树。。
//Author: dream_maker
#include<bits/stdc++.h>
using namespace std;
//----------------------------------------------
typedef pair<int, int> pi;
typedef long long ll;
typedef double db;
#define fi first
#define se second
#define fu(a, b, c) for (int a = b; a <= c; ++a)
#define fd(a, b, c) for (int a = b; a >= c; --a)
#define fv(a, b) for (int a = 0; a < (signed)b.size(); ++a)
const int INF_of_int = 1e9;
const ll INF_of_ll = 1e18;
template <typename T>
void Read(T &x) {
bool w = 1;x = 0;
char c = getchar();
while (!isdigit(c) && c != '-') c = getchar();
if (c == '-') w = 0, c = getchar();
while (isdigit(c)) {
x = (x<<1) + (x<<3) + c -'0';
c = getchar();
}
if (!w) x = -x;
}
template <typename T>
void Write(T x) {
if (x < 0) {
putchar('-');
x = -x;
}
if (x > 9) Write(x / 10);
putchar(x % 10 + '0');
}
//----------------------------------------------
const int N = 1e5 + 10;
int n, k, c[N], d[N], s[N], w[N];
int dp[2][N];
struct Node {
int pre, nxt, id;
} p[N];
bool operator < (const Node a, const Node b) {
return a.nxt < b.nxt;
}
#define LD (t << 1)
#define RD (t << 1 | 1)
int val[N << 2], tag[N << 2];
void pushnow(int t, int vl) {
val[t] += vl;
tag[t] += vl;
}
void pushup(int t) {
val[t] = min(val[LD], val[RD]);
}
void pushdown(int t) {
if (tag[t]) {
pushnow(LD, tag[t]);
pushnow(RD, tag[t]);
tag[t] = 0;
}
}
void build(int t, int l, int r) {
tag[t] = 0;
val[t] = INF_of_int;
if (l == r) return;
int mid = (l + r) >> 1;
build(LD, l, mid);
build(RD, mid + 1, r);
}
void modify(int t, int l, int r, int pos, int vl) {
if (l == r) {
val[t] = vl;
return;
}
pushdown(t);
int mid = (l + r) >> 1;
if (pos <= mid) modify(LD, l, mid, pos, vl);
else modify(RD, mid + 1, r, pos, vl);
pushup(t);
}
void modify(int t, int l, int r, int ql, int qr, int vl) {
if (ql > qr) return;
if (ql <= l && r <= qr) {
pushnow(t, vl);
return;
}
int mid = (l + r) >> 1;
pushdown(t);
if (qr <= mid) modify(LD, l, mid, ql, qr, vl);
else if (ql > mid) modify(RD, mid + 1, r, ql, qr, vl);
else {
modify(LD, l, mid, ql, mid, vl);
modify(RD, mid + 1, r, mid + 1, qr, vl);
}
pushup(t);
}
int find_pre(int x) {
int l = 1, r = x, res = 0;
while (l <= r) {
int mid = (l + r) >> 1;
if (s[x] >= d[x] - d[mid]) res = mid, r = mid - 1;
else l = mid + 1;
}
return res;
}
int find_nxt(int x) {
int l = x, r = n, res = 0;
while (l <= r) {
int mid = (l + r) >> 1;
if (s[x] >= d[mid] - d[x]) res = mid, l = mid + 1;
else r = mid - 1;
}
return res;
}
int main() {
#ifdef dream_maker
freopen("input.txt", "r", stdin);
freopen("output.txt", "w", stdout);
#endif
memset(dp, 0, sizeof(dp));
Read(n), Read(k);
fu(i, 2, n) Read(d[i]);
fu(i, 1, n) Read(c[i]);
fu(i, 1, n) Read(s[i]);
fu(i, 1, n) Read(w[i]);
fu(i, 1, n) {
p[i].pre = find_pre(i);
p[i].nxt = find_nxt(i);
p[i].id = i;
}
sort(p + 1, p + n + 1);
int ans = INF_of_int, now = 1, sumw = 0;
int ind = 0;
fu(i, 1, n) {
dp[ind][i] = sumw + c[i];
while (now <= n && p[now].nxt <= i)
sumw += w[p[now].id], ++now;
}
fu(j, 2, k + 1) {
build(1, 1, n);
ind ^= 1, now = 1;
fu(i, 1, n) {
dp[ind][i] = val[1] + c[i];
while (now <= n && p[now].nxt <= i)
modify(1, 1, n, 1, p[now].pre - 1, w[p[now].id]), ++now;
modify(1, 1, n, i, dp[ind ^ 1][i]);
}
ans = min(ans, val[1]);
}
Write(ans);
return 0;
}
BZOJ1835: [ZJOI2010]base 基站选址【线段树优化DP】的更多相关文章
- BZOJ1835: [ZJOI2010]base 基站选址(线段树优化Dp)
Description 有N个村庄坐落在一条直线上,第i(i>1)个村庄距离第1个村庄的距离为Di.需要在这些村庄中建立不超过K个通讯基站,在第i个村庄建立基站的费用为Ci.如果在距离第i个村庄 ...
- 【BZOJ1835】[ZJOI2010]base 基站选址 线段树+DP
[BZOJ1835][ZJOI2010]base 基站选址 Description 有N个村庄坐落在一条直线上,第i(i>1)个村庄距离第1个村庄的距离为Di.需要在这些村庄中建立不超过K个通讯 ...
- 洛谷$P2605\ [ZJOI2010]$基站选址 线段树优化$dp$
正解:线段树优化$dp$ 解题报告: 传送门$QwQ$ 难受阿,,,本来想做考试题的,我还造了个精妙无比的题面,然后今天讲$dp$的时候被讲到了$kk$ 先考虑暴力$dp$?就设$f_{i,j}$表示 ...
- [ZJOI2010]基站选址,线段树优化DP
G. base 基站选址 内存限制:128 MiB 时间限制:2000 ms 标准输入输出 题目类型:传统 评测方式:文本比较 题目描述 有N个村庄坐落在一条直线上,第i(i>1)个村庄距离 ...
- luogu P2605 [ZJOI2010]基站选址 线段树优化dp
LINK:基站选址 md气死我了l达成1结果一直调 显然一个点只建立一个基站 然后可以从左到右进行dp. \(f_{i,j}\)表示强制在i处建立第j个基站的最小值. 暴力枚举转移 复杂度\(n\cd ...
- BZOJ 1835 [ZJOI2010]基站选址 (线段树优化DP)
题目大意:略 洛谷题面传送门 BZOJ题面传送门 注意题目的描述,是村庄在一个范围内去覆盖基站,而不是基站覆盖村庄,别理解错了 定义$f[i][k]$表示只考虑前i个村庄,一共建了$k$个基站,最后一 ...
- luogu2605 基站选址 (线段树优化dp)
设f[i][j]表示在第i个村庄建第j个基站的花费 那么有$f[i][j]=min\{f[k][j-1]+w[k,i]\}$,其中w[k,i]表示在k,i建基站,k,i中间的不能被满足的村庄的赔偿金之 ...
- bzoj1835[ZJOI2010]base基站选址
据说正解是什么线段树优化DP,但是作为脑子有坑选手,我们需要5k的做法: 主席树+决策单调性..... F[m][i]表示已经放置了m个基站,第m个基站放置在第i个村庄,第i个村庄及之前的村庄的总最少 ...
- 2018.11.06 bzoj1835: [ZJOI2010]base 基站选址(线段树优化dp)
传送门 二分出每个点不需要付www贡献的范围,然后可以推出转移式子: f[i][j]=f[i−1][k]+value(k+1,j)+c[i]f[i][j]=f[i-1][k]+value(k+1,j) ...
- bzoj1835: [ZJOI2010]base 基站选址
新的一年新的开始.结果第一题就用了几乎一周.而且感觉很不好. 先检讨自己.最近写的各种数据结构模板基本没打过出来,各种细节崩盘,这题线段树都居然被lazy标记没清零卡挂. DP还是博大精深,这东西感觉 ...
随机推荐
- C++中的常量定义
本篇笔记总结自一次代码检视. 一般来说,使用C语言编程时我们都习惯在代码当中使用C当中的宏定义来定义一个数值常量: #define MY_CONST 7 在C++开发项目时,也会经常存在沿袭C当中常量 ...
- js判断浏览器的类型和获得浏览器的版本
<!DOCTYPE html><html> <head> <meta charset="UTF-8"> ...
- VS2010/MFC编程入门之十七(对话框:文件对话框)
上一讲鸡啄米介绍的是消息对话框,本节讲解文件对话框.文件对话框也是很常用的一类对话框. 文件对话框的分类 文件对话框分为打开文件对话框和保存文件对话框,相信大家在Windows系统中经常见 ...
- 2018 Multi-University Training Contest 9 Solution
A - Rikka with Nash Equilibrium 题意:构造一个$n * m$的矩阵,使得$[1, n * m]$ 中每个数只出现一次,并且纳什均衡只出现一次. 思路:从大到小的放置,每 ...
- ZW云推客即将登场
ZW云推客即将登场 即"4K云字库"框架文件发布后,ZW团队即将发布一全功能的:ZW云推客系统. ZW云推客系统,是原zBrow"百万社区营销系统".&qu ...
- 20145318《网络对抗》逆向及Bof基础
实践目标 本次实践的对象是一个名为pwn1的linux可执行文件. 该程序正常执行流程是:main调用foo函数,foo函数会简单回显任何用户输入的字符串. 该程序同时包含另一个代码片段,getShe ...
- 【异常记录(六)】vs文件乱码:文件加载,使用Unicode(UTF-8)编码加载文件xxx时,有些字节已用Unicode替换字符替换。保存该文件将不会保留原始文件内容。
VS2013偶遇这种情况,页面汉字编码出现乱码. .... 按照网上查到的: 工具>选项>文本编辑器> 勾选了 然并卵,还是乱码... 其实炒鸡简单 用记事本打开另存为,选择 ...
- python 去除不可见的控制字符
尤其是在json load的时候,字符串中的不可见控制字符可能会导致错误,应该先对字符串进行控制字符过滤. 对网页文本同样适用,最好在处理网页文本时先进性控制字符清洗. Replace null by ...
- shell 变量赋值运算
.变量赋值: name=lbg 等号前后不能有空格 name="Lebron James" 变量值中有空格要用双引号 echo ${name} 用${}更保险 shopt -s - ...
- python 解压zip压缩包
在当前路径解压zip压缩包,生成同名文件夹,内部目录结构与压缩包一致 import zipfile import os def un_zip(file_name): """ ...