bzoj1911: [Apio2010]特别行动队
2017-07-09 21:59
176 查看
传送门
首先得到转移方程:
f[i]=max(f[j]+a*(sum[i]-sum[j-1])^2+b*(sum[i]-sum[j-1])+c)
然后我们斜率优化一发:
如果j>k且j比k更优
f[j]-f[k]+a*sum[j]^2-a*sum[k]^2+b*(sum[k]-sum[j])>2*a*(sum[j]-sum[k])*sum[i]
然后……….就没有了
首先得到转移方程:
f[i]=max(f[j]+a*(sum[i]-sum[j-1])^2+b*(sum[i]-sum[j-1])+c)
然后我们斜率优化一发:
如果j>k且j比k更优
f[j]-f[k]+a*sum[j]^2-a*sum[k]^2+b*(sum[k]-sum[j])>2*a*(sum[j]-sum[k])*sum[i]
然后……….就没有了
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cstdlib> #include<cmath> #define ll long long #define N 1000005 using namespace std; ll sum ,f ,n,a,b,c,x; int q ,l,r; inline int read(){ int k=0; char ch=getchar(); for (;ch<'0'||ch>'9';ch=getchar()); for (;ch>='0'&&ch<='9';ch=getchar()) k=k*10+ch-48; return k; } ll sqr(ll x){return x*x;} double slop(int k,int j){ return (double)(f[j]-f[k]+a*(sqr(sum[j])-sqr(sum[k]))+b*(sum[k]-sum[j]))/(double)(2*a*(sum[j]-sum[k])); } int main(){ n=read(); scanf("%lld%lld%lld",&a,&b,&c); for (int i=1;i<=n;i++){ x=read(); sum[i]=sum[i-1]+x; } for (int i=1;i<=n;i++){ while (l<r&&slop(q[l],q[l+1])<sum[i]) l++; int t=q[l]; f[i]=f[t]+a*sqr(sum[i]-sum[t])+b*(sum[i]-sum[t])+c; while (l<r&&slop(q[r-1],q[r])>slop(q[r],i)) r--; q[++r]=i; } printf("%lld",f ); }
相关文章推荐
- 【BZOJ 1911】【APIO 2010】特别行动队
- 【bzoj 1911】 [Apio2010]特别行动队
- [BZOJ1911][Apio2010]特别行动队
- bzoj1911【APIO2010】特别行动队
- bzoj 1911 [Apio2010]特别行动队
- 斜率优化专题4——bzoj 1911: [Apio2010] 特别行动队 题解
- BZOJ 1911: [Apio2010]特别行动队
- 【斜率优化】BZOJ1911(Apio2010)[特别行动队]题解
- [斜率优化]BZOJ 1911—— [Apio2010]特别行动队
- [BZOJ1911][Apio2010]特别行动队(斜率优化dp)
- [斜率DP优化]BZOJ 1911: [Apio2010]特别行动队 题解
- bzoj1911 [Apio2010]特别行动队
- Bzoj1911 [Apio2010]特别行动队
- BZOJ 1911 [Apio2010]特别行动队(单调队列优化DP)
- BZOJ 1911: [Apio2010]特别行动队
- BZOJ 1911: [Apio2010]特别行动队 斜率优化dp
- [BZOJ 1911][Apio2010]特别行动队:DP斜率优化
- [BZOJ1010][HNOI2008]玩具装箱[BZOJ1911][APIO2010]特别行动队[BZOJ1492][NOI2007]货币兑换 斜率优化
- 【BZOJ1911】【APIO2010】特别行动队(斜率优化DP)
- BZOJ1911 Apio2010 特别行动队