1. 程式人生 > >codeforces 703D(樹狀陣列+離線處理)

codeforces 703D(樹狀陣列+離線處理)

題意: 給出n個數,m次詢問,每次詢問區間 [l,r]中出現次數為偶數次的數字的異或和

思路:

如果直接求某段區間的異或和,求出來的結果實際上是出現奇數次的數字的異或和。為了求偶數次,我們可以將出現奇數次的數字變為偶數次,將出現偶數次的數字變為奇數次。所以,找出該區間裡所有不同的數字,想象一個操作,將這些數字全部扔進原來的區間中,就可以將出現偶數次的數字變為出現奇數次。

如: 1 3 4 5 3 現在詢問區間[2,5]
區間[2,5]所有數字的異或和=3^4^5^3=1
區間[2,5]所有不同數字的異或和=3^4^5=2
區間[2,5]出現偶數次數字的異或和=1^2=3

1.求所有數字的異或和,可以維護字首異或和實現。

2.求不同數字的異或和,通過維護樹狀陣列實現。

具體實現:

離線處理:將所有詢問按照右端點從小到大排序,左端點無所謂。

記錄每個數字最後出現的位置。對於每個詢問[l,r],將前面r個數字全部更新到樹狀陣列中。如果該數字未出現過,則直接更新,如果該數字出現過,則將上次出現的數字置為0(異或該數字自己變為0),並更新這次的數字及其位置。

如:區間[2,5]={3,4,5,3} (以下指a陣列,非樹狀陣列)
①直接更新3 –>3 0 0 0
②直接更新4 –>3 4 0 0
③直接更新5 –>3 4 5 0
④3在之前已經出現過,將之前的3置為0 –> 0 4 5 3

因為已經按照右端點從小到大排序,所以如果有用到3,必定是用到最後出現的那個3,左邊的那個3被置為0不會對答案有任何影響,並且解決了重複的問題。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
#include<cstring>
#include<map>
#include<vector>
#include<set>
using namespace std;
#define Max(a,b) (a>b?a:b)
#define Min(a,b) (a<b?a:b)
#define INF ((1<<30)-1)
#define in(x) scanf("%d", &x) #define lin(x) scanf("%I64d", &x) #define out(x) printf("%d\n", x) #define lout(x) printf("%I64d\n", x) #define rep(i,a,b) for(int i=a; i<=(b); i++) #define ll __int64 #define mem(a,b) memset(a, b, sizeof(a)) #define bug(x) cout<<#x<<" = "<<x<<"; " #define N 1000005 #define nn printf("\n") int n, q; ll a[N], sum[N]; ll tree[N<<2]; map<ll,int>mp; struct A { int l, r; int no; ll ok; }ask[N]; bool cmp(A x, A y) { return x.r<y.r; } int lowbit(int k) { return k&(-k); } void update(int xb, ll num) { while(xb<=n) { tree[xb]^=num; xb+=lowbit(xb); } } //詢問前k個數的字首異或和 ll query(int k) { ll ans=0; while(k) { ans^=tree[k]; k-=lowbit(k); } return ans; } bool cmp2(A x, A y) { return x.no<y.no; } int main() { in(n); rep(i,1,n) lin(a[i]); in(q); rep(i,1,q) in(ask[i].l), in(ask[i].r), ask[i].no=i; sum[1]=a[1]; rep(i,2,n) sum[i]=sum[i-1]^a[i];//處理有重複數字的字首異或和 sort(ask+1,ask+1+q,cmp); mp.clear(); mem(tree,0); int p=1; rep(i,1,q) { int l=ask[i].l, r=ask[i].r; while(p<=r) { if(mp[a[p]]==0)//該數字未出現過 { mp[a[p]]=p;//記錄該數字當前位置 update(p,a[p]);//更新該位置為該數字 } else { update(mp[a[p]],a[p]);//更新上次位置為0 update(p,a[p]);//更新現在位置為該數字 mp[a[p]]=p;//記錄該數字當前位置 } p++; } ask[i].ok=sum[r]^sum[l-1]^query(r)^query(l-1); } sort(ask+1,ask+1+q,cmp2); rep(i,1,q) lout(ask[i].ok); return 0; }