JZOJ 5947. 【NOIP2018模拟11.02】初音未来（miku）

Description

Hercier作为一位喜爱Hatsune Miku的OIer，痛下决心，将Vocaloid买回了家。打开之后，你发现界面是一个长为n的序列，代表音调，并形成了全排列。你看不懂日语，经过多次尝试，你只会用一个按钮：将一段区间按升序排序。不理解音乐的Hercier决定写一个脚本，进行m次操作，每次对一段区间进行操作。可惜Hercier不会写脚本，他找到了在机房里的你，请你模拟出最后的结果。

Input

Output

Sample Input

5 2 1 5
5 1 4 2 3
1 3
2 4

样例2
见下发文件

Sample Output

1 2 4 5 3

样例解释：

Data Constraint

Solution

• 这题只有升序排序，那么序列肯定会趋于有序。

• 长度为 $n$ 的序列逆序对个数最多为 $\frac{n(n-1)}{2}$ ，排序后逆序对个数不会变多。

• 于是我们用线段树记录一个区间是否还存在逆序对，

• 那么对区间进行升序操作时若该区间已经无逆序对，说明区间已经有序，就不用再操作了。

• 于是只要区间还有逆序对，我们就暴力递归下去修改，最多有 $O(n^2)$ 个修改。

• 所以复杂度 $O((n^2+m)logn)$ 。

Code

#include<cstdio>
#include<algorithm>
#include<cctype>
using namespace std;
const int N=1505;
int n,qx,qy,qz;
bool pd;
int a[N];
bool f[N<<2];
inline int read()
{
    int X=0,w=0; char ch=0;
    while(!isdigit(ch)) w|=ch=='-',ch=getchar();
    while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return w?-X:X;
}
void write(int x)
{
	if(x>9) write(x/10);
	putchar(x%10+'0');
}
void make(int v,int l,int r)
{
	if(l==r)
	{
		f[v]=a[l-1]>a[l];
		return;
	}
	int mid=l+r>>1;
	make(v<<1,l,mid);
	make(v<<1|1,mid+1,r);
	f[v]=f[v<<1]|f[v<<1|1];
}
void update(int v,int l,int r)
{
	if(l==r)
	{
		f[v]=a[l-1]>a[l];
		return;
	}
	int mid=l+r>>1;
	if(qz<=mid) update(v<<1,l,mid); else update(v<<1|1,mid+1,r);
	f[v]=f[v<<1]|f[v<<1|1];
}
void modify(int v,int l,int r)
{
	if(!f[v]) return;
	if(l==r)
	{
		pd=true;
		bool pl=a[l-2]>a[l-1],pr=a[l]>a[l+1];
		swap(a[l],a[l-1]);
		if((a[l-2]>a[l-1])^pl)
		{
			qz=l-1;
			update(1,1,n);
		}
		if(l<n && (a[l]>a[l+1])^pr)
		{
			qz=l+1;
			update(1,1,n);
		}
		f[v]=false;
		return;
	}
	int mid=l+r>>1;
	modify(v<<1,l,mid);
	if(pd) return;
	modify(v<<1|1,mid+1,r);
	f[v]=f[v<<1]|f[v<<1|1];
}
void change(int v,int l,int r)
{
	if(!f[v]) return;
	if(qx<=l && r<=qy)
	{
		modify(v,l,r);
		return;
	}
	int mid=l+r>>1;
	if(qx<=mid) change(v<<1,l,mid);
	if(pd) return;
	if(qy>mid) change(v<<1|1,mid+1,r);
	f[v]=f[v<<1]|f[v<<1|1];
}
int main()
{
	freopen("miku.in","r",stdin);
	freopen("miku.out","w",stdout);
	n=read();
	int m=read(),L=read(),R=read();
	for(int i=1;i<=n;i++) a[i]=read();
	make(1,1,n);
	while(m--)
	{
		qx=read()+1,qy=read();
		if(qx>qy) continue;
		pd=false;
		change(1,1,n);
		while(pd)
		{
			pd=false;
			change(1,1,n);
		}
	}
	for(int i=L;i<=R;i++) write(a[i]),putchar(' ');
	return 0;
}