bzoj4810: [Ynoi2017]由乃的玉米田

莫队。然后学了下bitset ORZ OZY巨佬

先莫队，然后乘法暴力枚举因子，减的话让bsa自己左移(右移应该也可以)x位然后就和原来的自己and一下，这样如果有1的位说明当前位和当前位+左移的位数是有的，也就是可以减出x。

然后加法可以转换为减法，具体的做法就是定上限为n+1，再定义一个bitset bsb和bsa恰好相反，那么让bsb移动（n+1-x）位和bsa and，如果有1的位说明(n+1-x)是可以得到的，又因为bsb和bsa相反，bsb的（n+1-i）位就是bsa的第i位，所以说x实质上就是得到了。

update 之前n写了m，但是好像bzoj的数据n和m是一样的。。让我狗过去了。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<bitset>
using namespace std;

int n,a[110000];

struct query
{
    int l,r,wi,id,k;
}q[110000];int m;

int block,st[110000];
bool cmp(query n1,query n2)
{
    if(st[n1.l]<st[n2.l]||
        (st[n1.l]==st[n2.l]&&st[n1.r]<st[n2.r]))
            return true;
    return false;
}

int as[110000];
int cnt[110000];
bitset<110000> bsa,bsb;
void change(int x,int w)
{
    if(w==1)
    {
        cnt[x]++;
        if(cnt[x]==1) bsa[x]=1, bsb[m+1-x]=1;
    }
    else 
    {
        cnt[x]--;
        if(cnt[x]==0) bsa[x]=0, bsb[m+1-x]=0;
    }
}
int red(int x)
{
    if( ((bsa>>x)&(bsa)).any()==true )return 1;
    else return 0;
}
int plu(int x)
{
    if( ((bsb>>(m+1-x))&(bsa)).any()==true )return 1;
    else return 0;
}
int mul(int k)
{
    for(int i=1;i*i<=k;i++)
        if(k%i==0)
            if(bsa[i]==1&&bsa[k/i]==1)
                return 1;
    return 0;
}
int solve()
{
    int l=1,r=0;
    memset(cnt,0,sizeof(cnt));
    for(int i=1;i<=m;i++)
    {
        while(l<q[i].l) change(a[l],-1), l++;
        while(l>q[i].l) l--, change(a[l],1);
        while(r<q[i].r) r++, change(a[r],1);
        while(r>q[i].r) change(a[r],-1), r--;
        
             if(q[i].wi==1)as[q[i].id]=red(q[i].k);
        else if(q[i].wi==2)as[q[i].id]=plu(q[i].k);
        else if(q[i].wi==3)as[q[i].id]=mul(q[i].k);
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    
    for(int i=1;i<=m;i++)
        scanf("%d%d%d%d",&q[i].wi,&q[i].l,&q[i].r,&q[i].k), q[i].id=i;
        
        
    block=sqrt(n);
    for(int i=1;i<=n;i++)st[i]=(i-1)/block+1;
    sort(q+1,q+m+1,cmp);
    
    solve();
    
    for(int i=1;i<=m;i++)
        if(as[i]==1)printf("yuno\n");
        else         printf("yumi\n");
    
    return 0;
}