http://uoj.ac/problem/29

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

``````#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上是交互题的形式

## 随机推荐

1. Windows8笔记

一.WIN8 硬盘变为可弹出设备 如果电脑使用的是nforce芯片组的主板和串口硬盘.在系统中安装好所有的硬件驱动程序后,会出现这种情况.这是nforce芯片组在安装了IDE-SW主板驱动以后,系统会 ...

2. IT人学习方法论（三）：高效学习

一些有关“怎么学”的建议 首先需要明确方向,否则即使学习方法再高效,也不免南辕北辙,背离自己的目标.关于学习方向的讨论,请参见之前的一篇文章 .下面我来重点说一说有关“怎么学”的建议. IT技术,不是 ...

3. Android 图片轮播（最简单的）

布局文件 <RelativeLayout xmlns:android="http://schemas.android.com/apk/res/android" android ...

#import "RootViewController.h" @interface RootViewController () @end @implementation RootV ...

5. String类详解，StringBuffer

先说一下String类的equals()方法. 下面我们先看一段代码: 这段代码输出的结果为: ture true -------------- false 咋看之下貌似Object类比较特别,那么我 ...

6. wordpress常用插件汇总

WordPress之所以能成为目前最具人气的独立博客程序,除了无数爱好者为它开发的主题外,源源不断的插件支持也是重要的原因之一.wordpress的强大,也是在于无数爱好者源源不断的主题和插件. wo ...

7. JAVAEE——spring01：介绍、搭建、概念、配置详解、属性注入和应用到项目

一.spring介绍 1.三层架构中spring位置 2.spring一站式框架 正是因为spring框架性质是属于容器性质的. 容器中装什么对象就有什么功能.所以可以一站式. 不仅不排斥其他框架,还 ...

8. Springboot 系列（四）Spring Boot 日志框架

注意:本 Spring Boot 系列文章基于 Spring Boot 版本 v2.1.1.RELEASE 进行学习分析,版本不同可能会有细微差别. 前言 Spring 框架选择使用了 JCL 作为默 ...

9. usermod - linux修改用户帐户信息

usermod - 修改用户帐户信息 modify a user account usermod [options] user_name usermod 命令修改系统帐户文件来反映通过命令行指定的变化 ...

10. Asp.Net Core 发布异常 502.5 [The Application process failed to Start]

出现这个问题大部分时间都是因为发布时,少打包了一些文件.. 只打包了.Net Core的运行时库,没有打包Asp.Net Core 运行时.. 需要在打包指导文件中加入以下节点 <Propert ...