[BZOJ4810][Ynoi2017]由乃的玉米田

Time Limit: 30 Sec
Memory Limit: 256 MB

由乃在自己的农田边散步，她突然发现田里的一排玉米非常的不美。这排玉米一共有N株，它们的高度参差不齐。
由乃认为玉米田不美，所以她决定出个数据结构题

这个题是这样的：
给你一个序列a，长度为n，有m次操作，每次询问一个区间是否可以选出两个数它们的差为x，或者询问一个区间是否可以选出两个数它们的和为x，或者询问一个区间是否可以选出两个数它们的乘积为x ，这三个操作分别为操作1,2,3选出的这两个数可以是同一个位置的数

Input

第一行两个数n,m
后面一行n个数表示$a_i$
后面m行每行四个数opt, l, r, x
opt表示这个是第几种操作，l,r表示操作的区间，x表示这次操作的x
定义c为每次的x和ai中的最大值，$a_i>=0$，每次的$x>=n,m,c<=100000$

Output

对于每个询问，如果可以，输出yuno，否则输出yumi

Sample Input

5 5
1 1 2 3 4
2 1 1 2
1 1 2 2
3 1 1 1
3 5 5 16
1 2 3 4

Sample Output

yuno
yumi
yuno
yuno
yumi

题解：
没想到是分块暴力吧23333，莫队分块，对于每个块用两个bitset分别存储查询区间的元素$a[i]$的存在性(tp)和$100000-a[i]$的存在性(rtp)。对于询问3，我们直接$sqrt(x)$暴力（注意x==0的情况），对于询问2，我们询问集合tp和(rtp>>(100000-x-1))[即$x-a[i]\in [L,R]$的集合]是否有交集，对于询问2，我们询问集合tp和(tp>>x) [即$a[i]-x\in [L,R]$的集合] 是否有交集。

复杂度O(n^2)但是常数很小的哦。

#include<bits/stdc++.h>
#define LiangJiaJun main
#define N 100001
using namespace std;
int n,m;
bitset<N>rtp,tp,now;
int cnt[100004];
int belong[100004];
struct ask{
    int op,l,r,x,sit;
}q[100004];
int a[100004];
bool ans[100004];
inline bool dex(ask A,ask B){
    return belong[A.l]==belong[B.l]?A.r<B.r:belong[A.l]<belong[B.l];
}
int solve(){
    sort(q+1,q+m+1,dex);
    tp.reset();
    rtp.reset();
    int l=0,r=0;
    for(int i=1;i<=m;i++){
        while(r<q[i].r){
            ++r;
            if(!cnt[a[r]]){
                tp.set(a[r]);
                rtp.set(N-a[r]-1);
            }
            cnt[a[r]]++;
        }
        while(l>q[i].l){
            --l;
            if(!cnt[a[l]]){
                tp.set(a[l]);
                rtp.set(N-a[l]-1);
            }
            cnt[a[l]]++;
        }
        while(r>q[i].r){
            if(cnt[a[r]]==1){
                tp.reset(a[r]);
                rtp.reset(N-a[r]-1);
            }
            cnt[a[r]]--;
            r--;
        }
        while(l<q[i].l){
            if(cnt[a[l]]==1){
                tp.reset(a[l]);
                rtp.reset(N-a[l]-1);
            }
            cnt[a[l]]--;
            l++;
        }

        if(q[i].op==3){
            if(q[i].x==0&&tp[0]){
                ans[q[i].sit]=1;
                continue;
            }
            for(int j=1;j<=sqrt(q[i].x);j++){
                if(q[i].x%j!=0)continue;
                if(tp[j]&&tp[q[i].x/j]){
                    ans[q[i].sit]=1;
                    break;
                }
            }
        }
        if(q[i].op==2){
            now=(tp&(rtp>>(N-q[i].x-1)));
            if(now.any())ans[q[i].sit]=1;
        }
        if(q[i].op==1){
            now=(tp&(tp>>q[i].x));
            if(now.any())ans[q[i].sit]=1;
        }
    }
    return 0;
}
int w33ha(){
    memset(cnt,0,sizeof(cnt));
    memset(ans,0,sizeof(ans));
    int block=(int)sqrt(n);
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)belong[i]=(i-1)/block+1;
    for(int i=1;i<=m;i++){
        scanf("%d%d%d%d",&q[i].op,&q[i].l,&q[i].r,&q[i].x);
        q[i].sit=i;
    }
    solve();
    for(int i=1;i<=m;i++){
        if(!ans[i])puts("yumi");
        else puts("yuno");
    }
    return 0;
}
int LiangJiaJun(){
    while(scanf("%d%d",&n,&m)!=EOF)w33ha();
    return 0;
}