【BZOJ4810】[Ynoi2017]由乃的玉米田 bitset+莫队
【BZOJ4810】[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
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
yumi
yuno
yuno
yumi
题解:区间询问当然是用莫队来水啊~
对于差为x,直接在bitset上搞就行了;对于和为x,我们维护一个翻转的bitset,然后转化成差为x的问题来搞;对于乘积为x,我们直接开个桶,用sqrt(x)的复杂度枚举x的约数就行了。
小号交题一次提交顺利成为status倒数第一,大号提交先TLE了两次,然后加了一堆优化后顺利成为status倒数第二,试问还有谁~~~所以千万不要拿我的代码当做标程来拍。
#include <cstdio> #include <iostream> #include <cstring> #include <bitset> #include <algorithm> #include <cmath> using namespace std; const int maxn=100010; bitset<maxn*2> b1,b2; int n,m,B; int v[maxn],st[maxn*2],ans[maxn]; struct QUERY { int qa,qb,qk,qx,org; }q[maxn]; bool cmp(QUERY a,QUERY b) { return (a.qa/B==b.qa/B)?(a.qb<b.qb):(a.qa/B<b.qa/B); } void ins(int x) { if(!st[x]) b1[x]=b2[100000-x]=1; st[x]++; } void del(int x) { st[x]--; if(!st[x]) b1[x]=b2[100000-x]=0; } int rd() { int ret=0,f=1; char gc=getchar(); while(gc<'0'||gc>'9') {if(gc=='-')f=-f; gc=getchar();} while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar(); return ret*f; } int main() { n=rd(),m=rd(); B=(int)sqrt((double)n); int i,j,l=1,r=0; for(i=1;i<=n;i++) v[i]=rd(); for(i=1;i<=m;i++) q[i].qk=rd(),q[i].qa=rd(),q[i].qb=rd(),q[i].qx=rd(),q[i].org=i; sort(q+1,q+m+1,cmp); for(i=1;i<=m;i++) { while(l>q[i].qa) ins(v[--l]); while(r<q[i].qb) ins(v[++r]); while(l<q[i].qa) del(v[l++]); while(r>q[i].qb) del(v[r--]); if(q[i].qk==1) ans[q[i].org]=((b1<<q[i].qx)&b1).any(); if(q[i].qk==2) ans[q[i].org]=((b1<<100000)&(b2<<q[i].qx)).any(); if(q[i].qk==3) for(j=1;j*j<=q[i].qx&&!ans[q[i].org];j++) if(q[i].qx%j==0&&st[j]&&st[q[i].qx/j]) ans[q[i].org]=1; } for(i=1;i<=m;i++) { if(ans[i]) printf("yuno\n"); else printf("yumi\n"); } return 0; }
| 欢迎来原网站坐坐! >原文链接<