[HNOI2008]玩具装箱

OJ题号:
  BZOJ1010

思路:
  斜率优化动态规划。
  由题意得状态转移方程为$f_i=\displaystyle{\min_{j=0}^{i-1}}\{f_j+\left(i-j-1+\displaystyle{\sum_{k=j+1}^i}c_k-L\right)^2\}$。
  用$a_i$表示$c_i$的前缀和,则原式为$f_i=\displaystyle{\min_{j=0}^{i-1}}\{f_j+\left(i-j-1+a_i-a_j-L\right)^2\}$。
  考虑现在有两个状态$j$和$k$都可以转移到$i$。
  假设$j$比$k$更优,则有:$f_j+\left(i-j-1+a_i-a_j-L\right)^2<f_k+\left(i-k-1+=a_i-a_k-L\right)^2$。
  将与$i$有关的项提取出来,设$x=i-1+a_i-L$。
  则原式变为$f_j+\left(x-j-a_j\right)^2<f_k+\left(x-k-a_k\right)^2$。
  化简得$f_j+\left(j+a_j\right)^2-f_k-\left(k+a_k\right)^2<2x(j+a_j-k-a_k)$。
  即$\frac{f_j+\left(j+a_j\right)^2-f_k-\left(k+a_k\right)^2}{(j+a_j-k-a_k)}<2x$。
  对于状态$j<k<l$,若要使$k$为一个有用的状态,则有$\frac{f_k+(k+a_k)^2-f_j-(j+a_j)^2}{2(k+a_k-j-a_j)}<x\leq\frac{f_l+(l+a_l)^2-f_k-(k+a_k)^2}{2(l+a_l-k-a_k)}$。
  然后我们可以维护一个单调队列,使队列中的相邻元素的斜率单调递增。
  每当插入一个元素时,我们比较队列前端两个元素的斜率是否小于$x$,如果是,则将第一个元素弹出队列。
  这时候队列前端的元素一定是最优的一个状态。
  然后尝试将这个元素加入队列,为了保证队列中相邻元素之间的斜率单调递增,每次比较队列后端两个元素的斜率$x1$和队列最末端元素与当前元素$i$的斜率$x2$。
  如果$x1>x2$,即新加入元素后不满足单调性,则将队列末端元素弹出。
  由于每个元素最多只会进队一次,最后的时间复杂度是$O(n)$的。

 1 #include<cstdio>
 2 #include<cctype>
 3 inline int getint() {
 4     register char ch;
 5     while(!isdigit(ch=getchar()));
 6     register int x=ch^‘0‘;
 7     while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^‘0‘);
 8     return x;
 9 }
10 const int N=50001;
11 int n,l,a[N]={0},q[N]={0},h=0,t=0;
12 long long f[N]={0};
13 inline long long sqr(const long long x) {
14     return x*x;
15 }
16 inline double slope(const int &j,const int &k) {
17     return double(f[j]+sqr(j+a[j])-f[k]-sqr(k+a[k]))/double(j+a[j]-k-a[k]);
18 }
19 inline bool check(const int &i,const int &j,const int &k) {
20     return slope(k,j)>2*(i+a[i]-l-1);
21 }
22 inline bool check2(const int &i,const int &j,const int &k) {
23     return slope(j,i)<slope(k,j);
24 }
25 int main() {
26     n=getint(),l=getint();
27     for(register int i=1;i<=n;i++) {
28         a[i]=a[i-1]+getint();
29         while(h<t&&!check(i,q[h],q[h+1])) h++;
30         const int &j=q[h];
31         f[i]=f[j]+sqr(i-j-1+a[i]-a[j]-l);
32         while(h<t&&!check2(q[t-1],q[t],i)) t--;
33         q[++t]=i;
34     }
35     printf("%lld\n",f[n]);
36     return 0;
37 }
时间: 2024-10-20 09:18:37

[HNOI2008]玩具装箱的相关文章

BZOJ 1010: [HNOI2008]玩具装箱toy [DP 斜率优化]

1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 9812  Solved: 3978[Submit][Status][Discuss] Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P

BZOJ 1010 [HNOI2008]玩具装箱toy

1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 7184  Solved: 2724[Submit][Status][Discuss] Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P

BZOJ 1010: [HNOI2008]玩具装箱toy(DP+斜率优化)

[HNOI2008]玩具装箱toy Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P教授要求在一个一维容器中的玩具编号是连续的.同时如果一个一维容器中有多个玩具,那么两件玩具之间要加入一个单位长度的填充物,形式地说如果将第i件玩具到第j个玩具放到一个容器中,那么容器的长度将为

bzoj 1010 [HNOI2008]玩具装箱toy(DP的斜率优化)

1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 7874  Solved: 3047[Submit][Status][Discuss] Description P 教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维 容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理

HNOI2008 玩具装箱toy (BZOJ1010,斜率dp)

传送门 1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MB Description P 教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维 容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P教授要求在一个一维容器中的玩具编号是连续的. 同时如果一个一维容器中有多个玩具,那么两件玩

BZOJ 1010: [HNOI2008]玩具装箱toy 斜率优化DP

1010: [HNOI2008]玩具装箱toy Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P教授要求在一个一维容器中的玩具编号是连续的.同时如果一个一维容器中有多个玩具,那么两件玩具之间要加入一个单位长度的填充物,形式地说如果将第i件玩具到第j个玩具放到一个容器中,那么容器

bzoj 1010: [HNOI2008]玩具装箱toy 2011-12-27

1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 2254  Solved: 744[Submit][Status][Discuss] DescriptionP 教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维 容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P

HNOI2008玩具装箱 (斜率优化)

总算A了,心情好激动-- 如果会了一类斜率优化,基本上这类题就成了套模版了-- 只是k函数不同 1 var n,l,x,tail,head,m:int64; 2 i,j:longint; 3 dp,q,s:array[0..100000] of int64; 4 function k(x,y:longint):double; 5 begin 6 k:=1.0*((dp[x]+s[x]*s[x]-dp[y]-s[y]*s[y])/(s[x]-s[y])); 7 end; 8 procedure m

斜率优化专题2——bzoj 1010 [HNOI2008]玩具装箱toy 题解

[原题] 1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MB Submit: 5434  Solved: 1969 [Submit][Status] Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具,第i件玩具经过压缩后变成一维长度为Ci.为了方便整理,P教授

BZOJ 1010 HNOI2008 玩具装箱

1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec Memory Limit: 162 MB Submit: 11703 Solved: 4960 [Submit][Status][Discuss] Description P教授要去看奥运,但是他舍不下他的玩具,于是他决定把所有的玩具运到北京.他使用自己的压缩器进行压 缩,其可以将任意物品变成一堆,再放到一种特殊的一维容器中.P教授有编号为1...N的N件玩具,第i件玩具经过 压缩后变成一维长度为Ci.为了方便整