BZOJ4543 [POI2014]Hotel加强版
Hotel加强版
有一个树形结构,每条边的长度相同,任意两个节点可以相互到达。选3个点。两两距离相等。有多少种方案?
数据范围:n<=100000
yyb的题解
我们先考虑一个\(O(n^2)\)的dp,也就是原题的做法。
我们考虑一下,三个点两两的距离相同是什么情况,
存在一个三个点公共的LCA,所以我们在LCA统计答案即可。
存在一个点,使得这个点到另外两个子树中距离它为d的点以及这个点的d次祖先。
所以,设计DP状态为
\(f[i][j]\)表示以\(i\)为根的子树中,距离当前点为\(j\)的点数。
\(g[i][j]\)表示以\(i\)为根的子树中,两个点到LCA的距离为\(d\),并且他们的LCA到\(i\)的距离为\(d−j\)的点对数,简单来说就是\(i\)往其他地方走\(j\)步就能找到一组解。
考虑合并的时候的转移:
\]
转移的正确性比较显然,不在多讲了,并不是这里的重点。这样子的复杂度是\(O(n^2)\)的。
我们观察一下转移的时候有这样两步:
\]
如果我们钦定一个儿子的话,那么这个数组是可以直接赋值的,并不需要再重复计算。
所以我们用指针来写,也就是:\(f[i]=f[son]−1,g[i]=g[son]+1\)。
如果整棵树是链我们发现复杂度可以做到O(n),既然如此,我们推广到树。我们进行长链剖分,每次钦定从重儿子直接转移,那么我们还需要从轻儿子进行转移。不难证明所有轻儿子都是一条重链的顶部,转移时的复杂度是重链长度。
那么,复杂度拆分成两个部分:直接从重儿子转移\(O(1)\),从轻儿子转移\(O(∑len)\)。发现每个点有且仅有一个父亲,因此一条重链算且仅被一个点暴力转移,而每次转移复杂度是链长。所以全局复杂度是∑链长,也就是\(O(n)\),因此总复杂度就是\(O(n)\)。
这样子写下来,发现长链剖分之后,我们的复杂度变为了线性。但是注意到复杂度证明中的一点:转移和链长相关。而链长和什么相关呢?深度。所以说对于这一类与深度相关的、可以快速合并的信息,使用长链剖分可以优化到一个非常完美的复杂度。如果需要维护的与深度无关的信息的话,或许dsu on tree是一个更好的选择。
代码
DP是的for是在用相对深度,比较简单的实现方法是之前统计重儿子的时候用高度代替深度。
然后tmp必须开到4倍是因为g数组指针给儿子的时候在前移。
co int N=1e5+1;
int n,head[N],to[N*2],nx[N*2],tot;
void add(int x,int y){to[++tot]=y,nx[tot]=head[x],head[x]=tot;}
int dep[N],md[N],son[N];
void dfs1(int x,int fa){
for(int i=head[x];i;i=nx[i]){
int y=to[i];if(y==fa) continue;
dfs1(y,x),md[x]=std::max(md[x],md[y]);
if(md[y]>md[son[x]]) son[x]=y;
}
md[x]=md[son[x]]+1;
}
ll*f[N],*g[N],tmp[N*4],*id=tmp,ans;
void dfs2(int x,int fa){
if(son[x]) f[son[x]]=f[x]+1,g[son[x]]=g[x]-1,dfs2(son[x],x);
f[x][0]=1,ans+=g[x][0];
for(int i=head[x];i;i=nx[i]){
int y=to[i]; if(y==fa||y==son[x]) continue;
f[y]=id,id+=md[y]*2,g[y]=id,id+=md[y]*2;
dfs2(y,x);
for(int j=0;j<md[y];++j){
if(j)ans+=f[x][j-1]*g[y][j];
ans+=g[x][j+1]*f[y][j];
}
for(int j=0;j<md[y];++j){
g[x][j+1]+=f[x][j+1]*f[y][j];
if(j)g[x][j-1]+=g[y][j];
f[x][j+1]+=f[y][j];
}
}
}
int main(){
// freopen("BZOJ4543.in","r",stdin);
// freopen(".out","w",stdout);
read(n);
for(int i=1,x,y;i<n;++i){
read(x),read(y);
add(x,y),add(y,x);
}
dfs1(1,0);
f[1]=id,id+=md[1]*2,g[1]=id,id+=md[1]*2;
dfs2(1,0);
printf("%lld\n",ans);
return 0;
}
BZOJ4543 [POI2014]Hotel加强版的更多相关文章
- BZOJ4543 POI2014 Hotel加强版 【长链剖分】【DP】*
BZOJ4543 POI2014 Hotel加强版 Description 同OJ3522 数据范围:n<=100000 Sample Input 7 1 2 5 7 2 5 2 3 5 6 4 ...
- bzoj4543 [POI2014]Hotel加强版 长链剖分+树形DP
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4543 题解 这道题的弱化版 bzoj3522 [POI2014]Hotel 的做法有好几种吧. ...
- BZOJ4543[POI2014]Hotel加强版——长链剖分+树形DP
题意参见BZOJ3522 n<=100000 数据范围增强了,显然之前的转移方程不行了,那么不妨换一种. 因为不能枚举根来换根DP,那么我们描述的DP方程每个点要计算三个点都在这个点的子树内的方 ...
- 2019.01.08 bzoj4543: [POI2014]Hotel加强版(长链剖分+dp)
传送门 代码: 长链剖分好题. 题意:给你一棵树,问树上选三个互不相同的节点,使得这个三个点两两之间距离相等的方案数. 思路: 先考虑dpdpdp. fi,jf_{i,j}fi,j表示iii子树中离 ...
- 【BZOJ4543】[POI2014]Hotel加强版 长链剖分+DP
[BZOJ4543][POI2014]Hotel加强版 Description 同OJ3522数据范围:n<=100000 Sample Input 7 1 2 5 7 2 5 2 3 5 6 ...
- 【BZOJ4543】Hotel加强版(长链剖分)
[BZOJ4543]Hotel加强版(长链剖分) 题面 BZOJ,没有题面 洛谷,只是普通版本 题解 原来我们的\(O(n^2)\)做法是设\(f[i][j]\)表示以\(i\)为根的子树中,距离\( ...
- 4543: [POI2014]Hotel加强版
4543: [POI2014]Hotel加强版 链接 分析: f[u][i]表示子树u内,距离u为i的点的个数,g[u][i]表示在子树u内,已经选了两个深度一样的点,还需要在距离u为i的一个点作为第 ...
- 【BZOJ4543】Hotel加强版
[BZOJ4543]Hotel加强版 题面 bzoj 洛谷 $ps:$在洛谷看题在bzoj交... 题解 我们分析一下这个问题,要怎么样的点才满足三点距离两两相等呢? 1.存在三个点有共同的$LCA$ ...
- bzoj4543[POI2014]Hotel
题目链接 bzoj4543 [POI2014]Hotel 题解 这不是裸地点分嘛 ,我真傻,真的 n^2 这不是是sb题,~滑稽 ~ 枚举点转换为无根树,暴力子树中点的深度 计数转移 令a b c d ...
随机推荐
- 《深入理解mybatis原理6》 MyBatis的一级缓存实现详解 及使用注意事项
<深入理解mybatis原理> MyBatis的一级缓存实现详解 及使用注意事项 0.写在前面 MyBatis是一个简单,小巧但功能非常强大的ORM开源框架,它的功能强大也体现在它的缓 ...
- python 连接操作 各类数据库
转载自MySQL Loners 一,python 操作 MySQL:详情见:这里 #!/bin/env python # -*- encoding: utf-8 -*- #-------------- ...
- 初识async函数
为什么会出现async函数 首先从大的方面来说,出现async函数时为了解决JS编程中的异步操作,再往具体说就是为了对以往异步编程方法的一种改进,也有人说仅仅只是Generator 函数的语法糖,这个 ...
- 杀死dialog
先 pkill -9 normal.sh 和 pkill -9 terminal_ui.sh 然后在pkill -9 dialog
- python 语法错误记录
1 Missing parameter end_time in docstring less... (Ctrl+F1) 参数位置错误 注意:只有在形参表末尾的那些参数可以有默认参数值,即你不能在声明 ...
- laravel中上传图片之后图片的处理
$file=Input::file('file'); if ($file->isValid()){ $entension=$file->getClientOriginalExtension ...
- gitblit-部署
什么是 Gitblit Gitblit是一个开源的用于管理,查看和提供Git仓库. 它主要设计为希望托管集中存储库的小工作组的工具. Gitblit有什么特点 Gitblit部署示例1 日常维护添加步 ...
- 使用 create-react-app 构建 react应用程序
原文 http://blog.csdn.net/github_squad/article/details/57452333#
- Mysql双机热备--预备知识
1.双机热备 对于双机热备这一概念,我搜索了很多资料,最后,还是按照大多数资料所讲分成广义与狭义两种意义来说. 从广义上讲,就是对于重要的服务,使用两台服务器,互相备份,共同执行同一服务.当一台服务器 ...
- 修改myelipse中部署路径deploy location内容的方法
在new web project 中的project name等内容,可以打开.mymetadata文件进行修改 <?xml version="1.0" encoding=& ...