简介

这是我自己的一点理解,可能写的不好

点分治都学过吧。。

点分治每次找重心把树重新按重心的深度重建成了一棵新的树,称为分治树

这个树最多有log层。。。

动态点分治:记录下每个重心的上一层重心,这棵分治树就确定了

修改就暴力在分治树中向上改,反正是log的

至于为什么叫动态点分治我不知道。。。我觉得就是点分治

做题时最主要的难点不在点分治,在于维护什么和怎样维护

例题

Bzoj1095: [ZJOI2007]Hide 捉迷藏

先搞出这个分治树,然后基本和点分治无关了

以下基于分治树

每个重心开两个堆

第一个堆记录子树中所有节点到重心的距离

第二个堆记录所有子节点的第一个堆的堆顶

那么一个节点的第二个堆堆中的最大值和次大值加起来就是子树中经过这个节点的最长链

然后开一个全局的堆,记录所有第二个堆中最大值和次大值之和

堆顶就是答案

修改就分治树中暴跳重心,大力讨论一番

代码和思路来自hzwer

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
const int _(2e5 + 10); IL ll Read(){
RG ll x = 0, z = 1; RG char c = getchar();
for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
return x * z;
} int n, nxt[_], to[_], fst[_], cnt, deep[_], size[_], frt[_], rt, mx[_], sz;
int eul[20][_], pw[20] = {1}, lg[_], id[_], tot;
bool cls[_], vis[_];
struct Heap{
priority_queue <int> A, B;
IL void Push(RG int x){ A.push(x); } IL void Del(RG int x){ B.push(x); } IL void Pop(){ while(!B.empty() && A.top() == B.top()) A.pop(), B.pop(); A.pop(); } IL int Top(){ while(!B.empty() && A.top() == B.top()) A.pop(), B.pop(); return !A.empty() ? A.top() : 0; } IL int Size(){ return A.size() - B.size(); } IL int _Top(){ if(Size() < 2) return 0; RG int x = Top(); Pop(); RG int y = Top(); Push(x); return y; }
} A, B[_], C[_]; IL void Add(RG int u, RG int v){ to[cnt] = v; nxt[cnt] = fst[u]; fst[u] = cnt++; }
/***********************************建立分治树***************************************/
IL void Getroot(RG int u, RG int ff){
size[u] = 1; mx[u] = 0;
for(RG int e = fst[u]; e != -1; e = nxt[e]){
if(to[e] == ff || vis[to[e]]) continue;
Getroot(to[e], u);
size[u] += size[to[e]];
mx[u] = max(mx[u], size[to[e]]);
}
mx[u] = max(mx[u], sz - size[u]);
if(mx[u] < mx[rt]) rt = u;
} IL void Create(RG int u, RG int ff){
frt[u] = ff; /*记录上层重心*/ vis[u] = 1;
for(RG int e = fst[u]; e != -1; e = nxt[e]){
if(vis[to[e]]) continue;
rt = 0; sz = size[to[e]];
Getroot(to[e], 0);
Create(rt, u);
}
}
/***********************************************************************************/
IL void Dfs(RG int u, RG int ff){
eul[0][++cnt] = deep[u]; id[u] = cnt;
for(RG int e = fst[u]; e != -1; e = nxt[e]){
if(to[e] == ff) continue;
deep[to[e]] = deep[u] + 1;
Dfs(to[e], u);
eul[0][++cnt] = deep[u];
}
} IL int Query(RG int x, RG int y){
x = id[x]; y = id[y]; if(x > y) swap(x, y); RG int len = y - x + 1;
return min(eul[lg[len]][x], eul[lg[len]][y - pw[lg[len]] + 1]);
} IL int MaxDis(RG int x, RG int y){ return deep[x] + deep[y] - 2 * Query(x, y); } IL void Close(RG int x, RG int y){
if(x == y){
B[x].Push(0);
if(B[x].Size() == 2) A.Push(B[x].Top());
}
if(!frt[x]) return;
RG int ff = frt[x], dis = MaxDis(ff, y), tmp = C[x].Top(); C[x].Push(dis);
if(dis > tmp){
RG int mmx = B[ff].Top() + B[ff]._Top(), sszz = B[ff].Size();
if(tmp) B[ff].Del(tmp); B[ff].Push(dis);
RG int _mmx = B[ff].Top() + B[ff]._Top();
if(_mmx > mmx){
if(sszz >= 2) A.Del(mmx);
if(B[ff].Size() >= 2) A.Push(_mmx);
}
}
Close(ff, y);
} IL void Open(RG int x, RG int y){
if(x == y){
if(B[x].Size() == 2) A.Del(B[x].Top());
B[x].Del(0);
}
if(!frt[x]) return;
RG int ff = frt[x], dis = MaxDis(ff, y), tmp = C[x].Top(); C[x].Del(dis);
if(dis == tmp){
RG int mmx = B[ff].Top() + B[ff]._Top(), sszz = B[ff].Size();
B[ff].Del(dis); if(C[x].Top()) B[ff].Push(C[x].Top());
RG int _mmx = B[ff].Top() + B[ff]._Top();
if(_mmx < mmx){
if(sszz >= 2) A.Del(mmx);
if(B[ff].Size() >= 2) A.Push(_mmx);
}
}
Open(ff, y);
} int main(RG int argc, RG char* argv[]){
Fill(fst, -1); tot = sz = n = Read(); mx[0] = 2e9;
for(RG int i = 2; i < _; ++i) lg[i] = lg[i >> 1] + 1;
for(RG int i = 1; i < 20; ++i) pw[i] = pw[i - 1] << 1;
for(RG int i = 1, u, v; i < n; i++) u = Read(), v = Read(), Add(u, v), Add(v, u);
cnt = 0; Dfs(1, 0);
for(RG int i = 1; i <= lg[cnt]; ++i)
for(RG int j = 1; j + pw[i] - 1 <= cnt; ++j)
eul[i][j] = min(eul[i - 1][j], eul[i - 1][j + pw[i - 1]]);
Getroot(1, 0); Create(rt, 0);
for(RG int i = 1; i <= n; ++i) C[i].Push(0), cls[i] = 1, Close(i, i);
for(RG int Q = Read(); Q; --Q){
RG char op; RG int x; scanf(" %c", &op);
if(op == 'G'){
if(tot <= 1) printf("%d\n", tot - 1);
else printf("%d\n", A.Top());
}
else{
x = Read();
if(cls[x]) Open(x, x), --tot, cls[x] = 0;
else Close(x, x), ++tot, cls[x] = 1;
}
}
return 0;
}

