http://uoj.ac/problem/29

cdq四次处理出一直向左, 一直向右, 向左后回到起点, 向右后回到起点的dp数组,最后统计答案。

举例:\(fi\)表示一直向右走i天能参观的最多景点数。

其中有一个很重要的条件\(fi≤fi+1fi≤fi+1\),这个条件是分治的前提。

关于这个条件的证明,我想了好久才想出来,用反证法证明一下就行。

分治时需要用主席树维护路径上的前k大和。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include"holiday.h"
using namespace std;
typedef long long ll;
const int N = 100003;
const int M = 250003; struct node {
int l, r, s;
ll sum;
node() {l = r = s = sum = 0;}
} T[N * 20];
int a[N], H[N], cnt = 0, root[N], top, st;
ll f[M], g[M], f1[M], g1[M]; void update(int &pos, int l, int r, int key) {
T[++cnt] = T[pos]; pos = cnt;
++T[pos].s; T[pos].sum += H[key];
if (l == r) return;
int mid = (l + r) >> 1;
if (key <= mid) update(T[pos].l, l, mid, key);
else update(T[pos].r, mid + 1, r, key);
} ll query(int tl, int tr, int l, int r, int num) {
if (l == r) {return min(T[tr].sum - T[tl].sum, 1ll * H[l] * num);}
int mid = (l + r) >> 1, s = T[T[tr].r].s - T[T[tl].r].s;
if (s >= num)
return query(T[tl].r, T[tr].r, mid + 1, r, num);
else
return T[T[tr].r].sum - T[T[tl].r].sum + query(T[tl].l, T[tr].l, l, mid, num - s);
} void cdq_f(int l, int r, int tmp_l, int tmp_r) {
if (l > r) return;
int mid = (l + r) >> 1, pos = tmp_l;
ll t;
for(int i = tmp_l; i - st <= mid && i <= tmp_r; ++i)
if ((t = query(root[st - 1], root[i], 1, top, mid - i + st)) > f[mid])
f[mid] = t, pos = i;
cdq_f(l, mid - 1, tmp_l, pos);
cdq_f(mid + 1, r, pos, tmp_r);
} void cdq_f1(int l, int r, int tmp_l, int tmp_r) {
if (l > r || tmp_l > tmp_r) return;
int mid = (l + r) >> 1, pos = tmp_l;
ll t;
for(int i = tmp_l; ((i - st) << 1) <= mid && i <= tmp_r; ++i)
if ((t = query(root[st], root[i], 1, top, mid - ((i - st) << 1))) > f1[mid])
f1[mid] = t, pos = i;
cdq_f1(l, mid - 1, tmp_l, pos);
cdq_f1(mid + 1, r, pos, tmp_r);
} void cdq_g(int l, int r, int tmp_l, int tmp_r) {
if (l > r) return;
int mid = (l + r) >> 1, pos = tmp_r;
ll t;
for(int i = tmp_r; st - i <= mid && i >= tmp_l; --i) {
if ((t = query(root[i - 1], root[st], 1, top, mid - st + i)) > g[mid])
g[mid] = t, pos = i;
}
cdq_g(l, mid - 1, pos, tmp_r);
cdq_g(mid + 1, r, tmp_l, pos);
} void cdq_g1(int l, int r, int tmp_l, int tmp_r) {
if (l > r || tmp_l > tmp_r) return;
int mid = (l + r) >> 1, pos = tmp_r;
ll t;
for(int i = tmp_r; ((st - i) << 1) <= mid && i >= tmp_l; --i) {
if ((t = query(root[i - 1], root[st - 1], 1, top, mid - ((st - i) << 1))) > g1[mid])
g1[mid] = t, pos = i;
}
cdq_g1(l, mid - 1, pos, tmp_r);
cdq_g1(mid + 1, r, tmp_l, pos);
} ll findMaxAttraction(int n, int start, int d, int attraction[]) {
for(int i = 0; i < n; ++i) H[++cnt] = attraction[i];
sort(H + 1, H + cnt + 1);
cnt = unique(H + 1, H + cnt + 1) - H;
for(int i = 1; i <= n; ++i) a[i] = lower_bound(H + 1, H + cnt, attraction[i - 1]) - H;
top = cnt - 1; cnt = 0; for(int i = 1; i <= n; ++i) {
root[i] = root[i - 1];
update(root[i], 1, top, a[i]);
} st = start + 1;
cdq_f(0, d, st, n);
cdq_f1(0, d, st + 1, n);
cdq_g(0, d, 1, st);
cdq_g1(0, d, 1, st - 1); ll ans = max(f[d], g[d]);
for(int i = 1; i <= d; ++i)
ans = max(ans, max(f[i] + g1[d - i], g[i] + f1[d - i]));
return ans;
}

UOJ上是交互题的形式

