【斜率优化】特别行动队

原创 2012年03月23日 21:29:41

特别行动队
【问题描述】
你有一支由n名预备役士兵组成的部队,士兵从1n编号,要将他们拆分
成若干特别行动队调入战场。出于默契的考虑,同一支特别行动队中队员的编号
应该连续,即为形如(i,i + 1, …, i + k)的序列。
编号为i的士兵的初始战斗力为xi,一支特别行动队的初始战斗力x为队内
士兵初始战斗力之和,即x= xi + xi+1 + … + xi+k
通过长期的观察,你总结出一支特别行动队的初始战斗力x将按如下经验公
式修正为x'x'= ax
2
+bx + c
,其中a,b, c是已知的系数(a< 0)。
作为部队统帅,现在你要为这支部队进行编队,使得所有特别行动队修正后
战斗力之和最大。试求出这个最大和。
例如,你有4名士兵,x1= 2, x2 = 2, x3 = 3, x4 = 4。经验公式中的参数为a= –1,
b= 10, c = –20
。此时,最佳方案是将士兵组成3个特别行动队:第一队包含士兵
1
和士兵2,第二队包含士兵3,第三队包含士兵4。特别行动队的初始战斗力分
别为4,3, 4,修正后的战斗力分别为4,1, 4。修正后的战斗力和为9,没有其它
方案能使修正后的战斗力和更大。
【输入格式】
输入由三行组成。第一行包含一个整数n,表示士兵的总数。第二行包含三

个整数a,b, c,经验公式中各项的系数。第三行包含n个用空格分隔的整数x1,

x2,…, xn,分别表示编号为1,2, …, n的士兵的初始战斗力。

【输出格式】
输出一个整数,表示所有特别行动队修正后战斗力之和的最大值。
【样例输入】
4
-110 -20
2 2 3 4
【样例输出】
9
【数据范围】
20%
的数据中,n≤ 1000
50%
的数据中,n≤ 10,000
100%
的数据中,1≤ n ≤ 1,000,000,–5≤ a ≤ –1|b|≤ 10,000,000|c|≤
10,000,000
1≤ x ≤ 100



【题解】

f[i]表示从1i划分为若干组能取得的最大战斗力。

