由题开始==

例题:求在一棵有权树上,是否存在一条路径满足权值和为K

解法:以每个点为根一次,看在他的子树间是否存在两段,其和为K;O(==)

和例题一样,对于树上问题,求某些要求的路径(数量或者存在性等),

往往可以先对一条经过根节点的路径操作,

后再删去这个根,对他的子树们同样的操作

显然是一个分治过程

原理就是,一条路径,要么是由一个点经过根节点,与其他子树内的节点形成;

要么就是只在这棵子树内形成路径

大概图示意思(红绿为两条上述路径)



<

当我们的树比较平衡时,每个点被路径计算是$ logn $ 的,但是当树是一条链的时候,就退化成$ n^2 $ 了

为了避免这种情况,可以用树的重心代替,成为新的根。此时总复杂度为$ O (nlogn) 。$

原因就是,一棵树,怎么为根都还是一颗树,但是以重心为根的时候,这颗树是最好看的最平衡的

rt,将链按箭头提起来:





这看向去更像是一棵树

所以

总结一下

一棵树先确定他的重心,以重心为根,确定经过根节点的路径;再把根节点删掉,对于删掉他的子树们,按上述同样操作;

已证$ O(nlogn) $。

luogu模板题代码仅供参考,不解释。

#include<bits/stdc++.h>
using namespace std;
const int N = 1e4+50;
const int K = 1e7+50;
int n,m;
struct node{int next,to,dis;}edge[N<<1];
int head[N],cnt;
inline void add(int from,int to,int dis) {
edge[++cnt].to=to,edge[cnt].dis=dis,edge[cnt].next=head[from],head[from]=cnt;
}
int q[N],ans[N],maxp[N],size[N],visited[N],tmp[N],dis[N],judge[K];
int rt,sum,tot;
void getrt(int u,int f) {
size[u]=1,maxp[u]=0;
for(int i=head[u];i;i=edge[i].next) {
int v=edge[i].to;if(v==f||visited[v]) continue;
getrt(v,u);
size[u]+=size[v];maxp[u]=max(maxp[u],size[v]);
}
maxp[u]=max(maxp[u],sum-maxp[u]);
if(maxp[u]<maxp[rt]) rt=u;//要求最大的最小
}
void getdis(int u,int f) {
tmp[++tot]=dis[u];
for(int i=head[u];i;i=edge[i].next) {
int v=edge[i].to;if(v==f||visited[v]) continue;
dis[v]=dis[u]+edge[i].dis; getdis(v,u);
}
}
queue<int> que;
void solve(int u) {
for(int i=head[u];i;i=edge[i].next) {
int v=edge[i].to;if(visited[v]) continue;
dis[v]=edge[i].dis;
tot=0;getdis(v,u);
for(int j=1;j<=tot;j++)
for(int k=1;k<=m;k++)
if(q[k]>=tmp[j])
ans[k]|=judge[q[k]-tmp[j]];
for(int j=1;j<=tot;j++) que.push(tmp[j]),judge[tmp[j]]=1;
}
while(!que.empty()) judge[que.front()]=0,que.pop();//数组过大,memset超时
}
void divide(int u) {
judge[0]=visited[u]=1;solve(u);
for(int i=head[u];i;i=edge[i].next) {
int v=edge[i].to;if(visited[v]) continue;
maxp[rt=0]=sum=size[v];
getrt(v,0),getrt(rt,0);
divide(rt);
}
}
int main() {
cin>>n>>m;
for(int i=1;i<n;i++) { int a,b,c;scanf("%d%d%d",&a,&b,&c);add(a,b,c),add(b,a,c);}
for(int i=1;i<=m;i++) scanf("%d",q+i);
maxp[rt=0]=sum=n;//初始化
getrt(1,0),getrt(rt,0);//找重心
divide(rt);//点分治
for(int i=1;i<=m;i++) if(ans[i]) puts("AYE");else puts("NAY");
return 0;
}

