5380: Function 斜率优化

题解：

斜率优化真的是差的一塌糊涂……这题搞了一个下午，晚上才清醒点AC了。
题目显然就是求从$y$位置开始往右走，中途每走到一个格子可以至少获得一次它的权值，获得$x$次权值的最小值。那么策略显然就是走到一个比较优秀的格子停下来，假设走到$i$停下来，那么$v=sum_y-sum_i+a_i\times i+a_i\times(x-y)$，移项得$a_i\times(y-x)+v=sum_y-sum_i+a_i\times i$，那么这个显然是可以斜率优化的，因为要令截距$v$最小，所以维护一个下凸包，但是有个问题，这些点的横坐标不是单调的，不过我们发现，如果前面的$a_i$比后面的$a_j$大，那么$a_i$是没用的，利用这个可以使点的横坐标单调。求答案的时候由于给出的斜率不单调，所以要在凸包上二分求解。注意求解时并非整个凸包上的点都是可以用的，因为从$y$开始最多只能走$x$步，可以二分出合法左端点。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=500010;
const int inf=2147483647;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
int n,q,a[Maxn],sta[Maxn],top=0;
LL sum[Maxn],ans[Maxn],b[Maxn];
struct Query{int x,y,id;}Q[Maxn];
bool cmp(Query a,Query b){return a.y<b.y;}
double K(int x,int y){return (double)(b[x]-b[y])/(a[x]-a[y]);}
int binary_search(int pos)
{
    int l=1,r=top;
    while(l<=r)
    {
        int mid=l+r>>1;
        if(sta[mid]>=pos)r=mid-1;
        else l=mid+1;
    }
    return r+1;
}
int main()
{
    n=read();sum[0]=0;
    for(int i=1;i<=n;i++)
    {
        a[i]=read();
        sum[i]=sum[i-1]+a[i];
        b[i]=(LL)i*a[i]-sum[i];
    }
    q=read();
    for(int i=1;i<=q;i++)Q[i].x=read(),Q[i].y=read(),Q[i].id=i;
    sort(Q+1,Q+1+q,cmp);
    int now=0;
    for(int i=1;i<=q;i++)
    {
        while(now<Q[i].y)
        {
            now++;
            while(top&&a[now]<=a[sta[top]])top--;
            while(top>1&&K(now,sta[top])<K(sta[top],sta[top-1]))top--;
            sta[++top]=now;
        }
        int l=binary_search(Q[i].y-Q[i].x+1)+1,r=top;
        int k=Q[i].y-Q[i].x;
        while(l<=r)
        {
            int mid=l+r>>1;
            if(k>=K(sta[mid],sta[mid-1]))l=mid+1;
            else r=mid-1;
        }
        l--;
        ans[Q[i].id]=sum[Q[i].y]-sum[sta[l]]+(LL)a[sta[l]]*(Q[i].x-Q[i].y+sta[l]);
    }
    for(int i=1;i<=q;i++)printf("%lld\n",ans[i]);
}