BZOJ 1010 斜率优化DP

题意：

http://www.lydsy.com/JudgeOnline/problem.php?id=1010

思路：

DP方程不难想：
DP[I]=MIN(DP[J]+(SUM[I]-SUM[J]+I-J-1-L)^2) ,J<I
其中DP[I]表示处理到第I个玩具时最少的费用。

 

直接来显然TLE，一维方程一般用斜率优化。

先化简下方程：
令F[I]=SUM[I]+I,C=1+L,则化为
DP[I]=MIN(DP[J]+(F[I]-F[J]-C)^2)

 

 

以下是用这方法的步骤：
1.证明较优决策点对后续状态影响的持续性
假设在状态I处,有两个决策点J,K(J<K),且决策K比决策J好,即
   DP[K]+(F[I]-F[K]-C)^2<=DP[J]+(F[I]-F[J]-C)^2
则,对于状态I之后的某状态T,  F[T]=F[I]+V
要证
DP[K]+(F[T]-F[K]-C)^2<=DP[J]+(F[T]-F[J]-C)^2
只需证
DP[K]+(F[I]+V-F[K]-C)^2<=DP[J]+(F[I]+V-F[J]-C)^2
只需证
DP[K]+(F[I]-F[K]-C)^2+2*(F[I]-F[K]-C)*V+V^2<=DP[J]+(F[I]-F[J]-C)^2+2*(F[I]-F[J]-C)*V+V^2
联系假设只需要证:
2*(F[I]-F[K]-C)*V<=2*(F[I]-F[J]-C)*V
即证:
F[K]>=F[J]
F[I]显然为增函数,且J<K,所以上式得证.

 

 

2.求斜率方程:一般化为左边是J，K，右边是I的形式
由DP[K]+(F[I]-F[K]-C)^2<=DP[J]+(F[I]-F[J]-C)^2
展开得:
DP[K]+F[I]^2-2*F[I]*(F[K]+C)+(F[K]+C)^2<=DP[J]+F[I]^2-2*F[I]*(F[J]+C)+(F[J]+C)^2
即:
DP[K]+(F[K]+C)^2-DP[J]-(F[J]+C)^2<=2*F[I]*(F[K]-F[J])
即
(DP[K]+(F[K]+C)^2-DP[J]-(F[J]+C)^2)/(2*(F[K]-F[J]))<=F[I]
令
G(K,J)=DP[K]+(F[K]+C)^2-DP[J]-(F[J]+C)^2
S(K,J)=2*(F[K]-F[J])
X(K,J)=G(K,J)/S(K,J)
所以斜率方程为X(K,J)<=F[I]

 

3.规定队列的维护规则
队首维护:
    假设A,B(A<B)是队首元素,若X(B,A)<=F[I],则B比A好,删除A,否则不需维护.

队尾维护:
    假设A,B,C(A<B<C)是队尾元素
a.若X(B,A)<=F[I],且X(C,B)<=F[I],则C比B好,B比A好
b.若X(B,A)<=F[I],且X(C,B)>F[I],则B比C好,B比A好,B为极大值
c.若X(B,A)>F[I],A比B好

a,c情况直接删掉B,b情况保留.b情况可改为X(B,A)<X(C,B)

 

以上内容转自：http://blog.sina.com.cn/s/blog_5f5353cc0100jx41.html

 

自我总结：

斜率优化的题前提还是很多的。。

正如风炎大神说的：

1.证明较优决策点对后续状态影响的持续性

2.求斜率方程:一般化为左边是J，K，右边是I的形式

3.规定队列的维护规则（这个规则实质上是在维护一个凸性函数，此题为下凸函数，所有转移的点都是凸性函数上的点，而这个凸性函数上的点又符合“较优决策点对后续状态影响的持续性”）

4.此题满足决策单调，转移时0(1)的，对于一些决策不单调的，可以用数据结构做到0(logn)的

 

View Code

 1 #include <iostream>

 2 #include <cstring>

 3 #include <cstdio>

 4 #include <cstdlib>

 5 #include <algorithm>

 6 

 7 #define N 55000

 8 

 9 using namespace std;

10 

11 int n;

12 long long c,sum[N],f[N],dp[N];

13 int q[N];

14 

15 inline void read()

16 {

17     scanf("%d%lld",&n,&c); c+=1;

18     for(int i=1;i<=n;i++)

19     {

20         scanf("%lld",&sum[i]);

21         sum[i]+=sum[i-1];

22         f[i]=sum[i]+(long long)i;

23     }

24 }

25 

26 inline long long S(int x,int y)

27 {

28     return 2*(f[x]-f[y]);

29 }

30 

31 inline long long G(int x,int y)

32 {

33     return dp[x]-dp[y]+(f[x]+c)*(f[x]+c)-(f[y]+c)*(f[y]+c);

34 }

35 

36 inline void go()

37 {

38     int h=1,t=2;

39     q[1]=0;

40     for(int i=1;i<=n;i++)

41     {

42         while(h<t-1&&G(q[h+1],q[h])<=S(q[h+1],q[h])*f[i]) h++;

43         dp[i]=dp[q[h]]+(f[i]-f[q[h]]-c)*(f[i]-f[q[h]]-c);

44         while(h<t-1&&G(i,q[t-1])*S(q[t-1],q[t-2])<=G(q[t-1],q[t-2])*S(i,q[t-1])) t--;

45         q[t++]=i;

46     }

47     printf("%lld\n",dp[n]);

48 }

49 

50 int main()

51 {

52     read();

53     go();

54     return 0;

55 }

 

 

PS：BZOJ 好似不能用__int64，一用就CE