【分块+BIT】BZOJ3787 Gty的文艺妹子序列

【题目】
BZOJ
正整数序列$A$，支持单点修改，区间逆序对数，所有数字不超过$5\times 10^4$。强制在线。

【解题思路】
毒瘤数据结构题不妨考虑分块，设分块大小为$S$。

维护$su{m}_{i,j}$表示前$i$块中元素$j$的个数，每次修改只会影响$O(S)$个值，同样用树状数组维护可以做到$O(S\log n)$。
由于我们询问时需要对第二维进行区间查询，所以第二维要用树状数组维护起来。

维护$an{s}_{i,j}$表示第$i$块和第$j$块之间元素组成的逆序对个数，每次修改只会影响$O(\frac{n}{S})$个值，利用上面维护的$sum$可以快速进行更新，做到$O(\frac{n}{S}\log n)$。
由于我们询问时需要对第二维进行区间查询，所以第二维要用树状数组维护起来。

最后维护$nu{m}_i$表示第$i$块内逆序对个数，每次修改可以$O(S)$暴力算。

询问时分成四部分：块内统计$nu{m}_i$可以做到$O(\frac{n}{S})$，块间通过维护$an{s}_{i,j}$第二维的树状数组可以做到$O(\frac{n}{S}\log S)$，散块和整块之间可以通过维护$su{m}_{i,j}$第二维的树状数组可以做到$O(S\log S)$，散块部分可以暴力维护做到$O(S\log n)$

当$S=\sqrt{n}$时可以取到最优复杂度$O(n\sqrt{n}\log n)$

想清楚了很好写。所以一点也不毒瘤。
【参考代码】

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

typedef long long ll;
const int N=5e4+10,L=233;

namespace IO
{
	int read()
	{
		int ret=0;char c=getchar();
		while(!isdigit(c)) c=getchar();
		while(isdigit(c)) ret=ret*10+(c^48),c=getchar();
		return ret;
	}
	void write(int x){if(x>9)write(x/10);putchar(x%10^48);}
	void writeln(int x){write(x);putchar('\n');}
}
using namespace IO;

namespace BITisgood
{
#define lowbit(x) (x&(-x))
	struct BIT
	{
		vector<int>c;int lim;
		void init(int x){c.clear();c.resize(x+1,0);lim=x;}
		void update(int x,int v){for(;x<=lim;x+=lowbit(x))c[x]+=v;}
		int qry(int x){int res=0;for(;x;x-=lowbit(x))res+=c[x];return res;}
		int query(int l,int r){return l>r?0:qry(r)-qry(l-1);}
	};
#undef lowbit
}
using namespace BITisgood;

namespace DreamLolita
{
	int n,Q,cnt,a[N],bl[N];
	ll fans;
	struct Block
	{
		int num,st,ed;
	}B[N/L+10];
	BIT sum[N/L+10],ans[L+10],tmp;//wrong because ans[L],wrong memory!
	ll Query(int l,int r)
	{
		ll res=0;
		if(bl[l]==bl[r]) 
		{
			for(int i=l;i<=r;++i) res+=tmp.query(a[i]+1,n),tmp.update(a[i],1);
			for(int i=l;i<=r;++i) tmp.update(a[i],-1);
			return res;
		}
		for(int i=bl[l]+1;i<bl[r];++i) 
		{
			res+=B[i].num;
			res+=ans[i].query(i+1,bl[r]-1);
		}
		if(bl[l]+1<=bl[r]-1)
		{
			for(int i=l;i<=B[bl[l]].ed;++i) res+=sum[bl[r]-1].query(1,a[i]-1)-sum[bl[l]].query(1,a[i]-1);
			for(int i=B[bl[r]].st;i<=r;++i) res+=sum[bl[r]-1].query(a[i]+1,n)-sum[bl[l]].query(a[i]+1,n);
		}
		for(int i=l;i<=B[bl[l]].ed;++i) res+=tmp.query(a[i]+1,n),tmp.update(a[i],1);
		for(int i=B[bl[r]].st;i<=r;++i) res+=tmp.query(a[i]+1,n),tmp.update(a[i],1);
		for(int i=l;i<=B[bl[l]].ed;++i) tmp.update(a[i],-1);
		for(int i=B[bl[r]].st;i<=r;++i) tmp.update(a[i],-1);
		return res;
	}
	void Update(int p,int v)
	{
		for(int i=1;i<bl[p];++i) ans[i].update(bl[p],-(sum[i].query(a[p]+1,n)-sum[i-1].query(a[p]+1,n)));
		for(int i=bl[p]+1;i<=cnt;++i) ans[bl[p]].update(i,-(sum[i].query(1,a[p]-1)-sum[i-1].query(1,a[p]-1)));
		for(int i=bl[p];i<=cnt;++i) sum[i].update(a[p],-1),sum[i].update(v,1); a[p]=v;
		for(int i=1;i<bl[p];++i) ans[i].update(bl[p],(sum[i].query(a[p]+1,n)-sum[i-1].query(a[p]+1,n)));
		for(int i=bl[p]+1;i<=cnt;++i) ans[bl[p]].update(i,(sum[i].query(1,a[p]-1)-sum[i-1].query(1,a[p]-1)));
		B[bl[p]].num=0;
		for(int i=B[bl[p]].st;i<=B[bl[p]].ed;++i) B[bl[p]].num+=tmp.query(a[i]+1,n),tmp.update(a[i],1);
		for(int i=B[bl[p]].st;i<=B[bl[p]].ed;++i) tmp.update(a[i],-1);
	}
	void solution()
	{
		n=read();cnt=(n-1)/L+1;
		for(int i=1;i<=n;++i) a[i]=read(),bl[i]=(i-1)/L+1,B[bl[i]].ed=i;
		for(int i=n;i;--i) B[bl[i]].st=i;
		for(int i=0;i<=cnt;++i) sum[i].init(n),ans[i].init(cnt);
		tmp.init(n);
		for(int i=1;i<=cnt;++i) 
		{
			for(int j=0;j<=n;++j) sum[i].c[j]=sum[i-1].c[j];
			for(int j=B[i].st;j<=B[i].ed;++j) 
			{
				sum[i].update(a[j],1);
				for(int k=1;k<i;++k) ans[k].update(i,sum[k].query(a[j]+1,n)-sum[k-1].query(a[j]+1,n));
				for(int k=B[i].st;k<j;++k) if(a[k]>a[j]) B[i].num++;
			}
		}
		Q=read();
		while(Q--)
		{
			int op=read(),x=read()^fans,y=read()^fans;
			if(op) Update(x,y);
			else printf("%lld\n",(fans=Query(x,y)));
		}
	}
}

int main()
{
#ifndef ONLINE_JUDGE
	freopen("BZOJ3787.in","r",stdin);
	freopen("BZOJ3787.out","w",stdout);
#endif
	DreamLolita::solution();
	return 0;
}