码迷,mamicode.com
首页 > Windows程序 > 详细

[BZOJ1911][Apio2010]特别行动队

时间:2016-04-20 21:43:31      阅读:276      评论:0      收藏:0      [点我收藏+]

标签:

[BZOJ1911][Apio2010]特别行动队

试题描述

技术分享

输入

技术分享

输出

技术分享

输入示例

4
-1 10 -20
2 2 3 4

输出示例

9

数据规模及约定

技术分享

题解

DP,令 f(i) 表示前 i 个士兵可以达到的最大修正后的战斗力,则有 f(i) = max{ f(j) + d(j+1 ~ i) | 0 < j < i },其中d(i ~ j)表示第 i 个士兵到第 j 个士兵组成一队修正之后的战斗力。展开之后得到

f(i) - a·Si2 - b·Si - c = f(j) + a·Sj2 - b·Sj - 2·a·Si · Sj

发现蓝色部分只与i相关,黄色部分只与j相关,绿色部分随 i 的增加是单调递减的,橙色部分只与j相关,把它们依次换元成b, y, k, x,得到

b = y - k · x,变形得 y = kx + b,斜率k单调递增,我们要求的是b的最大值,所以维护所有(x, y)点构成的上凸壳即可。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <stack>
#include <vector>
#include <queue>
#include <cstring>
#include <string>
#include <map>
#include <set>
using namespace std;

const int BufferSize = 1 << 16;
char buffer[BufferSize], *Head, *tail;
inline char Getchar() {
    if(Head == tail) {
        int l = fread(buffer, 1, BufferSize, stdin);
        tail = (Head = buffer) + l;
    }
    return *Head++;
}
int read() {
    int x = 0, f = 1; char c = Getchar();
    while(!isdigit(c)){ if(c == ‘-‘) f = -1; c = Getchar(); }
    while(isdigit(c)){ x = x * 10 + c - ‘0‘; c = Getchar(); }
    return x * f;
}

#define maxn 1000010
#define LL long long
int n, a, b, c, S[maxn], q[maxn], l, r;
LL f[maxn];

double slop(int i, int j) {
	double yi = i ? f[i] + (double)a * S[i] * S[i] - (double)b * S[i] : 0.0, yj = f[j] + (double)a * S[j] * S[j] - (double)b * S[j];
	double xi = i ? S[i] : 0, xj = S[j];
	return (yi - yj) / (xi - xj);
}

int main() {
	n = read(); a = read(); b = read(); c = read();
	for(int i = 1; i <= n; i++) S[i] = S[i-1] + read();
	
	f[0] = 0;
	l = r = 1; q[1] = 0;
	for(int i = 1; i <= n; i++) {
		while(l < r && slop(q[l], q[l+1]) > 2.0 * a * S[i]) l++;
		int j = q[l];
		f[i] = f[j] + (LL)a * S[j] * S[j] - (LL)b * S[j] - 2ll * a * S[i] * S[j] + (LL)a * S[i] * S[i] + (LL)b * S[i] + c;
		while(l < r && slop(q[r], i) >= slop(q[r-1], q[r])) r--;
		q[++r] = i;
	}
	
	printf("%lld\n", f[n]);
	
	return 0;
}

 

[BZOJ1911][Apio2010]特别行动队

标签:

原文地址:http://www.cnblogs.com/xiao-ju-ruo-xjr/p/5414479.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!