題目大意：給你一個長度為n的序列${a_1....a_n}$，有$m$次操作

每次操作有兩種情況：修改$a_i$的值，詢問$[l,r]$中所有子區間的異或和。

數據範圍:$n,m≤10^5$，$a_i≤1000$。

對於序列$a$，我們對每一個二進制位開一個線段樹，對於每個節點，我們存儲六個值：

$sum$：該區間內所有位的異或和。

$ans$：該區間內所有子區間異或和為1的數量。

$l_0$：該區間內以區間左端點為起點的所有區間中，異或和為0的區間數量。

$l_1$：該區間內以區間左端點為起點的所有區間中，異或和為1的區間數量。

$r_0$：該區間內以區間右端點為起點的所有區間中，異或和為0的區間數量。

$r_1$：該區間內以區間右端點為起點的所有區間中，異或和為1的區間數量。

關於pushup的過程，可以參考代碼。

然後隨便搞一搞就沒了。

 1 #include<bits/stdc++.h>
 2 #define M 400005
 3 #define L long long
 4 #define MOD 100000007
 5 using namespace std;
 6  
 7 struct node{
 8     int l[2],r[2];L sum,ans;
 9     node(){l[0]=l[1]=r[0]=r[1]=sum=ans=0;}
10     node(int
 
 x){
11         l[1]=r[1]=sum=ans=x;
12         l[0]=r[0]=x^1;
13     }
14     friend node operator +(node a,node b){
15         node c;
16         c.sum=a.sum^b.sum;
17         c.l[0]=a.l[0]+b.l[a.sum];
18         c.l[1]=a.l[1]+b.l[a.sum^1];
19         c.r[0]=a.r[b.sum]+b.r[0];
20         c.r[1]=a.r[b.sum^1
 
]+b.r[1];
21         c.ans=a.ans+b.ans+1LL*a.r[0]*b.l[1]+1LL*a.r[1]*b.l[0];
22         return c;
23     }
24 };
25 struct seg{
26     node a[M];
27     void updata(int x,int lc,int rc,int k,int op){
28         if(lc==rc) return void(a[x]=node(op));
29         int mid=(lc+rc)>>1;
30         if(k<=mid) updata(x<<1,lc,mid,k,op);
31         else updata(x<<1|1,mid+1,rc,k,op);
32         a[x]=a[x<<1]+a[x<<1|1];
33     }
34     node query(int x,int lc,int rc,int ll,int rr){
35         if(ll<=lc&&rc<=rr) return a[x];
36         int mid=(lc+rc)>>1;
37         node res;
38         if(ll<=mid) res=res+query(x<<1,lc,mid,ll,rr);
39         if(mid<rr) res=res+query(x<<1|1,mid+1,rc,ll,rr);
40         return res;
41     }
42 }p[10];
43 int n,m; 
44 int main(){
45     scanf("%d",&n);
46     for(int i=1;i<=n;i++){
47         int x; scanf("%d",&x);
48         for(int j=0;j<10;j++)
49         p[j].updata(1,1,n,i,(x>>j)&1);
50     }
51     scanf("%d",&m);
52     while(m--){
53         char op[10]; int x,y; 
54         scanf("%s%d%d",op,&x,&y);
55         if(op[0]==‘Q‘){
56             L ans=0;
57             for(int j=0;j<10;j++){
58                 node res=p[j].query(1,1,n,x,y);
59                 ans+=res.ans<<j;
60             }
61             printf("%lld\n",ans%MOD);
62         }else{
63             for(int j=0;j<10;j++)
64             p[j].updata(1,1,n,x,(y>>j)&1);
65         }
66     }
67 }

【GDKOI2018】 魔卡少女 線段樹