[bzoj3289]Mato的文件管理_莫队_树状数组

Posted 2021-01-06 shurak

Mato的文件管理 bzoj-3289

题目大意：给定一个n个数的序列。m次询问：一段区间中的逆序对个数。

注释：$1le n,mle 5cdot 10^4$。

---

想法：

开始想这个题的大佬们，给您点儿提示吧：$O(nlognsqrt(n))$可过哦！

所以这个题就是莫队的裸题了。

我们的莫队上的区间在动的时候随时更新树状数组上的信息即可。、

然后碰见了一整块区间，我们就直接求逆序对即可，

最后，附上丑陋的代码... ...

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 50010 
using namespace std;
struct Query {int l,r,id;}q[N]; int L[N],R[N],maxval,tree1[N<<1],tree2[N<<1],blg[N],a[N],ans[N];
inline bool cmp_dispose(const Query &x,const Query &y) {return blg[x.l]==blg[y.l]?x.r<y.r:blg[x.l]<blg[y.l];}
inline char nc() {static char *p1,*p2,buf[100000]; return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}
int rd() {int x=0; char c=nc(); while(!isdigit(c)) c=nc(); while(isdigit(c)) x=(x<<3)+(x<<1)+(c^48),c=nc(); return x;}
inline int lowbit(int x) {return x&(-x);}
void update1(int x,int val) {for(int i=x;i>=1;i-=lowbit(i)) tree1[i]+=val;}
void update2(int x,int val) {for(int i=x;i<=maxval;i+=lowbit(i)) tree2[i]+=val;}
inline void update(int x,int val) {update1(x,val); update2(x,val);}
int query1(int x) {int ans=0; for(int i=x;i<=maxval;i+=lowbit(i)) ans+=tree1[i]; return ans;}
int query2(int x) {int ans=0; for(int i=x;i>=1;i-=lowbit(i)) ans+=tree2[i]; return ans;}
int main()
{
	int n=rd(); int t=sqrt(n); int blck=n/t;
	for(int i=1;i<=blck;i++)
	{
		int bfr=(i-1)*t;
		L[i]=bfr+1; R[i]=bfr+t;
		for(int j=1;j<=t;j++)
		{
			blg[bfr+j]=i;
			a[bfr+j]=rd();
			maxval=max(maxval,a[bfr+j]);
		}
	}
	if(blck*t<n)
	{
		L[blck+1]=blck*t; R[blck+1]=n; blck++;
		for(int j=blck*t+1;j<=n;j++)
		{
			blg[j]=blck; a[j]=rd(); maxval=max(maxval,a[j]);
		}
	}
	int m=rd();
	for(int i=1;i<=m;i++) q[i].l=rd(),q[i].r=rd(),q[i].id=i;
	sort(q+1,q+m+1,cmp_dispose);
	int point_l=q[1].l,point_r=q[1].r; for(int i=point_l;i<=point_r;i++)
	{
		ans[q[1].id]+=query1(a[i]+1); update(a[i],1);
	}
	for(int i=2;i<=m;i++)
	{
		int idx=q[i].id;
		while(point_l<q[i].l)
		{
			update(a[point_l],-1);
			ans[idx]-=query2(a[point_l]-1);
			point_l++;
		}
		while(point_l>q[i].l)
		{
			point_l--;
			update(a[point_l],1);
			ans[idx]+=query2(a[point_l]-1);
		}
		while(point_r<q[i].r)
		{
			point_r++;
			update(a[point_r],1);
			ans[idx]+=query1(a[point_r]+1);
		}
		while(point_r>q[i].r)
		{
			update(a[point_r],-1);
			ans[idx]-=query1(a[point_r]+1);
			point_r--;
		}
	}
	for(int i=1;i<=m;i++) printf("%d
",ans[i]);
	return 0;
}

小结：莫队真的强..