Description
韩父有N个儿子,分别是韩一,韩二…韩N。由于韩家演技功底深厚,加上他们间的密切配合,演出获得了巨大成功,票房甚至高达2000万。舟子是名很有威望的公知,可是他表面上两袖清风实则内心阴暗,看到韩家红红火火,嫉妒心遂起,便发微薄调侃韩二们站成一列时身高参差不齐。由于舟子的影响力,随口一句便会造成韩家的巨大损失,具体亏损是这样计算的,韩一,韩二…韩N站成一排,损失即为C*(韩i与韩i+1的高度差(1<=i< N))之和,搞不好连女儿都赔了.韩父苦苦思索,决定给韩子们内增高(注意韩子们变矮是不科学的只能增高或什么也不做),增高1cm是很容易的,可是增高10cm花费就很大了,对任意韩i,增高Hcm的花费是H^2.请你帮助韩父让韩家损失最小。
Input
有若干组数据,一直处理到文件结束。
每组数据第一行为两个整数:韩子数量N(1<=N<=50000)和舟子系数C(1<=C<=100)
接下来N行分别是韩i的高度(1<=hi<=100)。
Output
对每组测试数据用一行输出韩家的最小损失。
Sample Input
5 2
2
3
5
1
4
Sample Output
15
#include<stdio.h>
#include<algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
int dp[2][101],q[101],h,t,cur,x,n,c,ans,k;
inline const int read()
{
register int x=0,f=1;
register char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
int main()
{
while(scanf("%d%d",&n,&c)!=EOF)
{ ans=inf;
x=read();
cur=0;//cur 是为了滚动数组的
for(int i=0;i<x;i++) dp[cur][i]=inf;//不能变矮,赋为超大值
for(int i=x;i<=100;i++) dp[cur][i]=(i-x)*(i-x);//增高的花费
for(register int i=1;i<n;i++)
{
x=read(),cur=cur^1,h=1,t=0;//单调队列初始化
for(int j=0;j<=100;j++)//前一个比当前高
{
k=dp[cur^1][j]-j*c;
while(h<=t&&q[t]>k) t--;
q[++t]=k;
if(j<x) dp[cur][j]=inf;
else dp[cur][j]=q[h]+j*c+(x-j)*(x-j);
}
h=1,t=0;
for(int j=100;j>=0;j--)//前一个比当前矮
{
k=dp[cur^1][j]+j*c;
while(h<=t&&q[t]>k) t--;
q[++t]=k;
if(j>=x) dp[cur][j]=min(dp[cur][j],q[h]-j*c+(x-j)*(x-j));
}
}
for(int i=0;i<=100;i++) ans=min(ans,dp[cur][i]);
printf("%d\n",ans);
}
return 0;
}
其实保留一个最小值就够了
#include <cstdlib>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define inf 0x3fffffff
using namespace std;
int n,c,h[50005],dp[2][105],cur,limit,Min;
inline const int read()
{
register int x=0,f=1;
register char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
int main()
{
while (scanf("%d%d",&n,&c)==2) {
cur=1;
limit=-inf;
Min=inf;
for(register int i=1;i<=n;++i){h[i]=read();limit=max(limit,h[i]);}
for(register int i=0;i<=1;++i) memset(dp[i],0x3f,sizeof(dp[i]));
for(int i=h[1];i<=limit;++i) dp[cur][i]=(i-h[1])*(i-h[1]);
for(register int i=2;i<=n;++i) {
cur=cur^1;
Min=inf;
for(int j=h[i-1];j<=limit;++j) {
Min=min(Min,dp[cur^1][j]-c*j);
if (j>=h[i]) dp[cur][j]=Min+c*j+(j-h[i])*(j-h[i]);
}
Min=inf;
for(int j=limit;j>=h[i];--j) {
Min=min(Min,dp[cur^1][j]+c*j);
if(j>=h[i]) dp[cur][j]=min(dp[cur][j],Min-c*j+(j-h[i])*(j-h[i]));
}
memset(dp[cur^1],0x3f,sizeof(dp[cur^1]));
}
for (int i=h[n];i<=limit;++i) {
Min=min(Min,dp[cur][i]);
}
printf("%d\n",Min);
}
return 0;
}