标签：p1 val int mn lft CF765F 505

分块。

\(f[i][j]\)：\(i\) 一直到第 \(i\) 所在块 \(x\) 尾端，对 \(x+1\sim j\) 块造成贡献/\(i\) 一直到 \(i\) 所在块 \(x\) 开头，对 \(j\sim x-1\) 块造成贡献。

\(mn[i][j]\)：块 \(i\sim j\) 的答案。

预处理：块内元素排序，双指针求 \(f\) 后得到 \(mn\)。

查询：整块以及散块对整块的贡献查询 \(mn\) 和 \(f\)，散块暴力归并。

时间复杂度 \(O(n\sqrt{n}))\)。

分块题做得少，好像这是个很经典的分块（？）

#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;
}

标签：p1,val,int,mn,lft,CF765F,505
来源： https://www.cnblogs.com/stinger/p/16637066.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。