动态点分治:Bzoj1095: [ZJOI2007]Hide 捉迷藏的更多相关文章

  1. [bzoj1095][ZJOI2007]Hide 捉迷藏 点分树,动态点分治

    [bzoj1095][ZJOI2007]Hide 捉迷藏 2015年4月20日7,8876 Description 捉迷藏 Jiajia和Wind是一对恩爱的夫妻,并且他们有很多孩子.某天,Jiaji ...

  2. BZOJ1095: [ZJOI2007]Hide 捉迷藏【动态点分治】

    Description 捉迷藏 Jiajia和Wind是一对恩爱的夫妻,并且他们有很多孩子.某天,Jiajia.Wind和孩子们决定在家里玩 捉迷藏游戏.他们的家很大且构造很奇特,由N个屋子和N-1条 ...

  3. BZOJ1095 [ZJOI2007]Hide 捉迷藏 动态点分治 堆

    原文链接https://www.cnblogs.com/zhouzhendong/p/BZOJ1095.html 题目传送门 - BZOJ1095 题意 有 N 个点,每一个点是黑色或者白色,一开始所 ...

  4. BZOJ1095 [ZJOI2007]Hide 捉迷藏 【动态点分治 + 堆】

    题目链接 BZOJ1095 题解 传说中的动态点分治,一直不敢碰 今日一会,感觉其实并不艰涩难懂 考虑没有修改,如果不用树形dp的话,就得点分治 对于每个重心,我们会考虑其分治的子树内所有点到它的距离 ...

  5. 2019.01.10 bzoj1095: [ZJOI2007]Hide 捉迷藏(动态点分治)

    传送门 蒟蒻真正意义上做的第一道动态点分治! 题意:给一棵最开始所有点都是黑点的树,支持把点的颜色变成从黑/白色变成白/黑色,问当前状态树上两个最远黑点的距离. 思路: 首先考虑不带修改一次点分治怎么 ...

  6. bzoj1095: [ZJOI2007]Hide 捉迷藏 动态点分治学习

    好迷啊...感觉动态点分治就是个玄学,蜜汁把树的深度缩到logn (静态)点分治大概是递归的时候分类讨论: 1.答案经过当前点,暴力(雾)算 2.答案不经过当前点,继续递归 由于原树可以长的奇形怪状( ...

  7. BZOJ1095: [ZJOI2007]Hide 捉迷藏【线段树维护括号序列】【思维好题】

    Description 捉迷藏 Jiajia和Wind是一对恩爱的夫妻,并且他们有很多孩子.某天,Jiajia.Wind和孩子们决定在家里玩 捉迷藏游戏.他们的家很大且构造很奇特,由N个屋子和N-1条 ...

  8. bzoj千题计划245:bzoj1095: [ZJOI2007]Hide 捉迷藏

    http://www.lydsy.com/JudgeOnline/problem.php?id=1095 查询最远点对,带修改 显然可以用动态点分治 对于每个点,维护两个堆 堆q1[x] 维护 点分树 ...

  9. BZOJ1095 [ZJOI2007]Hide 捉迷藏

    动态树分治,用三个set分别维护每个重心到每一个子树的距离种类.每个重心所有子树的最大值和次大值.全局答案的最大值.复杂度O(nlogn^2) 代码 #include<cstdio> #i ...

随机推荐

  1. Java并发编程Semaphore

    信号量 信号量类Semaphore,用来保护对唯一共享资源的访问.一个简单的打印队列,并发任务进行打印,加入信号量同时之能有一个线程进行打印任务 . import java.util.concurre ...

  2. [Python Study Notes]进程信息(丁丁软件监控进程,http-post)

    ''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' ...

  3. php生成文字水印和图片水印

    生成文字水印   //文字水印 /*打开图片*/ //1.配置图片路径 $src = "4.jpg"; //2.获取图片的信息(得到图片的基本信息) $info = getimag ...

  4. word在页眉中插入页码

    编辑页眉时,插入-页码-当前位置-普通数字

  5. 中小研发团队架构实践之应用监控Metrics

    一.Metrics简介        应用监控系统Metrics由Metrics.NET+InfluxDB+Grafana组合而成,通过客户端Metrics.NET在业务代码中埋点,Metrics.N ...

  6. 项目构建工具Maven

  7. POJ - 2698 贪心

    经典面替换算法,每次选择最远的那个碟片请求进行替换. AC代码 #include <cstdio> #include <cmath> #include <algorith ...

  8. Centos搭建mysql/Hadoop/Hive/Hbase/Sqoop/Pig

    目录: 准备工作 Centos安装 mysql Centos安装Hadoop Centos安装hive JDBC远程连接Hive Hbase和hive整合 Centos安装Hbase 准备工作: 配置 ...

  9. java中的mvc和三层结构究竟是什么关系

    一件事,要知其然往往很简单,要知其所以然通常不是那么容易,就如最近重新巩固spring的过程中,就觉得还有许多问题其实并不是十分明了. 屈指一算,手头上做过的正式项目也有了四五六七个了,不管用的数据库 ...

  10. UBOOT添加命令的执行流程

    BootLoader(引导装载程序)是嵌入式系统软件开发的第一个环节,它把操作系统和硬件平台衔接在一起,对于嵌入式系统的后续软件开发十分重要,在整个开发中也占有相当大的比例.U-BOOT是当前比较流行 ...