bzoj4810 [Ynoi2017]由乃的玉米田 bitset优化+暴力+莫队

[Ynoi2017]由乃的玉米田

Time Limit: 30 Sec  Memory Limit: 256 MB
Submit: 917  Solved: 447
[Submit][Status][Discuss]

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

HINT

 

Source

By 佚名提供

题解：一开始以为是数据结构题，然后想了好久，题目也有一个坑点，可以选两个同位置的数，就是可以选相同的数。

　　　然后对于乘法的话，根号枚举即可，对于加法减法怎么办，可以用bitset优化，因为值的数据范围不大，

　　　所以我们可以x-y=z,那么x这个bitset右移z位，判断一下，和原来有没有交即可，对于区间的话还是需要莫队一下的

　　　对于加法x+y=z，貌似可以向左移，但是这里转换了以下，换成c-y，然后，x右移 c-(c-y)=+y，所以就可以了，

　　　这样复杂度是呢n^2/32的，莫队是并列的不是嵌套。

 #include<cstring>
 #include<cmath>
 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<bitset>

 #define N 400007
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(!isdigit(ch)){if(ch=='-')f=-;ch=getchar();}
     while(isdigit(ch)){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return x*f;
 }

 int n,m,c,blo;
 int a[N],b[N],bel[N],ans[N],num[N];
 struct Node
 {
     int opt,l,r,x,id;
     friend bool operator<(Node x,Node y)
     {
         if (bel[x.l]==bel[y.l]) return x.r<y.r;
         else return bel[x.l]<bel[y.l];
     }
 }q[N];
 bitset<>f,g,h;

 void add(int x)
 {
     num[a[x]]++;
     if (num[a[x]]==)
     {
         f[a[x]]=;
         g[b[x]]=;
     }
 }
 void del(int x)
 {
     num[a[x]]--;
     if (num[a[x]]==)
     {
         f[a[x]]=;
         g[b[x]]=;
     }
 }
 bool query(Node w)
 {
     if (w.opt==)
     {
         int up=sqrt(w.x);
         for (int i=;i<=up;i++)
             if (w.x%i==) if (num[i]&&num[w.x/i]) return true;
     }
     else if (w.opt==)
     {
         h=g;
         h>>=(c-w.x);
         h&=f;
         if (h.count()) return true;
     }
     else
     {
         h=f;
         h>>=w.x;
         h&=f;
         if (h.count()) return true;
     }
     return false;
 }
 void solve_modui()
 {
     int l=,r=;
     for (int i=;i<=m;i++)
     {
         for (;r<q[i].r;r++)add(r+);
         for (;r>q[i].r;r--)del(r);
         for (;l<q[i].l;l++)del(l);
         for (;l>q[i].l;l--)add(l-);
         ans[q[i].id]=query(q[i]);
     }
     for (int i=;i<=m;i++)
         if (ans[i]) puts("yuno");
         else puts("yumi");
 }
 int main()
 {
     n=read(),m=read(),blo=*sqrt(n),c=;
     for (int i=;i<=n;i++)
         a[i]=read(),b[i]=c-a[i],bel[i]=(i-)/blo+;
     for (int i=;i<=m;i++)
         q[i].opt=read(),q[i].l=read(),q[i].r=read(),q[i].x=read(),q[i].id=i;
     sort(q+,q+m+);
     solve_modui();
 }