思路：

看完這道題根本沒有思路啊....

然后我就膜拜了一波題解...

這神tm亂搞思路

維護(hù)兩個(gè)bitset

第一個(gè)bitset代表當(dāng)前區(qū)間哪些數(shù)出現(xiàn)過

第二個(gè)bitset是 maxp-p出現(xiàn)過

差為x的時(shí)候 ?就用第一個(gè)bitset與一下它右移x就好了

和為x的時(shí)候 就第一個(gè)bitset與一下第二個(gè)bitset右移maxp-x

乘積為x的時(shí)候 就枚舉約數(shù),, 暴力判斷一下

復(fù)雜度是O(nsqrt(n)+n^2/32)的?

?

?

?

//By SiriusRen #include <cmath> #include <bitset> #include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int N=100050,M=100000; int n,m,a[N],Block,block[N],cnt[N],ans[N]; bitset<N>f,g; struct Node{int op,l,r,x,id;}node[N]; bool cmp(Node a,Node b){if(block[a.l]==block[b.l])return a.r<b.r;return a.l<b.l; } int main(){scanf("%d%d",&n,&m),Block=sqrt(n);for(int i=1;i<=n;i++)scanf("%d",&a[i]),block[i]=(i-1)/Block+1;for(int i=1;i<=m;i++)scanf("%d%d%d%d",&node[i].op,&node[i].l,&node[i].r,&node[i].x),node[i].id=i;sort(node+1,node+1+m,cmp);for(int i=1,l=1,r=0;i<=m;i++){for(;r<node[i].r;r++)cnt[a[r+1]]++,f.set(a[r+1]),g.set(M-a[r+1]);for(;l>node[i].l;l--)cnt[a[l-1]]++,f.set(a[l-1]),g.set(M-a[l-1]);for(;l<node[i].l;l++){cnt[a[l]]--;if(!cnt[a[l]])f.reset(a[l]),g.reset(M-a[l]);}for(;r>node[i].r;r--){cnt[a[r]]--;if(!cnt[a[r]])f.reset(a[r]),g.reset(M-a[r]);}if(node[i].op==1){if((f&(f>>node[i].x)).any())ans[node[i].id]=1;}else if(node[i].op==2){if((f&(g>>(M-node[i].x))).any())ans[node[i].id]=1;}else{for(int j=1;j*j<=node[i].x;j++)if(node[i].x%j==0)if(f[j]&&f[node[i].x/j]){ans[node[i].id]=1;break;}}}for(int i=1;i<=m;i++)puts(ans[i]?"yuno":"yumi"); }

?

轉(zhuǎn)載于:https://www.cnblogs.com/SiriusRen/p/6685505.html

總結(jié)

以上是生活随笔為你收集整理的BZOJ 4810 莫队+bitset的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。