bzoj3437

Posted 123456

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj3437相关的知识,希望对你有一定的参考价值。

斜率优化

斜率优化是指对于dp[i]=max/min(dp[j]+a[i]*b[j]+c[j])这样的方程的优化

-a[i]*b[j]+dp[i]=dp[j]+c[j]

把b[j]看成x,dp[j]+c[j]看成y

-a[i]=k dp[i]=b

这就是一次函数的形式,可以看成用斜率为-a[i]的直线交(x,y),dp就是截距

对于最大值,我们自然希望截距大,那么就是上凸壳,否则就是下凸壳

当(x,y)中x单调且-a[i]单调的时候我们可以用队列维护凸壳,因为x单调所以我们可以直接在队尾添加点,又因为-a[i]单调所以每次选的决策肯定是单调的

如果-a[i]不单调那么我们就得在凸壳上二分

如果x不单调,那么我们就得用cdq或者splay

这道题的dp式是dp[i]=dp[j]+(s2[i]-s2[j])*i-(s1[i]-s1[j])+a[i]

s2是b[i]的前缀和,s1[i]是b[i]*i的前缀和

那么把只和i有关的去掉

s2[j]*i+dp[i]=dp[j]+s1[j]

s2[j]递增 i递增

由于要最小值,维护下凸壳

技术分享图片
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e6 + 5;
int n;
ll a[N], b[N], s1[N], s2[N], dp[N];
int q[N];
double slope(int i, int j) 
{
    return (double)(dp[i] - dp[j] + s1[i] - s1[j]) / (double)(s2[i] - s2[j]);
}
int main()
{
    scanf("%d", &n);
    for(int i = 1; i <= n; ++i) scanf("%lld", &a[i]);
    for(int i = 1; i <= n; ++i) scanf("%lld", &b[i]), s1[i] = s1[i - 1] + b[i] * i, s2[i] = s2[i - 1] + b[i];
    int l = 1, r = 1;
    for(int i = 1; i <= n; ++i) 
    {
        while(l < r && slope(q[l], q[l + 1]) <= i) ++l;
        dp[i] = dp[q[l]] + (s2[i] - s2[q[l]]) * i - (s1[i] - s1[q[l]]) + a[i];
        while(l < r && slope(q[r], i) <= slope(q[r], q[r - 1])) --r;
        q[++r] = i;
    }
    printf("%lld\n", dp[n]);
    return 0;
}
View Code

 

以上是关于bzoj3437的主要内容,如果未能解决你的问题,请参考以下文章

bzoj 3437: 小P的牧场 -- 斜率优化

BZOJ3437 小P的牧场(斜率优化dp)

bzoj3437

BZOJ 3437: 小P的牧场 斜率优化DP

bzoj 3437小P的牧场

3437: 小P的牧场