poj 3612

DP+简单优化，状态转移方程容易看出，优化看了http://hi.baidu.com/forverlin1204/blog/item/df9c62dc0fb16bd38d1029d8.html。哎，太粗心了，初始化时忘了得平方，看了1个多小时才发现

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int inf=2000000000;
const int maxn=100000+10;
int f[100+10],a[maxn],high[maxn],low[maxn];
int n,c;
int main()
{
    while(~scanf("%d%d",&n,&c))
    {
    int i;
    int maxv=0;
    for(i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        maxv=max(maxv,a[i]);
    }
    int j;
    for(i=0;i<a[1];i++) f[i]=inf;
    for(i=a[1];i<=maxv;i++) f[i]=(i-a[1])*(i-a[1]);
    for(i=2;i<=n;i++)
    {
        low[0]=inf;
        for(j=1;j<=maxv;j++)
            low[j]=min(f[j]-j*c,low[j-1]);
        high[maxv+1]=inf;
        for(j=maxv;j>=1;j--)
        {
            high[j]=min(f[j]+j*c,high[j+1]);
            f[j]=inf;
        }
        for(j=a[i];j<=maxv;j++)
            f[j]=min(low[j]+j*c+(j-a[i])*(j-a[i]),high[j]-j*c+(j-a[i])*(j-a[i]));

    }
    int ans=inf;
    for(j=a[n];j<=maxv;j++) ans=min(ans,f[j]);
    printf("%d\n",ans);
    }
    return 0;
}