【UOJ #29】【IOI 2014】holiday的更多相关文章

  1. JAVA 基础编程练习题29 【程序 29 求矩阵对角线之和】

    29 [程序 29 求矩阵对角线之和] 题目:求一个 3*3 矩阵对角线元素之和 程序分析:利用双重 for 循环控制输入二维数组,再将 a[i][i]累加后输出. package cskaoyan; ...

  2. uoj 41 【清华集训2014】矩阵变换 婚姻稳定问题

    [清华集训2014]矩阵变换 Time Limit: 20 Sec  Memory Limit: 256 MB 题目连接 http://uoj.ac/problem/41 Description 给出 ...

  3. AC日记——【清华集训2014】奇数国 uoj 38

    #38. [清华集训2014]奇数国 思路: 题目中的number与product不想冲: 即为number与product互素: 所以,求phi(product)即可: 除一个数等同于在模的意义下乘 ...

  4. [官方软件] Easy Sysprep v4.3.29.602 【系统封装部署利器】(2016.01.22)--skyfree大神

    [官方软件] Easy Sysprep v4.3.29.602 [系统封装部署利器](2016.01.22) Skyfree 发表于 2016-1-22 13:55:55 https://www.it ...

  5. 【UOJ】67 新年的毒瘤 &【BZOJ】1123 BLO

    [UOJ 67] 题目链接: 传送门 题解: 第一眼很懵逼……这什么鬼. 思考什么点复合条件……(o(>﹏<)o 1.树,也就是说还剩n-2条边,等价于要删去一个度数为m-n+2的点. 2 ...

  6. 【UOJ#236】[IOI2016]railroad(欧拉回路,最小生成树)

    [UOJ#236][IOI2016]railroad(欧拉回路,最小生成树) 题面 UOJ 题解 把速度看成点,给定的路段看成边,那么现在就有了若干边,然后现在要补上若干边,以及一条\([inf,\) ...

  7. 【UOJ#177】欧拉回路

    [UOJ#177]欧拉回路 题面 UOJ 题解 首先图不连通就没啥好搞的了. 对于无向图而言,每个点度数为偶数. 对于有向图而言,每个点入度等于出度. 然后就是一本通上有的做法,直接\(dfs\)一遍 ...

  8. 【UOJ#311】【UNR #2】积劳成疾(动态规划)

    [UOJ#311][UNR #2]积劳成疾(动态规划) UOJ Solution 考虑最大值分治解决问题.每次枚举最大值所在的位置,强制不能跨过最大值,左右此时不会影响,可以分开考虑. 那么设\(f[ ...

  9. 【UOJ#450】【集训队作业2018】复读机(生成函数,单位根反演)

    [UOJ#450][集训队作业2018]复读机(生成函数,单位根反演) 题面 UOJ 题解 似乎是\(\mbox{Anson}\)爷的题. \(d=1\)的时候,随便怎么都行,答案就是\(k^n\). ...

随机推荐

  1. java 解决中文乱码

    //1.代码解决 public class LuanMa { public static String getNewString(String luanma){ String result = &qu ...

  2. 手机开启HDR后拍照有什么不同?

    转自http://www.leiphone.com/news/201406/hdr-pic.html HDR这词大家可能并不陌生,但大部分人只是听过而已,在脑子并没有太多的概念,可能只是简单的认为HD ...

  3. 为什么 Java 不提供无符号类型呢?

    网上查资料,无意中找到一个java写的开源论坛,用的人还挺多 http://jforum.net/ 查MD5,了解到 Java getBytes方法详解(字符集问题) http://liushilan ...

  4. Eclipse启动时布局不合理调整

    1. 关掉 启动页 2. 关掉InstSearch页 3.修正InSearch布局 3.1 默认InstSearch不合理,影响使用. 3.2 Inst 搜索一次,然后最小化InstSearch框,再 ...

  5. Android保存ArrayList至SharedPreferences

    保存ArrayList至SharedPreferences 其中ArrayList中每个元素为String List<String> environmentList = new Array ...

  6. webpack里的module选项配置

    webpack里的module 有loaders选项和noParse选项,noParse选项配置不需要解析的目录和文件 module:{ loaders:[], noParse: [ path.joi ...

  7. webapp开发需要注意的浏览器内核知识

    Web App:1.开发成本较低使用web开发技术就可以轻松的完成web app的开发2.升级较简单升级不需要通知用户,在服务端更新文件即可,用户完全没有感觉3.维护比较轻松和一般的web一样,维护比 ...

  8. jdbc连接数据库总结

    jdbc支持多种数据库,比如说oracle, mysql, mssql,现在总结一下连接各种数据库的相关知识 1,mysql连接,代码如下 Class.forName("com.mysql. ...

  9. 内网机(无网络安装 .NET Core win开发环境

    1.安装 vs2015 update3 2.按顺序安装以下包 DotNetCore.1.0.0-SDK.Preview2-x64.exe aspnetcoremodule_x64_en_rc2_14. ...

  10. Castle.ActiveRecord 多对多关系 引发的错误处理

    在Castle.ActiveRecord 实体类中,如果两个对象有 “多对多” 关系,一般的做法是将其分解为 两个“一对多”关系,但有时引发了 “您要删除 或 引用 的对象#2在数据库中不存在”的异常 ...