題目鏈接：點(diǎn)擊查看

題目大意：給出一個(gè)長(zhǎng)度為 $n$ 的序列，現(xiàn)在需要執(zhí)行 $m$ 次操作，每次操作分為下列兩種類(lèi)型：

0 l r x：區(qū)間 $[l,r]$ 內(nèi)的數(shù)都異或上 $x$

1 l r x：區(qū)間 $[l,r]$ 內(nèi)的數(shù)異或上一個(gè)首項(xiàng)為 $x$ 公差為 $1$ 的等差序列，具體來(lái)說(shuō)就是 $a[l]\oplus x,a[l+1]\oplus (x+1),...,a[r]\oplus (x+r?l)$

最后需要輸出最終的序列

題目分析：第一個(gè)操作通過(guò)差分很容易維護(hù)，關(guān)鍵是第二個(gè)操作的等差數(shù)列該如何維護(hù)，因?yàn)楫惢虿僮鲗?duì)于加法操作不具有分配律，所以我們需要將區(qū)間拆解

具體而言就是，設(shè) $k=log(lowbit(x))$，對(duì)于值域區(qū)間 $[x,x+2^k)$ 來(lái)說(shuō)，因?yàn)?$x$ 轉(zhuǎn)換成二進(jìn)制后的低 $i$ 位都是 $0$，所以原本需要異或的 $x+(i?l)$ 就可以拆成 $x\oplus (i?l)$，所以可以將每次需要進(jìn)行第二種操作的區(qū)間都拆成 $log$ 個(gè)類(lèi)似的區(qū)間進(jìn)行維護(hù)，更具體的，維護(hù)一個(gè)倍增差分數(shù)組 $c[i][j]$，代表區(qū)間 $[i,i+(1<<j))$ 需要分別異或上等差數(shù)列 $\{1,2,3,...,2^j?1\}$，每次更新的時(shí)候同時(shí)維護(hù)原差分?jǐn)?shù)組和倍增差分?jǐn)?shù)組就可以了

還原的時(shí)候類(lèi)似于線段樹(shù)的pushdown，也類(lèi)似于倍增的下壓，$O(nlogn)$ 處理一下就可以得到正確的原差分?jǐn)?shù)列了

代碼：

// Problem: Kuriyama Mirai and Exclusive Or // Contest: NowCoder // URL: https://ac.nowcoder.com/acm/contest/11254/I // Memory Limit: 131072 MB // Time Limit: 6000 ms // // Powered by CP Editor (https://cpeditor.org)// #pragma GCC optimize(2) // #pragma GCC optimize("Ofast","inline","-ffast-math") // #pragma GCC target("avx,sse2,sse3,sse4,mmx") #include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cmath> #include<cstring> #include<algorithm> #include<stack> #include<climits> #include<queue> #include<map> #include<set> #include<sstream> #include<cassert> #include<bitset> #include<list> #include<unordered_map> #define lowbit(x) (x&-x) using namespace std; typedef long long LL; typedef unsigned long long ull; template<typename T> inline void read(T &x) {T f=1;x=0;char ch=getchar();while(0==isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}while(0!=isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();x*=f; } template<typename T> inline void write(T x) {if(x<0){x=~(x-1);putchar('-');}if(x>9)write(x/10);putchar(x%10+'0'); } const int inf=0x3f3f3f3f; const int N=1e6+100; int a[N],b[N],c[N][20],n,m; void add(int l,int r,int x) {if(l>r) {return;}if(l<=n) {b[l]^=x;}if(r+1<=n) {b[r+1]^=x;} } void add_tag(int x,int k) {c[x][k]^=1; } int main() { #ifndef ONLINE_JUDGE // freopen("data.in.txt","r",stdin); // freopen("data.out.txt","w",stdout); #endif // ios::sync_with_stdio(false);read(n),read(m);for(int i=1;i<=n;i++) {read(a[i]);}while(m--) {int op,l,r,x;read(op),read(l),read(r),read(x);if(op==0) {add(l,r,x);} else if(op==1) {for(int i=0;i<=19;i++) {if(((x>>i)&1)&&l+(1<<i)-1<=r) {int ll=l,rr=l+(1<<i)-1;add(ll,rr,(x>>i)<<i);add_tag(l,i);x+=(1<<i);l+=(1<<i);}}for(int i=19;i>=0;i--) {if(l+(1<<i)-1<=r) {int ll=l,rr=l+(1<<i)-1;add(ll,rr,(x>>i)<<i);add_tag(l,i);l+=(1<<i);x+=(1<<i);}}}}for(int i=19;i>=1;i--) {for(int j=1;j<=n;j++) {if(!c[j][i]) {continue;}int l=j+(1<<(i-1)),r=j+(1<<i)-1;c[j][i-1]^=1;if(l<=n) {c[l][i-1]^=1;add(l,r,1<<(i-1));}}}for(int i=1;i<=n;i++) {b[i]^=b[i-1];printf("%d ",a[i]^b[i]);}return 0; }

總結(jié)

以上是生活随笔為你收集整理的2021牛客多校3 - Kuriyama Mirai and Exclusive Or(差分+倍增)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

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