bzoj4810 [Ynoi2017]由乃的玉米田 莫队+bitset

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

HINT

分析：这题的想法挺好，不是很难但是也不水quq。
我们用bitset存储权值数组，每一位为0表示这个下标的权值不存在，否则存在。
那么我们对于询问，用莫队离线处理以后，更新的时候，如果是减法，即a-b=x，那么我们移项，a=b+x。我们已知x，那么对于这个区间内的bitset，我们向后移动x位，然后和原来的and一下就好了，加法把bitset反过来就行了，乘法随便yy一波就好。。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<iostream>
#include<bitset>
#include<queue>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define N 100005
#define M N

using namespace std;
int n,m,k;
typedef long long ll;
bitset<N> g,f,rev;
int pos[N],mx,t[N],a[N];
struct node
{
    int id,l,r,op,x;
    int ans;
}q[N];
inline bool cmp(node a,node b)
{
    return pos[a.l]<pos[b.l]||pos[a.l]==pos[b.l]&&a.r<b.r;
}
inline bool cmpid(node a,node b)
{
    return a.id<b.id;
}
inline void updata(int x,int y)
{
    if (y==1)
    {
        if (!t[x])f[x]=1,rev[mx-x]=1;
        t[x]++;
    }
    else 
    {
        t[x]--;
        if (!t[x])f[x]=0,rev[mx-x]=0;
    }
}
inline int query(int op,int x)
{
    if (op==1)
    {
        g=f;
        g>>=x;
        g&=f;
        if (g.count())return 1;
        else return 0;
    }
    else if (op==2)
    {
        g=rev;
        g>>=(mx-x);
        g&=f;
        if (g.count())return 1;
        else return 0;
    }
    else
    {
        int w=sqrt(x);
        fo(i,1,w)
        if (x%i==0&&f[i]&&f[x/i])return 1;
        return 0;
    }
}
inline void solve()
{
    //updata(a[1],1);
    int l=1,r=0;
    fo(i,1,m)
    {
        for(;r<q[i].r;r++)updata(a[r+1],1);
        for(;l>q[i].l;l--)updata(a[l-1],1);
        for(;r>q[i].r;r--)updata(a[r],-1);
        for(;l<q[i].l;l++)updata(a[l],-1);
        q[i].ans=query(q[i].op,q[i].x);
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    int block=int(sqrt(n));
    fo(i,1,n)scanf("%d",&a[i]),pos[i]=(i+block-1)/block,mx=max(mx,a[i]);
    fo(i,1,m)
    {
        scanf("%d%d%d%d",&q[i].op,&q[i].l,&q[i].r,&q[i].x);
        mx=max(q[i].x,mx);
        q[i].id=i;
    }
    sort(q+1,q+1+m,cmp);
    solve();
    sort(q+1,q+1+m,cmpid);
    fo(i,1,m)
    if (q[i].ans)printf("yuno\n");
    else printf("yumi\n");
    return 0;
}