BZOJ 4810([Ynoi2017]由乃的玉米田-莫队)

Description

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

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

Input

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

用bitset+莫队。

#include<bits/stdc++.h> 
using namespace std;
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
#define
                        For(j,m-1) cout<<a[i][j]<<' ';\
                        cout<<a[i][m]<<endl; \
                        } 
#pragma
#define
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return ((a-b)%F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
    return x*f;
} 
#define
#define
bitset<100000> f,g;
int b[MAXN],res[MAXN],belong[MAXN],cnt[MAXN]={};
struct node{
    int op,l,r,x,id;
    friend bool operator<(node a,node b) {
        if (belong[a.l] ^ belong[b.l] )
            return belong[a.l] < belong[b.l];
        return a.r<b.r; 
    }
}a[MAXN];
int main()
{
//  freopen("bzoj4810.in","r",stdin);
//  freopen(".out","w",stdout);
    int n=read(),m=read();
    int BS=sqrt(n);
    For(i,n)
        b[i]=read();
    For(i,m){
        a[i]=node{read(),read(),read(),read(),i};
    }
    For(i,n)
        belong[i]=(i-1)/BS+1;
    sort(a+1,a+m+1);
    int l=0,r=0;
    For(i,m){
        while(l>a[i].l){
            l--,cnt[b[l]]++;f[b[l]]=1;g[100000-b[l]]=1;
        }while(r<a[i].r){
            r++,cnt[b[r]]++;f[b[r]]=1;g[100000-b[r]]=1;
        }while(l<a[i].l){
            cnt[b[l]]--;if(!cnt[b[l]])f[b[l]]=0,g[100000-b[l]]=0;l++;
        }while(r>a[i].r){
            cnt[b[r]]--;if(!cnt[b[r]])f[b[r]]=0,g[100000-b[r]]=0;r--;
        }
        if(a[i].op==1){
            res[a[i].id] = ((f>>a[i].x)&f).any();
        }else if(a[i].op==2){
            res[a[i].id] = ((g>>(100000-a[i].x))&f).any();
        }else{
            for(int j=1;j*j<=a[i].x;j++) if(a[i].x%j==0){
                if (f[j]==1&&f[a[i].x/j]==1) {res[a[i].id]=1;break;}
            }
            if(a[i].x==0&&f[0]) res[a[i].id]=1;
        }
    }
    For(i,m)
        puts(res[i]?"yuno":"yumi");
    return 0;
}