P2202 题解
前言
提供一个平衡树做法,虽然和 std::set
一个道理就是了。(那你为啥不写 set!!!!)
前置知识
如何判断两个点对应的正方形相交?正方形的边长是
类似地,重叠部分的面积就是
思路
很容易想到对
接着加入
很容易想到平衡树。(std::set
:????)以
显然
代码
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
const int N = 2e5 + 5;
int n, k;
struct Poi {int x, y;} a[N];
bool cmp(Poi p, Poi q) {return p.x < q.x;}
bool chk(Poi p, Poi q) {return abs(p.x - q.x) < k && abs(p.y - q.y) < k;}
long long ans;
void update(Poi p, Poi q)
{
if (!chk(p, q)) return;
if (ans) puts("-1"), exit(0); //有两个了,寄
ans = 1ll * (k - abs(p.x - q.x)) * (k - abs(p.y - q.y));
}
bool operator < (Poi p, Poi q) {if (p.y != q.y) return p.y < q.y; return p.x < q.x;}
bool operator == (Poi p, Poi q) {return p.x == q.x && p.y == q.y;}
bool operator <= (Poi p, Poi q) {return p < q || p == q;}
struct ScapegoatTree { //替罪羊树板子
#define ls tr[pos].l
#define rs tr[pos].r
#define tmp (tr[ls].siz + (!tr[pos].del))
const double T = 0.7; const int inf = 2147483647;
struct Node {int l, r, tot, siz; Poi val; bool del; void clear() {l = r = del = 0, siz = tot = 1;}} tr[N];
int root, stk[N], top, order[N], cnt;
ScapegoatTree() {for (int i = N - 1; i; i--) stk[++top] = i;}
void pushup(int pos) {tr[pos].siz = tr[ls].siz + tr[rs].siz + 1, tr[pos].tot = tr[ls].tot + tr[rs].tot + 1;}
void inorder(int pos) {if (!pos) return; inorder(ls); tr[pos].del ? stk[++top] = pos : order[++cnt] = pos; inorder(rs);}
void build(int l, int r, int &pos) {if (l > r) return pos = 0, void(); int mid = (l + r) >> 1; pos = order[mid]; if (l == r) return tr[pos].clear(); build(l, mid - 1, ls), build(mid + 1, r, rs); pushup(pos);}
void rebuild(int &pos) {cnt = 0, inorder(pos); build(1, cnt, pos);}
int getrnk(int pos, Poi val) {if (!pos) return 0; if (val <= tr[pos].val) return getrnk(ls, val); else return tmp + getrnk(rs, val);} int getrnk(Poi val) {return getrnk(root, val) + 1;}
Poi getval(int pos, int k) {if (!pos) return (Poi){-inf, -inf}; if (!tr[pos].del && tmp == k) return tr[pos].val; if (k <= tmp) return getval(ls, k); else return getval(rs, k - tmp);} Poi getval(int k) {return getval(root, k);}
void insert(int &pos, Poi val) {if (!pos) return pos = stk[top--], tr[pos].val = val, tr[pos].clear(); tr[pos].siz++, tr[pos].tot++; if (val <= tr[pos].val) insert(ls, val); else insert(rs, val); if (T * tr[pos].siz <= max(tr[ls].siz, tr[rs].siz)) rebuild(pos);} void insert(Poi val) {insert(root, val);}
void delett(int &pos, int k) {tr[pos].siz--; if (!tr[pos].del && tmp == k) return tr[pos].del = true, void(); if (k <= tmp) delett(ls, k); else delett(rs, k - tmp);} void delett(Poi val) {delett(root, getrnk(val)); if (T * tr[root].tot >= tr[root].siz) rebuild(root);}
Poi getpre(Poi val) {Poi t = getval(getrnk(val) - 1); return (t.x == -inf ? (Poi){-inf, -inf} : t);} Poi getsuf(Poi val) {Poi t = getval(getrnk((Poi){val.x + 1, val.y})); return (t.x == -inf ? (Poi){inf, inf} : t);}
#undef ls
#undef rs
} T;
int main()
{
scanf("%d%d", &n, &k);
for (int i = 1; i <= n; i++) scanf("%d%d", &a[i].x, &a[i].y);
sort(a + 1, a + n + 1, cmp);
for (int l = 1, r = 1; r <= n; r++)
{
while (l < r && a[r].x - a[l].x >= k) T.delett(a[l++]);
Poi pre = T.getpre(a[r]), suf = T.getsuf(a[r]);
if (pre.x != -T.inf) update(a[r], pre);
if (suf.y != T.inf) update(a[r], suf);
T.insert(a[r]);
}
cout << ans << '\n';
return 0;
}
希望能帮助到大家!
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 在鹅厂做java开发是什么体验
· 百万级群聊的设计实践
· WPF到Web的无缝过渡:英雄联盟客户端的OpenSilver迁移实战
· 永远不要相信用户的输入:从 SQL 注入攻防看输入验证的重要性
· 全网最简单!3分钟用满血DeepSeek R1开发一款AI智能客服,零代码轻松接入微信、公众号、小程