[NOIp2018提高组]赛道修建

题目大意:

给你一棵\(n(n\le5\times10^4)\)个结点的树,从中找出\(m\)个没有公共边的路径,使得第\(m\)长的路径最长。问第\(m\)长的路径最长可以是多少。

思路:

二分答案+树形DP。\(f[x]\)表示以\(x\)为根的子树中最多能找出几个长度\(\ge k\)的路径。\(g[x]\)表示去掉已经满足的路径,从\(x\)子树内往上连的最长的路径有多长。

转移时将所有子结点的贡献\(g[y]+w\)排序。若贡献已经\(\ge k\),那么就直接计入答案。否则从小到大枚举每一个贡献,找到能与其配对的最小的贡献,计入答案。如果找不到能与之配对的贡献,那么就用它来更新\(g[x]\)。可以证明这样能够在保证\(f[x]\)最大化的情况下,最大化\(g[x]\)。

时间复杂度\(\mathcal O(n\log n\log\)值域\()\)。

源代码:

  1. #include<set>
  2. #include<cstdio>
  3. #include<cctype>
  4. #include<vector>
  5. #include<climits>
  6. #include<algorithm>
  7. inline int getint() {
  8. register char ch;
  9. while(!isdigit(ch=getchar()));
  10. register int x=ch^'0';
  11. while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
  12. return x;
  13. }
  14. const int N=5e4+1;
  15. struct Edge {
  16. int to,w;
  17. };
  18. std::vector<Edge> e[N];
  19. inline void add_edge(const int &u,const int &v,const int &w) {
  20. e[u].push_back((Edge){v,w});
  21. e[v].push_back((Edge){u,w});
  22. }
  23. std::multiset<int> t;
  24. int f[N],g[N],len;
  25. void dfs(const int &x,const int &par) {
  26. f[x]=0;
  27. for(auto &j:e[x]) {
  28. const int &y=j.to;
  29. if(y==par) continue;
  30. dfs(y,x);
  31. f[x]+=f[y];
  32. }
  33. for(auto &j:e[x]) {
  34. const int &y=j.to;
  35. if(y==par) continue;
  36. t.insert(g[y]+j.w);
  37. }
  38. while(!t.empty()) {
  39. const int u=*t.rbegin();
  40. if(u>=len) {
  41. f[x]++;
  42. t.erase(t.find(u));
  43. } else {
  44. break;
  45. }
  46. }
  47. g[x]=0;
  48. while(!t.empty()) {
  49. const int u=*t.begin();
  50. t.erase(t.begin());
  51. auto p=t.lower_bound(len-u);
  52. if(p==t.end()) {
  53. g[x]=u;
  54. } else {
  55. t.erase(p);
  56. f[x]++;
  57. }
  58. }
  59. t.clear();
  60. }
  61. inline int calc(const int &k) {
  62. len=k;
  63. dfs(1,0);
  64. return f[1];
  65. }
  66. int main() {
  67. const int n=getint(),m=getint();
  68. int l=INT_MAX,r=0;
  69. for(register int i=1;i<n;i++) {
  70. const int u=getint(),v=getint(),w=getint();
  71. add_edge(u,v,w);
  72. l=std::min(l,w);
  73. r+=w;
  74. }
  75. while(l<=r) {
  76. const int mid=(l+r)>>1;
  77. if(calc(mid)>=m) {
  78. l=mid+1;
  79. } else {
  80. r=mid-1;
  81. }
  82. }
  83. printf("%d\n",l-1);
  84. return 0;
  85. }

