luogu3674 小清新人渣的本愿 (bitset+莫队)

对于加减，用bitset维护当前每个数有没有

对于乘，暴力枚举约数

然后莫队

复杂度$O(m(\sqrt{n}+\frac{c}{64}))$

 #include<bits/stdc++.h>
 #define pa pair<ll,ll>
 #define CLR(a,x) memset(a,x,sizeof(a))
 #define MP make_pair
 using namespace std;
 typedef long long ll;
 const int maxn=1e5+;

 inline char gc(){
     return getchar();
     static const int maxs=<<;static char buf[maxs],*p1=buf,*p2=buf;
     return p1==p2&&(p2=(p1=buf)+fread(buf,,maxs,stdin),p1==p2)?EOF:*p1++;
 }
 inline ll rd(){
     ll x=;char c=gc();bool neg=;
     while(c<''||c>''){if(c=='-') neg=;c=gc();}
     while(c>=''&&c<='') x=(x<<)+(x<<)+c-'',c=gc();
     return neg?(~x+):x;
 }

 int N,V,NN,M,v[maxn];
 struct Node{
     int o,l,r,x,i;
 }q[maxn];
 bool ans[maxn];
 int cnt[maxn];
 bitset<> b1,b2;

 inline bool cmp(Node a,Node b){
     return a.l/NN==b.l/NN?((a.l/NN)&?a.r<b.r:a.r>b.r):a.l<b.l;
 }

 inline void solve(int p,int d){
     if(!cnt[v[p]]) b1[v[p]]=,b2[V-v[p]]=;
     cnt[v[p]]+=d;
     if(!cnt[v[p]]) b1[v[p]]=,b2[V-v[p]]=;
 }

 int main(){
     //freopen("","r",stdin);
     int i,j,k;
     N=rd(),NN=sqrt(N),M=rd();
     for(i=;i<=N;i++) v[i]=rd(),V=max(V,v[i]);
     for(i=;i<=M;i++){
         q[i].o=rd(),q[i].l=rd(),q[i].r=rd(),q[i].x=rd(),q[i].i=i;
     }sort(q+,q+M+,cmp);
     int l=,r=;
     for(i=;i<=M;i++){
         while(r<q[i].r) solve(++r,);
         while(r>q[i].r) solve(r--,-);
         while(l<q[i].l) solve(l++,-);
         while(l>q[i].l) solve(--l,);
         if(q[i].o==){
             ans[q[i].i]=(b1&(b1<<q[i].x)).count();
         }else if(q[i].o==){
             ans[q[i].i]=(b1&(q[i].x>V?(b2<<(q[i].x-V)):(b2>>(V-q[i].x)))).count();
         }else{
             for(j=;j*j<=q[i].x;j++){
                 if(q[i].x%j==&&b1[j]&&b1[q[i].x/j]){
                     ans[q[i].i]=;break;
                 }
             }
         }
     }
     for(i=;i<=M;i++){
         if(ans[i]) printf("hana\n");
         else printf("bi\n");
     }
     return ;
 }