NOIP模拟（11.06）T2 序列操作

序列操作

题目背景：

11.06 NOIP模拟T2

分析：线段树 or 贪心

 

垃圾数据范围······请问你10⁵的数据范围有意思吗······过得了10⁵的有谁过不了10⁶考场上全部在想O(n)结果连treap的nlogn都能水过去（手动再见······

讲题，考虑每一次你一定会从当前最大的c[i]个数上面减一，这样可以保证进行尽量多的操作，考虑如何维护这个过程，如果一开始我们的序列是排好序了的，那么我们将最后的c[i]个数全部减一，如果第c[i]个数和第c[i] + 1个数是不相同的，那么这样的减法，不会影响原来序列的单调不降，而如果第c[i]个数和第c[i] + 1个数是相同的，假设这个数为x，那么如果前c[i]个数中有p个x，而一共有q个，我们可以选择将对前p个x的操作改成对q个x后面p个的操作，那么这样的减法依然不会对原序列的单调存在影响，考虑如何维护这种操作，我们可以利用线段树，我们先找到第c[i]个数的大小，记做x，我们可以在线段树上二分找到，大于x个最小的数字在哪里，显然大于x的数全部都会被减一，然后如果当前大于x的数一共有了a个，那么我们还需要更改c[i] - a个x，显然我们可以在线段树上二分找到第一个x的位置，然后改掉钱c[i] - a个x就好了。

听说题解可以贪心，然而这种正确性都没有的算法，谁特喵知道是啥玩意儿

Source:

/*
	created by scarlyw
*/
#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <string>
#include <cstring>
#include <cctype>
#include <vector>
#include <queue>
#include <set>
#include <ctime> 

const int MAXN = 1000000 + 10;
int n, m, x, pos, cur;
int h[MAXN], c[MAXN], bucket[MAXN];

struct node {
	int min, tag, max;
} tree[MAXN << 2];

inline void modify(int k, int x) {
	tree[k].min -= x, tree[k].max -= x, tree[k].tag += x;
}

inline void push_down(int k) {
	if (tree[k].tag != 0) {
		modify(k << 1, tree[k].tag), modify(k << 1 | 1, tree[k].tag);
		tree[k].tag = 0;
	}
}

inline void update(int k) {
	tree[k].min = std::min(tree[k << 1].min, tree[k << 1 | 1].min);
	tree[k].max = std::max(tree[k << 1].max, tree[k << 1 | 1].max);
}

inline void get(int k, int l, int r) {
	if (l == r) {
		std::cout << tree[k].max << '\n';
		return ;
	}
	push_down(k);
	int mid = l + r >> 1;
	get(k << 1, l, mid), get(k << 1 | 1, mid + 1, r);
	update(k);
}

inline void build(int k, int l, int r) {
	if (l == r) return (void)(tree[k].min = tree[k].max = h[l]);
	int mid = l + r >> 1;
	build(k << 1, l, mid), build(k << 1 | 1, mid + 1, r);
	update(k);
}

inline int query(int k, int l, int r, int pos) {
	if (l == r) return tree[k].min;
	int mid = l + r >> 1;
	push_down(k);
	if (pos <= mid) return query(k << 1, l, mid, pos);
	else return query(k << 1 | 1, mid + 1, r, pos); 
}

inline int find(int k, int l, int r, int x) {
	if (l == r) return l;
	int mid = l + r >> 1;
	push_down(k);
	if (tree[k << 1].max >= x) return find(k << 1, l, mid, x);
	else return find(k << 1 | 1, mid + 1, r, x);
}

inline void modify(int k, int l, int r, int ql, int qr) {
	if (ql <= l && r <= qr) return modify(k, 1);
	int mid = l + r >> 1;
	push_down(k);
	if (ql <= mid) modify(k << 1, l, mid, ql, qr);
	if (qr > mid) modify(k << 1 | 1, mid + 1, r, ql, qr);
	update(k);
}

inline void solve() {
	scanf("%d%d", &n, &m); 
	for (int i = 1; i <= n; ++i) scanf("%d", &h[i]);
	std::sort(h + 1, h + n + 1), build(1, 1, n);
	for (int i = 1; i <= m; ++i) {
		scanf("%d", &x), pos = n - x + 1, cur = query(1, 1, n, pos);
		if (cur == 0) {
			std::cout << i - 1;
			return ;
		}
		int f1 = find(1, 1, n, cur), f2;
		if (cur != tree[1].max) f2 = find(1, 1, n, cur + 1);
		else f2 = n + 1;
		x -= (n - f2 + 1);
		if (f2 != n + 1) modify(1, 1, n, f2, n);
		modify(1, 1, n, f1, f1 + x - 1);
	}
	std::cout << m;
}

int main() {
	clock_t start = clock();
//	freopen("sequence.in", "r", stdin);
//	freopen("sequence.out", "w", stdout);
	solve(); 
	std::cerr << clock() - start << '\n';
	return 0;
}