[BZOJ1492][NOI2007]货币兑换Cash && CDQ分治+斜率优化

46 篇文章 0 订阅
4 篇文章 0 订阅

这种分治思想我也是醉了Orz

首先对于这道题 我们可以发现 如果某一天你要买进或卖出 那一定是尽可能的买或卖

那么我们就可以用一个状态f[i]表示某一天的最大收益

那么就有 f[i] = (rate[j] * f[j] * a[i] + f[j] * b[i]) / (rate[i] * a[i] + b[i]

那么朴素转移就很容易了 怎样优化呢

令 

y[j] = f[j] / (rate[i] * a[i] + b[i])

x[j] = rate[j] * f[j] / (rate[i] * a[i] + b[i]) = rate[j] * y[j];

则有 f[i] = a[i] * x[j] + b[i] * y[j];

假设决策点j 优于 决策点k 并且 j < k 当且仅当

(y[k] - y[j]) / (x[k] - x[j]) < -a[i] / b[i]

我们可以将所有 -a[i] / b[i] 从大到小排序 那么就维护一个斜率单调递减的凸包


然后就是CDQ分治了 步骤如下

1. 把区间分成两段 得到mid

2. 在原序列中小于mid的部分 进行递归

3. 当L == R 时返回答案

4. 得到答案后 维护刚刚求解部分加入后的凸包 并用他们更新mid+1到R部分的f值

5. 递归处理右子区间


话虽这么说我还是感觉很难啊QAQ

 
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<queue>
#include<cmath>
#define SF scanf
#define PF printf
using namespace std;
typedef long long LL;
const double eps = 1e-9;
const double INF = 1e20;
const int MAXN = 100000;
int n, sta[MAXN+10], top, S;
double f[MAXN+10];
struct Node {
    double x, y, a, b, rate, k;
    int id;
    bool operator < (const Node &t) const {
        return k > t.k;
    }
} A[MAXN+10], tmp[MAXN+10];
double slope(int i, int j) {
    if(!j) return -INF;
    if(fabs(A[i].x - A[j].x) < eps) return INF;
    return (A[j].y - A[i].y) / (A[j].x - A[i].x);
}
void CDQ(int L, int R) {
    if(L == R) {
        f[L] = max(f[L], f[L-1]);
        A[L].y = f[L] / (A[L].a * A[L].rate + A[L].b);
        A[L].x = A[L].rate * A[L].y;
        return ;
    }
     
    int mid = (L + R) >> 1, l1 = L, l2 = mid+1;
    for(int i = L; i <= R; i++) {
        if(A[i].id <= mid) tmp[l1++] = A[i];
        else tmp[l2++] = A[i];
    }
    for(int i = L; i <= R; i++) A[i] = tmp[i];
    CDQ(L, mid);
    top = 0;
    for(int i = L; i <= mid; i++) {
        while(top > 1 && slope(sta[top-1], sta[top]) < slope(sta[top-1], i) + eps) top--;
        sta[++top] = i;
    }
    sta[++top] = 0;
    int j = 1;
    for(int i = mid+1; i <= R; i++) {
        while(j < top && slope(sta[j], sta[j+1]) + eps > A[i].k) j++;
        f[A[i].id] = max(f[A[i].id], A[sta[j]].x * A[i].a + A[sta[j]].y * A[i].b);
    }
    CDQ(mid+1, R);
    l1 = L; l2 = mid+1;
    for(int i = L; i <= R; i++) {
        if(l1 <= mid && (l2 > R || A[l1].x < A[l2].x || (fabs(A[l1].x - A[l2].x) < eps && A[l1].y < A[l2].y)))
            tmp[i] = A[l1++];
        else tmp[i] = A[l2++];
    }
    for(int i = L; i <= R; i++) A[i] = tmp[i];
}
int main() {
    SF("%d%d", &n, &S);
    f[0] = S;
    for(int i = 1; i <= n; i++) {
        SF("%lf%lf%lf", &A[i].a, &A[i].b, &A[i].rate);
        A[i].k = -A[i].a / A[i].b;
        A[i].id = i;
    }
    sort(A+1, A+1+n);
    CDQ(1, n);
    PF("%.3f", f[n]);
    return 0;
}


### 回答1: bzoj[1597][usaco2008 mar]土地购买 斜率优化 这道题是一道经典的斜率优化题目,需要用到单调队列的思想。 首先,我们可以将题目中的式子进行变形,得到: f[i] = f[j] + (sum[i] - sum[j] - m) ^ 2 + k 其中,sum[i] 表示前缀和,m 和 k 都是常数。 我们可以将式子中的 sum[i] 和 k 看作常数,那么我们需要优化的就是 (sum[i] - sum[j] - m) ^ 2 这一项。 我们可以将其展开,得到: (sum[i] - sum[j] - m) ^ 2 = sum[i] ^ 2 - 2 * sum[i] * (sum[j] + m) + (sum[j] + m) ^ 2 我们可以将其看作一个二次函数,其中 a = 1,b = -2 * (sum[j] + m),c = (sum[j] + m) ^ 2。 我们可以发现,当 j < k 时,如果 f[j] + a * sum[j] + b * sum[j] <= f[k] + a * sum[k] + b * sum[k],那么 j 就不可能是最优决策点,因为 k 比 j 更优。 因此,我们可以用单调队列来维护决策点。具体来说,我们可以维护一个单调递增的队列 q,其中 q[i] 表示第 i 个决策点的下标。每次加入一个新的决策点 i 时,我们可以将队列尾部的决策点 j 弹出,直到队列为空或者 f[j] + a * sum[j] + b * sum[j] <= f[i] + a * sum[i] + b * sum[i]。然后,我们将 i 加入队列尾部。 最后,队列头部的决策点就是最优决策点。我们可以用类似于双指针的方法来维护队列头部的决策点是否在当前区间内,如果不在,就弹出队列头部。 时间复杂度为 O(n)。 ### 回答2: 这道题目属于斜率优化的经典题目,难度较高,需要掌握一定的数学知识。 首先,我们可以将题目中的“最大利润”转化为“最小成本”,这样问题就变成了找到一个方案,使得购买土地的成本最小。 接着,我们考虑如何用斜率优化来解决这个问题。我们可以定义一个函数f(i),表示前i块土地的最小成本。 显然,f(1)=0,因为不需要购买任何土地。 对于f(i),它可以由f(j)+b(i)×a(j+1)得到,其中j<i,a(j+1)表示第j+1块土地的面积,b(i)表示第i块土地的价格。这个式子的含义是,我们现在要购买第i块土地,那么前面的土地(即前j块)就都要买,所以f(j)表示前j块土地的最小成本,b(i)×a(j+1)表示购买第i块土地的成本。 那么,我们可以得到递推公式: f(i)=min{f(j)+b(i)×a(j+1)},其中j<i。 这个公式看起来很简单,但是要注意的是,当b(i)×a(j+1)的斜率相同时,我们需要取其中面积较小的土地,因为它的价格更低。因此,我们需要对斜率进行排序,并在递推中用单调队列维护斜率相等的情况下面积最小的土地。 最终,f(n)就是题目所求的最小成本。 总之,这道题目需要深入理解斜率优化算法的原理和实现方式,并且需要注意细节处理,如果能够顺利地解决这个问题,那么对于斜率优化算法的掌握程度就有了很大的提升。 ### 回答3: 土地购买问题可以采用斜率优化算法来解决。这个问题可以转化为一个单调队列的问题。 首先,我们需要对土地价格按照边长从小到大排序。然后,对于每块土地,我们需要求出它的贡献。设 $f_i$ 表示前 $i$ 块土地连续的最小代价。 设当前处理到第 $i$ 块土地,已经求出了前 $j$ 块土地的最小代价 $f_j$。那么我们可以得到下面这个式子: $$ f_i=\min\limits_{j=1}^{i-1}\{f_j+(S_i-S_j)^2+P\} $$ 式子中,$S_i$ 表示前 $i$ 块土地的边长和,$P$ 表示额外购买土地的代价。首先,不考虑额外购买土地,我们可以使用动态规划来求出 $f_i$。但是,考虑到额外购买土地的代价 $P$ 是一个固定值,我们可以考虑将它与某一块土地的代价合并起来,这样就可以使用斜率优化技术来优化动态规划算法。 我们定义一个决策点 $j$,表示我们当前要处理第 $i$ 块土地时,已经处理过 $j$ 块土地,并将第 $j+1$ 块土地到第 $i$ 块土地购买,所需的最小代价。我们假设 $S_i>S_j$,则可以得到下面这个式子: $$ f_i=\min\limits_{j=1}^{i-1}\{f_j+(S_i-S_j)^2+P\} $$ 将它整理成斜率截距式可以得到: $$ y=kx+b $$ 其中 $k=(S_j)^2-2S_iS_j$,$b=f_j+(S_i)^2+P-S_j^2$,$x=S_j$,$y=f_j+(S_j-S_i)^2-S_j^2$。我们发现 $k$ 是一个单调递减的函数,因此我们可以使用一个单调队列来维护所有可能成为决策点的点。对于每个点,我们计算函数 $y$ 的值并将它们加入队列,然后取队头元素的值作为 $f_i$。 综上所述,我们可以使用斜率优化技术来解决土地购买问题,时间复杂度为 $O(n)$。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值