BZOJ 4552: [Tjoi2016&Heoi2016]排序

Description

在2016年，佳媛姐姐喜欢上了数字序列。因而他经常研究关于序列的一些奇奇怪怪的问题，现在他在研究一个难题

，需要你来帮助他。这个难题是这样子的：给出一个1到n的全排列，现在对这个全排列序列进行m次局部排序，排

序分为两种：1:(0,l,r)表示将区间[l,r]的数字升序排序2:(1,l,r)表示将区间[l,r]的数字降序排序最后询问第q

位置上的数字。

Input

输入数据的第一行为两个整数n和m。n表示序列的长度，m表示局部排序的次数。1 <= n, m <= 10^5第二行为n个整

数，表示1到n的一个全排列。接下来输入m行，每一行有三个整数op, l, r, op为0代表升序排序，op为1代表降序

排序, l, r 表示排序的区间。最后输入一个整数q，q表示排序完之后询问的位置, 1 <= q <= n。1 <= n <= 10^5

，1 <= m <= 10^5

Output

 输出数据仅有一行，一个整数，表示按照顺序将全部的部分排序结束后第q位置上的数字。

Sample Input

6 3
1 6 2 5 3 4
0 1 4
1 3 6
0 2 4
3

Sample Output

5

本想练习线段树合并的。。。。。

但是标算太过奇特，想象奇异，令人不能自拔。。。。。

如何用线段树多次区间排序？

用线段树查出每一个值在区间中的出现次数，然后从小到大从左到右插入.......

如果只有两种不同的值，那么这个排序方式是logn的，十分优秀。。。。。。

那么我们用二分答案，

设答案为a

如果我们能够判断a与mid的关系，那么这就是一个简单的二分答案

那么我们将小于等于mid的数设为1，其他设为0

用线段树排完序后如果我们要查询的位置上为1，那么a<=mid,否则a>mid

多简单的思路。。。。。。

#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#define maxn 100005
#define lc now<<1
#define rc now<<1|1
using namespace std;

int n,m,q,seq[maxn],a[maxn];
int sum[4*maxn],lazy[4*maxn],L[4*maxn],R[4*maxn];
inline void upd(int now){ sum[now]=sum[lc]+sum[rc]; }
inline void dcd(int now)
{
	if(lazy[now]!=-1)
	{
		lazy[lc]=lazy[rc]=lazy[now];
		sum[lc]=(R[lc]-L[lc]+1)*lazy[lc];
		sum[rc]=(R[rc]-L[rc]+1)*lazy[rc];
		lazy[now]=-1;
	}
}

void build(int now,int l,int r)
{
	L[now]=l,R[now]=r,lazy[now]=-1;
	if(l==r) 
	{
		sum[now]=seq[l];
		return;
	}
	int mid=(l+r)>>1;
	build(lc,l,mid);
	build(rc,mid+1,r);
	upd(now);
}

void insert(int now,int ql,int qr,int val)
{
	dcd(now);
	if(qr<ql || L[now]>qr || R[now]<ql) return;
	if(ql<=L[now] && R[now]<=qr) 
	{
		lazy[now]=val;
		sum[now]=(R[now]-L[now]+1)*lazy[now];
		return;
	}
	insert(lc,ql,qr,val);
	insert(rc,ql,qr,val);
	upd(now);
}

int Query(int now,int ql,int qr)
{
	dcd(now);
	if(qr<ql || L[now]>qr || R[now]<ql) return 0;
	if(ql<=L[now] && R[now]<=qr)
		return sum[now];
	int tmp=Query(lc,ql,qr)+Query(rc,ql,qr);
	upd(now);
	return tmp;
}

int x[maxn],y[maxn],op[maxn];

bool check(int num)
{
	for(int i=1;i<=n;i++) seq[i]=(a[i]<=num);
	build(1,1,n);
	for(int i=1,sum1,sum2;i<=m;i++)
	{
		sum1=Query(1,x[i],y[i]);
		sum2=(y[i]-x[i]+1)-sum1;
		if(op[i]==0) insert(1,x[i],x[i]+sum1-1,1),insert(1,x[i]+sum1,y[i],0);
		else insert(1,x[i],x[i]+sum2-1,0),insert(1,x[i]+sum2,y[i],1);
	}
	return Query(1,q,q)==1;
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=m;i++)
		scanf("%d%d%d",&op[i],&x[i],&y[i]);
	scanf("%d",&q);
	
	int l=1,r=n,mid;
	for(;l<r;)
	{
		mid=(l+r)>>1;
		if(check(mid)) r=mid;
		else l=mid+1;
	}
	printf("%d\n",l);
}