排列
题目背景:
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;
}