poj 3612
DP+简单优化，状态转移方程容易看出，优化看了http://hi.baidu.com/forverlin1204/blog/item/df9c62dc0fb16bd38d1029d8.html。哎，太粗心了，初始化时忘了得平方，看了1个多小时才发现
1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 using namespace std; 5 const int inf=2000000000; 6 const int maxn=100000+10; 7 int f[100+10],a[maxn],high[maxn],low[maxn]; 8 int n,c; 9 int main() 10 { 11 while(~scanf("%d%d",&n,&c)) 12 { 13 int i; 14 int maxv=0; 15 for(i=1;i<=n;i++) 16 { 17 scanf("%d",&a[i]); 18 maxv=max(maxv,a[i]); 19 } 20 int j; 21 for(i=0;i<a[1];i++) f[i]=inf; 22 for(i=a[1];i<=maxv;i++) f[i]=(i-a[1])*(i-a[1]); 23 for(i=2;i<=n;i++) 24 { 25 low[0]=inf; 26 for(j=1;j<=maxv;j++) 27 low[j]=min(f[j]-j*c,low[j-1]); 28 high[maxv+1]=inf; 29 for(j=maxv;j>=1;j--) 30 { 31 high[j]=min(f[j]+j*c,high[j+1]); 32 f[j]=inf; 33 } 34 for(j=a[i];j<=maxv;j++) 35 f[j]=min(low[j]+j*c+(j-a[i])*(j-a[i]),high[j]-j*c+(j-a[i])*(j-a[i])); 36 37 } 38 int ans=inf; 39 for(j=a[n];j<=maxv;j++) ans=min(ans,f[j]); 40 printf("%d\n",ans); 41 } 42 return 0; 43 }
