[BZOJ3156]防御准备（斜率优化dp）

题目：
我是超链接
题解：
woc！调了一大节课，还对拍样程（“完全没有问题啊，拍了一晚上都不带错的“），结果看好数量级之后强转longlong才A了
代码：
第一种解法：两个点之间的距离是（i-j）（i-j-1)/2
方程：f[i]=min(f[j]+（i-j）（i-j-1)/2)+a[i]

#include <cstdio>
#include <iostream>
#define LL long long
#define N 1000005
using namespace std;
int n,i;
LL f[N];LL a[N];int q[N];
LL K(int j){return -j;}
LL B(int j){return (LL)j*(j+1)/2+(LL)f[j];}
LL Y(int i,int j){return (LL)i*(LL)K(j)+(LL)B(j);}
bool fg(int x1,int x2,int x3)
{
    LL w1=(LL)(K(x1)-K(x2))*(LL)(B(x3)-B(x1));
    LL w2=(LL)(K(x1)-K(x3))*(LL)(B(x2)-B(x1));
    return w1>=w2;
}
int main()
{
    scanf("%d",&n);
    for (i=n;i>=1;i--)
      scanf("%lld",&a[i]);
    LL ans=a[1]+(LL)n*(LL)(n-1)/(LL)2; 
    int l,r;
    l=r=0;
    f[1]=a[1];q[r]=1;
    for (i=2;i<=n;i++)
    {
        while (l<r && Y(i,q[l])>=Y(i,q[l+1])) l++;
        f[i]=(LL)Y(i,q[l])+(LL)i*(i-1)/2+(LL)a[i];
        ans=min(ans,f[i]+(LL)(n-i+1)*(LL)(n-i)/(LL)2);  
        while (l<r && fg(i,q[r],q[r-1])) r--;
        q[++r]=i;
    }
    printf("%lld\n",ans);
}

第二种解法：fi表示在i点建守卫塔的费用总和，转移方程：f[i]=min(f[j]+s[i-1]-s[j]-d[j]*(c[i-1]-c[j]))+a[i]
其中di表示1～i的距离，ci表示1～i的点数，s表示的是d的前缀和。
注意这里第一个必须选，所以开始的队列不能为0，应该为1。并且最后一个可以不选，单独枚举。

#include <cstdio>
#include <iostream>
#define LL long long
#define N 1000005
using namespace std;
LL f[N],a[N],s[N];int c[N],d[N],q[N];
LL K(int j){return -d[j];}
LL B(int j){return (LL)d[j]*c[j]-s[j]+f[j];}
LL Y(int i,int j){return K(j)*c[i-1]+(LL)B(j);}
bool cover(int x1,int x2,int x3)
{
    LL w1=(LL)(K(x1)-K(x2))*(LL)(B(x3)-B(x1));
    LL w2=(LL)(K(x1)-K(x3))*(LL)(B(x2)-B(x1));
    return w1<=w2;
}
int main()
{
    int n,i;
    scanf("%d",&n);
    for (i=n;i>=1;i--) scanf("%lld",&a[i]);
    d[1]=0; for (i=2;i<=n;i++) d[i]=d[i-1]+1;
    for (i=1;i<=n;i++) c[i]=c[i-1]+1;
    for (i=1;i<=n;i++) s[i]=s[i-1]+d[i];
    int l,r;l=r=0;
    LL ans=a[1]+(LL)n*(LL)(n-1)/2;
    q[r]=1;
    f[1]=a[1];
    for (i=2;i<=n;i++)
    {
        while (l<r && Y(i,q[l])>=Y(i,q[l+1])) l++;
        f[i]=(LL)Y(i,q[l])+a[i]+(LL)s[i-1];
        ans=min(ans,f[i]+(LL)(n-i)*(LL)(n-i+1)/2);
        while (l<r && cover(i,q[r-1],q[r])) r--;
        q[++r]=i;
    }
    printf("%lld",ans);
}