树的合并(ybtoj-树上dp)

文章目录

  • 题目描述
  • 前言
  • 解析
  • 代码
  • thanks for reading!

题目描述

在这里插入图片描述

前言

全网唯一AC!!!
妙啊
而且还是完全自己想出来的做法
开心
(APIO还是没白听)

但是思路出来后代码实现十分坎坷
建两个图分别dfs3次那个地方就算是写不明白了
代码能力需要加强

解析

定义数组dis:

dis[i]表示i在自己那棵树上所能连出的最长简单路径

不难看出:
对于将两棵树上的边x,y相连时
新的直径要么是两棵树原来一棵的直径,要么就是经过了x和y
也就是:

max3(直径A,直径B,dis[x]+dis[y]+1)

通过对两棵树分别进行dfs,我们是可以预处理出两棵树的直径和dis数组的
那么:

X=max(直径A,直径B)

就成了定值

那么我们现在考虑如何找出所有的x,y,使得二者相连的的贡献大于X
也就是:

dis[x]+dis[y]>X

我们发现:
可以把A树的dis降序排序一下
然后再枚举B的dis值
这样在A的dis序列中与B当前的dis符合条件的点分布就是单调的
可以通过二分求解
假设对于B的dis1,A的dis序列的1-p的点i满足:

dis1+disA[i]>X

那么通过dis1得到的贡献之和就是

∑(disA[i]+dis1+1)(1<=i<=p)+X*(n-p)

把西格玛拆开:

∑disA[i] (1<=i<=p)+(dis1+1)p+X(n-p)

前面的西格玛可以用前缀和预处理处来
这样我们就可以在二分求出p后O1求出dis1的贡献了
总时间复杂度:
(n+m)logn
(从复杂度可以看出来,还可以在m<n时将两棵树交换以优化复杂度,不过本题并不需要了)

代码

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+100;
#define ll long long
int n,m,root,mod;struct node{int to,nxt;
}p[N*4];
int fi[N],cnt=-1;
void addline(int x,int y){p[++cnt]=(node){y,fi[x]};fi[x]=cnt;
//	printf("x=%d y=%d id=%d\n",x,y,cnt);
}
int fi2[N];
void addline2(int x,int y){p[++cnt]=(node){y,fi2[x]};fi2[x]=cnt;
}int dis2A[N],disA[N],dis2B[N],disB[N];
int r=-1;
int xa,ya,xb,yb;
ll ans;
void dfsA(int x,int f){
//	printf("f=%d x=%d dis=%d\n",f,x,disA[x]);for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;
//		printf("id=%d to=%d\n",i,to);if(to==f) continue;disA[to]=disA[x]+1;if(r==-1||disA[r]<disA[to]) r=to;dfsA(to,x);}
}
void dfs2A(int x,int f){
//	printf("f=%d x=%d dis=%d\n",f,x,disA[x]);for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;
//		printf("id=%d to=%d\n",i,to);if(to==f) continue;dis2A[to]=dis2A[x]+1;dfs2A(to,x);}
}
void dfsB(int x,int f){for(int i=fi2[x];~i;i=p[i].nxt){int to=p[i].to;if(to==f) continue;disB[to]=disB[x]+1;if(r==-1||disB[r]<disB[to]) r=to;dfsB(to,x);}
}
void dfs2B(int x,int f){for(int i=fi2[x];~i;i=p[i].nxt){int to=p[i].to;if(to==f) continue;dis2B[to]=dis2B[x]+1;dfs2B(to,x);}
}int X;
void solve(){dfsA(1,0);xa=r;r=-1;memset(disA,0,sizeof(disA));dfsA(xa,0);ya=r;
//	for(int i=1;i<=n;i++) printf("i=%d dis=%d\n",i,disA[i]);dfs2A(ya,0);for(int i=1;i<=n;i++) disA[i]=max(disA[i],dis2A[i]);X=disA[xa];
//	for(int i=1;i<=n;i++) printf("i=%d dis=%d\n",i,disA[i]);dfsB(1,0);xb=r;r=-1;memset(disB,0,sizeof(disB));dfsB(xb,0);yb=r;dfs2B(yb,0);for(int i=1;i<=m;i++) disB[i]=max(disB[i],dis2B[i]);
//	for(int i=1;i<=m;i++) printf("i=%d dis=%d\n",i,disB[i]);X=max(X,disB[xb]);//	printf("%d %d %d %d X=%d\n",xa,ya,xb,yb,X);
}bool cmp(int x,int y){return x>y;}
ll sum[N];
void compute(){sort(disB+1,disB+m+1,cmp);for(int i=1;i<=m;i++){sum[i]=sum[i-1]+disB[i];}for(int i=1;i<=n;i++){int now=disA[i];int st=0,ed=m;while(st<ed){int mid=(st+ed+1)>>1;if(disB[mid]+now+1>=X) st=mid;else ed=mid-1;}ans+=(ll)sum[st]+(ll)st*now+st+(ll)(m-st)*X;}
}
int main(){memset(fi,-1,sizeof(fi));memset(fi2,-1,sizeof(fi2));int a,b,c,d;scanf("%d%d",&n,&m);for(int i=1;i<n;i++){scanf("%d%d",&a,&b);addline(a,b);addline(b,a);}for(int i=1;i<m;i++){scanf("%d%d",&c,&d);addline2(c,d);addline2(d,c);}
//	printf("dkesf=%d\n",p[0].to);solve();compute();printf("%lld",ans);
}
/*
3 5
1 2
1 3
1 2
1 3
3 4
4 5
*/

