题解:

神题。。。我看到的时候直接吓懵了。。。

这是一道STL题。。。否则可能要写可持久化ETT或者可持久化Toptree?

用bitset来维护每个蘑菇上哪里有杂草,那么

对于操作1和操作2:可以预处理每个点为跟的bitset;

对于操作3和操作4:预处理两个点到根这条链上的bitset,先异或再或两个点的lca的bitset;

对于操作5和操作6:直接暴力区间清零即可;

询问直接贪心,由于随机所以复杂度O(玄学)(实际上应该是$O(\frac{n}{32}+\frac{n}{ln200})$?)

用STL的bitset会T掉一个点。。。所以我只能看着dalao们的手写bitset瑟瑟发抖(然而也并没有完全看懂)

代码:

 #include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#include<map>
#define inf 2147483647
#define eps 1e-9
using namespace std;
typedef long long ll;
typedef unsigned int uint;
struct edge{
int v,next;
}a[];
int n,m,u,v,op,L,now,cnt=,tot=,lg[],dep[],head[],rts[],fa[][];
uint *bt[],s1[][],s2[][];
map<uint*,int>mp;
void New(int k){
bt[++cnt]=bt[k];
mp[bt[cnt]]++;
}
void clr(int l,int r){
if(r-l<){
for(int i=l;i<=r;i++){
bt[now][i>>]&=~((uint)<<(i&));
}
return;
}
int L=(l>>)+,R=(r>>)-;
for(int i=L;i<=R;i++)bt[now][i]=;
for(int i=l;(i>>)!=L;i++)bt[now][i>>]&=~((uint)<<(i&));
for(int i=r;(i>>)!=R;i--)bt[now][i>>]&=~((uint)<<(i&));
}
void rec(int k){
if(mp[bt[k]]>){
uint *bb=bt[k];
bt[k]=(uint*)malloc(L*sizeof(uint));
mp[bt[k]]++;
for(int i=;i<L;i++)bt[k][i]=bb[i];
mp[bt[k]]--;
}
}
void Or(uint a[],uint b[]){for(int i=;i<L;i++)a[i]|=b[i];}
void And(uint a[],uint b[]){for(int i=;i<L;i++)a[i]&=b[i];}
void Nand(uint a[],uint b[]){for(int i=;i<L;i++)a[i]&=~b[i];}
void AND(uint a[],uint b[],uint c[],int k){
uint t=a[k>>]>>(k&)&;
for(int i=;i<L;i++)a[i]&=b[i]^c[i];
a[k>>]|=t<<(k&);
}
void NAND(uint a[],uint b[],uint c[],int k){
for(int i=;i<L;i++)a[i]&=~(b[i]^c[i]);
a[k>>]&=~((uint)<<(k&));
}
void add(int u,int v){
a[++tot].v=v;
a[tot].next=head[u];
head[u]=tot;
}
void dfs(int u,int ff,int dpt){
dep[u]=dpt;
fa[u][]=ff;
for(int i=;i<=;i++)fa[u][i]=fa[fa[u][i-]][i-];
for(int i=;i<L;i++)s1[u][i]=s1[ff][i];
s1[u][u>>]^=(uint)<<(u&);
for(int tmp=head[u];tmp!=-;tmp=a[tmp].next){
int v=a[tmp].v;
if(v!=ff){
dfs(v,u,dpt+);
}
}
s2[u][u>>]^=(uint)<<(u&);
for(int i=;i<L;i++)s2[ff][i]|=s2[u][i];
}
int lca(int u,int v){
if(dep[u]<dep[v])swap(u,v);
int l=dep[u]-dep[v];
for(int i=;i>=;i--){
if((<<i)&l)u=fa[u][i];
}
if(u==v)return u;
for(int i=;i>=;i--){
if(fa[u][i]!=fa[v][i]){
u=fa[u][i],v=fa[v][i];
}
}
return fa[u][];
}
int query(int k){
int ans=;
for(int i=;i<=n;){
if(bt[now][i>>]>>(i&)){
uint t=(bt[now][i>>]>>(i&(uint)))<<(i&),tt=t&-t;
//printf("%lld %lld %lld\n",(ll)bt[now][i>>5],(ll)t,(ll)tt);
if(tt>)tt=lg[tt>>]+;
else tt=lg[tt];
ans++;
i=((i&(~(uint)))|tt)+k+;
}else i=((i>>)+)<<;
}
return ans;
}
int main(){
memset(head,-,sizeof(head));
scanf("%d",&n);
for(int i=;i<n;i++){
scanf("%d%d",&u,&v);
add(u,v);
add(v,u);
}
L=(n+)>>;
for(int i=;i<=;i++)lg[<<i]=i;
bt[]=(uint*)malloc(L*sizeof(uint));
bt[]=(uint*)malloc(L*sizeof(uint));
for(int i=;i<L;i++)bt[][i]=bt[][i]=;
for(int i=;i<=n;i++)bt[][i>>]|=(uint)<<(i&);
dfs(,,);
scanf("%d",&m);
for(int i=;i<=m;i++){
scanf("%d%d",&op,&now);
switch(op){
case :{New(now);break;}
case :{
scanf("%d",&u);
rec(now);
Or(bt[now],bt[u]);
break;
}
case :{
scanf("%d",&u);
rec(now);
Nand(bt[now],s2[u]);
break;
}
case :{
scanf("%d",&u);
rec(now);
And(bt[now],s2[u]);
break;
}
case :{
scanf("%d%d",&u,&v);
rec(now);
NAND(bt[now],s1[u],s1[v],lca(u,v));
break;
}
case :{
scanf("%d%d",&u,&v);
rec(now);
AND(bt[now],s1[u],s1[v],lca(u,v));
break;
}
case :{
scanf("%d%d",&u,&v);
rec(now);
clr(u,v);
break;
}
case :{
scanf("%d%d",&u,&v);
rec(now);
clr(,u-),clr(v+,n);
break;
}
case :{
scanf("%d",&u);
printf("%d\n",query(u));
break;
}
}
}
return ;
}

