P1502 視窗的星星

纯粹的發表於2024-07-20

原題連結

題解

1.首先要解決的問題是,怎樣的視窗能包裹住儘可能多的星星?

這裡有一個很巧妙的思維點,那就是我們構造一個以星星為左下角的矩形,矩形重疊的部分可以構造一個右上角在其中的視窗,這樣這個視窗就能覆蓋矩陣重合的星星

2.這樣,問題就轉換成了求最大並矩陣權值和(表示視窗的右上角在某個矩陣內),再轉換一下,就變成了某條豎線上,區間重合最多的值
我們可以離散化+掃描線的方法來做,即構建一個y軸的線段樹,然後把矩形分割開來(有點類似於差分),離散化是為了方便用線段樹操作

code

#include <bits/stdc++.h>
#define ll long long
using namespace std;

vector<int> tree(800005, 0);        // 線段樹陣列
vector<int> lazytag(800005, 0);     // 延遲標記陣列

struct node {
    int x, y1, y2, val;             // 節點結構體
    bool operator<(const node &b) const {
        if (x != b.x) return x < b.x;
        else return val < b.val;
    }
};

void build(int node, int l, int r) {
    tree[node] = 0;
    lazytag[node] = 0;
    if (l != r) {
        int mid = (l + r) / 2;
        build(node * 2, l, mid);
        build(node * 2 + 1, mid + 1, r);
    }
}

void pushdown(int node, int l, int r) {
    if (!lazytag[node]) return;
    tree[node] += lazytag[node];
    if (l != r) {
        lazytag[node * 2] += lazytag[node];
        lazytag[node * 2 + 1] += lazytag[node];
    }
    lazytag[node] = 0;
}

void update(int node, int l, int r, int x, int y, int val) {
    pushdown(node, l, r);
    if (l > y || r < x) return;

    if (l >= x && r <= y) {
        lazytag[node] += val;
        pushdown(node, l, r);
        return;
    }

    int mid = (l + r) / 2;
    update(node * 2, l, mid, x, y, val);
    update(node * 2 + 1, mid + 1, r, x, y, val);
    tree[node] = max(tree[node * 2], tree[node * 2 + 1]);
}

void solve() {
    int n, w, h;
    cin >> n >> w >> h;

    vector<node> seg;
    set<int> haxi;

    for (int i = 1; i <= n; i++) {
        int x, y, v;
        cin >> x >> y >> v;

        int endx = x + w;
        int y1 = y, y2 = y + h - 1;

        haxi.insert(y1);
        haxi.insert(y2);

        seg.push_back({x, y1, y2, v});
        seg.push_back({endx, y1, y2, -v});
    }

    vector<int> haxi_vec(haxi.begin(), haxi.end());
    int len = haxi_vec.size();

    build(1, 1, len);

    sort(seg.begin(), seg.end());

    int ans = 0;
    for (auto it : seg) {
        auto [x, y1, y2, val] = it;

        int pos1 = lower_bound(haxi_vec.begin(), haxi_vec.end(), y1) - haxi_vec.begin() + 1;
        int pos2 = lower_bound(haxi_vec.begin(), haxi_vec.end(), y2) - haxi_vec.begin() + 1;

        update(1, 1, len, pos1, pos2, val);
        ans = max(ans, tree[1]);
    }
    cout << ans << '\n';
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t = 1;
    cin >> t;
    while (t--) solve();
    return 0;
}

相關文章