原文出处:http://txhwind.blog.163.com/blog/static/2035241792012112285146817/(版权为原作者所有,本博文无营利目的,不构成侵权)
题目大意:(POJ 3017 也是这个)
#include<stdio.h>
#include<set> using std::multiset; const int MAXN=+; int h[MAXN],s[MAXN],w[MAXN];
long long f[MAXN];
int main()
{
freopen("bookshelf.in","r",stdin);
freopen("bookshelf.out","w",stdout);
int n,l,i;
int sum=,head=,tail=;
multiset<long long> val;
scanf("%d%d",&n,&l);
for(i=;i<=n;++i)
{
scanf("%d%d",h+i,w+i); //s记录每段开头
for(s[++tail]=i-;head<tail && h[i]>=h[s[tail]];--tail)
val.erase(val.find(f[s[tail-]]+h[s[tail]]));
val.insert(f[s[tail]]+h[i]);
s[tail+]=i;
for(sum+=w[i];sum>l;sum-=w[s[head]])
{
val.erase(val.find(f[s[head]]+h[s[head+]]));
if(++s[head]==s[head+])
++head;
else val.insert(f[s[head]]+h[s[head+]]);
}
f[i]=*val.begin();
}
printf("%lld\n",f[n]);
return ;
}
以下是我的同学的解题报告:不涉及任何侵权问题。。同学的权,随便侵。。。。
这道题的原型是POJ的3017,只是W[I]换成数列中的权值而已;
对题目进行分析,不难得到动规方程f[i]=min(f[j]+max(h[j+1],h[j+2],h[j+3]....h[i]))其中w[j+1]+w[j+2]+...w[j]<=l;
但是这只是O(n^2)的做法,必须想出优化办法;
首先我们会发现一个i状态会与许多j的子状态相关,这种情况大多都是与单调队列的优化相关;但是这道题并不这么明显;
首先,我们向求max的h[j]值方向思考,不难维护出一个单调递减,队头到i的w[i]之和小于l的最大值单调队列;但是这样就将f的值给忽略掉了;
所以我们需要想出f数组的特点,不难看出这是单调递增的,
所以我们猜想:每一个单调队列中的决策点相对于后面的一段都有可能是最优决策点,但是如何证明呢?
如果这个不是最优决策点,那么他后面的一段区间的最大值都与该点相等,而f非递增,所以决策都没有在单调队列中的决策好;
既然我们发现f与单调队列也有关系,那么此题单调队列的方向就已经确定;
维护上述的一个单调队列,在计算f[i]时在单调队列中的所有点都有可能成为决策点,对于这种情况可以用平衡树进行维护,但代码还有许多细节处理
#include<cstdio>
#include<set>
#include<iostream>
using namespace std;
const int mm=;
int a[mm],q[mm],w[mm];
long long f[mm],sum,tmp,m;
multiset<long long> sbt;
int i,l,r,p,n;
int main()
{
scanf("%d%lld",&n,&m);
sbt.clear();
sum=l=,f[n]=r=-;
for(p=i=;i<=n;++i)
{
scanf("%d%d",&a[i],&w[i]);
sum+=w[i];
while(sum>m)sum-=w[p++];
while(l<=r&&a[i]>=a[q[r]])
{
if(l<r)sbt.erase(f[q[r-]]+a[q[r]]);
--r;
}
q[++r]=i;
if(l<r)sbt.insert(f[q[r-]]+a[q[r]]); //手动模拟下就会发现对于边界外决策点,其待定决策值为f[q[i]-1]+a[q[i]];
while(q[l]<p)
{
if(l<r)sbt.erase(f[q[l]]+a[q[l+]]);
++l;
}
f[i]=f[p-]+a[q[l]]; //在i为当前最大值时用p数组更新 tmp=*sbt.begin();
if(f[i]>tmp)f[i]=tmp;
}
cout<<f[n];
return ;
}
摘自另一地方。http://www.cnblogs.com/kedebug/archive/2013/03/01/2939428.html
题意:
长度为 n 的数列,要求把这个数列划分为任意块,每块的元素和小于 m,使得所有块的最大值的和最小。
思路:
1. 很明显的一个转移方程是: dp[i] = max(dp[j] + max(a[j+1], a[j+2], ..., a[i])); 其中满足 sum[i] - sum[j] <= M 的 j 都算,这是一个 O(N * N) 的解法。
2. 如果不去优化枚举 j 的这个过程,就会直接 TLE,如何优化?联系到求最大值,可以想到用最大值的单调队列,队列里面存放的是以 i 为结尾,窗口里面的元素和
恰好不大于 M 的 j 为左边界,即恰好有 sum[i] - sum[j-1] <= M.
3. 如果是这样的话,还不能够达到目标,因为 1 中最优值点的 j 的选择可能是窗口中的任意一个值,这时候遇到一个困难:确定 j 的位置,使情况最优。
4. 考虑到 dp[i] 是一个单调增的数组,而窗口中的元素是关于最大值递减的,先让 M 足够大,从简单的角度来考虑:对于 a[1] = 8, a[2] = 5, a[3] = 7,对于 dp[3]
则有 dp[3] = min(dp[0] + 8, dp[1] + 7, dp[2] + 7),显然 dp[1] <= dp[2],则进一步优化有 dp[3] = min(dp[0] + 8, dp[1] + 7);
5. 从 4 中我们可以看出点端倪,单调队列中所指向的元素也恰好是 8, 7,所以得出的结论有:对于 dp[i] 只要选择单调队列中所指向的位置的 j 即可保证结果是最优的。
6. 如何维护单调队列中的 dp[j] + max(a[j+1], a[j+2], ..., a[i]),则考虑使用一个平衡树,这样可以把时间复杂度降到最低。
#include <iostream>
#include <set>
#include <algorithm>
using namespace std; #define LL long long int const int MAXN = 100010; int num[MAXN], deq[MAXN];
LL dp[MAXN];
multiset<LL> optset; int main()
{
LL N, M; scanf("%lld %lld", &N, &M);
for (int i = 1; i <= N; ++i)
scanf("%d", &num[i]); LL sum = 0;
int s = 0, e = -1, p = 1;
bool flag = false; optset.clear();
for (int i = 1; i <= N; ++i)
{
if (num[i] > M)
{
flag = true;
break;
} sum += num[i];
while (sum > M)
sum -= num[p++]; while (s <= e && num[i] >= num[deq[e]])
{
if (s < e)
optset.erase(dp[deq[e-1]] + num[deq[e]]);
--e;
} deq[++e] = i;
if (s < e)
optset.insert(dp[deq[e-1]] + num[deq[e]]); while (deq[s] < p)
{
if (s < e)
optset.erase(dp[deq[s]] + num[deq[s+1]]);
++s;
} dp[i] = dp[p-1] + num[deq[s]];
if (s < e && dp[i] > *optset.begin())
dp[i] = *optset.begin();
}
if (flag)
dp[N] = -1;
printf("%lld\n", dp[N]);
return 0;
}