P4198 楼房重建 题解

一道线段树题目，思路很巧妙。

首先先转化一下题意，发现如果后面的楼房能够被前面的挡住，一定是后面楼房的斜率比前面楼房小，斜率就是 /(/dfrac{H_i}{i}/)，于是这道题变成了单点修改，全局查询哪几个点斜率是前缀最大值中的最大值。

维护两个值 /(ans,Maxn/) 分别表示只存在这个区间时的答案和这个区间的斜率最大值，发现不需要建树，单点修改也好写，查询更好写，但是这个 Update(Pushup) 不好写，因为前面的区间是会干扰到后面的区间的。

考虑我们现在要 Update 一个节点 /(p/)，右儿子的左儿子 /(ls/)，右儿子的右儿子 /(rs/)，根据当前点斜率最大值 /(Maxn(p)/) 和左儿子斜率最大值 /(Maxn(ls)/) 分类讨论一下：

• 如果 /(Maxn(p) /geq Maxn(ls)/)，说明左儿子的区间对答案的贡献是 0（全部被挡住了），直接递归查右儿子的贡献。
• 如果 /(Maxn(p) < Maxn(ls)/)，说明此时左儿子实际上是有贡献的（有的没有被挡住），那么考虑往左儿子递归查贡献。
  那么右儿子呢？注意到此情况下，左儿子对右儿子的约束一定是比 /(Maxn(p)/) 对右儿子的约束要强的，换言之如果左儿子挡不住则 /(p/) 一定也挡不住，所以右儿子是都有贡献的，然后这里有个坑，右儿子贡献是 /(ans(p 的右儿子)-ans(ls)/) 而不是 /(ans(rs)/)，因为右儿子的部分点是会被左儿子挡住的所以要这么写。

讨论完毕，上述两种情况在往下递归时传下的斜率最大值依然是 /(Maxn(p)/)，上面这段的左儿子与右儿子指的是 /(p/) 的右儿子的左右儿子。

这段话会有点绕，我觉得还是要看一下代码里面的 Update 函数会更好理解。

Update 解决了这道题也就解决了。

GitHub：CodeBase-of-Plozia

Code：

/*
========= Plozia =========
    Author:Plozia
    Problem:P4198 楼房重建
    Date:2022/1/2
========= Plozia =========
*/

#include <bits/stdc++.h>

typedef long long LL;
const int MAXN = 1e5 + 5;
int n, m;
struct node
{
    int l, r, ans;
    double Maxn;
    #define l(p) tree[p].l
    #define r(p) tree[p].r
    #define ans(p) tree[p].ans
    #define Maxn(p) tree[p].Maxn
}tree[MAXN << 2];

int Read()
{
    int sum = 0, fh = 1; char ch = getchar();
    for (; ch < '0' || ch > '9'; ch = getchar()) fh -= (ch == '-') << 1;
    for (; ch >= '0' && ch <= '9'; ch = getchar()) sum = sum * 10 + (ch ^ 48);
    return sum * fh;
}
double Max(double fir, double sec) { return (fir > sec) ? fir : sec; }
double Min(double fir, double sec) { return (fir < sec) ? fir : sec; }

void Build(int p, int l, int r)
{
    l(p) = l, r(p) = r; if (l == r) return ;
    int mid = (l(p) + r(p)) >> 1; Build(p << 1, l, mid); Build(p << 1 | 1, mid + 1, r);
}

int Ask(int p, double k)
{
    if (Maxn(p) <= k) return 0;
    if (l(p) == r(p)) return (Maxn(p) > k);
    if (Maxn(p << 1) <= k) return Ask(p << 1 | 1, k);
    return ans(p) - ans(p << 1) + Ask(p << 1, k);
}

void Change(int p, int x, int k)
{
    if (l(p) == r(p) && l(p) == x) { ans(p) = 1; Maxn(p) = (double) k / x; return ; }
    int mid = (l(p) + r(p)) >> 1;
    if (x <= mid) Change(p << 1, x, k); else Change(p << 1 | 1, x, k);
    Maxn(p) = Max(Maxn(p << 1), Maxn(p << 1 | 1));
    ans(p) = ans(p << 1) + Ask(p << 1 | 1, Maxn(p << 1));
}

int main()
{
    n = Read(), m = Read(); Build(1, 1, n);
    for (int i = 1; i <= m; ++i)
    {
        int x = Read(), y = Read();
        Change(1, x, y);
        printf("%d/n", ans(1));
    }
    return 0;
}