thanks for reading!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318558.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

牛客题霸 [矩阵元素查找] C++题解/答案

牛客题霸 [矩阵元素查找] C题解/答案 题目描述 已知int一个有序矩阵mat&#xff0c;同时给定矩阵的大小n和m以及需要查找的元素x&#xff0c;且矩阵的行和列都是从小到大有序的。设计查找算法返回所查找元素的二元数组&#xff0c;代表该元素的行号和列号(均从零开始)。保证元…

51nod1821-最优集合【贪心】

正题 题目链接:http://www.51nod.com/Challenge/Problem.html#problemId1821 题目大意 nnn个可重集合&#xff0c;TTT次询问&#xff0c;询问将集合SxS_xSx​和一个选出大小为kkk的SyS_ySy​的子集合并成的独立集S′SS′满足能够找到一个最大的zzz使得对于任意的正整数i∈[1,z…

.net core中的高效动态内存管理方案

.net core在新增的System.Buffers中引入了一大堆高效内存管理的类&#xff0c;如span和memory、内存池。本文今天这里介绍一个高效动态内存访问方案。ReadOnlySequenceSegment<T>在我们读取数据的过程&#xff0c;很多时候会出现如下场景&#xff1a;不知道数据实际大小一…

C. Sum of Log(数位dp)

C. Sum of Log Code1 暴力记的状态 TLE #include<bits/stdc.h> using namespace std; using lllong long; template <class Tint> T rd() {T res0;T fg1;char chgetchar();while(!isdigit(ch)) {if(ch-) fg-1;chgetchar();}while( isdigit(ch)) res(res<<…

APIO游记(5.19-5.23)

概要 本来以为钱基本是白花了 但惊喜的发现其实没有 APIO之后做了几道题&#xff0c;感觉对题目的分析能力增强了 大概虽然不会任何前置算法&#xff0c;但其中的思路多多少少还是吸收到一些 还有一枚Cu诶 针不戳 课程 5.19决策单调性与四边形不等式 第一天 见识到了天高地…

.net core 上 K8S(三)Yaml文件运行.netcore程序

正文上一章我们通过kubectl run简单运行了一个.netcore网站&#xff0c;但实际的开发中&#xff0c;我们都是通过yaml来实现的。1.编写yaml文件关于yaml文件的格式在此就不多描述了&#xff0c;不熟悉的可以去网上搜一下示例。2.运行yamlkubectl create -f netcore.yaml 我们可…

51nod1676-无向图同构【乱搞】

正题 题目连接:http://www.51nod.com/Challenge/Problem.html#problemId1676 题目大意 给出两张nnn个点mmm条边的无向图&#xff0c;求这两张图是否同构。 1≤n≤200,1≤m≤4000,1≤T≤201\leq n\leq 200,1\leq m\leq 4000,1\leq T\leq 201≤n≤200,1≤m≤4000,1≤T≤20 解题…

Jozky模板

文章目录字符串处理后缀数组manacherhashKMP最大最小表达法数论约瑟夫环欧拉函数莫比乌斯反演逆序对归并排序求逆序对素数线性筛欧几里得与扩展欧几里得欧几里得算法&#xff1a;扩展欧几里得算法&#xff1a;逆元扩展欧几里得费马小定理欧拉定理递推求逆元__int128高精度运算唯…

Visual Studio 2017 15.9 版本发布:推出全新的导入 / 导出配置功能

