P4396 [AHOI2013]作业

最新推荐文章于 2024-10-08 11:39:04 发布

crazy morning

最新推荐文章于 2024-10-08 11:39:04 发布

阅读量88

点赞数

分类专栏：分块&莫队

本文链接：https://blog.csdn.net/weixin_49843717/article/details/119490709

版权

分块&莫队专栏收录该内容

7 篇文章 0 订阅

订阅专栏

本文介绍了如何结合莫队算法和值域分块技术来解决具有多组查询且不需在线处理的区间操作问题。通过实例展示了代码实现，包括分块的处理和数据结构维护，以及在排序询问中的应用。强调了莫队算法和输入序列分块、值域分块的区别，并提供了完整的C++代码示例。

摘要由CSDN通过智能技术生成

题目

思路

看见这种多测+无修改+不强制在线的就可以莫队。
然后这里询问的东西和数值有关，所以是值域分块
所以这破题是莫队套分块？
莫队就是大家耳熟能详的莫队，分块就是大家耳熟能详的分块……
莫队是没什么好说的，就老老实实搞
值域分块的时候~~也是老老实实搞~~ 就是要存3个：单点个数，块个数，块数值个数
另外排序询问用的分块是输入序列分块，值域分块是值域分块，不要搞混了~~话说这应该就我一个错~~
code:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<iostream>
#include<cmath>
using namespace std;
inline int read()
{
    int ret,c,f=1;
    while (((c=getchar())> '9'||c< '0')&&c!='-');
    if (c=='-') f=-1,ret=0;
    else ret=c-'0';
    while ((c=getchar())>='0'&&c<='9') ret=ret*10+c-'0';
    return ret*f;
}
int n,m,a[100006],bcnt,cnt,v,l,r,blo[100006],ans[100006],ans2[100006],nublo[100006],vs[100006],vb[100006],vv[100006];
struct f{
	int l,r,a,b,id;
} q[100006];
bool cmp(f a,f b)
{
	if (blo[a.l]==blo[b.l])
	{
		return a.r<b.r;
	}
	return blo[a.l]<blo[b.l];
}
int get_ans(int l,int r)
{
	int ans=0;
	r=min(r,v);
	if (l>r) return 0;
	if (nublo[l]==nublo[r])
	{
		for (int i=l;i<=r;i++) ans+=vs[i];
		return ans;
	}
	for (int i=nublo[l]+1;i<nublo[r];i++)
	{
		ans+=vb[i];
	}
	for (int i=l;nublo[i]==nublo[l]&&i<=v;i++) ans+=vs[i];
	for (int i=r;nublo[i]==nublo[r]&&i>=0;i--) ans+=vs[i];
	return ans;
}
int get_ans2(int l,int r)
{
	int ans=0;
	r=min(r,v);
	if (l>r) return 0;
	if (nublo[l]==nublo[r])
	{
		for (int i=l;i<=r;i++) ans+=(vs[i]>0);
		return ans;
	}
	for (int i=nublo[l]+1;i<nublo[r];i++)
	{
		ans+=vv[i];
	}
	for (int i=l;nublo[i]==nublo[l]&&i<=v;i++) ans+=(vs[i]>0);
	for (int i=r;nublo[i]==nublo[r]&&i>=0;i--) ans+=(vs[i]>0);
	return ans;
}
int main()
{
	n=read(),m=read();
	for (int i=1;i<=n;i++) a[i]=read(),v=max(a[i],v);
	for (int i=1;i<=m;i++)
	{
		q[i].id=i;
		q[i].l=read(),q[i].r=read(),q[i].a=read(),q[i].b=read();
	}
	bcnt=sqrt(v),cnt=sqrt(n);
	for (int i=1;i<=n;i++) blo[i]=(i-1)/cnt;
	for (int i=0;i<=v;i++) nublo[i]=i/bcnt;
	sort(q+1,q+1+m,cmp);
	int l=1,r=0;
	for (int i=1;i<=m;i++)
	{
		while (l<q[i].l)
		{
			vs[a[l]]--;
			vb[nublo[a[l]]]--;
			if (vs[a[l]]==0) vv[nublo[a[l]]]--;
			l++;
		}
		while (l>q[i].l)
		{
			l--;
			vs[a[l]]++;
			vb[nublo[a[l]]]++;
			if (vs[a[l]]==1) vv[nublo[a[l]]]++;
		}
		while (r<q[i].r)
		{
			r++;
			vs[a[r]]++;
			vb[nublo[a[r]]]++;
			if (vs[a[r]]==1) vv[nublo[a[r]]]++;
		}
		while (r>q[i].r)
		{
			vs[a[r]]--;
			vb[nublo[a[r]]]--;
			if (vs[a[r]]==0) vv[nublo[a[r]]]--;
			r--;
		}
		ans[q[i].id]=get_ans(q[i].a,q[i].b);
		ans2[q[i].id]=get_ans2(q[i].a,q[i].b);
	}
	for (int i=1;i<=m;i++) printf("%d %d\n",ans[i],ans2[i]);
	return 0;
}