BZOJ4810: [Ynoi2017]由乃的玉米田 莫队 bitset

题意：一个序列（长度100000），三种询问：区间内是否存在两个数和为x、差为x、积为x
莫队，用bitset维护当前存在哪些数，还得维护一个反转的bitset。和为x用反转bitset右移后和正bitset取交判定，差为x用正bitset右移后与正bitset取交判定，积为x直接暴力分解因数检查是否存在，总复杂度n*sqrt(n)+n*n/32
这复杂度怎么过去的。。。。BZ上各种不科学复杂度的题，总在考试时给人一种某某做法能卡过去的错觉

#include<cstdio>
#include<cmath>
#include<bitset>
#include<algorithm>
#define gm 100004
using namespace std;
struct Istream
{
    static const size_t str=1<<16;
    char buf[str],*s,*t;
    Istream():buf(),s(),t(){}
    inline char get()
    {
        return (s==t)?(t=buf+fread(s=buf,1,str,stdin),*s++):(*s++);
    }
    inline Istream& operator>>(int &x)
    {
        register char c;
        do c=get(); while(c<'0'||c>'9');
        x=0;
        while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-'0',c=get();
        return *this;
    }
}cin;
int n,m;
bitset<gm> f,g;
int a[gm];
int size;
struct query
{
    int opt,l,r,x;
    bool *ptr;
    friend Istream& operator>>(Istream &in,query &v)
    {
        return in>>v.opt>>v.l>>v.r>>v.x;
    }
    bool operator < (const query &b) const
    {
        return (l-1)/size!=(b.l-1)/size?(l-1)/size<(b.l-1)/size:r<b.r;
    }
}q[gm];
int cnt[gm];
bool ans[gm];
int l=1,r=0;
inline void push(int x)
{
    if(++cnt[x]==1) f.set(x),g.set(100000-x);
}
inline void pop(int x)
{
    if(--cnt[x]==0) f.reset(x),g.reset(100000-x);
}
int main()
{
    cin>>n>>m;
    size=ceil(sqrt(n));
    for(int i=1;i<=n;++i)
    cin>>a[i];
    for(int i=1;i<=m;++i)
    {
        cin>>q[i];
        q[i].ptr=ans+i;
    }
    sort(q+1,q+m+1);
    for(int i=1;i<=m;++i)
    {
        int u=q[i].l,v=q[i].r;
        while(u<l) push(a[--l]);
        while(r<v) push(a[++r]);
        while(l<u) pop(a[l++]);
        while(v<r) pop(a[r--]);
        switch(q[i].opt)
        {
            case 1:
                *q[i].ptr=((f>>q[i].x)&f).any();
                break;
            case 2:
                *q[i].ptr=((g>>100000-q[i].x)&f).any();
                break;
            default:
                int x=q[i].x,m=sqrt(x);
                for(int j=1;j<=m;++j)
                if(!(x%j)&&(f.test(j)&&f.test(x/j)))
                {
                    *q[i].ptr=1;
                    break;
                }
        }
    }
    for(int i=1;i<=m;++i)
    puts(ans[i]?"yuno":"yumi");
    return 0;
}