bzoj 2038 [2009国家集训队]小Z的袜子(hose)

这题听说就是莫队的超级例题,所以来做一做。

大佬%d%a%l%a%o

这个博客写得不错。
我们要先学会分块。
我们可以将n个数分成根号n个块。
我们按照l所属的块为第一关键字,r为第二关键字来从小到大排序。
这样子每次l最多移动根号n,r最多移动n。但O(n*m)是不可能达到的。
所以时间复杂度可以过。
然后莫队要记得的是:

删去一个位置的时候要先更新ans再l++/r- -

添加一个位置的时候要先l- -/r++再更新ans

上标:

#include<cstdio>
#include<algorithm>
#include<cmath>
#define N 50010
#define ll long long
using namespace std;
struct node{int l,r,fr;}b[N];
int n,m,a[N],st,bl[N];
ll tot[N],X[N],Y[N],ans=0,gcd;

int cmp(node x,node y) {return bl[x.l]==bl[y.l] ? x.r<y.r : bl[x.l]<bl[y.l];}

inline int read()
{
	int x=0; char c=getchar();
	while (c<'0' || c>'9') c=getchar();
	while (c>='0' && c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x;
}

ll C(ll x) {return x*(x-1)/2;}

void update(int x,int add)
{
	ans-=C(tot[a[x]]);
	tot[a[x]]+=add;
	ans+=C(tot[a[x]]);
}

ll GCD(ll x,ll y) {return !y ? x:GCD(y,x%y);}

int main()
{
	freopen("2038.in","r",stdin);
	freopen("2038.out","w",stdout);
	n=read(),m=read();st=sqrt(n);
	for (int i=1;i<=n;i++)
		a[i]=read(),bl[i]=(i-1)/st+1;
	for (int i=1;i<=m;i++)
		b[i].l=read(),b[i].r=read(),b[i].fr=i;
	sort(b+1,b+m+1,cmp);
	int l=1,r=0;
	for (int i=1,fr;i<=m;i++)
	{
		while (l<b[i].l) update(l,-1),l++;
		while (l>b[i].l) l--,update(l,1);
		while (r<b[i].r) r++,update(r,1);
		while (r>b[i].r) update(r,-1),r--;
		fr=b[i].fr;
		if (l==r) {X[fr]=0,Y[fr]=1; continue;}
		X[fr]=ans,Y[fr]=C(r-l+1);
		gcd=GCD(X[fr],Y[fr]);
		X[fr]/=gcd,Y[fr]/=gcd;
	}
	for (int i=1;i<=m;i++)
		printf("%lld/%lld
",X[i],Y[i]);
	return 0;
}
转载需注明出处。
原文地址:https://www.cnblogs.com/jz929/p/11817582.html