CF765F
分块。
: 一直到第 所在块 尾端,对 块造成贡献/ 一直到 所在块 开头,对 块造成贡献。
:块 的答案。
预处理:块内元素排序,双指针求 后得到 。
查询:整块以及散块对整块的贡献查询 和 ,散块暴力归并。
时间复杂度 。
分块题做得少没见过,好像这是个很经典的分块(?),基本思路就是考虑整块与散块之间互相的贡献与内部贡献吧。
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#define gc (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 100000, stdin), p1 == p2) ? EOF : *p1 ++)
const int S = 200;
inline void cmn(int &x, const int y) {if (x > y) x = y;}
char buf[100000], *p1, *p2;
inline int read() {
char ch; int x = 0; while ((ch = gc) < 48);
do x = x * 10 + ch - 48; while ((ch = gc) >= 48); return x;
}
struct node {
int id, val;
inline bool operator < (const node a) const {return val < a.val;}
} b[100005];
int f[200000][505], mn[505][505], a[100005], lft[505], rgt[505], c[405], n, q;
bool mark[100005];
void calc(int l, int r) {
for (int i = lft[l], j = lft[r]; i <= rgt[l]; ++ i) {
while (j < rgt[r] && b[j + 1].val <= b[i].val) ++ j;
if (b[j].val > b[i].val) f[b[i].id][r] = b[j].val - b[i].val;
else {
f[b[i].id][r] = b[i].val - b[j].val;
if (j < rgt[r]) cmn(f[b[i].id][r], b[j + 1].val - b[i].val);
}
}
}
void build() {
memset(f, 0x3f, sizeof f), memset(mn, 0x3f, sizeof mn);
for (int i = 1; i <= (n + S - 1) / S; ++ i) {
lft[i] = (i - 1) * S + 1, rgt[i] = std::min(i * S, n);
std::sort(b + lft[i], b + rgt[i] + 1);
for (int j = lft[i]; j < rgt[i]; ++ j) cmn(mn[i][i], b[j + 1].val - b[j].val);
}
for (int l = 1; l <= (n + S - 1) / S; ++ l) {
for (int r = l - 1; r; -- r)
for (int i = (calc(l, r), cmn(f[lft[l]][r], f[lft[l]][r + 1]), lft[l] + 1); i <= rgt[l]; ++ i)
cmn(f[i][r], std::min(f[i - 1][r], f[i][r + 1]));
for (int r = l + 1; r <= (n + S - 1) / S; ++ r)
for (int i = (calc(l, r), cmn(f[rgt[l]][r], f[rgt[l]][r - 1]), rgt[l] - 1); i >= lft[l]; -- i)
cmn(f[i][r], std::min(f[i + 1][r], f[i][r - 1]));
}
for (int i = 1; i <= (n + S - 1) / S; ++ i)
for (int j = i + 1; j <= (n + S - 1) / S; ++ j) mn[i][j] = std::min({mn[i][j - 1], mn[j][j], f[rgt[j]][i]});
}
int getmin(int l, int r) {
int i = (l - 1) / S + 1, ans = 1e9, cnt = 0;
for (int j = lft[i]; j <= rgt[i]; ++ j) if (l <= b[j].id && b[j].id <= r) c[++ cnt] = b[j].val;
for (int i = 1; i < cnt; ++ i) cmn(ans, c[i + 1] - c[i]);
return ans;
}
int merge_sort(int l1, int r1, int l2, int r2) {
int ans = 1e9, cnt = 0, mid = 0, i = (l1 - 1) / S + 1, j = (l2 - 1) / S + 1;
for (int k = lft[i]; k <= rgt[i]; ++ k) if (l1 <= b[k].id && b[k].id <= r1) c[++ cnt] = b[k].val;
mid = cnt;
for (int k = lft[j]; k <= rgt[j]; ++ k) if (l2 <= b[k].id && b[k].id <= r2) c[++ cnt] = b[k].val;
std::inplace_merge(c + 1, c + mid + 1, c + cnt + 1);
for (int k = 1; k < cnt; ++ k) cmn(ans, c[k + 1] - c[k]);
return ans;
}
int query(int l, int r) {
int i = (l - 1) / S + 1, j = (r - 1) / S + 1;
if (i == j) return getmin(l, r);
return std::min({mn[i + 1][j - 1], f[l][j - 1], f[r][i + 1], merge_sort(l, i * S, (j - 1) * S + 1, r)});
}
int main() {
n = read();
for (int i = 1; i <= n; ++ i) a[i] = b[i].val = read(), b[i].id = i;
build(), q = read();
for (int i = 1, l, r; i <= q; ++ i) l = read(), r = read(), printf("%d\n", query(l, r));
return 0;
}
本文作者:zqs2020
本文链接:https://www.cnblogs.com/stinger/p/16637066.html
版权声明:本作品采用知识共享署名-非商业性使用-禁止演绎 2.5 中国大陆许可协议进行许可。
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步