Kattis - cupid莫队

最新推荐文章于 2024-08-07 09:16:33 发布

原创最新推荐文章于 2024-08-07 09:16:33 发布 · 245 阅读

0 ·

CC 4.0 BY-SA版权

本文通过一个具体的题目实例，详细介绍了裸莫队算法的实现细节，包括如何通过预处理和区间更新来解决区间查询问题，并提供了完整的代码示例。

题目链接点这里

裸莫队,,然而少加一个return..re2小时

#include<algorithm>
#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
using namespace std;
#define mem(x,y) memset(x,y,sizeof(x))
#define FIN freopen("input.txt","r",stdin)
#define fuck(x) cout<<x<<endl
const int MX=1000100;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
int n,m,k;
int cnta[MX],cntb[MX],sum[MX];
int a[MX],b[MX],block;
struct Nod
{
	int l,r,id,ans;
	bool operator < ( const Nod &a)const
	{
		if(l/block!=a.l/block) return l/block<a.l/block;
		else return r<a.r;
	}
} p[MX];
bool cmp(Nod a,Nod b)
{
	return a.id<b.id;
}
void del(int x,int &ans)
{
	if(cnta[a[x]]<=cntb[a[x]]) ans--;
	cnta[a[x]]--;
	if(cntb[b[x]]<=cnta[b[x]]) ans--;
	cntb[b[x]]--;
}
void add(int x,int &ans)
{
	if(cnta[a[x]]<cntb[a[x]]) ans++;
	cnta[a[x]]++;
	if(cntb[b[x]]<cnta[b[x]]) ans++;
	cntb[b[x]]++;
}
int main()
{
	FIN;
	scanf("%d%d%d",&n,&m,&k);
	block=sqrt(n+0.5);
	for(int i=0; i<n; i++)scanf("%d",&a[i]);
	for(int i=0; i<n; i++)scanf("%d",&b[i]);
	for(int i=0; i<m; i++)
	{
		scanf("%d%d",&p[i].l,&p[i].r);
		p[i].id=i;
	}
	sort(p,p+m);
	mem(cntb,0);
	mem(cnta,0);
	int l=0,r=0;
	cnta[a[l]]++;
	cntb[b[l]]++;
	int ans;
	if(a[l]==b[l]) ans=1;
	else ans=0;
	for(int i=0; i<m; i++)
	{

		while(r>p[i].r)	del(r--,ans);
		while(r<p[i].r) add(++r,ans);
		while(l<p[i].l) del(l++,ans);
		while(l>p[i].l) add(--l,ans);
		sum[p[i].id]=ans;
	}
	for(int i=0; i<m; i++) printf("%d\n",sum[i]);
	return 0;
}