bzoj4810 [ YNOI2017 ] -- 莫队+bitset

最新推荐文章于 2020-11-19 20:08:17 发布

原创最新推荐文章于 2020-11-19 20:08:17 发布 · 401 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

莫队同时被 2 个专栏收录

11 篇文章

订阅专栏

bitset

3 篇文章

订阅专栏

本文介绍了一种基于莫队算法的实现方案，利用bitset进行高效的状态维护与查询操作。适用于处理区间查询问题，特别是在数据范围较大时表现优异。文章详细解析了算法原理，并提供了完整的代码示例。

莫队。维护 $2$ 个 $bitset$ ，第一个记录数字是否出现，第二个是第一个翻转后的结果。那么查询差时将第一个左移 $x$ 位再与原来的 $and$ 一下；查询和时将第一个与第二个左移 $max-x$ 的结果 $and$ 一下；查询积时枚举每个因数就可以了。
时间复杂度 $O(n\sqrt n + {n^2 \over {32}})$

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<bitset>
#include<cmath>
using namespace std;
#define N 100010
bitset<N>f,g,t;
struct Node{
    int l,r,p,f,x;
}c[N];
int i,j,k,n,m,x,a[N],b[N],s,l,r,Cnt[N],M;
bool Ans[N];
inline bool Cmp(Node x,Node y){
    return b[x.l]<b[y.l]||(b[x.l]==b[y.l]&&x.r<y.r);
}
inline void Update1(int x){
    if(++Cnt[x]==1)f[x]=g[M-x]=1;
}
inline void Update2(int x){
    if(!(--Cnt[x]))f[x]=g[M-x]=0;
}
inline bool Query(int x,int y){
    if(x==1){
        t=f&(f>>y);
        return t.count();
    }
    if(x==2){
        t=f&(g>>(M-y));
        return t.count();
    }
    s=sqrt((double)y);
    for(int i=1;i<=s;i++)
    if(!(y%i)&&f[i]&&f[y/i])return 1;
    return 0;
}
inline int Max(int x,int y){
    return x<y?y:x;
}
int main(){
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)scanf("%d",&a[i]),M=Max(M,a[i]);
    for(i=1;i<=m;i++)scanf("%d%d%d%d",&c[i].p,&c[i].l,&c[i].r,&c[i].x),c[i].f=i,M=Max(M,c[i].x);
    s=sqrt((double)n);
    for(i=1;i<=n;i++)b[i]=(i-1)/s+1;
    sort(c+1,c+n+1,Cmp);
    for(l=i=1;i<=n;i++){
        while(c[i].l<l)Update1(a[--l]);
        while(c[i].r>r)Update1(a[++r]);
        while(c[i].l>l)Update2(a[l++]);
        while(c[i].r<r)Update2(a[r--]);
        Ans[c[i].f]=Query(c[i].p,c[i].x);
    }
    for(i=1;i<=n;i++)if(Ans[i])printf("yuno\n");else printf("yumi\n");
    return 0;
}