bzoj千题计划245:bzoj1095: [ZJOI2007]Hide 捉迷藏
http://www.lydsy.com/JudgeOnline/problem.php?id=1095
查询最远点对,带修改
显然可以用动态点分治
对于每个点,维护两个堆
堆q1[x] 维护 点分树x的子树中,所有黑点到x的点分树中父节点的距离
堆q2[x]维护 点分树x的子节点的堆q1的堆顶,即若y是x在点分树中的子节点,则q2[x].push(q1[y].top())
再来维护一个全局的堆Q,维护所有q2的堆顶,即Q.push(q2[x].top())
#include<cmath>
#include<queue>
#include<cstdio>
#include<iostream> using namespace std; #define N 100001 #define Swap(a,b) ( (a)^=(b),(b)^=(a),(a)^=(b) ) int n; int front[N],to[N<<],nxt[N<<],tot; int fa[N]; bool light[N];
int sum; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} void add(int u,int v)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot;
to[++tot]=u; nxt[tot]=front[v]; front[v]=tot;
} namespace LCA
{
int fa[N][],dep[N],bin[],lim; void dfs(int x)
{
for(int i=front[x];i;i=nxt[i])
if(to[i]!=fa[x][])
{
dep[to[i]]=dep[x]+;
fa[to[i]][]=x;
dfs(to[i]);
}
} int query(int x,int y)
{
if(dep[x]<dep[y]) swap(x,y);
for(int i=lim;i>=;--i)
if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
if(x==y) return x;
for(int i=lim;i>=;--i)
if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
return fa[x][];
} int dist(int x,int y)
{
return dep[x]+dep[y]-(dep[query(x,y)]<<);
} void main()
{
lim=log(n)/log();
bin[]=;
for(int i=;i<;++i) bin[i]=bin[i-]<<;
dep[]=;
dfs();
for(int i=;i<=lim;++i)
for(int j=;j<=n;++j)
fa[j][i]=fa[fa[j][i-]][i-];
}
} struct HEAP
{
priority_queue<int>heap,trash; void Pop(int x) { trash.push(x); } void Push(int x) { heap.push(x); } int Size() { return heap.size()-trash.size(); } int Top()
{
if(Size())
{
while(!trash.empty() && heap.top()==trash.top()) heap.pop(),trash.pop();
return heap.top();
}
return ;
} int Top_Sec()
{
if(Size()>=)
{
int first,second;
while(!trash.empty() && heap.top()==trash.top()) heap.pop(),trash.pop();
first=heap.top(); heap.pop();
while(!trash.empty() && heap.top()==trash.top()) heap.pop(),trash.pop();
second=heap.top(); heap.push(first);
return first+second;
}
else return Top();
} }q1[N],q2[N],Q; namespace Point_divide
{
int siz[N],mx[N];
bool vis[N]; int root,min_size; void get_dist(int x,int pa,int fa)
{
q1[root].Push(LCA::dist(pa,x));
for(int i=front[x];i;i=nxt[i])
if(to[i]!=fa && !vis[to[i]]) get_dist(to[i],pa,x);
} void get_size(int x,int fa)
{
siz[x]=; mx[x]=;
for(int i=front[x];i;i=nxt[i])
if(!vis[to[i]] && to[i]!=fa)
{
get_size(to[i],x);
siz[x]+=siz[to[i]];
if(siz[to[i]]>mx[x]) mx[x]=siz[to[i]];
}
} void get_root(int x,int pa,int fa)
{
mx[x]=max(mx[x],siz[pa]-siz[x]);
if(mx[x]<min_size) min_size=mx[x],root=x;
for(int i=front[x];i;i=nxt[i])
if(to[i]!=fa && !vis[to[i]]) get_root(to[i],pa,x);
} void work(int x,int pa)
{
min_size=n+;
get_size(x,);
get_root(x,x,);
fa[root]=pa;
vis[root]=true;
q2[root].Push();
q1[root].Push(LCA::dist(pa,root));
for(int i=front[root];i;i=nxt[i])
if(!vis[to[i]]) get_dist(to[i],pa,root);
q2[pa].Push(q1[root].Top());
int rt=root;
for(int i=front[root];i;i=nxt[i])
if(!vis[to[i]]) work(to[i],rt);
if(q2[rt].Size()>=) Q.Push(q2[rt].Top_Sec());
} void main()
{
work(,);
} } void turn_off(int x)
{
if(q2[x].Size()>=) Q.Pop(q2[x].Top_Sec());
q2[x].Push();
if(q2[x].Size()>=) Q.Push(q2[x].Top_Sec());
for(int u=x;fa[u];u=fa[u])
{
if(q2[fa[u]].Size()>=) Q.Pop(q2[fa[u]].Top_Sec());
if(q1[u].Size()) q2[fa[u]].Pop(q1[u].Top());
q1[u].Push(LCA::dist(fa[u],x));
q2[fa[u]].Push(q1[u].Top());
if(q2[fa[u]].Size()>=) Q.Push(q2[fa[u]].Top_Sec());
}
} void turn_on(int x)
{
if(q2[x].Size()>=) Q.Pop(q2[x].Top_Sec());
q2[x].Pop();
if(q2[x].Size()>=) Q.Push(q2[x].Top_Sec());
for(int u=x;fa[u];u=fa[u])
{
if(q2[fa[u]].Size()>=) Q.Pop(q2[fa[u]].Top_Sec());
q2[fa[u]].Pop(q1[u].Top());
q1[u].Pop(LCA::dist(x,fa[u]));
if(q1[u].Size()) q2[fa[u]].Push(q1[u].Top());
if(q2[fa[u]].Size()>=) Q.Push(q2[fa[u]].Top_Sec());
}
} int main()
{
//freopen("hide.in","r",stdin);
//freopen("hide.out","w",stdout);
read(n);
int u,v;
for(int i=;i<n;++i)
{
read(u); read(v);
add(u,v);
}
LCA::main();
Point_divide::main();
sum=n;
int m; char s[];
read(m);
// printf("%d\n",Q.Top());
while(m--)
{
scanf("%s",s);
if(s[]=='G')
{
if(sum>=) printf("%d\n",Q.Top());
else printf("%d\n",sum-);
}
else
{
read(u);
if(light[u])turn_off(u),sum++;
else turn_on(u),sum--;
light[u]^=;
//printf("%d\n",Q.Top());
}
}
return ;
}
bzoj千题计划245:bzoj1095: [ZJOI2007]Hide 捉迷藏的更多相关文章
- 动态点分治:Bzoj1095: [ZJOI2007]Hide 捉迷藏
简介 这是我自己的一点理解,可能写的不好 点分治都学过吧.. 点分治每次找重心把树重新按重心的深度重建成了一棵新的树,称为分治树 这个树最多有log层... 动态点分治:记录下每个重心的上一层重心,这 ...
- [bzoj1095][ZJOI2007]Hide 捉迷藏 点分树,动态点分治
[bzoj1095][ZJOI2007]Hide 捉迷藏 2015年4月20日7,8876 Description 捉迷藏 Jiajia和Wind是一对恩爱的夫妻,并且他们有很多孩子.某天,Jiaji ...
- bzoj千题计划300:bzoj4823: [Cqoi2017]老C的方块
http://www.lydsy.com/JudgeOnline/problem.php?id=4823 讨厌的形状就是四联通图 且左右各连一个方块 那么破坏所有满足条件的四联通就好了 按上图方式染色 ...
- bzoj千题计划252:bzoj1095: [ZJOI2007]Hide 捉迷藏
http://www.lydsy.com/JudgeOnline/problem.php?id=1095 点分树+堆 请去看 http://www.cnblogs.com/TheRoadToTheGo ...
- bzoj千题计划163:bzoj1060: [ZJOI2007]时态同步
http://www.lydsy.com/JudgeOnline/problem.php?id=1060 以激发器所在节点为根 终止节点一定是叶节点 记录点的子树内最深的终止节点 然后从根往下使用道具 ...
- bzoj千题计划196:bzoj4826: [Hnoi2017]影魔
http://www.lydsy.com/JudgeOnline/problem.php?id=4826 吐槽一下bzoj这道题的排版是真丑... 我还是粘洛谷的题面吧... 提供p1的攻击力:i,j ...
- bzoj千题计划280:bzoj4592: [Shoi2015]脑洞治疗仪
http://www.lydsy.com/JudgeOnline/problem.php?id=4592 注意操作1 先挖再补,就是补的范围可以包含挖的范围 SHOI2015 的题 略水啊(逃) #i ...
- bzoj千题计划177:bzoj1858: [Scoi2010]序列操作
http://www.lydsy.com/JudgeOnline/problem.php?id=1858 2018 自己写的第1题,一遍过 ^_^ 元旦快乐 #include<cstdio> ...
- bzoj千题计划317:bzoj4650: [Noi2016]优秀的拆分(后缀数组+差分)
https://www.lydsy.com/JudgeOnline/problem.php?id=4650 如果能够预处理出 suf[i] 以i结尾的形式为AA的子串个数 pre[i] 以i开头的形式 ...
随机推荐
- 【CF245H】Queries for Number of Palindromes(回文树)
[CF245H]Queries for Number of Palindromes(回文树) 题面 洛谷 题解 回文树,很类似原来一道后缀自动机的题目 后缀自动机那道题 看到\(n\)的范围很小,但是 ...
- 【Luogu1876】开灯(数论)
[Luogu1876]开灯(数论) 题面 题目描述 首先所有的灯都是关的(注意是关!),编号为1的人走过来,把是一的倍数的灯全部打开,编号为二的的把是二的倍数的灯全部关上,编号为3的人又把是三的倍数的 ...
- [Luogu3066][USACO12DEC]逃跑的BarnRunning Away From…
题面 题目描述 给出以1号点为根的一棵有根树,问每个点的子树中与它距离小于等于l的点有多少个. 输入格式: Line 1: 2 integers, N and L (1 <= N <= 2 ...
- [HNOI2010]CHORUS 合唱队
题面 Bzoj Sol 简单\(DP\) # include <bits/stdc++.h> # define IL inline # define RG register # defin ...
- 使用 vscode将本地项目上传到github以及删除github上的某个文件夹
安装Git后,可以看到windows环境下有两个命令输入窗口Git CMD 和Git Bash Git GUI是可视化图形界面 Git中的Bash是基于CMD的,在CMD的基础上增添一些新的命令与功能 ...
- 关系型数据库工作原理-数据库查询器(翻译自Coding-Geek文章)
本文翻译自Coding-Geek文章:< How does a relational database work>.原文链接:http://coding-geek.com/how-data ...
- 学习Javascript闭包(Closure)及几个经典面试题理解
今天遇到一个面试题,结果让我百思不得其解.后来在查阅了各种文档后,理清了来龙去脉.让我们先来看看这道题: function Foo( ){ var i = 0; return function( ){ ...
- C# Redis实战(七)
七.修改数据 在上一篇 C# Redis实战(六)中介绍了如何查询Redis中数据,本篇将介绍如何修改Redis中相关数据.大家都知道Redis是key-value型存储系统,所以应该可以修改key, ...
- 灵感&小思路
1.被处理的数据量级比较大.没有明显规律可循,可以分割.或者提取文件,采用分步和保存思想去解决. 比如在使用python的正则re模块时候,根据被处理对象提取属性抽象一个基类,方法实现分步,并且把中间 ...
- redis缓存的应用详解
在现在的很多项目,基本上都需要引入缓存机制,那么缓存到底是什么呢? 缓存 也就是数据交互的缓冲区 Cache 在java-web项目中实现缓存,也就是需要首先把数据库需要用到的数据备份一份作为副本 ...