bzoj4810: [Ynoi2017]由乃的玉米田

链接

　　http://www.lydsy.com/JudgeOnline/problem.php?id=4810

题解

　　$1A$二连击。我好像从来就没三连击过。
　　$bitset$题，$YNOI$这种直接抄原题的行为我是不想说啥。
　　大体说下思想，首先莫队，然后对于每个询问：
　　$a-b=x$可以化为$a=b+x$，就是对于区间内每个数你要判定下是否存在一个数等于它加$x$，可以直接用$bitset$，把$bitset$左移$x$位，和原来的$bitset$按位与以下，然后看看结果只要非零就说明这个查询的答案是$yes$。
　　对于$a+b=x$，就是$a=-b+x$，就是对于区间内每个数的相反数，是否存在一个数等于它加上$x$，那就维护一个相反数的$bitset$，同理搞搞。
　　对于$ab=x$，这个东西可以直接根号暴力。

代码

//bitset
#include <cstdio>
#include <algorithm>
#include <bitset>
#include <cmath>
#define maxn 100000
using namespace std;
int N, M, size, cnt[maxn+10], L[maxn+10], R[maxn+10], num[maxn+10], a[maxn+10],
    type[maxn+10], ans[maxn+10], x[maxn+10];
bitset<maxn+1> f, g;
inline int read(int x=0)
{
    char c=getchar();
    while(c<48 or c>57)c=getchar();
    while(c>=48 and c<=57)x=(x<<1)+(x<<3)+c-48, c=getchar();
    return x;
}
inline bool cmp(int a, int b)
{return L[a]/size==L[b]/size?R[a]<R[b]:L[a]/size<L[b]/size;}
void init()
{
    int i;
    N=read(), M=read(); size=sqrt(N);
    for(i=1;i<=N;i++)a[i]=read();
    for(i=1;i<=M;i++)type[i]=read(), L[i]=read(), R[i]=read(), x[i]=read(), num[i]=i;
    sort(num+1,num+M+1,cmp);
}
inline void add(int &x){if(++cnt[x]==1)f[x]=g[-x+maxn]=1;}
inline void del(int &x){if(--cnt[x]==0)f[x]=g[-x+maxn]=0;}
inline int solve(int n)
{
    int opt=type[n];
    if(opt==1)return (((f<<x[n])&f).any());
    if(opt==2)return (((g>>maxn-x[n])&f).any());
    for(int i=1;i*i<=x[n];i++)if(cnt[i] and x[n]%i==0 and cnt[x[n]/i])return true;
    return false;
}
int modui()
{
    int l=0, r=2, i, ql, qr;
    add(a[1]);
    for(i=1;i<=M;i++)
    {
        ql=L[num[i]], qr=R[num[i]];
        while(l>=ql)add(a[l]),l--;
        while(r<=qr)add(a[r]),r++;
        while(l<ql-1)l++,del(a[l]);
        while(r>qr+1)r--,del(a[r]);
        ans[num[i]]=solve(num[i]);
    }
}
int main()
{
    init();
    modui();
    for(int i=1;i<=M;i++)printf(ans[i]?"yuno\n":"yumi\n");
    return 0;
}