f[i]=max{f[j]+G(s[i]-s[j])}。(0<=j<i

其中s[i]1i的战斗力之和,G(x)=A*x^2+B*x+C;直接动规O(n^2),预计40~50分。

这种1D/1D的方程可以用斜率优化。

设有决策j,k。其中ji的最优决策,k为除j外任意决策。

f[j]+A*(s[i]-s[j])^2+B*(s[i]-s[j])+C>=f[k]+A*(s[i]-s[k])^2+B*(s[i]-s[k])+C

=>f[j]-f[k]+A*sj*sj-A*sk*sk-B*sj+B*sk>=2A*si*(sj-sk)

j<k时,(f[j]-f[k]+A*sj*sj-A*sk*sk-B*sj+B*sk/(sj-sk)<=2A*si

H[j,k]=f[j]-f[k]+A*sj*sj-A*sk*sk-B*sj+B*sk/(s[j]-s[k])

依据此函数与2A*s[i]的比较就可知两个决策的优劣。

然后维护一个决策值序列d,其中d1<d2<d3<……<dk。满足H[d1,d2]>=H[d2,d3]>=……>=H[dk-1,dk]

每次在队首剔除H[dl,dl+1]>2*A*si的决策,直至H[dl,dl+1]恰好小于等于2A*si,此时的dl就是最优决策。

只要多做此类的题就可以了。

除法的效率极低,耗时为乘法的几十倍,建议将一切除法改为乘法。

代码照标程序打(没法,太菜了) 


#include <iostream>
#include <cstdio>

using namespace std;

const  int maxn=1000000+5;
long long f[maxn],s[maxn];
int q[maxn];
int n;
long long a,b,c;

long long calc(int j,int k)
{
     return f[j]-f[k]+a*(s[j]*s[j]-s[k]*s[k])-b*(s[j]-s[k]);
}

int main()
{
    freopen("commando.in","r",stdin);
    freopen("commando.out","w",stdout);
    scanf("%d\n",&n);
    scanf("%I64d %I64d %I64d",&a,&b,&c);
    for (int i=1;i<=n;++i)
    {
        scanf("%d",&s[i]);
        s[i]+=s[i-1];
    }
    long long x;
    int y,l=0,r=0;
    for (int i=1;i<=n;++i)
    {
        while (l<r&&calc(q[l],q[l+1])<=2*a*s[i]*(s[q[l]]-s[q[l+1]]))l++;
        y=q[l];
        x=s[i]-s[y];
        f[i]=f[y]+a*x*x+b*x+c;
        while (l<r&&calc(q[r-1],q[r])*(s[q[r]]-s[i])<=(s[q[r-1]]-s[q[r]])*calc(q[r],i))r--;
        q[++r]=i;
    }
    printf("%I64d",f[n]);
    return 0;
}

BZOJ 1911: [Apio2010]特别行动队(斜率优化)

BZOJ 1911: [Apio2010]特别行动队题目概述:有n个士兵,其战斗力为x,某一段连续的士兵的修正战斗力为 题目分析:1.先思考一下朴素的DP, ans[i]表示到i点的最大修正战...

bzoj1911: [Apio2010]特别行动队(斜率优化+Dp)

题目传送门 好开心又独立做了一道斜率优化原来斜率优化也不是辣么难嘛。。解法: 这道题Dp方程还是蛮好想的。 s[i]表示1到i的和。 f[i]=min(f[i],f[j]+a*(s[i]-s[...
  • Hanks_o
  • Hanks_o
  • 2017年10月26日 15:48
  • 46

[APIO2010]特别行动队(斜率优化dp)

【题解】 设s[i]=x[1]+x[2]+……+x[i] 则:f[0]=0     f[i]=max{ f[j]+zdl(s[i]-s[j]) } ( i>0 , 0      若 j不比k...
  • cjk_cjk
  • cjk_cjk
  • 2015年05月06日 02:22
  • 303

[斜率DP优化]BZOJ 1911: [Apio2010]特别行动队 题解

(传送门)题目大意给出一个有n个数的数组aia_i,要求将a分割成若干个区间,每个区间的价值是 令这个区间所有数的和为x,则该区间价值为A*x*x+B*x+C 一个分割方案ans=所有分割的区间的...

【APIO2010】特别行动队(斜率优化dp)

一路走来形单影只 转身回望大家已经远去 即使如此我依然前行 这正是我的坚强...

BZOJ 1911([Apio2010]特别行动队-斜率优化第一题)

1911: [Apio2010]特别行动队 Time Limit: 4 Sec  Memory Limit: 64 MB Submit: 1684  Solved: 706 [Submit][S...

BZOJ 1911-特别行动队(DP+斜率优化)

1911: [Apio2010]特别行动队 Time Limit: 4 Sec  Memory Limit: 64 MB Submit: 4662  Solved: 2253 [Submit][Sta...

Bzoj 1911: [Apio2010]特别行动队(斜率优化)

1911: [Apio2010]特别行动队 Time Limit: 4 Sec Memory Limit: 64 MB Description Input O...

bzoj 1911 [Apio2010]特别行动队斜率优化+DP

bzoj 1911 特别行动队 传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1911 Description Input ...
  • DJS_K_D
  • DJS_K_D
  • 2015年02月09日 22:42
  • 274

[省选前题目整理][BZOJ 1911][APIO 2010]特别行动队(斜率优化DP)

题目链接http://www.lydsy.com/JudgeOnline/problem.php?id=1911思路 注:单调递减的队列里,相邻两个点的连线的斜率是单调递增的,即维护一个下凸壳代码#...
  • qpswwww
  • qpswwww
  • 2015年03月31日 20:35
  • 469
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【斜率优化】特别行动队
举报原因:
原因补充:

(最多只允许输入30个字)