(2016北京集训十三)【xsy1533】mushroom - bitset的更多相关文章

  1. (2016北京集训十三)【xsy1531】魔法游戏 - Nim游戏

    题解: 好题!我的结论很接近正解了... 把一个数化成二进制,每次至少要拿走一位,最多全拿走,不能不拿.那么这就是一个经典的Nim问题了,子树异或起来就是根节点的答案,随便递推一下就行了. 代码: # ...

  2. (2016北京集训十三)【xsy1532】网络战争 - 最小割树+树上倍增+KD树

    题解: 好题!! 这题似乎能上我代码长度记录的前五? 调试时间长度应该也能上前五QAQ 首先题目要求的明显就是最小割,当然在整个森林上求Q次最小割肯定是会GG的,所以我们需要一个能快速求最小割的算法— ...

  3. (2016北京集训十)【xsy1528】azelso - 概率期望dp

    北京集训的题都是好题啊~~(于是我爆0了) 注意到一个重要的性质就是期望是线性的,也就是说每一段的期望步数可以直接加起来,那么dp求出每一段的期望就行了... 设$f_i$表示从$i$出发不回到$i$ ...

  4. 【2016北京集训】Mushroom

    Portal --> broken qwq Description 一开始有个蘑菇,蘑菇里面有\(n\)个房间,是一棵有根树,\(1\)号是根,每个房间里面都有杂草,现在要支持以下操作:将某个指 ...

  5. 2016北京集训测试赛(十三) Problem B: 网络战争

    Solution KD tree + 最小割树

  6. [2016北京集训试题6]mushroom-[bitset]

    Description Solution bitset是个好东西啊..强行压位什么的真是够orz. 由于所有的蘑菇上房间的长相是一样的,我们针对每个房间,算出它到根节点的bitset和以它为根的子树的 ...

  7. 【2016北京集训测试赛(十)】 Azelso (期望DP)

    Time Limit: 1000 ms   Memory Limit: 256 MB Description 题解 状态表示: 这题的状态表示有点难想...... 设$f_i$表示第$i$个事件经过之 ...

  8. 【2016北京集训测试赛(二)】 thr (树形DP)

    Description 题解 (这可是一道很早就碰到的练习题然后我不会做不想做,没想到在Contest碰到欲哭无泪......) 题目大意是寻找三点对的个数,使得其中的三个点两两距离都为d. 问题在于 ...

  9. 【2016北京集训测试赛(八)】 crash的数列 (思考题)

    Description 题解 题目说这是一个具有神奇特性的数列!这句话是非常有用的因为我们发现,如果套着这个数列的定义再从原数列引出一个新数列,它居然还是一样的...... 于是我们就想到了能不能用多 ...

