【GDKOI2016Day1T1-魔卡少女】【拆位】线段树维护区间内所有连续子区间的异或和
题意:给出N个数,M个操作。操作有修改和询问两种,每次修改将一个数改成另一个数,每次询问一个区间的所有连续子区间的异或和。n,m<=100000,ai<=1000
题解:
当年(其实也就是今年)做不出来的题。。D1T1啊。。。
因为ai<=1000,我们可以拆位处理。拆成10个二进制位,每位开1棵线段树。
对于每个节点,维护:
d:这段区间的异或和
L[0],L[1]:子区间一定从左端点开始,异或和为0,1的子区间分别有多少个
R[0],R[1]:子区间一定从右端点开始,异或和为0,1的子区间分别有多少个
s[0],s[1]:异或和为0,1的子区间分别有多少个
然后重点就是合并啦。
node upd(int ind,int tmp,node lc,node rc)
{
int dl=lc.d,dr=rc.d;
node x;
if(tmp!=) x=t[ind][tmp];
x.d=lc.d^rc.d;
x.L[]=(lc.L[]+rc.L[(dl==) ? :])%mod;
x.L[]=(lc.L[]+rc.L[(dl==) ? :])%mod;
x.R[]=(rc.R[]+lc.R[(dr==) ? :])%mod;
x.R[]=(rc.R[]+lc.R[(dr==) ? :])%mod;
x.s[]=(lc.s[]+rc.s[]+(lc.R[]*rc.L[])%mod+(lc.R[]*rc.L[])%mod)%mod;
x.s[]=(lc.s[]+rc.s[]+(lc.R[]*rc.L[])%mod+(lc.R[]*rc.L[])%mod)%mod;
return x;
}
我打成node形式。。因为最后查询的时候有多个区间也要合并。。
代码:
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std; typedef long long LL;
const int N=;
const LL mod=;
struct node{
int l,r,lc,rc,d;
LL L[],R[],s[];
//L:从左开始
//R:从右开始
//s:总答案
}t[][*N];
char c[];
int n,m,tl,a[N][];
LL bit[]; node upd(int ind,int tmp,node lc,node rc)
{
int dl=lc.d,dr=rc.d;
node x;
if(tmp!=) x=t[ind][tmp];
x.d=lc.d^rc.d;
x.L[]=(lc.L[]+rc.L[(dl==) ? :])%mod;
x.L[]=(lc.L[]+rc.L[(dl==) ? :])%mod;
x.R[]=(rc.R[]+lc.R[(dr==) ? :])%mod;
x.R[]=(rc.R[]+lc.R[(dr==) ? :])%mod;
x.s[]=(lc.s[]+rc.s[]+(lc.R[]*rc.L[])%mod+(lc.R[]*rc.L[])%mod)%mod;
x.s[]=(lc.s[]+rc.s[]+(lc.R[]*rc.L[])%mod+(lc.R[]*rc.L[])%mod)%mod;
return x;
} int bt(int ind,int l,int r)
{
int x=++tl;
t[ind][x].l=l;t[ind][x].r=r;
t[ind][x].lc=t[ind][x].rc=;
t[ind][x].d=;
memset(t[ind][x].L,,sizeof(t[ind][x].L));
memset(t[ind][x].R,,sizeof(t[ind][x].R));
memset(t[ind][x].s,,sizeof(t[ind][x].s));
if(l<r)
{
int mid=(l+r)/;
t[ind][x].lc=bt(ind,l,mid);
t[ind][x].rc=bt(ind,mid+,r);
int lc=t[ind][x].lc,rc=t[ind][x].rc;
t[ind][x]=upd(ind,x,t[ind][lc],t[ind][rc]);
}
else
{
int d=a[l][ind];
t[ind][x].d=d;
t[ind][x].L[d]=t[ind][x].R[d]=t[ind][x].s[d]=;
}
return x;
} void change(int ind,int x,int p,int d)
{
if(t[ind][x].l==t[ind][x].r)
{
t[ind][x].d=d;
t[ind][x].L[d]=t[ind][x].R[d]=t[ind][x].s[d]=;
t[ind][x].L[d^]=t[ind][x].R[d^]=t[ind][x].s[d^]=;
return ;
}
int lc=t[ind][x].lc,rc=t[ind][x].rc,mid=(t[ind][x].l+t[ind][x].r)/;
if(p<=mid) change(ind,lc,p,d);
else change(ind,rc,p,d);
t[ind][x]=upd(ind,x,t[ind][lc],t[ind][rc]);
} node query(int ind,int x,int l,int r)
{
if(t[ind][x].l==l && t[ind][x].r==r) return t[ind][x];
int lc=t[ind][x].lc,rc=t[ind][x].rc,mid=(t[ind][x].l+t[ind][x].r)/;
if(r<=mid) return query(ind,lc,l,r);
else if(l>mid) return query(ind,rc,l,r);
else
{
node a0=query(ind,lc,l,mid);
node a1=query(ind,rc,mid+,r);
return upd(,,a0,a1);
}
} void output(int ind,int x)
{
int lc=t[ind][x].lc,rc=t[ind][x].rc;
printf("l=%d r=%d d=%d l0=%lld l1=%lld r0=%lld r1=%lld s0=%lld s1=%lld\n",t[ind][x].l,t[ind][x].r,t[ind][x].d,t[ind][x].L[],t[ind][x].L[],t[ind][x].R[],t[ind][x].R[],t[ind][x].s[],t[ind][x].s[]);
if(lc) output(ind,lc);
if(rc) output(ind,rc);
} int main()
{
freopen("a.in","r",stdin);
freopen("me.out","w",stdout);
// freopen("cardcaptor.in","r",stdin);
// freopen("cardcaptor.out","w",stdout);
scanf("%d",&n);
int x,ind;node now;
bit[]=;
for(int i=;i<=;i++) bit[i]=bit[i-]*;
memset(a,,sizeof(a));
for(int i=;i<=n;i++)
{
scanf("%d",&x);
ind=;
while(x)
{
a[i][ind]=x%;
x/=;
ind++;
}
}
scanf("%d",&m);
for(int i=;i<;i++) {tl=;bt(i,,n);}
for(int i=;i<=m;i++)
{
scanf("%s",c);
if(c[]=='Q')
{
int l,r;LL ans=;
scanf("%d%d",&l,&r);
for(int j=;j<;j++)
{
now=query(j,,l,r);
ans=(ans+(bit[j]*now.s[])%mod)%mod;
}
printf("%lld\n",ans);
}
else
{
int ind=,p,d;
scanf("%d%d",&p,&d);
while(d)
{
change(ind,,p,d%);
d/=;
ind++;
}
for(int j=ind;j<;j++) change(j,,p,);
}
}
return ;
}
【GDKOI2016Day1T1-魔卡少女】【拆位】线段树维护区间内所有连续子区间的异或和的更多相关文章
- POJ.2763 Housewife Wind ( 边权树链剖分 线段树维护区间和 )
POJ.2763 Housewife Wind ( 边权树链剖分 线段树维护区间和 ) 题意分析 给出n个点,m个询问,和当前位置pos. 先给出n-1条边,u->v以及边权w. 然后有m个询问 ...
- Can you answer these queries V SPOJ - GSS5 (分类讨论+线段树维护区间最大子段和)
recursion有一个整数序列a[n].现在recursion有m次询问,每次她想知道Max { A[i]+A[i+1]+...+A[j] ; x1 <= i <= y1 , x2 &l ...
- 线段树维护区间前k小
线段树维护区间前k小 $ solution: $ 觉得超级钢琴太麻烦?在这里线段树提供一条龙服务 . 咳咳,开始讲正题!这道题我们有一个和超级钢琴复杂度一样 $ ~O(~\sum x\times lo ...
- CodeForces - 587E[线段树+线性基+差分] ->(线段树维护区间合并线性基)
题意:给你一个数组,有两种操作,一种区间xor一个值,一个是查询区间xor的结果的种类数 做法一:对于一个给定的区间,我们可以通过求解线性基的方式求出结果的种类数,而现在只不过将其放在线树上维护区间线 ...
- hdu_5726_GCD(线段树维护区间+预处理)
题目链接:hdu_5726_GCD 题意: 给你n个数(n<=1e5)然后m个询问(m<=1e5),每个询问一个区间,问你这个区间的GCD是多少,并且输出从1到n有多少个区间的GCD和这个 ...
- FJUT3568 中二病也要敲代码(线段树维护区间连续最值)题解
题意:有一个环,有1~N编号,m次操作,将a位置的值改为b,问你这个环当前最小连续和多少(不能全取也不能不取) 思路:用线段树维护一个区间最值连续和.我们设出两个变量Lmin,Rmin,Mmin表示区 ...
- 滑动窗口(poj,线段树维护区间最值)
题目描述 现在有一堆数字共N个数字(N<=10^6),以及一个大小为k的窗口.现在这个从左边开始向右滑动,每次滑动一个单位,求出每次滑动后窗口中的最大值和最小值. 例如: The array i ...
- 51nod 1376【线段树维护区间最大值】
引自:wonter巨巨的博客 定义 dp[i] := 以数字 i(不是下标 i)为结尾的最长上升长度 然后用线段树维护 dp[i]: 每个节点维护 2 个信息,一个是当前区间的最大上升长度,一个是最大 ...
- [CSP-S模拟测试]:椎(线段树维护区间最值和单调栈)
题目描述 虽不能至,心向往之. $Treap=Tree+Heap$ 椎$=$树$+$堆 小$\pi$学习了计算机科学中的数据结构$Treap$. 小$\pi$知道$Treap$指的是一种树. 小$\p ...
随机推荐
- databales详解(一)
学习可参考:http://www.guoxk.com/node/jquery-datatables http://yuemeiqing2008-163-com.iteye.com/blog/20069 ...
- JSON解析与序列化
JSON之所以流行,拥有与JavaScript类似的语法并不是全部原因.更重要的一个原因是,可以把JSON数据结构解析为有用的 JavaScript对象.与XML数据结构要解析成DOM文档而且从中提取 ...
- MFC消息处理
1.MFC窗口如何与AfxWndProc建立联系. 当一个新的CWnd派生类创建时,在调用CWnd::CreateEx()过程中,MFC都会安装AfxCbtFilterHook().这个Hook将拦截 ...
- JDK源码分析 – ArrayList
ArrayList类的申明 ArrayList是一个支持泛型的,底层通过数组实现的一个可以存任意类型的数据结构,源码中的定义如下: public class ArrayList<E> ex ...
- C语言的世界
大家好,我是一名大一的学生,我叫陈由钧,我来自计算机系,一开始选择这门专业的时候,是出于对计算机的热爱,我喜欢计算机,喜欢没事琢磨琢磨计算的各种程序,各种软件,所以我选择学习计算机这门专业,第一周我就 ...
- 关于已部署的WCF服务升级的问题
在日常的开发过程中,我们会经常迭代发布不同的版本,所以WCF服务的接口也会经常处于变动的状态,比如在传递实体类中新加一个字段.修改参数名称等等关于服务升级的问题.但是我们不可能让已发布的版本重新引用新 ...
- ICE checkbox 用法
Hello everybody, I have a datable which contain multiple lines gotten from database, in the header o ...
- WPF对某控件添加右键属性
代码创建右键属性 ContextMenu cm = new ContextMenu(); MenuItem mi = new MenuItem(); mi.Header = "打开此文件所有 ...
- HTML5拖拽表格中单元格间的数据库
效果图:
- arc076 F - Exhausted? (霍尔定理学习)
题目链接 Problem Statement There are M chairs arranged in a line. The coordinate of the i-th chair ($$$1 ...