【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(n\sqrt{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=,q=; char c=getchar();
     while((c<'' || c>'') && c!='-') c=getchar();
     if (c=='-')  q=, c=getchar(); while (c>='' && c<='') w=w*+c-'', c=getchar();
     return q ? -w : w;
 }

 bitset<MA+>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=;i<=n;i++) a[i]=getint();
     for (llg i=;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+,ask+m+,cmp);
 }

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

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

         }

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

         if (ask[i].ty==)
         {
             for (llg k=;k*k<=ask[i].x;k++)
                 if (ask[i].x%k==)
                 {
                     if (f[k] && f[ask[i].x/k]) {ans[ask[i].num]=; break;}
                 }
         }
         if (ask[i].x== && f[]) ans[ask[i].num]=;
     }
     for (llg i=;i<=m;i++) puts(ans[i]?"yuno":"yumi");
     return ;
 }