玩具装箱 bzoj1010 斜率优化

斜率优化的题好像都是这样的方程：左边关于j，k的一个(...)/(...)的式子，右边是个只与i有关的可算的数字；

然后把它放到二维坐标轴上，用单调队列维护一个凸壳，O（n）的复杂度；

这道题但是我发现我wrong了，找了程序看了一下，才发现斜率优化还有一点没理解；才明白上午T2能A是由于数据太水，出题人万岁！

#include<iostream>
#include<cstdio>
using namespace std;
#define LL long long
const int maxn=50500;
LL n,L,a[maxn],g[maxn],q[maxn],tail=0,head=1,f[maxn];
void init(){
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
    scanf("%d%d",&n,&L);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        g[i]=g[i-1]+a[i];
    }
}
double p(LL k,LL j){
    return double(f[j]+g[j]*g[j]-f[k]-g[k]*g[k])/(double)(2*(g[j]-g[k]));
}
LL squ(LL x){return x*x;}
void work(){
    for(int i=1;i<=n;i++)g[i]+=i;
    q[++tail]=0;
    for(int i=1;i<=n;i++){
        while(head<tail&&p(q[head],q[head+1])<g[i]-1-L)head++;
        cout<<q[head]<<endl;
        f[i]=f[q[head]]+squ(g[i]-g[q[head]]-L-1);
        while(head<tail&&p(q[tail],i)<p(q[tail],q[tail-1]))tail--;
        q[++tail]=i;
    }
    cout<<f[n]<<endl;
}
int main(){
    init();
    work();
}