poj3709题目传送
看到这名字第一反应是那个面具
斜率优化的动态规划,第一次做斜率优化的dp题。。。而且自己也不会,是看着书慢慢弄懂的。书是《挑战》p341 ,这一面有很多错误,刚开始自己写完后照书上一对,一模一样后发现答案很不对,后来才发现书上的状态转移方程存在错误(在下面的公式里你自己对一遍就知道了) 总之就是,书上的不一定对,自己领悟并理解才最重要。
首先不难得出, a0 是最小的值,为了使操作的次数较小,我们不可能把最小值拿着减。所以我们发现,对于 a0 这样的较小的项,我们可以选择从小到大来选择需要减少到 a0 的项。设 dpi 为只考虑前i项的前提下最少的操作次数,则有 dpi =min{ dpj +( a[j+1] - a[j] )+…+( a[i−1] - a[j] )},其中 0≤j≤i−k 。
直接计算自然没问题,但是 O3 的复杂度不太好的样子,于是考虑优化:
例如我们记一个前缀和,
有 dpi =min{ dpj + S[i] - S[j+1] - aj *(i-j-1)}
观察方程我们可以发现,若我们设一个函数(用x替换i) fj(x)=−aj∗x+dp[j]−S[j+1]+aj∗(j+1) ,则有 dp[i]=S[i]+minfj(i) ,于是计算dp[n]就变为了在i-k+1条直线中寻找x=i的最小值就好了。
于是我们考虑用双端队列去维护所有可能成为最小值的直线的集合,并以从小到大的顺序排列。那么我们如何判断要去掉某条直线呢?
设:
f1(x)=a1∗x+b1
f2(x)=a2∗x+b2
f3(x)=a3∗x+b3
我们有
(a2−a1)∗(b3−b2)≥(b2−b1)∗(a3−a2)
时,
f2
为要排除的直线,为什么呢?
简证:当
f1
与
f3
的交点在
f2
以下时,联立
f1
与
f3
得
x0=(b1−b3a3−a1)
然后对于
f1
,
f3
,有
f(x0)=a1(b1−b3a3−a1)+b1=(b1∗a3−b3∗a1a3−a1)
而
f2(x0)=a2∗(b1−b3a3−a1)+b2
当
f2(x0)≤f1(x0)和f3(x0)
时,有
a2∗(b1−b3a3−a1)+b2≤(b1−b3a3−a1)
拿纸来算吧 化简就有
(a2−a1)∗(b3−b2)≥(b2−b1)∗(a3−a2)
了。
所以我们就可以根据上述方法来做双端队列,详见标程:
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
typedef long long ll;
const ll maxn=505000;
ll cas=0;
ll n,k;
ll a[maxn];
ll deq[maxn];//虽然是双端队列 但存的是值
ll dp[maxn*10];
ll sh[maxn];//前缀和
bool check(ll f1,ll f2,ll f3)
{
ll a1=-a[f1];
ll b1=dp[f1]-sh[f1+1]+a[f1]*(f1+1);
ll a2=-a[f2];
ll b2=dp[f2]-sh[f2+1]+a[f2]*(f2+1);
ll a3=-a[f3];
ll b3=dp[f3]-sh[f3+1]+a[f3]*(f3+1);
return (a2-a1)*(b3-b2)>=(b2-b1)*(a3-a2);
}
ll f(ll j,ll x)
{
return -a[j]*x+dp[j]-sh[j+1]+a[j]*(j+1);
}
int main()
{
//freopen("1.in","r",stdin);
//freopen("2.out","w",stdout);
scanf("%lld",&cas);
for(ll z=0;z<cas;z++)
{
memset(sh,0,sizeof(sh));
memset(a,0,sizeof(a));
memset(dp,0,sizeof(dp));
memset(deq,0,sizeof(deq));
scanf("%lld%lld",&n,&k);
for(ll i=0;i<n;i++)
{
scanf("%lld",&a[i]);
}
for(ll i=0;i<n;i++)
{
sh[i+1]=sh[i]+a[i];
}
ll s=0;
ll t=1;
deq[s]=0;
dp[0]=0;
for(ll i=k;i<=n;i++)
{
if(i-k>=k)
{
while(s+1<t && check(deq[t-2],deq[t-1],i-k))
{
t--;
}
deq[t]=i-k;
t++;
}
while(s+1<t && f(deq[s],i)>=f(deq[s+1],i))
{
s++;
}
dp[i]=sh[i]+f(deq[s],i);
}
printf("%lld\n",dp[n]);
}
return 0;
}
书上没看懂得就来看一下这个吧QAQ