No.7 线段树

crf 出生的第三秒

1 Description
crf 是一个天才。
他出生的第三秒，偶然看到了一张简谱，可能是医院的工作人员闲的没事的时候创作的。
当时的他并不能理解简谱和音乐的关系（但这并不妨碍他今后成为世界著名小提琴演奏家），他只是单纯的看着这张写满1-7 的纸，觉得非常缺乏美感。他认为有序的数列才是最美的。
但毕竟他才出生第三秒，他的mogic 也只是初步显现，他可以改变一段连续的音符的顺序，但他的mogic 不足以一次改掉整张简谱。
这是crf 这辈子第一次犯蠢，也是唯一一次犯蠢。他尝试着将一段一段的音符调成升序，但他并没有意识到就算这样调了也不能让整个序列变得有序。
crf 的mogic 可以选出一段连续的音符，将它们重新排为升序之后再放回，他尝试了一些这样的操作之后发现并没有将序列变得有序就放弃了。
后来的程序员为了纪念伟大的crf，将他唯一犯蠢的这段经历重新变成了一道算法题。
虽然这道题大部分题面都是对crf 的赞美，但可能你们不敢兴趣所以就直接上题吧。
给定一个长度为n 的只有1-7 的数字的序列，有q 次操作，每次操作是将一段连续的数字拿出来，排成升序或者降序后放回。求最终序列。

2 Input
输入的第一行为两个整数n; q，表示序列的长度和操作的数量。
接下来一行有n 个整数，表示crf 的数字序列。
接下来q 行，每行有三个整数l; r; op，l; r 含义如题目中所示，op 为0 表示升序，为1 表示降序。

3 Output
输出共1 行，包含n 个整数，表示最终的序列，中间用空格隔开。

Sample Input
5 2
3 1 2 5 4
1 3 0
2 5 1
Sample Output
1 5 4 3 2

4 Hint
对于30% 的数据，保证1 <= n, q <= 1000。
对于100% 的数据，保证1 <= n , q <= 100000。

思路：
线段树，因为只有7种数字，用线段树维护每一段的7种数字分别有多少，每次操作即为一次查询和七次覆盖。

#include <cstdio>#include <cmath>#include <algorithm>#include <cstring>using namespace std;const int maxn = 100010;int a[maxn], n, T, cnt[20], sum[20];inline const int read(){    register int f = 1, x = 0;    register char ch = getchar();    while(ch < '0' || ch > '9') ch = getchar();    while(ch >= '0' && ch <= '9') x = (x<<3) + (x<<1) + ch - '0', ch = getchar();    return x;}struct Node{    Node *ls, *rs;    int cnt[8], flag;    Node(){ flag = 0 ; memset( cnt , 0 , sizeof( cnt )) ;}    Node operator+( const Node &a ) const{        Node tmp = Node();        for(int i=1; i<=7; ++i)            tmp.cnt[i] = cnt[i] + a.cnt[i] ;        return tmp;    }    void pushdown(int lf, int rg){        if( flag ){            int mid = (lf + rg) >> 1;            ls->cnt[flag] = mid - lf + 1;            ls->flag = flag;            rs->cnt[flag] = rg - mid;            rs->flag= flag;            for(int i=1; i<=7; ++i)                 if(i != flag) ls->cnt[i] = rs->cnt[i] = 0;            flag = 0;        }    }    void update(){        for(int i=1; i<=7; ++i)            cnt[i] = ls->cnt[i] + rs->cnt[i];    }}pool[maxn*20], *tail = pool, *rot;inline Node* build(int lf, int rg){    Node *nd = ++tail;    if(lf == rg) {        for(int i=1; i<=7; ++i)             nd->cnt[i] = 0;        nd->cnt[a[lf]]++;         return nd;    }    int mid = (lf + rg) >> 1;    nd->ls = build(lf, mid);    nd->rs = build(mid+1, rg);    nd->update();    return nd;}inline void modify(Node *nd, int lf, int rg, int L, int R, int data){    if(L <= lf && rg <= R){        nd->cnt[data] = rg - lf + 1;         nd->flag = data;        for(int i=1; i<=7; ++i)             if(i != data) nd->cnt[i] = 0;        return;    }     nd->pushdown(lf, rg);    int mid = (lf + rg) >> 1;    if(L <= mid) modify(nd->ls, lf, mid, L, R, data);    if(R > mid) modify(nd->rs, mid+1, rg, L, R, data);    nd->update();    return;}/*inline int query(Node *nd, int lf, int rg, int L, int R, int data){    if(L <= lf && rg <= R)         return nd->cnt[data];    nd->pushdown(lf, rg);    int mid = (lf + rg) >> 1, rt = 0;    if(L <= mid) rt += query(nd->ls, lf, mid, L, R, data);    if(R > mid) rt += query(nd->rs, mid+1, rg, L, R, data);    nd->update();    return rt;}*/inline Node query(Node *nd, int lf, int rg , int L , int R ){    if( L <= lf && rg <= R )        return *nd;    Node rt = Node() ;    nd->pushdown(lf, rg) ;    int mid = ( lf + rg ) >> 1 ;    if( L <= mid ) rt = rt + query( nd->ls, lf , mid , L , R ) ;    if( R >  mid ) rt = rt + query( nd->rs, mid+1, rg , L , R ) ;    return rt ;}//返回一个结构体（内有cnt[1~7]），就只需要run一次，而不是七次了。 inline int query_point(Node *nd, int lf, int rg, int pos){    if(lf == rg){        for(int i=1; i<=7; ++i)            if(nd->cnt[i]) return i;    }    nd->pushdown(lf, rg);    int mid = (lf + rg) >> 1, rt = 0;    if(pos <= mid) rt = query_point(nd->ls, lf, mid, pos);    else rt = query_point(nd->rs, mid+1, rg, pos);    nd->update();    return rt;}int main(){    freopen("third.in", "r", stdin);    freopen("third.out", "w", stdout);    n = read(), T = read();    for(register int i=1; i<=n; ++i) a[i] = read();    rot = build(1, n);    for(int i=1; i<=T; ++i){        int l = read(), r = read(), opt = read();        Node t = query(rot, 1, n, l, r);//答案结构体         for(int j=1; j<=7; ++j) sum[j] = t.cnt[j];        int st = l;        if( !opt ){            for(int j=1; j<=7; ++j)                if( sum[j] ) modify(rot, 1, n, st, st+sum[j]-1, j), st += sum[j];        }        else{            for(int j=7; j>=1; --j)                if( sum[j] ) modify(rot, 1, n, st, st+sum[j]-1, j), st += sum[j];        }    }    for(register int j=1; j<=n; ++j) printf("%d ", query_point(rot, 1, n, j));    return 0;}