skkyk:点分治的更多相关文章

  1. [bzoj2152][聪聪和可可] (点分治+概率)

    Description 聪聪和可可是兄弟俩,他们俩经常为了一些琐事打起来,例如家中只剩下最后一根冰棍而两人都想吃.两个人都想玩儿电脑(可是他们家只有一台电脑)……遇到这种问题,一般情况下石头剪刀布就好 ...

  2. POJ 2965. The Pilots Brothers' refrigerator 枚举or爆搜or分治

    The Pilots Brothers' refrigerator Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 22286 ...

  3. [poj1741][tree] (树/点分治)

    Description Give a tree with n vertices,each edge has a length(positive integer less than 1001). Def ...

  4. 【教程】简易CDQ分治教程&学习笔记

    前言 辣鸡蒟蒻__stdcall终于会CDQ分治啦!       CDQ分治是我们处理各类问题的重要武器.它的优势在于可以顶替复杂的高级数据结构,而且常数比较小:缺点在于必须离线操作. CDQ分治的基 ...

  5. BZOJ 3262 陌上花开 ——CDQ分治

    [题目分析] 多维问题,我们可以按照其中一维排序,然后把这一维抽象的改为时间. 然后剩下两维,就像简单题那样,排序一维,树状数组一维,按照时间分治即可. 挺有套路的一种算法. 时间的抽象很巧妙. 同种 ...

  6. BZOJ 1176 [Balkan2007]Mokia ——CDQ分治

    [题目分析] 同BZOJ2683,只需要提前处理s对结果的影响即可. CDQ的思路还是很清晰的. 排序解决一维, 分治时间, 树状数组解决一维. 复杂度是两个log [代码] #include < ...

  7. BZOJ 2683 简单题 ——CDQ分治

    [题目分析] 感觉CDQ分治和整体二分有着很本质的区别. 为什么还有许多人把他们放在一起,也许是因为代码很像吧. CDQ分治最重要的是加入了时间对答案的影响,x,y,t三个条件. 排序解决了x ,分治 ...

  8. HDU5977 Garden of Eden(树的点分治)

    题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=5977 Description When God made the first man, he ...

  9. Tsinsen A1493 城市规划(DP + CDQ分治 + NTT)

    题目 Source http://www.tsinsen.com/A1493 Description 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了. 刚才说过, 阿狸的国家有n个城市, 现在 ...

随机推荐

  1. C#中转换运算符explicit、implicit、operator、volatile研究

    C#中的这个几个关键字:explicit.implicit与operator,估计好多人的用不上,什么情况,这是什么?字面解释:explicit:清楚明白的;易于理解的;(说话)清晰的,明确的;直言的 ...

  2. 对python函数后面有多个括号的理解?

    一般而言,函数后面只有一个括号.如果看见括号后还有一个括号,说明第一个函数返回了一个函数,如果后面还有括号,说明前面那个也返回了一个函数.以此类推. 比如fun()() def fun(): prin ...

  3. PHP 日期之间所有日期

    /** * 获取起止日期之间所有日期 * @param $sdate * @param $edate * @return array */ function get_dates($sdate, $ed ...

  4. SQL Server优化技巧——如何避免查询条件OR引起的性能问题

    之前写过一篇博客"SQL SERVER中关于OR会导致索引扫描或全表扫描的浅析",里面介绍了OR可能会引起全表扫描或索引扫描的各种案例,以及如何优化查询条件中含有OR的SQL语句的 ...

  5. Cocos2d-x3.0网络通信学习(一)

    配置:win7+Cocos2d-x.3.0+VS2012 摘要:建立基本的http通信并得到返回信息. 一.添加项目与编译库 1.添加头文件 在需要用到Http网络相关类的文件中加入头文件 #incl ...

  6. 逆向学习周记-C语言空函数

    实验环境:WIN7虚拟机 软件:VC6 首先在VC6里面写一个空函数Fun(): F7编译运行一下,没有出错,接着在函数处使用F9下断点,使程序运行到Fun函数时停下. 接着F5开始运行这个程序 程序 ...

  7. Qt的诞生和本质

    有没有发现,在这个主函数里面没有了消息循环了.消息循环到哪里去了呢? 我们知道,每个GUI应用程序都有消息循环,一般都需要将消息循环while那个函数写到main的最后位置.那既然如此,为了代码复用, ...

  8. springboot 多环境

    springboot 多环境 --spring.profiles.active=dev 查看 Ioc 容器 PostProcessorRegistrationDelegate

  9. Git终端命令行的常用操作

    一.git源代码管理的优点 方便多人协同开发.工作 降低代码的管理成本 良好的分支管理机制 二.结构分析 服务端和客户端都有版本控制能力,都能进行代码的提交.合并 结构一: 结构二: 三.工作区的创建 ...

  10. Element-ui 下拉列表 选项过多时通过自定义搜索来解决卡顿问题

    当使用Select选择器时,如果下拉列表的数据量太多,会有一个明显的卡顿体验,例如: <!DOCTYPE html> <html lang="en"> &l ...