P6186 [NOI Online 提高组]冒泡排序

P6186 [NOI Online 提高组]冒泡排序

题目描述

分析

对于这道题，我们可以先求出原序列中的逆序对数量，然后执行两个操作

为了方便进行操作，我们先将原序列离散化到1~n，然后维护一个数组con[i]表示序列的第i个数前有con[i]个数比它大，那么逆序对数量cnt=(sumlimits_{i=1}^ncon[i])

对于第一个操作，判断交换的两个数的大小关系，然后改变cnt和con[i]即可

对于第二个操作，由引理1可得：k轮冒泡排序后减少的逆序对数量是(sumlimits_{con[i]leq k}con[i]+sumlimits_{con[i]>k}k)分别维护两个树状数组即可

引理1

在一个序列进行一轮冒泡排序之后，(con[i]=max(con[i+1]-1,0))，特别的(con[n]=max(con[1]-1,0))

证明：

第k轮冒泡排序时，对于序列中第i个数((iin[2,k-2]))，要么置换到(i-1)的位置，要么置换到(i+1)的位置，前者有con[i-1]=con[i]-1，后者有con[i]=con[i+1]-1

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;

const int N=1e6+5;
struct Node
{
	ll val,id;
}a[N];
struct node
{
	ll x,y;
}que[N];
ll n,m,c[4*N],sq[N],con[N],cnt,c1[N],c2[N];

bool cmp(Node x,Node y)
{
	return x.val < y.val;
}

ll lowbit(ll x)
{
	return x & (-x);
}

void add(ll x,ll val,ll c[N])
{
	if( x==0 ) x=1;
	while( x <= n )
	{
		c[x]+=val;
		x+=lowbit(x);
	}
	return;
}

ll sum(ll x,ll c[N])
{
	ll tem=0;
	while( x )
	{
		tem+=c[x];
		x-=lowbit(x);
	}
	return tem;
}

void Input()
{
	scanf("%lld%lld",&n,&m);
	for(ll i=1;i<=n;i++) scanf("%lld",&a[i].val),a[i].id=i;
	for(ll i=1;i<=m;i++) scanf("%lld%lld",&que[i].x,&que[i].y);
	return;
}

void Prepare()
{
	sort(a+1,a+n+1,cmp);
	for(ll i=1;i<=n;i++) sq[a[i].id]=i;
	for(ll i=1;i<=n;i++) add(sq[i],1,c1),con[i]=i-sum(sq[i],c1),cnt+=con[i];
	return;
}

void deal1(ll id)
{
	add(con[id],-1,c1); add(con[id],-con[id],c2);
	add(con[id+1],-1,c1); add(con[id+1],-con[id+1],c2);
	if( sq[id] < sq[id+1] ) con[id]++,cnt++;
	else con[id+1]--,cnt--;
	swap(con[id],con[id+1]);
	swap(sq[id],sq[id+1]);
	add(con[id],1,c1); add(con[id],con[id],c2);
	add(con[id+1],1,c1); add(con[id+1],con[id+1],c2);
	return;
}

void deal2(ll k)
{
	ll tem;
	if( k >= n ) k=n-1;
	tem=sum(k,c2);
	tem+=k*(n-sum(k,c1));
	printf("%lld
",cnt-tem);
	return;
}

void work()
{
	memset(c1,0,sizeof(c1));
	for(ll i=1;i<=n;i++) add(con[i],1,c1),add(con[i],con[i],c2);
	for(ll i=1;i<=m;i++)
	{
		if( que[i].x==1 ) deal1(que[i].y);
		else deal2(que[i].y);
	}
	return;
}

int main()
{
	Input();
	Prepare();
	work();
	return 0;
}