[NOIp2018提高组]赛道修建的更多相关文章

  1. [NOIp2018提高组]旅行

    [NOIp2018提高组]旅行: 题目大意: 一个\(n(n\le5000)\)个点,\(m(m\le n)\)条边的连通图.可以从任意一个点出发,前往任意一个相邻的未访问的结点,或沿着第一次来这个点 ...

  2. [NOIp2018提高组]货币系统

    [NOIp2018提高组]货币系统 题目大意: 有\(n(n\le100)\)种不同的货币,每种货币的面额为\([1,25000]\)之间的一个整数.若两种货币系统能够组合出来的数是相同的的,那我们就 ...

  3. [NOIp2013提高组]积木大赛/[NOIp2018提高组]铺设道路

    [NOIp2013提高组]积木大赛/[NOIp2018提高组]铺设道路 题目大意: 对于长度为\(n(n\le10^5)\)的非负数列\(A\),每次可以选取一个区间\(-1\).问将数列清零至少需要 ...

  4. NOIP2018提高组省一冲奖班模测训练(六)

    NOIP2018提高组省一冲奖班模测训练(六) https://www.51nod.com/Contest/ContestDescription.html#!#contestId=80 20分钟AC掉 ...

  5. NOIP2018提高组省一冲奖班模测训练(五)

    NOIP2018提高组省一冲奖班模测训练(五) http://www.51nod.com/Contest/ContestDescription.html#!#contestId=79 今天有点浪…… ...

  6. NOIP2018提高组金牌训练营——动态规划专题

    NOIP2018提高组金牌训练营——动态规划专题 https://www.51nod.com/Live/LiveDescription.html#!#liveId=19 多重背包 二进制优化转化成01 ...

  7. NOIP2018提高组省一冲奖班模测训练(四)

    NOIP2018提高组省一冲奖班模测训练(四) 这次比赛只AC了第一题,而且花了40多分钟,貌似是A掉第一题里面最晚的 而且还有一个半小时我就放弃了…… 下次即使想不出也要坚持到最后 第二题没思路 第 ...

  8. NOIP2018提高组省一冲奖班模测训练(三)

    NOIP2018提高组省一冲奖班模测训练(三) 自己按照noip的方式考,只在最后一两分钟交了一次 第一题过了,对拍拍到尾. 第二题不会.考试时往组合计数的方向想,推公式,推了一个多小时,大脑爆炸,还 ...

  9. NOIP2018提高组金牌训练营——搜索专题

    NOIP2018提高组金牌训练营——搜索专题 1416 两点 福克斯在玩一款手机解迷游戏,这个游戏叫做”两点”.基础级别的时候是在一个n×m单元上玩的.像这样: 每一个单元有包含一个有色点.我们将用不 ...

随机推荐

  1. C++ Primer 笔记——控制内存分配

    1.当我们使用一条new表达式时,实际执行了三步操作. new表达式调用一个名为operator new ( 或者operator new[] ) 的标准库函数.该函数分配一块足够大的,原始的,未命名 ...

  2. 微信如何获取unionid 并且打通微信公众号和小程序

    准备 1.微信公众号 2.微信小程序 3.微信开发平台帐号 没有在开发平台绑定的小程序和公众号是没法获取unionid的 只需要在开发平台绑定小程序和公众号,便可以获取unionid 其中对于小程序和 ...

  3. 论文阅读笔记十八:ENet: A Deep Neural Network Architecture for Real-Time Semantic Segmentation(CVPR2016)

    论文源址:https://arxiv.org/abs/1606.02147 tensorflow github: https://github.com/kwotsin/TensorFlow-ENet ...

  4. python函数之各种器

    一: 装饰器 1:装饰器模板 def wrapper(func): def inner(*args,**kwargs): ret =func(*args,**kwargs) return ret re ...

  5. python基础面试题(一)

    1.   简述Python代码的运行机制 1.把原始代码编译成字节码         编译后的字节码是特定于Python的一种表现形式,它不是二进制的机器码,需要进一步编译才能被机器执行. 2.把编译 ...

  6. mybatis 遍历map;

    mybatis 遍历map; 参考http://blog.csdn.net/hj7jay/article/details/78652050 ps: ${m[key]}这是显示 打印的key读value ...

  7. C# float与UInt16互转

    //float拆分成两个UInt16 public static UInt16 FloatToTwoUInt16(float f) { byte[] bs = BitConvert.GetBytes( ...

  8. [转] 一张图理解prototype、proto和constructor的三角关系

    前面的话 javascript里的关系又多又乱.作用域链是一种单向的链式关系,还算简单清晰:this机制的调用关系,稍微有些复杂:而关于原型,则是prototype.proto和constructor ...

  9. 【CF446D】DZY Loves Games

    题解: 不错的题目 首先要求的黑点个数非常多 比较容易想到矩阵乘法 于是我们可以求出从某个黑点出发到任意一个黑点之间的概率 发现不同出发点带来的变化只有常数项 于是我们可以预处理出从每个方程转移的系数 ...

  10. 【AtCoder】ARC075

    ARC075 在省选前一天听说正式选手线画到省二,有了别的女选手,慌的一批,然后刷了一个ARC来稍微找回一点代码感觉 最后还是挂分了,不开心 果然水平退化老年加重啊 原题链接 C - Bugged 直 ...