码迷,mamicode.com
首页 > 其他好文 > 详细

hdu 4258 Covered Walkway

时间:2015-04-09 23:14:33      阅读:162      评论:0      收藏:0      [点我收藏+]

标签:

题目大意:

一个N个点的序列,要将他们全部覆盖,求总最少费用;费用计算: c+(x-y)2

分析:

斜率优化DP

我们假设k<j<i。如果在j的时候决策要比在k的时候决策好,那么也是就是dp[j]+M+(sum[i]-sum[j])^2<dp[k]+M+(sum[i]-sum[k])^2。(因为是最小花费嘛,所以优就是小于)

两边移项一下,得到:(dp[j]+num[j]^2-(dp[k]+num[k]^2))/(2*(num[j]-num[k]))<sum[i]。我们把dp[j]-num[j]^2看做是yj,把2*num[j]看成是xj。

那么不就是yj-yk/xj-xk<sum[i]么?   左边是不是斜率的表示? 

那么yj-yk/xj-xk<sum[i]说明了什么呢?  我们前面是不是假设j的决策比k的决策要好才得到这个表示的? 如果是的话,那么就说明g[j,k]=yj-jk/xj-xk<sum[i]代表这j的决策比k的决策要更优。

 

关键的来了:现在从左到右,还是设k<j<i,如果g[i,j]<g[j,k],那么j点便永远不可能成为最优解,可以直接将它踢出我们的最优解集。为什么呢?

我们假设g[i,j]<sum[i],那么就是说i点要比j点优,排除j点。

如果g[i,j]>=sum[i],那么j点此时是比i点要更优,但是同时g[j,k]>g[i,j]>sum[i]。这说明还有k点会比j点更优,同样排除j点。

排除多余的点,这便是一种优化!

接下来看看如何找最优解。

设k<j<i。

由于我们排除了g[i,j]<g[j,k]的情况,所以整个有效点集呈现一种上凸性质,即k j的斜率要大于j i的斜率。

这样,从左到右,斜率之间就是单调递减的了。当我们的最优解取得在j点的时候,那么k点不可能再取得比j点更优的解了,于是k点也可以排除。换句话说,j点之前的点全部不可能再比j点更优了,可以全部从解集中排除。

 

于是对于这题我们对于斜率优化做法可以总结如下:

1,用一个单调队列来维护解集。

2,假设队列中从头到尾已经有元素a b c。那么当d要入队的时候,我们维护队列的上凸性质,即如果g[d,c]<g[c,b],那么就将c点删除。直到找到g[d,x]>=g[x,y]为止,并将d点加入在该位置中。

3,求解时候,从队头开始,如果已有元素a b c,当i点要求解时,如果g[b,a]<sum[i],那么说明b点比a点更优,a点可以排除,于是a出队。最后dp[i]=getDp(q[head])。

 

原文链接  http://www.cnblogs.com/ka200812/archive/2012/08/03/2621345.html

 

对于这道题: f[i] = min(f[j] + (a[i]-a[j+1]) ^ 2+c) (N^2)

公式变形+数形结合:f[i] = min{f[j] + a[j+1]^2 - 2*a[i]*a[j+1] + a[i]^2 +C}

令x = a[j+1],y = f[j] + a[j+1]^2;

f[i] = y - 2*a[i]*x + a[i]^2 + C;

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<iostream>
 5 #include<algorithm>
 6 #include<queue>
 7 #define maxn 1000005
 8 using namespace std;
 9 typedef long long LL;
10 LL n,a[maxn];
11 LL dp[maxn];
12 int que[maxn];
13 double judge_k(int x,int y)
14 {
15     return (double)(dp[x]+a[x+1]*a[x+1]-(dp[y]+a[y+1]*a[y+1]))/((double)(2.0*(a[x+1]-a[y+1])));
16 }
17 int main()
18 {
19     int c;
20     while(scanf("%I64d %d",&n,&c)!=EOF)
21     {
22         if(n==0&&c==0)
23             break;
24         for(int i=1;i<=n;i++)
25             scanf("%I64d",&a[i]);
26         dp[0]=0;dp[1]=c;
27         que[0]=0;
28         que[1]=1;
29         int l=0,r=1;
30         for(int i=2;i<=n;i++)
31         {
32             while(l<r&&judge_k(que[l],que[l+1])<=a[i]) l++;
33             dp[i]=dp[que[l]]+c+(a[i]-a[que[l]+1])*(a[i]-a[que[l]+1]);
34             while(l<r&&judge_k(que[r],i)<=judge_k(que[r],que[r-1])) r--;
35             que[++r]=i;
36         }
37         printf("%I64d\n",dp[n]);
38     }
39     return 0;
40 }

 

hdu 4258 Covered Walkway

标签:

原文地址:http://www.cnblogs.com/tsw123/p/4412195.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!