随机推荐

  1. 修改properties文件后系统运行异常

    今天修改了项目的properties配置文件以后,运行会报异常,即使将内容改回,异常仍然存在.中间还会出现项目报错等问题,现将解决方法整理出来. 1.修改properties的打开方式,将打开方式从p ...

  2. Unity坐标系 左手坐标系 图

    x轴:从左指向右 y轴:从下指向上 z轴:指向屏幕里的是左手坐标系,指向屏幕外的是右手坐标系 记忆小技巧:都是X轴朝右,Y轴向上,跟平时画坐标一模一样,区别只是Z的朝向.你用手试一下就知道了,当大拇指 ...

  3. BZOJ 4229: 选择 LCT_独创方法_边双

    考虑如果两点在一个环中,那么这两点一定可以构出双联通分量. 考虑环和环镶嵌,那么两个环中的点一定都互为双联通分量. 由此,我们想到一个算法: 将删边转为反向加边,用LCT维护图. 当我们连接两个点时, ...

  4. 路飞学城Python-Day30

    11-僵尸进程与孤儿进程 现象:运行程序会产生父进程,在父进程中开子进程,这两个进程公用一个打印终端,运行的时候就只运行父进程,父进程虽然自己结束了,但是要等子进程结束完才会结束. 父进程可以开多个子 ...

  5. Use gdb attach pid and debug it

  6. BZOJ 4016 [FJOI2014]最短路径树问题 (贪心+点分治)

    题目大意:略 传送门 硬是把两个题拼到了一起= = $dijkstra$搜出单源最短路,然后$dfs$建树,如果$dis_{v}=dis_{u}+e.val$,说明这条边在最短路图内,然后像$NOIP ...

  7. ExcelToHtmlTable转换算法:将Excel转换成Html表格并展示(项目源码+详细注释+项目截图)

    功能概述 Excel2HtmlTable的主要功能就是把Excel的内容以表格的方式,展现在页面中.Excel的多个Sheet对应页面的多个Tab选项卡.转换算法的难点在于,如何处理行列合并,将Exc ...

  8. Golang-import-introduce

    本文主要讲解golang中import关键字的用法 import( "fmt" ) //然后在代码里面可以通过如下的方式调用 fmt.Println("hello wor ...

  9. STM32CubeMx的使用分享

    1. 新建立工程(以F103ZET6为例) 2. 配置引脚(以PA0为例)   3. 配置外设(以串口为例) 4. 配置时钟 5. 外设.GPIO.中断初始化 6. 生成工程 7. 添加自己的代码 8 ...

  10. 【BZOJ 1305】[CQOI2009]dance跳舞

    [链接] 我是链接,点我呀:) [题意] 在这里输入题意 [题解] 男生和女生每个人都分身成两个节点 即x[1],x[2]和y[1],y[2] 然后如果i和j不相互喜欢 那么add(x[i][2],y ...