BZOJ4810 Ynoi2017由乃的玉米田（莫队+bitset）

　　多组询问不强制在线，那么考虑莫队。bitset维护当前区间出现了哪些数，数组记录每个数的出现次数以维护bitset。对于乘法，显然应有一个根号范围内的因子，暴力枚举即可。对于减法，a[i]-a[j]=x移项得a[i]-x=a[j]，可以让bitset大力右移取and。对于加法，a[i]+a[j]=x移项得a[i]=x-a[j]，维护一个翻转的bitset大力右移取and。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<bitset>
using namespace std;
#define ll long long
#define N 100010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<''||c>'')) c=getchar();return c;}
int gcd(int n,int m){return m==?n:gcd(m,n%m);}
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
bitset<N> f,g;
int n,m,a[N],cnt[N];
bool flag[N];
struct data
{
    int op,l,r,x,k,i;
    bool operator <(const data&a) const
    {
        return k<a.k||k==a.k&&(k&?r>a.r:r<a.r);
    }
}q[N];
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj4810.in","r",stdin);
    freopen("bzoj4810.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    for (int i=;i<=n;i++) a[i]=read();
    int block=sqrt(n);
    for (int i=;i<=m;i++) q[i].op=read(),q[i].l=read(),q[i].r=read(),q[i].x=read(),q[i].k=q[i].l/block,q[i].i=i;
    sort(q+,q+m+);
    int l=,r=;
    for (int i=;i<=m;i++)
    {
        while (r<q[i].r) if ((cnt[a[++r]]++)==) f[a[r]]=,g[-a[r]]=;
        while (r>q[i].r) if ((--cnt[a[r--]])==) f[a[r+]]=,g[-a[r+]]=;
        while (l>q[i].l) if ((cnt[a[--l]]++)==) f[a[l]]=,g[-a[l]]=;
        while (l<q[i].l) if ((--cnt[a[l++]])==) f[a[l-]]=,g[-a[l-]]=;
        if (q[i].op==)
        {
            for (int j=;j*j<=q[i].x;j++)
            if (q[i].x%j==&&cnt[j]&&cnt[q[i].x/j]) {flag[q[i].i]=;break;}
        }
        else if (q[i].op==) flag[q[i].i]=(f&(f>>q[i].x)).count();
        else flag[q[i].i]=(f&(g>>-q[i].x)).count();
    }
    for (int i=;i<=m;i++)
    if (flag[i]) printf("yuno\n");
    else printf("yumi\n");
    return ;
}