【BZOJ】 4810: [Ynoi2017]由乃的玉米田

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=4810

---

思路还是比较显然，第一反应应该就是莫队。

考虑怎么维护三个询问，想到了要维护每一个数字是否出现就想到了要用bitset。那么这道题就做完了。

莫队相当于是在利用bitset维护一个01序列$f$，表示每一个数字是否出现在这一段区间。

1.如果差为$x$，将这个$f$左移$x$然后判断是否和原来的$f$集合有交集即可。

2.如果和为$x$，将这个$f$翻转并左移$100000-x$然后判断是否和原来的$f$集合有交集即可。

3.如果积为$x$，直接${O(sqrt{x})}$的枚举判断即可。

复杂度${O(nsqrt{n}+frac{n^{2}}{64})}$

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#include<map>
#include<bitset>
using namespace std;
#define llg int
#define maxn 100010
#define MA 100000
#define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);

inline llg getint()
{
    llg w=0,q=0; char c=getchar();
    while((c<'0' || c>'9') && c!='-') c=getchar();
    if (c=='-')  q=1, c=getchar(); while (c>='0' && c<='9') w=w*10+c-'0', c=getchar();
    return q ? -w : w;
}

bitset<MA+1>f,g;
llg ans[maxn],n,m,a[maxn],KUAI,cnt[maxn];

struct node {llg ty,l,r,x,num;}ask[maxn];

bool cmp(const node&a,const node&b)
{
    if (a.l/KUAI==b.l/KUAI) return a.r<b.r;
    else return a.l/KUAI<b.l/KUAI;
}

void init()
{
    cin>>n>>m;
    KUAI=sqrt(n);
    for (llg i=1;i<=n;i++) a[i]=getint();
    for (llg i=1;i<=m;i++) 
        ask[i].ty=getint(),ask[i].l=getint(),ask[i].r=getint(),ask[i].x=getint(),ask[i].num=i;
    sort(ask+1,ask+m+1,cmp);
}

int main()
{
    yyj("mo");
    init();
    llg l=0,r=0;
    for (llg i=1;i<=m;i++)
    {
        while (l>ask[i].l) l--,cnt[a[l]]++,f[a[l]]=1,g[MA-a[l]]=1;
        while (r<ask[i].r) r++,cnt[a[r]]++,f[a[r]]=1,g[MA-a[r]]=1;
        while (l<ask[i].l) {cnt[a[l]]--; if (!cnt[a[l]]) f[a[l]]=0,g[MA-a[l]]=0; l++;}
        while (r>ask[i].r) {cnt[a[r]]--; if (!cnt[a[r]]) f[a[r]]=0,g[MA-a[r]]=0; r--;}

        if (ask[i].ty==1)
        {
            if (((f>>ask[i].x)&f).any()) ans[ask[i].num]=1; else ans[ask[i].num]=0;
        
        }

        if (ask[i].ty==2)
        {
            if (((g>>(MA-ask[i].x))&f).any()) ans[ask[i].num]=1; else ans[ask[i].num]=0;
        }
        
        if (ask[i].ty==3)
        {
            for (llg k=1;k*k<=ask[i].x;k++) 
                if (ask[i].x%k==0)
                {
                    if (f[k] && f[ask[i].x/k]) {ans[ask[i].num]=1; break;}
                }
        }
        if (ask[i].x==0 && f[0]) ans[ask[i].num]=1;
    }
    for (llg i=1;i<=m;i++) puts(ans[i]?"yuno":"yumi");
    return 0;
}