bzoj#4722-由乃【倍增,抽屉原理,bitset】
生活随笔
收集整理的這篇文章主要介紹了
bzoj#4722-由乃【倍增,抽屉原理,bitset】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://darkbzoj.tk/problem/4722
題目大意
給出一個長度為nnn的序列值域為[0,v)[0,v)[0,v),要求支持操作
1≤n≤105,1≤v≤10001\leq n\leq 10^5,1\leq v\leq 10001≤n≤105,1≤v≤1000
解題思路
考慮如果我們選出了兩個有交的集合相等,那么我們把交的部分去掉就變成無交的了,所以無需考慮有沒有交。
然后根據抽屜原理對于nnn個元素所有集合總共能表示出2n2^n2n個和,但是和最大只有n×vn\times vn×v所以如果2n>n×v2^n>n\times v2n>n×v時就肯定有解,那么此時會發現當vvv最大時nnn超過131313就肯定有解。
如果元素個數小于或等于131313時我們可以先預處理出一個倍增數組加上一個樹狀數組來統計每個數最終被修改了多少次就可以得到每個數的具體值了。然后考慮dpdpdp,因為值域比較大可以用bitsetbitsetbitset優化就好了。
時間復雜度O(vlog?m+mlog?n+m13vω)O(v\log m+m\log n+m\frac{13v}{\omega})O(vlogm+mlogn+mω13v?)
code
#include<cstdio> #include<cstring> #include<algorithm> #include<bitset> #define lowbit(x) (x&-x) using namespace std; const int N=1e5+10; int n,m,v,f[1100][18],t[N],a[N]; bitset<13001>b; void Change(int x,int val){while(x<=n){t[x]+=val;x+=lowbit(x);}return; } int Ask(int x){int ans=0;while(x){ans+=t[x];x-=lowbit(x);}return ans; } int Step(int x,int b){for(int i=0;(1<<i)<=b;i++)if((b>>i)&1)x=f[x][i];return x; } int main() {scanf("%d%d%d",&n,&m,&v);for(int i=1;i<=n;i++)scanf("%d",&a[i]);for(int i=0;i<v;i++)f[i][0]=i*i*i%v;for(int j=1;(1<<j)<=n;j++)for(int i=0;i<v;i++)f[i][j]=f[f[i][j-1]][j-1];while(m--){int op,l,r;scanf("%d%d%d",&op,&l,&r);if(op==2)Change(l,1),Change(r+1,-1);else{if(r-l+1>13){puts("Yuno");continue;}b.reset();b[0]=1;bool flag=0;for(int i=l;i<=r;i++){int w=Step(a[i],Ask(i))+1;if((b&(b<<w)).any()){flag=1;break;}b=b|(b<<w);}if(flag)puts("Yuno");else puts("Yuki");}}return 0; }總結
以上是生活随笔為你收集整理的bzoj#4722-由乃【倍增,抽屉原理,bitset】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: CF835E-The penguin‘s
- 下一篇: 鸿海与 EXOLAUNCH 签署低轨卫星