码迷,mamicode.com
首页 > 其他好文 > 详细

HDU3507 Print Article(斜率优化+单调队列)

时间:2019-09-02 19:27:30      阅读:88      评论:0      收藏:0      [点我收藏+]

标签:i++   head   存在   for   def   时间复杂度   级别   out   end   

大佬的斜率优化博客,写的很好

https://www.cnblogs.com/orzzz/p/7885971.html

DP的斜率优化同四边形优化和数据结构优化类似,都是用一种方法快速找到最优的状态转移,以到达优化的目的。一次O(n)转移可以变成log或者常数级别的。

这个题斜率优化的步骤就是,比较不同的状态转移,写成斜率的形式,发现他们存在一个关系,只需要维护一个下凸的凸包结构,最优状态转移的点一定在这个凸包上,剔除了很多不可能是最优的状态转移点。这个题还存在一个单调性,维护单调队列即可,因为每次找到最优点前面的点都可在凸包上剔除。

这个题算斜率优化的模板,思路推导就见上面的博客吧,我也不会搞LaTeX。

1A代码:

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=5e5+10;
ll dp[N];
ll a[N];
ll pre[N];
ll que[N];
int n,m;
ll get_up(int i,int j){
    return dp[j]+pre[j]*pre[j]-dp[i]-pre[i]*pre[i];
}
ll get_down(int i,int j){
    return pre[j]-pre[i];
}
int main()
{
    while(~scanf("%d%d",&n,&m)){
        pre[0]=0;
        for(int i=1;i<=n;i++){
            scanf("%lld",&a[i]);
            pre[i]=pre[i-1]+a[i];
        }
        int head=0,tail=0;
        que[tail]=0;
        for(int i=1;i<=n;i++){
            while(head<tail&&get_up(que[head],que[head+1])<=2*pre[i]*get_down(que[head],que[head+1]))
                head++;
                //找到一个最优的转移,时间复杂度为为常数,保证队列非空
            dp[i]=dp[que[head]]+(pre[i]-pre[que[head]])*(pre[i]-pre[que[head]])+m;
            while(head<tail&&get_up(que[tail],i)*get_down(que[tail-1],que[tail])<=get_up(que[tail-1],que[tail])*get_down(que[tail],i))
                tail--;
            que[++tail]=i;

        }
        cout<<dp[n]<<endl;
    }
    return 0;
}

 

HDU3507 Print Article(斜率优化+单调队列)

标签:i++   head   存在   for   def   时间复杂度   级别   out   end   

原文地址:https://www.cnblogs.com/gzr2018/p/11448195.html

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