NOIP模拟（11.03）T2 排列

排列

题目背景：

11.03 NOIP模拟T2

分析：线段树

 

唯一一道可做题······

我们从n ~ 1来确定每个数的位置，显然，我们放入n的时候，他一定在一个逆序对数为0的位置，因为它的前面不可能有比它大的，其次，如果0位有多个，它一定在最后一个0位上，因为如果在前面的0位上，后面的0位上的数字一定会比它小，所以逆序对数不可能为0，那么我们现在确定了最大的n的位置，然后我们将n后面所有位置的逆序对数减一之后，那么去掉当前n的位置，我们就得到了一个新的逆序对数的排列，（因为n的影响相当于被完全消灭了，所以我们相当于重新讨论一个规模减小1的问题1 ~ n - 1）那么很显然的n - 1也应该放在当前序列的最后一个0位上，那么我们就可以一次递推下去，从大到小，每一次选择最后一个0位，然后将后面的数减一去掉贡献，那么我们就需要查询0位，以及区间减法，我们可以选择线段树维护区间减法和区间最小值，然后在线段树上二分找到最后一个0位然后进行区间减即可。

注意：因为选择了0位之后，你不能让当前位置影响答案，所以应该将当前位置赋值为INF，这样就不会在之后被选到了。

Source:

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

const int MAXN = 100000 + 10;
const int INF = 1000000000;

int n;
int a[MAXN], p[MAXN];
struct node {
	int min, tag;
} tree[MAXN << 2];

inline void modify(int k, int x) {
	tree[k].min -= 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);
}

inline void build(int k, int l, int r) {
	if (l == r) return (void)(tree[k].min = a[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) {
	if (l == r) return l;
	int mid = l + r >> 1;
	push_down(k);
	if (tree[k << 1 | 1].min == 0) return query(k << 1 | 1, mid + 1, r);
	else return query(k << 1, l, mid);
}

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 modify(int k, int l, int r, int pos) {
	if (l == r) return (void)(tree[k].min = INF);
	int mid = l + r >> 1;
	push_down(k);
	(pos <= mid) ? modify(k << 1, l, mid, pos) : 
		modify(k << 1 | 1, mid + 1, r, pos);
	update(k);
}

inline void read_in() {
	scanf("%d", &n);
	for (int i = 1; i <= n; ++i) scanf("%d", &a[i]);
	for (int i = n; i >= 1; --i) a[i] -= a[i - 1];
}

inline void solve() {
	build(1, 1, n);
	for (int i = n; i >= 1; --i) {
		int pos = query(1, 1, n);
		p[pos] = i, modify(1, 1, n, pos);
		if (pos != n) modify(1, 1, n, pos + 1, n);
	}
	for (int i = 1; i <= n; ++i) std::cout << p[i] << " ";
}

int main() {
//	freopen("premu.in", "r", stdin);
//	freopen("premu.out", "w", stdout);
	read_in();
	solve();
	return 0;
}