51nod 1488 帕斯卡小三角

前言

文化课弱鸡选手又来做题啦。。
(⊙o⊙)…尽量保证一周一或两题？（大雾）
可能过多几个星期文化课补上来就可以多做点题了吧。。

题意

· f(1,j)=a[j], 1≤j≤n.
· f(i,j)=min(f(i-1,j),f(i-1,j-1))+a[j], 2≤i≤n, i≤j≤n.
a是一个长度为n的数组。现在有若干个询问，输入x,y，求f(x,y)

题解

随意观察一下可以发现，这个东西其实就是你在第一行某一个找一个点x
然后你每一次可以往下走或者往右下走，问你走到这个点的代价最小
然后再胡乱观察一下，可以发现肯定是从起点一直往下走，然后剩下一段往右下走。。
然后问题就变成了要找一个s使得
c是前缀和
$c[y]-c[s]+a[s]*(y-x+s)$最小\
然后把式子整理一下
$c[y]+a[s]*(y-x)+(a[s]*s-c[s])$
很容易发现，如果一个a比前面的a要小，那么前面肯定是没有用的
这个就可以斜率优化了
维护一个下凸的凸包
每一次询问二分找最值就可以了
其实s的选择是有一个区间限制的
就是$[max(0,y-x),y]$
我们按y排序，每一次加入新的就可以了
然后听说不管左边界也可以AC。。然后我就偷工减料了。。
其实要管的话，就把凸包的那个单调栈换成set就可以了
时间复杂度还是一个$log$的，就是常数大一点点
一开始没有读优过不去啊
CODE:

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
using namespace std;
const int N=100005;
int a[N],c[N];
int ans[N];
struct qq
{
    int x,y,id;
}s[N];
int n,Q;
bool cmp (qq x,qq y)    {return x.y<y.y;}
int q[N];
int st,ed,now;
double Y (int x)    {return a[x]*x-c[x];}
double X (int x)    {return -a[x];}
double get_k (int x,int y)
{
    return (double)(Y(x)-Y(y))/(double)(X(x)-X(y));
}
int calc (int xx,int x,int y)//这个点行不行 
{
    return c[y]-c[xx]+a[xx]*(x-y+xx);
}
void solve ()
{
    st=1;ed=0;now=1;
    for (int u=1;u<=n;u++)
    {
        while (st<=ed&&a[q[ed]]>=a[u]) ed--;
        while (st<ed&&get_k(q[ed-1],q[ed])<get_k(q[ed],u))  ed--;   
        q[++ed]=u;
        while (now<=Q&&s[now].y==u)//可以用到这个了 
        {
            int l=st,r=ed,lalal;
            while (l<=r)
            {
                int mid=(l+r)>>1;
                if (mid==ed||calc(q[mid],s[now].x,s[now].y)<calc(q[mid+1],s[now].x,s[now].y)) {lalal=mid;r=mid-1;}
                else l=mid+1;
            }
            ans[s[now].id]=calc(q[lalal],s[now].x,s[now].y);
            now++;
        }
    }
}
int read()
{
    char ch=getchar();int x=0;
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=getchar();}
    return x;
}
int main()
{
    n=read();
    for (int u=1;u<=n;u++)  {a[u]=read();c[u]=c[u-1]+a[u];}
    Q=read();
    for (int u=1;u<=Q;u++)
    {
        s[u].x=read();s[u].y=read();
        s[u].id=u;
    }
    sort(s+1,s+1+Q,cmp);
    solve();
    for (int u=1;u<=Q;u++)  printf("%d\n",ans[u]);
    return 0;
}