标签：int sum hdu3507 斜率 出队 maxn 看做

不难得到递推式f[i]=min{f[j]+(sum[i]-sum[j])^2+m};

去掉min函数并展开：f[i]=f[j]+sum[i]^2-2*sum[i]*sum[j]+sum[j]^2+m

将含i的项与含j的项分离，并把单纯含j的项写在左边：f[j]+sum[j]^2=2*sum[i]*sum[j]+f[i]-sum[i]^2-m

现在把f[j]+sum[j]^2看做y，2*sum[i]看做k，sum[j]看做x，f[i]-sum[i]^2-m看做b。

对于每一个点(sum[j],f[j]+sum[j]^2)都是固定的，对于每一个i，斜率2*sum[i]也是固定的，而截距也就是b。

b越小，f[i]就越小。我们可以形象的理解为：一条斜率为2*sum[i]的线从下往上扫，扫到的第一个点就是答案。

我们构造一个队列，维护相邻两点间的斜率。如果能保证斜率单调递增，那么第一个满足斜率大于2*sum[i]的就是结果。

由于2*sum[i]也是单调递增的，所以不大于它的可以出队。时间复杂度N。

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=500000+10;
int n,m,x,s[maxn],f[maxn],q[maxn];
int yval(int x,int y){
    return f[y]-f[x]+s[y]*s[y]-s[x]*s[x];
}//计算y坐标的差
int xval(int x,int y){
    return s[y]-s[x];
}//计算x坐标的差
signed main(){
    while(cin>>n>>m){
        for(int i=1;i<=n;i++){
            scanf("%lld",&x);
            s[i]=s[i-1]+x;
        }
        memset(f,0x3f,sizeof(f));
        memset(q,0,sizeof(q));
        int l=1,r=1;
        q[l]=0,f[0]=0;
        for(int i=1;i<=n;i++){
            while(l<r&&yval(q[l],q[l+1])<=xval(q[l],q[l+1])*2*s[i])l++;
　　　　　　　//当队首不满足斜率大于当前斜率，则出队
            f[i]=f[q[l]]+(s[i]-s[q[l]])*(s[i]-s[q[l]])+m;
　　　　　　　//计算f[i]
            while(l<r&&yval(q[r-1],q[r])*xval(q[r],i)>=xval(q[r-1],q[r])*yval(q[r],i))r--;
　　　　　　　//若队尾不满足单调性，则出队
            q[++r]=i;
        }
        printf("%lld\n",f[n]);
    }
    return 0;
}

标签：int,sum,hdu3507,斜率,出队,maxn,看做
来源： https://www.cnblogs.com/syzf2222/p/12386788.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。