【JZOJ4744】同余

最新推荐文章于 2022-10-30 18:00:00 发布

原创最新推荐文章于 2022-10-30 18:00:00 发布 · 201 阅读

0 ·

CC 4.0 BY-SA版权

杂题同时被 3 个专栏收录

42 篇文章

订阅专栏

排序

34 篇文章

订阅专栏

分块

4 篇文章

订阅专栏

部署运行你感兴趣的模型镜像

description

analysis

对于 $p < = 100$ 的询问，可以设 $m o d [i] [j]$ 表示模 $i$ 余 $j$ 的数有几个
对于 $p > 100$ 的询问，设 $n u m [i]$ 表示 $i$ 有多少个
考虑把询问拆成两个， $l - 1$ 的删除和 $r$ 的插入操作
把询问排序，然后 $1$ 到 $n$ 扫过去，扫的中途更改 $m o d$ 和 $n u m$
当前位存在询问时，若 $p < = 100$ 直接询问 $m o d$ 数组得出答案
否则即答案为 $p k + q$ 的数的个数，枚举 $k$ 在 $n u m$ 数组里得出答案
最后把答案加到原来的询问里就好了

code

#pragma GCC optimize("O3")
#pragma G++ optimize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#define MAXN 200005
#define ll long long
#define reg register ll
#define fo(i,a,b) for (reg i=a;i<=b;++i)
#define fd(i,a,b) for (reg i=a;i>=b;--i)
#define O3 __attribute__((optimize("-O3")))

using namespace std;

ll a[MAXN],num[MAXN],mod[105][105],ans[MAXN];
ll n,m,tot,mx;

struct node
{
	ll x,id,p,q,val;
}f[MAXN*2];

O3 inline ll read()
{
	ll x=0,f=1;char ch=getchar();
	while (ch<'0' || '9'<ch){if (ch=='-')f=-1;ch=getchar();}
	while ('0'<=ch && ch<='9')x=x*10+ch-'0',ch=getchar();
	return x*f;
}
O3 inline bool cmp(node a,node b)
{
	return a.x<b.x;
}
O3 int main()
{
	//freopen("tongyu.in","r",stdin);
	n=read(),m=read();
	fo(i,1,n)mx=max(mx,a[i]=read());
	fo(i,1,m)
	{
		ll l=read(),r=read(),x=read(),y=read();
		f[++tot].x=l-1,f[tot].id=i,f[tot].p=x,f[tot].q=y,f[tot].val=-1;
		f[++tot].x=r,f[tot].id=i,f[tot].p=x,f[tot].q=y,f[tot].val=1;
	}
	sort(f+1,f+tot+1,cmp);
	ll i=1,j=1;
	while (f[j].x<1)++j;
	for (;j<=tot;++i)
	{
		++num[a[i]];
		fo(k,1,100)++mod[k][a[i]%k];
		while (j<=tot && f[j].x==i)
		{
			if (f[j].p<=100)ans[f[j].id]+=f[j].val*mod[f[j].p][f[j].q];
			else
			{
				for (ll k=f[j].q;k<=mx;k+=f[j].p)
				ans[f[j].id]+=f[j].val*num[k];
			}
			++j;
		}
	}
	fo(i,1,m)printf("%lld\n",ans[i]);
	return 0;
}