LuoguP4462 [CQOI2018]異或序列
阿新 • • 發佈:2018-04-24
href 一個數 reg SQ IT math num 極限 解析
https://zybuluo.com/ysner/note/1124952
題面
給你一個大小為\(n\)的序列,然後給你一個數字\(k\),再給出\(m\)組詢問,詢問給出一個區間,問這個區間裏面有多少個區間的異或結果為\(k\).
\(n,m\leq10^5\)
解析
莫隊裸題。
於是我交了份傻逼代碼。(於是RE成30分)struct line { int l,r,pos; bool operator < (const line &o) const {return (l/len)==(o.l/len)?(r<o.r):(l<o.l);} }a[N]; il void add(re int x,re int f,re int ysn) { re int tot=0; if(!f) fp(i,x,R) { tot^=p[i]; if(tot==k) now+=ysn; } else fq(i,x,L) { tot^=p[i]; if(tot==k) now+=ysn; } } int main() { n=gi();m=gi();k=gi();len=sqrt(n); fp(i,1,n) p[i]=gi(); fp(i,1,m) a[i].l=gi(),a[i].r=gi(),a[i].pos=i; sort(a+1,a+1+m); L=1,R=0; fp(i,1,m) { re int l=a[i].l,r=a[i].r; while(L>l) add(--L,0,1); while(R<r) add(++R,1,1); while(L<l) add(L++,0,-1); while(R>r) add(R--,1,-1); ans[a[i].pos]=now; } fp(i,1,m) printf("%d\n",ans[i]); return 0; }
餵餵餵,極限復雜度是\(O(n^2)\)呢。。。
於是用下腦子,發現可以存一下當前統計答案的區間中每種值的數目,每加上或減去\(x\)時,相應統計\(num[k\bigoplus x]\)的貢獻即可。// luogu-judger-enable-o2 #include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> #define re register #define il inline #define ll long long #define fp(i,a,b) for(re int i=a;i<=b;i++) #define fq(i,a,b) for(re int i=a;i>=b;i--) using namespace std; const int N=5e5+100; int n,m,k,p[N],len,ans[N],now,L,R,num[N]; struct line { int l,r,pos; bool operator < (const line &o) const {return (l/len)==(o.l/len)?(r<o.r):(l<o.l);} }a[N]; il int gi() { re char ch=getchar(); re int x=0,t=1; while(ch!=‘-‘&&(ch<‘0‘||ch>‘9‘)) ch=getchar(); if(ch==‘-‘) t=-1,ch=getchar(); while(ch>=‘0‘&&ch<=‘9‘) x=x*10+ch-48,ch=getchar(); return x*t; } il void add(re int x){now+=num[k^p[x]];++num[p[x]];} il void del(re int x){--num[p[x]];now-=num[k^p[x]];} int main() { n=gi();m=gi();k=gi();len=sqrt(n); fp(i,1,n) p[i]=gi()^p[i-1]; fp(i,1,m) a[i].l=gi()-1,a[i].r=gi(),a[i].pos=i; sort(a+1,a+1+m); L=0,R=-1; fp(i,1,m) { re int l=a[i].l,r=a[i].r; while(L>l) add(--L); while(R<r) add(++R); while(L<l) del(L++); while(R>r) del(R--); ans[a[i].pos]=now; } fp(i,1,m) printf("%d\n",ans[i]); return 0; }
LuoguP4462 [CQOI2018]異或序列