【DP】hdu3507

提议就是给出一个非负的序列，然后割成几段，每段的代价是那个公式。

然后dp式很好推也能一眼看出是斜率优化。

（d[i]-d[j]+sum[i]^2-sum[j]^2）/2(sum[i]-sum[j])<sum[k]   i>j

被玮神喷基本功不好了，

于是打算多写几个斜率优化练一下，

然后

1,起始项是q[1]=0不是q[1]=1;

2.传参的时候一定要对应，不要传值当做标号用

3.一定要注意队头出队时候写斜率的大小关系

4.不要手残把加号打成减号

5.样例坑爹

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
#include<set>
using namespace std;
long long f[701000],sum[701000];
int n,m,q[701000],a[700100];
bool xl(int x,int y,int z,int w)
{
	if (w==1) 
	{
		long long tmp1=(f[x]-f[y]+sum[x]*sum[x]-sum[y]*sum[y]);
		long long tmp2=(f[y]-f[z]+sum[y]*sum[y]-sum[z]*sum[z]);
		return (tmp1*(sum[y]-sum[z])>=tmp2*(sum[x]-sum[y]));
	}
	if (w==2) 
	{
		return f[y]-f[x]+sum[y]*sum[y]-sum[x]*sum[x]<sum[z]*2*(sum[y]-sum[x]);
	}
}
int main()
{
	freopen("test.in","r",stdin);
	int head,tail;
	while (~scanf("%d%d",&n,&m))
	{
		head=1;tail=1;
		memset(sum,0,sizeof(sum));
		memset(f,0,sizeof(f));
		memset(a,0,sizeof(a));
		memset(q,0,sizeof(q));
		for (int i=1;i<=n;i++)
		{
			scanf("%d",&a[i]);
			sum[i]=sum[i-1]+a[i];
		}
		q[1]=0;
		f[0]=0;
		sum[0]=0;
		for (int i=1;i<=n;i++)
		{
			while (tail>=head+1&&xl(q[head],q[head+1],i,2)) head++;
			f[i]=f[q[head]]+(sum[i]-sum[q[head]])*(sum[i]-sum[q[head]])+m;
			while (tail>=head+1&&xl(q[tail-1],q[tail],i,1)) tail--;
			q[++tail]=i;
		}
		printf("%I64d\n",f[n]);
	}
	return 0;
}