Microsoft 在开发 Visual Studio 2019 的同时&#xff0c;还在继续支持 VS2017 的用户。公司已经发布了 9 次更新&#xff0c;这展示了 Microsoft 在常规更新发布之后仍然会坚守继续支持 Visual Studio 的承诺。我们已经介绍过 15.9 版本中的一些新增内容&#xff0c;但是在最终…

P4780-Phi的反函数【dfs】

正题 题目链接:https://www.luogu.com.cn/problem/P4780 题目大意 给出nnn&#xff0c;求一个最小的xxx满足φ(x)n\varphi(x)nφ(x)n。 若不存在或者大于2312^{31}231则输出−1-1−1。 1≤n≤2311\leq n\leq 2^{31}1≤n≤231 解题思路 考虑用φ\varphiφ比较常用的公式&…

染色(树链剖分 洛谷-P2486)

文章目录题目描述解析代码thanks for reading&#xff01;传送门首先&#xff0c;对hash学姐对本题拔刀相助的debug行为表示衷心的感谢 题目描述 解析 用线段树维护颜色序列个数、最左颜色与最右颜色 合并时如果左儿子的最右颜色等于右儿子的最左颜色&#xff0c;就把加和-1 在…

第 45 届国际大学生程序设计竞赛(ICPC)亚洲区域赛(昆明) AC(带悔贪心)

下面两个题都是选择iii有个价值如果选择iii则不能选择i−1,i1i-1,i1i−1,i1&#xff0c;让价值最大或最小 P1792 [国家集训队]种树 野心qwq 的博客 #include<bits/stdc.h>using namespace std;const int N200010; int n,m; int a[N]; int fr[N],to[N]; priority_queue…

C# 8中的范围类型(Range Type)

C# 8.0中加入了一个新的范围类型(Range Type)。这里我们首先展示一些代码&#xff0c;并一步一步为代码添加一些不同的东西, 为大家展示一下范围类型的功能和用法。我们最原始的代码如下&#xff1a;这里我们显示的定义了我们查询数组的索引1-3, 并输出他们的值。毫无疑问&…

bzoj#4423-[AMPPZ2013]Bytehattan【并查集】

正题 题目链接:https://darkbzoj.tk/problem/4423 题目大意 给出一个n∗nn*nn∗n的网格图&#xff0c;然后四联通的点之间连接。每次删掉一条边求这条边的两个点是否连通。强制在线。 1≤n≤1500,1≤m≤2n(n−1)1\leq n\leq 1500,1\leq m\leq 2n(n-1)1≤n≤1500,1≤m≤2n(n−…

选数游戏(ybtoj-二叉堆)

文章目录题目描述解析代码thanks for reading&#xff01;题目描述 解析 一道很考验代码能力与思维的题 &#xff08;我不是在为自己的菜找理由&#xff09; 首先由于可以每一列都有类似于环的性质 所以我们可以忽略点的出入的纵坐标的位置&#xff0c;只考虑每列选几个 首先&…

CF741C Arpa’s overnight party and Mehrdad’s si

题目描述&#xff1a; 有n对情侣&#xff08;2n个人&#xff09;围成一圈坐在桌子边上&#xff0c;每个人占据一个位子&#xff0c;要求情侣不能吃同一 种食物&#xff0c;并且桌子上相邻的三个人的食物必须有两个人是不同的&#xff0c;只有两种食物&#xff08;1或者是2&…

开源库支付库Magicodes.Pay发布

Magicodes.Pay&#xff0c;是心莱科技团队提供的统一支付库&#xff0c;相关库均使用.NET标准库编写&#xff0c;支持.NET Framework以及.NET Core。目前支持以下支付方式和功能&#xff1a;支付宝APP支付支付宝Wap支付支付宝国际支付 支持分账微信小程序支付微信APP支付统一支…

火车载客(ybtoj-二叉堆)

文章目录题目描述解析我的思路代码题解思路题目描述 解析 我的思路 其实就是线段覆盖的一个变体 贪心的想&#xff1a; 把游客按右端点升序排序 后面的证明就和线段覆盖一样了 如果有两个游客冲突 我们应该选右端点靠右的 因为这样对以后继续在右边出现的游客来说肯定不会更差…

CF835E-The penguin‘s game【交互】

正题 题目链接:https://www.luogu.com.cn/problem/CF835E 题目大意 长度为nnn的序列中有两个yyy其他都是xxx&#xff0c;给出n,x,yn,x,yn,x,y。你每次可以询问一个下标集合的数字异或和&#xff0c;要求在191919次以内找到这两个yyy的位置。 1≤n≤1000,1≤x,y≤109,x≠y1\le…