【BZOJ4810】[Ynoi2017]由乃的玉米田 bitset+莫队

【BZOJ4810】[Ynoi2017]由乃的玉米田

Description

由乃在自己的农田边散步，她突然发现田里的一排玉米非常的不美。这排玉米一共有N株，它们的高度参差不齐。由乃认为玉米田不美，所以她决定出个数据结构题

这个题是这样的：

给你一个序列a，长度为n，有m次操作，每次询问一个区间是否可以选出两个数它们的差为x，或者询问一个区间是否可以选出两个数它们的和为x，或者询问一个区间是否可以选出两个数它们的乘积为x ，这三个操作分别为操作1,2,3选出的这两个数可以是同一个位置的数

Input

第一行两个数n,m

后面一行n个数表示ai

后面m行每行四个数opt l r x

opt表示这个是第几种操作，l,r表示操作的区间，x表示这次操作的x

定义c为每次的x和ai中的最大值，ai >= 0，每次的x>=2n,m,c <= 100000

Output

对于每个询问，如果可以，输出yuno，否则输出yumi

Sample Input

5 5
1 1 2 3 4
2 1 1 2
1 1 2 2
3 1 1 1
3 5 5 16
1 2 3 4

Sample Output

yuno
yumi
yuno
yuno
yumi

题解：区间询问当然是用莫队来水啊~

对于差为x，直接在bitset上搞就行了；对于和为x，我们维护一个翻转的bitset，然后转化成差为x的问题来搞；对于乘积为x，我们直接开个桶，用sqrt(x)的复杂度枚举x的约数就行了。

小号交题一次提交顺利成为status倒数第一，大号提交先TLE了两次，然后加了一堆优化后顺利成为status倒数第二，试问还有谁~~~所以千万不要拿我的代码当做标程来拍。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <bitset>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=100010;
bitset<maxn*2> b1,b2;
int n,m,B;
int v[maxn],st[maxn*2],ans[maxn];
struct QUERY
{
	int qa,qb,qk,qx,org;
}q[maxn];
bool cmp(QUERY a,QUERY b)
{
	return (a.qa/B==b.qa/B)?(a.qb<b.qb):(a.qa/B<b.qa/B);
}
void ins(int x)
{
	if(!st[x])	b1[x]=b2[100000-x]=1;
	st[x]++;
}
void del(int x)
{
	st[x]--;
	if(!st[x])	b1[x]=b2[100000-x]=0;
}
int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<'0'||gc>'9')	{if(gc=='-')f=-f;	gc=getchar();}
	while(gc>='0'&&gc<='9')	ret=ret*10+gc-'0',gc=getchar();
	return ret*f;
}
int main()
{
	n=rd(),m=rd();
	B=(int)sqrt((double)n);
	int i,j,l=1,r=0;
	for(i=1;i<=n;i++)	v[i]=rd();
	for(i=1;i<=m;i++)	q[i].qk=rd(),q[i].qa=rd(),q[i].qb=rd(),q[i].qx=rd(),q[i].org=i;
	sort(q+1,q+m+1,cmp);
	for(i=1;i<=m;i++)
	{
		while(l>q[i].qa)	ins(v[--l]);
		while(r<q[i].qb)	ins(v[++r]);
		while(l<q[i].qa)	del(v[l++]);
		while(r>q[i].qb)	del(v[r--]);
		if(q[i].qk==1)	ans[q[i].org]=((b1<<q[i].qx)&b1).any();
		if(q[i].qk==2)	ans[q[i].org]=((b1<<100000)&(b2<<q[i].qx)).any();
		if(q[i].qk==3)	for(j=1;j*j<=q[i].qx&&!ans[q[i].org];j++)	if(q[i].qx%j==0&&st[j]&&st[q[i].qx/j])	ans[q[i].org]=1;
	}
	for(i=1;i<=m;i++)
	{
		if(ans[i])	printf("yuno\n");
		else	printf("yumi\n");
	}
	return 0;
}